clear all;
clc;
N=1000;
u=0.01;
v=randn(1,N);
v1=sqrt(0.01)*randn(1,N);
a1=1;
a2=2;
w1=[];
w2=[];
d=[];
y=[];
e=[];
w1(1)=0;
w2(1)=0;
w1(2)=0;
w2(2)=0;
y(1)=0;
d(1)=a1*v(1)+v1(1);
e(1)=d(1)-y(1);
y(2)=0;
d(2)=a1*v(2)+a2*v(1)+v1(2);
e(2)=d(2)-y(2);
w1(3)=2*u*e(2)*v(2);
w2(3)=2*u*e(2)*v(1);
for n=3:N
    d(n)=a1*v(n)+a2*v(n-1)+v1(n);
    y(n)=w1(n)*v(n)+w2(n)*v(n-1);
    e(n)=d(n)-y(n);
    w1(n+1)=w1(n)+2*u*e(n)*v(n);
    w2(n+1)=w2(n)+2*u*e(n)*v(n-1);
end
figure(1);
subplot(2,1,1);
plot(w1);
title('權值w1(n)隨迭代次數n的變化曲線');
subplot(2,1,2);
plot(w2);
title('權值w2(n)隨迭代次數n的變化曲線');
figure(2);
plot(e.*e); 
title('迭代輸出誤差平方性能曲線e^2(n)');