语音共振峰提取

%从LPC谱中找出四个共振峰值
function [amp]=find_peak（SwFs） 
%共振峰频率p，带宽b和峰值amp
%求线性预测系数
p=14;
lpcc=lpc（Swp）;

 result=roots（lpcc）;%roots是求多项式的根
 result=result（imag（result）>0.001）; %取虚部
 theta=angle（result）;
 fre=abs（theta）*Fs/（2*pi）;
 a=polyval（lpccresult）;
 b1=-log（abs（result））*Fs/pi;
 [pt]=sort（fre）;
 
 for i=1:numel（t）
      b（i）=b1（t（i））;
 end
       
 for i=1:length（p）