filtre_lms.asv

来自「Very good matlab LMS filter program.」· ASV 代码 · 共 31 行

%%%%%%%%%%%%%%%%% [ UNIVERSITE DU QUEBEC A TROIS-RIVIERES  ]  %%%%%%%%%%
%%%%%%%%%%%%%DEPARTEMENT DE GENIE ELECTRIQUE ET DE GENIE INFORMATIQUE%%%%
%%%%%%%%%%%%%TRAITEMENT AVENCE DU SIGNAL ELECTRIQUE[AHMAD YEHYA FAWAZ]%%%
%%%%FILTRE LMS:
M=5;                      %%%%%%%%%%DIMENTION DU FILTRE%%%%%%%%%%
Y=0.2;                    %%%%%%%%%%PAS D'ADAPTATION%%%%%%%%%%%%%
W_LMS=zeros(1,M);         %%%%%%%%%%POIDS DU FILTRE%%%%%%%%%%%%%%

delai=0;                  %%%%%%%%%%(M+1)/2:%%%%%%%%%%%%%%%%%%%%%
for n=1:N
    %%%%%%%%%%%%%%%%%%%%%%%%%%%SORTIE DU FILTRE LMS%%%%%%%%%%%%%%
    V_estime(1,n)=0;
    for m=1:M
        if (n-m+1>0)
            V_estime(1,n)= V_estime(1,n)+W_LMS(1,m)*r(n-m+1,1);
        end
    end
    %%%%%%%%%%%%%%%%%%%%%%%%%%%ERREURE ESTIMI%%%%%%%%%%%%%%%%%%%%
    if(n>delai)
        ErreureEstime(1,n)=s(1,n-delai)-V_estime(1,n);
    else
        ErreureEstime(1,n)=0;
    end
    %%%%%%%%%%%%%%%%%%%%%%%%%%%ADAPTATION DE POIDS%%%%%%%%%%%%%%%
    for m=1:M
        if(n-m+1>0)
            W_LMS(1,m)=W_LMS(1,m)+Y*ErreureEstime(1,n)*r(n-m+1,1);
        end
    end 
end
    %%%%%%%%%%%%%%%%%%%%%%%%%% [DANIEL] %%%%%%%%%%%%%%%%%%%%%%%%%