%從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）
???