ICA特征提取

ICA特征提取matlab代码，同时结合支持向量机（SUV）对ECG信号进行特征处理

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