主动噪声控制FxAP算法

function [ eW ] = FxAP（xdJKLNugHsHs_m）
%多通道FXAP算法
%输入参数：
%p噪声源列矢量
%x参考信号列矢量
%其中p与x有相同长度
%J次级源数
%K误差数
%L滤波器长度
%N是投影阶数
%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 || N<=0 || u<=0）
    error（‘次级源数、误差数、滤波器长度、投影阶数、迭代步长必须大于0‘）;

elseif （length（Hs（1:））~=J*K）
    error（‘初级通路数量不对‘）;
end
num_iter=length（x）;%迭代点数
x=reshape（xnum_iter1）;

%控制器变量定义
M=length（Hs（:1））;%次级通路长度
MAX=ma