clc;
clear;
u=0.0001;        %  0<u<<1
var=1;
a1=-0.195;
a2=0.95;
N=5000;

v=0+sqrt(var)*randn(1,N);   %   零均值、方差為var的白噪聲

w1(3)=0;        %   設定各參數初值
w2(3)=0;
x(1)=v(1);
x(2)=-a1*x(1)+v(2);

for n=3:N
    x(n)=v(n)-a1*x(n-1)-a2*x(n-2);
    y(n)=w1(n)*x(n-1)+w2(n)*x(n-2);
    d(n)=x(n);
    e(n)=d(n)-y(n);
    w1(n+1)=w1(n)+2*u*x(n-1)*e(n);
    w2(n+1)=w2(n)+2*u*x(n-2)*e(n);
end

figure(1);
    plot(w1,'g');
    hold on
    plot(w2,'r');
    legend('w1','w2')
    xlabel('n');
    ylabel('w(n)');
    title('w(n)變化的曲線    u=0.0001, σ^2=1, a1=-0.195, a2=0.95');
    hold off

figure(2);
    plot(e.*e);
    xlabel('n');
    ylabel('e^2(n)');
    title('e^2(n)變化的曲線');