ica工具箱matlab

很好用的ica工具箱 matlab 程序

function [W Wefica Wwasobi ISRwa ISRef metoda]= combi（XAR_orderortho）
%
% combines EFICA WASOBI
%

if nargin<3
    ortho=false;
end

%COMMON PREPROCESSING
X = X-mean（X2）*ones（1size（X2））;
C = cov（X‘）;
CC = C^（-1/2）;
x = CC*X;


[Wefica ISRef]=efica（x）;
[Wwasobi WsobiISRwa]= ewasobi（x100.9）;
ISR1=sum（ISRef2）;
ISR2=sum（ISRwa2）;
imin=min（[ISR1 ISR2]）;
if imin（1）    hotove=find（ISR1    ostatni=find（ISR1>=imin（2））;
    W=Wefica（hotove:）;
    if ortho
        Wost=null（Wefica（hotove:））‘;
    else
        Wost=Wefica（ostatni:）;
    end
    y=Wost*x;
    metoda=ones（size（hotove‘））;
else
    hotove=find（ISR2    ostatni=find（ISR2>=imin（1））;
    W=Wwasobi（hotove:）;
    if ortho
        Wost=null（Wwasobi（hotove:））‘;
    else
        Wost=Wwasobi（ostatni:）;
    end
    y=Wost*x;
    metoda=2*ones（size（hotove‘））;
end

while length（ostatni）>1
    [Wefi ISRef0]=efica（y）;
    [WwaAOL_initISRwa0]= ewasobi（y100.9）;
    ISR1=sum（ISRef02）;
    ISR2=sum（ISRwa02）;
    imin=min（[ISR1 ISR2]）;
    if imin（1）        hotove=find（ISR1        ostatni=find（ISR1>=imin（2））;
        W=[W; Wefi（hotove:）*Wost];
        if ortho
            Wost=null（Wefi（hotove:））‘*Wost;
        else
            Wost=Wefi（ostatni:）*Wost;
        end
        y=Wost*x;
        metoda=[metoda ones（size（hotove‘））];
    else
        hotove=find（ISR2        ostatni=find（ISR2>=imin（1））;
        W=[W; Wwa（hotove:）*Wost];
        if ortho
            Wost=null（Wwa（hotove:））‘*Wost;
        else
            Wost=Wwa（ostatni:）*Wost;
        end
        y=Wost*x;
        metoda=[metoda 2*ones（size（hotove‘））];
    end
end
if ~isempty（ostatni）
    W=[W; Wost];
    metoda=[metoda 0];
end

W=W*CC;
Wefica=Wefica*CC;
Wwasobi=Wwasobi*CC;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       7166  2007-03-29 02:55  adv_opt_matrixH.p

     文件      37194  2007-03-29 03:05  adv_param.p

     文件       2014  2007-03-29 02:55  als_reg.p

     文件       3344  2007-03-29 02:55  amuse.p

     文件       1327  2007-03-29 02:55  amusepermuta.p

     文件       5260  2007-03-29 02:55  choseleng.p

     文件      19039  2007-03-29 02:55  chosemat.p

     文件       4995  2007-03-29 02:55  chosevar.p

     文件       1792  2007-02-08 17:49  combi.m

     文件       1531  2007-03-29 02:55  cost_off.p

     文件      27350  2007-03-29 02:55  display_chn.p

     文件      16585  2007-02-07 16:34  efica.m

     文件       5864  2004-02-23 13:27  erica.m

     文件       4174  2007-03-29 02:55  evd.p

     文件       8056  2007-03-29 02:55  evd24.p

     文件      16051  2007-02-13 18:31  ewasobi.m

     文件       5670  2007-03-29 02:55  f_x.p

     文件       6365  2007-03-29 02:55  fajdc4.p

     文件       2610  2007-03-29 02:55  fajdnew.p

     文件      12571  2007-03-29 11:25  fastcICA.p

     文件      11292  2007-03-29 02:55  fasticater.p

     文件       4013  2007-03-29 02:55  ffdiag.p

     文件      18080  2007-03-29 02:55  fjade.p

     文件       9114  2007-03-29 02:55  flexica.p

     文件       3655  2007-03-29 02:55  fobi.p

     文件       3340  2007-03-29 02:55  fobip.p

     文件       2086  2007-03-29 02:55  get_off.p

     文件        963  2007-03-29 02:55  getalgbyname.p

     文件     138301  2007-03-29 02:55  getalgotype.p

     文件       2281  2007-03-29 02:55  getval.p

............此处省略237个文件信息