eg__2.m

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

%  eg__2.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

v=randn(1000,1);
x=s+v;
figure(1);
subplot(2,2,2);
for i=1:100
    p=xcorr(s,i,'unbiased');  %ss
    P=p((i+1):(2*i+1));
    Psav=P(1);
    
    [R1,R]=corrmtx(x,i,'autocorrelation');
    R=inv(R);
    J=Psav-(P.')*R*P;
    plot(i,J,'*');
    hold on;
end
title('代价函数和滤波器的长度之间的关系(信噪比固定)');
xlabel('i'),ylabel('J');