filtre_lms.m

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

%%%%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
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%