k).^2)/N;Рfigure(2);Рplot(f,10*log10(Pxx));Рxlabel('f/Hz');ylabel('10log(PSD)');Рtitle('周期图法谱估计 N=128');РAR模型自相关法:Рclear;РN=1024;Рn=1:N;Рp=300;Рf0=0.200;f1=0.213;Рxn=20^(1/2)*sin(2*pi*f0*(n-1))+2^(1/2)*sin(2*pi*f1*(n-1)) +0.1*randn(1,N);Рfigure(1);Рplot(n,xn);Рxlabel('n');Рtitle('仿真信号');РRx=zeros(1,p);Rx0=sum(xn.^2)/N;Р%% 求自相关函数的有偏估计Рfor k=1:pР Rx(k)=xn(1:N-k)*(xn(k+1:N)');РendРRx=Rx./N;Р%%РA=zeros(p,p);Р%% 构造自相关矩阵Рfor i=1:p-1Р A(i,i:p)=([Rx0,Rx(1:p-i)]);Р i=i+1;РendРA=A+A'-Rx0*eye(p,p);Р%%РAa=-Rx';Рb=-Rx';Рa=A\b;РG=0;Рfor i=1:pР ro=a(i)*Rx(i);Р G=G+ro;РendРG=G+Rx0;РLs=([1,a']);Рfigure(2);Р[h,w]=freqz(1,Ls,128,1);РPxx=G*(abs(h)).^2;Рplot(w,10*log10(Pxx));Рxlabel('frequency/Hz');ylabel('10log10(PSD)');Рtitle('自相关算法功率谱估计N=1024 p=300');Р РAR模型Burg算法:Рclear all;Рsyms x y;РN=input('请输入采样点数N=');