x=sheng_500co2_2;

fs=62500;

u=filter([1 -.99],1,x);                     % 預(yù)加重

wlen=length(u);                             % 幀長(zhǎng)

p=12;                                      % LPC階數(shù)

a=lpc(u,p);                                 % 求出LPC系數(shù)

U=lpcar2pf(a,256);                          % 由LPC系數(shù)求出頻譜曲線

freq=(0:256)*fs/512;                        % 頻率刻度

df=fs/512;                                  % 頻率分辨率

U_log=10*log10(U);                          % 功率譜分貝值

subplot 211; plot(u,'k');                   % 作圖

axis([0 wlen -0.5 0.5]);

title('預(yù)加重波形');

xlabel('樣點(diǎn)數(shù)'); ylabel('幅值')

subplot 212; plot(freq,U,'k');

title('聲道傳遞函數(shù)功率譜曲線');

xlabel('頻率/Hz'); ylabel('幅值');

[Loc,Val]=findpeaks(U);                     % 在U中尋找峰值

ll=length(Loc);                             % 有幾個(gè)峰值

for k=1 : ll

   m=Val(k);                              % 設(shè)置m-1,m和m+1

   m1=m-1; m2=m+1;

   p=Loc(k);                              % 設(shè)置P(m-1),P(m)和P(m+1)

   p1=U(m1); p2=U(m2);

   aa=(p1+p2)/2-p;                        % 按式(9-3-4)計(jì)算

   bb=(p2-p1)/2;

   cc=p;

   dm=-bb/2/aa;                           % 按式(9-3-6)計(jì)算

   pp=-bb*bb/4/aa+cc;                     % 按式(9-3-8)計(jì)算

   m_new=m+dm;

    bf=-sqrt(bb*bb-4*aa*(cc-pp/2))/aa;      % 按式(9-3-13)計(jì)算

   F(k)=(m_new-1)*df;                     % 按式(9-3-7)計(jì)算

   Bw(k)=bf*df;                           % 按式(9-3-14)計(jì)算

   line([F(k) F(k)],[0 pp],'color','k','linestyle','-.');

end

fprintf('F =%5.2f   %5.2f  %5.2f   %5.2f\n',F)

<kbd id="odjak"><wbr id="odjak"></wbr></kbd>