eg__1.m

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

%  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

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