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