untitled5.m

来自「谱估计 (1)修正协方差法 (2)多重信号分类（MUSIC）算法 (3)E」· M 代码 · 共 39 行

clear all;clc;
M=2;
N=3000;
U=0.999;
a1=1.558;
a2=-.81;
V=randn(1,N);
x(1)=V(1);
x(2)=a1*x(1)+V(2);
v(2)=-2*x(2)*x(1);
u(2)=U*x(1)^2+x(2)^2+x(1)^2;
ef=zeros(N,M+1);
eb=zeros(N,M+1);
u=zeros(N,M+1);
v=zeros(N,M+1);
k1=[];
k2=[];
w1=[];
w2=[];
for i=3:N
    x(i)=a1*x(i-1)+a2*x(i-2)+V(i);
    ef(i,1)=x(i);
    eb(i,1)=x(i);
    for j=1:M
        v(i,j+1)=U*v(i-1,j+1)-2*ef(i,j)*eb(i-1,j);
        u(i,j+1)=U*u(i-1,j+1)+ef(i,j)^2+eb(i-1,j)^2;
        K(j)=v(i,j+1)/u(i,j+1);
        ef(i,j+1)=ef(i,j)+K(j)*eb(i-1,j);
        eb(i,j+1)=eb(i-1,j)+K(j)*ef(i,j);
    end
    k1(i)=K(1);
    k2(i)=K(2);
    w1(i)=-k1(i)-k1(i)*k2(i);
    w2(i)=-k2(i);
end
figure(3);
subplot(2,1,1);plot(w1);
subplot(2,1,2);plot(w2);