功率譜估計淺談-資料下載頁
2025-06-19 13:07本頁面
【正文】 % U2(k)是AR模型中的另一參數(shù)即誤差功率； FPE(k)=0。 U2(k+1)=0。endU2(1)=R(1)。 a(1,1)=R(2)/R(1)。 %由Levinson算法，計算一階模型參數(shù)a11；U2(2)=(1(abs(a(1,1)))^2)*U2(1)。 %由Levinson算法，計算一階模型參數(shù) 誤差功率；FPE(1)=U2(2)*(N+2)/N。 %計算一階模型的最終預(yù)測誤差準(zhǔn)則函數(shù)；S=0。for k=2:N %由Levinson梯歸公式計算K階模型參數(shù)和FPE函數(shù)； for l=1:k1 M=a(k1,l)*R(kl+1)。 S=S+M。 end a(k,k)=(R(k+1)+S)/U2(k)。 for i=1:k1 a(k,i)=a(k1,i)+a(k,k)*a(k1,ki)。 end U2(k+1)=(1(abs(a(k,k)))^2)*U2(k)。 FPE(k)=U2(k+1)*(N+k+1)/(Nk+1)。 S=0。end %梯歸計算完畢；%% 確定階數(shù)P min=FPE(1)。 %求出使得FPE函數(shù)取最小值的階數(shù)P；for k=2:N if FPE(k)min min=FPE(k)。 p=k。 endend %p=2 %為了調(diào)整效果可以在這里自行指定階數(shù)；disp(39。輸出模型參數(shù)a39。)。for k=1:p disp(a(p,k))。enddisp(39。輸出誤差功率39。)。disp(U2(p+1))。%%%%%%%%%%%%%%%%% AR模型參數(shù)確定后計算出功率譜Z=0。W=0::pi。 %功率譜以2*pi為周期，又信號為實信號，只需輸出0到PI即可；for k=1:p Z=Z+(a(p,k).*exp(j*k*W))。endPXX=U2(p+1)./((abs(1+Z)).^2)。 %得到功率譜函數(shù)；F=W*FS/(pi*2)。 %將角頻率坐標(biāo)換算成HZ坐標(biāo)，便于觀察；subplot(3,1,3)。plot(F,abs(PXX))。disp(39。輸出階數(shù)P39。)disp(p)。grid on。
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1