一道用Levinson-Durbin迭代法計算AR模型的系數并估計信號功率譜的題，下面是寫的程序：

N=256;? ?? ?? ?? ? %信號觀測樣本長度

Nf=1024;

p=16;? ?? ?? ?? ???%AR模型階數

noise=(randn(1,N)+j*randn(1,N))/sqrt(2);

%產生三個復正弦信號

f1=0.15;

f2=0.17;? ?? ?? ? %信號歸一化頻率

f3=0.26;

SNR1=30;

SNR2=30;? ?? ?? ? %信號的信噪比

SNR3=27;

A1=10^(SNR1/20);

A2=10^(SNR2/20);??%信號的幅度(電壓)

A3=10^(SNR3/20);

signal1=A1*exp(j*2*pi*f1*(0:N-1));

signal2=A2*exp(j*2*pi*f2*(0:N-1)); %生成復正弦信號

signal3=A3*exp(j*2*pi*f3*(0:N-1));

un=signal1+signal2+signal3+noise;??%生成觀察樣本u(n)

r0=xcorr(un,p,'biased');? ?? ?? ???%直接計算自相關函數

r=r0(p+1:2*p+1);? ?? ?? ?? ?? ?? ? %提取r(0),r(1),...,r(p)

%計算一階AR模型的系數與輸入方差

a(1,1)=-r(2)/r(1);? ?? ?? ?? ?? ???%1階AR模型的輸入方差

sigma(1)=r(1)-(abs(r(2)^2)/r(1));??%1階AR模型的輸入方差

for m=2:p

k(m)=-(r(m+1)+sum(a(m-1,1:m-1).*r(m:-1:2)))/sigma(m-1);%反射系數

a(m,m)=k(m);

for i=1:m-1

a(m,i)=a(m-1,i)+k(m)*conj(a(m-1,m-i));

end? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?%m階AR模型的系數

sigma(m)=sigma(m-1)*(1-abs(k(m))^2);? ?? ?? ?? ?? ?%m階AR模型的輸入方差

end

Par=sigma(p)./fftshift(abs(fft([1,a(p,:)],Nf)).^2);? ? %p階AR模型的功率譜

求解最后如何歸一化功率譜啊，要求圖像縱坐標是“歸一化功率譜”，橫坐標是“w/2pi”。

2015-10-6 21:35 上傳

題干，第三問

2015-10-6 21:36 上傳

最后結果圖是這樣的

2015-10-6 21:36 上傳

我不知道怎么歸一化的最后

總結

以上是生活随笔為你收集整理的matlab 归一化功率谱,求大神指点如何画归一化功率谱的图像啊的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。