3.m

来自「滤波器阶数固定时,代价函数和信噪比之间的关系 信噪比固定时,代价函数和滤波器长度」· M 代码 · 共 31 行

%  eg__1.m

a1=0.3;
a2=0.05;
s=zeros(1000,1);
w=randn(1002,1);
s_1=w(1001);
s_2=w(1002);
for i=1:1000
    s(i)=a1*s_1+a2*s_2+w(i);     %生成s信号
    s_2=s_1;
    s_1=s(i);
end

Psav＝xcorr(s,1,'unbiased')
%二阶前项一步预测, 比阶数大一
xr=xcorr(x,3,'unbiased'); % xx
px=p(5:7); %不包含均方值
Pxav=xr(4);
for i=0.1:0.01:2
    Pvav=i;
    v=Pvav^0.5*randn(1000,1);
    x=s+v;
    [R1,R]=corrmtx(x,2,'autocorrelation');% 阶数
    R=inv(R);
    snr=Psav/Pvav;
    J=Pxavx+(px')*(R*px);
    plot(snr,J,'*');
    hold on;
end