主动噪声控制FxLMS算法

function [ eWv] = FxLMS1（ xdJKLuHsHs_m ）
% %多通道FXLMS算法
% %输入参数：
% %p噪声源列矢量
% %x参考信号列矢量
% %其中p与x有相同长度
% %J次级源数
% %K误差数
% %L滤波器长度
% %u迭代步长
% %Hp为噪声源到误差初级通路，有K列[Hp1 Hp2 .. HpK]
% %Hs为次级源到误差的真实次级通路有J*K列[Hs11 Hs12 .. Hs1K Hs21 .. HsJK]
% %Hs_m为Hs的估计值
if （~isvector（x））
    error（‘参考信号必须是矢量‘）;
elseif （J<=0 || K<=0 || L<=0 || u<=0）
    error（‘次级源数、误差数、滤波器长度、迭代步长必须大于0‘）;

elseif （length（Hs（1:））~=J*K）
    error（‘初级通路数量不对‘）;
end
num_iter=length（x）;%迭代点数
M=length（Hs（:1））;%次级通路长度
MAX=max（ML）;
W=