oldsubstractsource.m

matlab的源程序

function [ContributionEst,Filter] = SubstractSource(Observations,Extracted,LFiltreSoustracteurAC,LFiltreSoustracteurC)
% SubstractSource identifie le filtre inverse permettant
% l'annulant des contributions de la source extraite
% des observations.

% On cherche un filtre de longueur 2L+1
% On utilise les moindres carr閑:
% min ||AX-B||^2 : X=(A'A)^-1A'B

% A contient le signal extrait
% B contient les observations

NObs = size(Observations,1);
LObs = size(Observations,2);

%LObs = min(LObs,2000);

u0 = 1;
% Longueur filtre anti-causal
LMax1 = LFiltreSoustracteurAC;
% Longueur filtre causal
LMax2 = LFiltreSoustracteurC;

% G閚閞ation de la matrice A:
for(ij=0:LObs-1) % Indice des lignes
    for (ii=LMax1:-1:-LMax2) % Indice de colonnes
        if (ij+u0+ii<1)
            A(ij+1,-ii+LMax1+1) = 0;
        elseif (ij+u0+ii>length(Extracted))
            A(ij+1,-ii+LMax1+1) = 0;
        else
            A(ij+1,-ii+LMax1+1) = Extracted(ij+u0+ii);
        end
    end
end

% La matrice B contient les observations
for (ii=1:NObs) % Indice de colonnes
    for(ij=0:LObs-1) % Indice des lignes
        B(ij+1,ii) = Observations(ii,u0+ij);
    end
end

Filter = inv(A'*A)*A'*B;

% On met le signal extrait dans une matrice pour le filtrage:
for(ij=1:length(Extracted)) % Indice des lignes
    for (ii=LMax1:-1:-LMax2) % Indice de colonnes
        u0 = ij;
        if (u0+ii<1)
            Ext(ij,-ii+LMax1+1) = 0;       
        elseif (u0+ii>length(Extracted))
            Ext(ij,-ii+LMax1+1) = 0;
        else
            Ext(ij,-ii+LMax1+1) = Extracted(u0+ii);
        end
    end
end

ContributionEst = Ext*Filter;