ex_lms1.asv

Very good matlab LMS filter program.

%%%%%%%%%%%%%%%%% [ UNIVERSITE DU QUEBEC A TROIS-RIVIERES  ]  %%%%%%%%%%
%%%%%%%%%%%%%DEPARTEMENT DE GENIE ELECTRIQUE ET DE GENIE INFORMATIQUE%%%%
%%%%%%%%%%%%%TRAITEMENT AVENCE DU SIGNAL ELECTRIQUE[AHMAD YEHYA FAWAZ]%%%
%%%%%EXEMPLE
clear all
close all
clc
%%%%%DONNIES
N=1000;
%%%signal x(n)
% for n=1:N
% x(n)=sin(0.05*pi*n);
% end

for n=1:N
x(n)=(1.1).^(-n);
end

n=1;
y(n,1)=u(n);
n=2;
y(n,1)=1.5*y(n-1)+u(n);
n=3;
y(n,1)=1.5*y(n-1)-0.7*y(n-2)+u(n)-u(n-2);
n=4;
y(n,1)=1.5*y(n-1)-0.7*y(n-2)+u(n)-u(n-2)+0.2*u(n-3);
n=5;
y(n,1)=1.5*y(n-1)-0.7*y(n-2)+u(n)-u(n-2)+0.2*u(n-3);
n=6;
y(n,1)=1.5*y(n-1)-0.7*y(n-2)+x(n-5)+u(n)-u(n-2)+0.2*u(n-3);
for n=7:N
    y(n,1)=1.5*y(n-1)-0.7*y(n-2)+x(n-5)+0.5*x(n-6)+u(n)-u(n-2)+0.2*u(n-3);
end

%%%%FILTRE LMS:%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%PAS D'ADAPTATION%%%%%%%%%%%%%
mu=0.01; 
%%%%%%%%%%POIDS DU FILTRE%%%%%%%%%%%%%%
w=[1.5,-0.7,1,0.5,1,-1,0.2]';
h=w;

%%%bruit
bruit;
%%%%%%%%%%%%APPLICATION DU FILTRE LMS
FILTRE_LMS;
%%%%FUGURE(1)%%%
figure(1);
plot(ErreureEstime);
hold on
plot(x,'-r')
erreure=0;
for j=1:20
    erreure=erreure+((x-(ErreureEstime)).^2);
    BRUIT;
    FILTRE_LMS;
end
erreure=erreure/20;
%%%%FIGURE(2)%%%
figure(2);
plot(s);
hold on
plot(x,'-r')
figure(3);  semilogy((mean(erreure,1)'),'r'); grid on;
            title('Convergence du RLS');
            ylabel('Erreur moyenne quadratique');
            xlabel('蒫hantillons');
                        %%%%%%%%%%%%%%%%[ DANIEL ]%%%%%%%%%%%