gexing.m

来自「自适应中关于格型算法和LMS的仿真程序」· M 代码 · 共 42 行

clc
clear all
a1=1.558;
a2=-0.81;
N=1000;
v=randn(1,N);
x(1)=0;
x(2)=0;
for n=3:N
    x(n)=a1*x(n-1)+a2*x(n-2)+v(n);
end
a=0.005;e0f=x;
e0b=x;
k1(1)=0;
k2(1)=0;
k1(2)=0;
k2(2)=0;
d1(1)=1;
d2(1)=1;
for n=2:N     
e1f(n)=e0f(n)+k1(n)*e0b(n-1);
e1b(n)=e0b(n-1)+k1(n)*e0f(n);
d1(n)=(1-a)*d1(n-1)+e1f(n)^2+e1b(n)^2;
beta1(n)=1/d1(n);
k1(n+1)=k1(n)-beta1(n)*[e1f(n)*e0b(n-1)+e1b(n)*e0f(n)];
e2f(n)=e1f(n)+k2(n)*e1b(n-1);
e2b(n)=e1b(n-1)+k2(n)*e1f(n);
d2(n)=(1-a)*d2(n-1)+e1f(n)^2+e1b(n)^2;
beta2(n)=1/d2(n);
k2(n+1)=k2(n)-beta2(n)*[e2f(n)*e1b(n-1)+e2b(n)*e1f(n)];
end
for n=1:N
a1(n)=-k1(n)*(1+k2(n));
a2(n)=-k2(n);
end
subplot(2,1,1);
plot(a1);
title('a1(n)-n');
subplot(2,1,2);
plot(a2);
title('a2(n)-n');