function?Z=ICA（X）
%-----------去均值---------
X=rand（220）;
[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

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

?????文件???????2227??2013-03-19?10:46??ICA\ICA.m

?????目錄??????????0??2013-03-26?15:51??ICA

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

?????????????????2227????????????????????2