功率谱估计的MATLAB实现

(M-1,n-1);Р eb(M,n)=eb(M-1,n-1)+xishu*ef(M-1,n);Р endР M=M+1;Р A=0;Р D=0;Р for n=M:NР A=A+ef(M-1,n)*eb(M-1,n-1);Р D=D+(abs(ef(M-1,n)))^2+(abs(eb(M-1,n-1)))^2;Р endР xishu=-2*A/D;Р a(M-1,M-1)=-2*A/D;Р P(M)=P(M-1)*(1-(abs(xishu))^2);Р FPE(M-1)=P(M)*(N+M)/(N-M);Р for m=1:M-2Р a(M-1,m)=a(M-2,m)+xishu*a(M-2,M-1-m);Р endРendРT=1/fs;sum1=0;Рf=0.01:0.01:0.5;Р for m=1:M-1;Р sum1=sum1+a(M-1,m)*exp(-j*2*pi*m*f*T);Р endРs1=(abs(1+sum1)).^2;Рs=P(M)*T./s1;Рplot(f,10*log10(s),'k');Рxlabel('f/fs');Рylabel('功率谱/dB');Р2、周期图谱估计Рclc;clear;РN=32;РSNR=30; Рfs=1;Рt=1:N;Рt=t/fs;Рy=sin(2*pi*0.175*t)+sin(2*pi*0.20*t);Рx = awgn(y,SNR);Рsum1=0;Рf=0.05:0.01:0.5;Р for m=1:NР sum1=sum1+x(m)*exp(-j*2*pi*m*f);Р endРs=(abs(sum1)).^2/N;Рplot(f,10*log10(s),'k');Рxlabel('f/fs');Рylabel('功率谱/dB');Р实验结果:Р1、最大熵法估计结果:Р2、周期图法估计结果: