wena.m

自适应滤波算法

function [out,h]=wena(direct,target,step)

cxx = (xcorr(direct,direct, step));
%ch2路反射波计算互相关矩阵
cxy = (xcorr(direct,target, step)); 
Cxx=[];
for i=1:step+1
    Cxx=[Cxx; cxx(step+2-i:2*step+2-i)];
end
Cxy = (cxy(step+1:-1:1)).';
h = inv(Cxx)*Cxy;
% LMS对消后的信号，ch2路反射波
out = target-filter(conj(h),1,direct);