titled62.m

LMS

clear;
close all;
clc;

N=1024;
a1=1.558;
a2=-0.81;

v=randn(1,N);
eb=zeros(3,N);
ef=zeros(3,N);
kb=zeros(3,N);
kf=zeros(3,N);
d=zeros(3,N);
r=zeros(3,N);
w1=zeros(1,N);
w2=zeros(1,N);

u=0.002;
w3=zeros(1,N+1);
w4=zeros(1,N+1);
e=zeros(1,N);
y=zeros(1,N);

x(1)=v(1);
x(2)=a1*x(1)+v(2);
for n=3:N
    x(n)=a1*x(n-1)+a2*x(n-2)+v(n);
    eb(1,n)=x(n);
    ef(1,n)=x(n);
    kb(1,n)=kf(1,n-1)+x(n)^2;
    kf(1,n)=kf(1,n-1)+x(n)^2;
end;
for n=1:N
    r(1,n)=1;
end

    wb=zeros(3,N+1);
    wf=zeros(3,N+1);
    
for m=1:2
    for n=5:N
        d(m+1,n)=d(m+1,n-1)+eb(m,n-1)*ef(m,n)/r(m,n-1);
        ef(m+1,n)=ef(m,n)-d(m+1,n)*eb(m,n-1)/kb(m,n-1);
        eb(m+1,n)=eb(m,n-1)-d(m+1,n)*ef(m,n)/kf(m,n);
        wb(m+1,n)=d(m+1,n)/kb(m,n-1);
        wf(m+1,n)=d(m+1,n)/kf(m,n);
        kf(m+1,n)=kf(m,n)-(d(m+1,n))^2/kb(m,n-1);
        kb(m+1,n)=kb(m,n-1)-(d(m+1,n))^2/kf(m,n);
        r(m+1,n-1)=r(m,n-1)-(eb(m,n-1))^2/kb(m,n-1);
    end
end
for n=2:N
    w1(1,n)=wb(2,n)-wf(2,n)*wb(3,n);
    w2(1,n)=wb(3,n);
end

%lms
w3(1)=4;
w4(1)=-10;
for n=3:N;
    y(n)=w3(n)*x(n-1)+w4(n)*x(n-2);
    e(n)=x(n)-y(n);
    w3(n+1)=w3(n)+2*u*e(n)*x(n-1);
    w4(n+1)=w4(n)+2*u*e(n)*x(n-2);
end;

%RLS
W=zeros(2,N+1);
C=1000*eye(2);
r=0.998;

for n=3:N
    er=x(n)-x(n-1:-1:n-2)*W(:,n);
    ur=x(n-1:-1:n-2)*C*x(n-1:-1:n-2)';
    g=C*x(n-1:-1:n-2)'/(r+ur);
    W(:,n+1)=W(:,n)+g*er;
    C=(C-g*x(n-1:-1:n-2)*C)/r;
end

%GAL
e0f=zeros(1,N);
e0b=zeros(1,N);
e1f=zeros(1,N);
e1b=zeros(1,N);
e2f=zeros(1,N);
e2b=zeros(1,N);
k1=zeros(1,N+1);
k2=zeros(1,N+1);
w5=zeros(1,N);
w6=zeros(1,N);
for n=1:N
    e0f(n)=x(n);
    e0b(n)=x(n);
end
for n=2:N
    e1f(n)=e0f(n)-k1(n)*e0b(n-1);
    e1b(n)=e0b(n-1)-k1(n)*e0f(n);
    k1(n+1)=k1(n)+2*u*[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);
    k2(n+1)=k2(n)+2*u*[e2f(n)*e1b(n-1)+e2b(n)*e1f(n)];
    
    w5(n)=k1(n)*[1-k2(n)];
    w6(n)=k2(n);
end


figure(1);
plot(w1,'y');
hold on;
plot(w2,'y');
hold on;
plot(w3,'r');
hold on;
plot(w4,'r');
hold on;
plot(w5,'g');
hold on;
plot(w6,'g');
hold on;
plot(W(1,:),'b');
hold on;
plot(W(2,:),'b');
axis([0,1200,-1,2]);
hold off;