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

end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1387??2014-05-29?11:11??ICA\ICA.m

?????文件????????909??2014-05-29?11:09??ICA\TEST_ICA.m

?????目錄??????????0??2014-05-29?11:20??ICA

-----------?---------??----------?-----??----

?????????????????2296????????????????????3