%下程序為ICA的調(diào)用函數(shù)，輸入為觀察的信號，輸出為解混后的信號
function?Z=ICA（X）
%-----------去均值---------
[MT]?=?size（X）;?%獲取輸入矩陣的行/列數(shù)，行數(shù)為觀測數(shù)據(jù)的數(shù)目，列數(shù)為采樣點(diǎn)數(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ù)
????????La

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

?????文件???????1497??2014-04-05?10:41??ICA\ICA.m

?????文件???????1094??2014-04-05?10:31??ICA\main.m

?????文件???????3306??2014-04-05?12:06??ICA\recognition.m

?????文件??????30966??2014-03-11?17:57??ICA\TestDatabase\1.bmp

?????文件??????30966??2014-03-11?17:57??ICA\TestDatabase\10.bmp

?????文件??????30966??2014-03-11?17:57??ICA\TestDatabase\2.bmp

?????文件??????30966??2014-03-11?17:57??ICA\TestDatabase\3.bmp

?????文件??????30966??2014-03-11?17:57??ICA\TestDatabase\4.bmp

?????文件??????30966??2014-03-11?17:57??ICA\TestDatabase\5.bmp

?????文件??????30966??2014-03-11?17:57??ICA\TestDatabase\6.bmp

?????文件??????30966??2014-03-11?17:57??ICA\TestDatabase\7.bmp

?????文件??????30966??2014-03-11?17:57??ICA\TestDatabase\8.bmp

?????文件??????30966??2014-03-11?17:57??ICA\TestDatabase\9.bmp

?????文件??????30966??2014-03-11?17:57??ICA\TrainDatabase\1.bmp

?????文件??????30966??2014-03-11?17:57??ICA\TrainDatabase\10.bmp

?????文件??????30966??2014-03-11?17:57??ICA\TrainDatabase\11.bmp

?????文件??????30966??2014-03-11?17:57??ICA\TrainDatabase\12.bmp

?????文件??????30966??2014-03-11?17:57??ICA\TrainDatabase\13.bmp

?????文件??????30966??2014-03-11?17:57??ICA\TrainDatabase\14.bmp

?????文件??????30966??2014-03-11?17:57??ICA\TrainDatabase\15.bmp

?????文件??????30966??2014-03-11?17:57??ICA\TrainDatabase\16.bmp

?????文件??????30966??2014-03-11?17:57??ICA\TrainDatabase\17.bmp

?????文件??????30966??2014-03-11?17:57??ICA\TrainDatabase\18.bmp

?????文件??????30966??2014-03-11?17:57??ICA\TrainDatabase\19.bmp

?????文件??????30966??2014-03-11?17:57??ICA\TrainDatabase\2.bmp

?????文件??????30966??2014-03-11?17:57??ICA\TrainDatabase\20.bmp

?????文件??????30966??2014-03-11?17:57??ICA\TrainDatabase\3.bmp

?????文件??????30966??2014-03-11?17:57??ICA\TrainDatabase\4.bmp

?????文件??????30966??2014-03-11?17:57??ICA\TrainDatabase\5.bmp

?????文件??????30966??2014-03-11?17:57??ICA\TrainDatabase\6.bmp

............此處省略10個文件信息