function?[ht1npm]?=?NMCFLMS（x?ro?lambda?M?L?tol?hr?h_hat）??

%?NMCFLMS:?Normalized?Multi-Channel?Frequency-Domain?LMS?
%??????????Adaptive?Algorithm?for?Blind?Channel?Identification
%
%?------
%?USAGE:?[h?en]?=?NMCFLMS（x?mu?lambda?M?L?tol?hr?hn）
%
%????????x??????:?multi-channel?outputs?（microphone?signals）;
%????????mu?????:?step?size;
%????????lambda?:?a?forgetting?factor;
%????????M??????:?the?number?of?channels;
%????????L??????:?the?number?of?taps?of?the?largest?channel?
%?????????????????impulse?response;
%????????tol????:?convergence?toleration?error;
%????????hr?????:?the?real?channel?impulse?responses;
%????????hn?????:?initial?estimate?of?the?channel?impulse?responses;
%????????h??????:?estimated?channel?impulse?responses;
%????????en?????:?the?norm?of?the?projection?error?vector.
%
%?%
%?------
%?Author:??Mohammad?Ariful?Haque
%??????????DSP?Lab?BUET
%?Last?updated?4/12/2010
%

[rows?cols]?=?size（x）;
hr?=?hr/norm（hr（:））;?
MBlocks?=?fix（rows?-?L）/L;
L2??=?2*L;
ulambda?=?1?-?lambda;
hf1?=?fft（h_hatL2）;
Cx?=?zeros（L2?M）;

for?m?=?1:MBlocks
????hf0?=?hf1;
????Cx??=?[Cx（L+1:L2?:）;?x（（m-1）*L+1:m*L:）];
????Dx??=?fft（Cx）;
????Dxc?=?conj（Dx）;

????Pxm?=?Dx.*Dxc;

????if?（m?==?1）
????????Px?=?Pxm;
????????%Px?=?sum（sum（Pxm））/L2*on