function?Z=ICA（X）
%-----------去均值---------
[MT]?=?size（X）;?%獲取輸入矩陣的行/列數，行數為觀測數據的數目，列數為采樣點數??????
average=?mean（X‘）‘;??%均值
for?i=1:M
????X（i:）=X（i:）-average（i）*ones（1T）;?
end
%---------白化/球化------
Cx?=?cov（X‘1）;????%計算協方差矩陣Cx
[eigvectoreigvalue]?=?eig（Cx）;?%計算Cx的特征值和特征向量
W=eigvalue^（-1/2）*eigvector‘;???%白化矩陣
Z=W*X;???%正交矩陣
?
%----------迭代-------
Maxcount=10000;????????%最大迭代次數
Critical=0.00001;???%判斷是否收斂
m=M;????????????????%需要估計的分量的個數
W=rand（m）;
for?n=1:m??
????WP=W（:n）;??%初始權矢量（任意）
%?????Y=WP‘*Z;
%?????G=Y.^3;%G為非線性函數，可取y^3等
%?????GG=3*Y.^2;??%G的導數
????count=0;
????LastWP=zeros（m1）;
????W（:n）=W（:n）/norm（W（:n））;
????while?abs（WP-LastWP）&abs（WP+LastWP）>Critical
????????count=count+1;???%迭代次數
????????LastWP=WP;??????%上次迭代的值
???????%?WP=1/T*Z*（（LastWP‘*Z）.^3）‘-3*LastWP;
????????for?i=1:m?????
WP（i）=mean（Z（i:）.*（tanh（（LastWP）‘*Z）））-（mean（1-（tanh（（LastWP））‘*Z）.^2））.*LastWP（i）;
????????end
????????WPP=zeros（m1）;
????????for?j=1:n-1
????????????WPP=WPP+（WP‘*W（:j））*W（