盲源分离/ICA经典算法——FastICA

这个是由芬兰人提出的FastICA算法，我编写了一个主函数，名字叫zhuhanshu.m,打开后，可以调用其他程序，查看盲源分离的结果

clear all;
clc;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%注意：本程序使用的测试信号是由功能函数生成的，并且是4路%%%%%%%%%%%%%%%%

[sigmixedsig]=demosig（）                     %调用功能函数demosig，生成了试验用的信号，得到源信号sig和混合信号mixedsig维数都为4*500

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%画出各路源信号的波形图%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
sig1=sig（1:）;   
sig2=sig（2:）;
sig3=sig（3:）;
sig4=sig（4:）;
figure（‘Name‘‘（1）Source Signal‘）
icaplot（‘classic‘sig500500‘lll‘）
subplot（411）plot（sig1）ylabel（‘正弦曲线‘）;
subplot（412）plot（sig2）ylabel（‘曲线‘）;
subplot（413）plot（sig3）ylabel（‘锯齿状‘）;
subplot（414）plot（sig4）ylabel（‘脉冲干扰‘）;
s=[sig1;sig2;sig3;sig4];

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%画出混合信号的波形图%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

mixedsig1=mixedsig（1:）;
mixedsig2=mixedsig（2:）;
mixedsig3=mixedsig（3:）;
mixedsig4=mixedsig（4:）;
figure（‘Name‘‘（2）Mixed Signal‘）
subplot（411）
plot（mixedsig1）
subplot（412）
plot（mixedsig2）
subplot（413）
plot（mixedsig3）
subplot（414）
plot（mixedsig4）
mixedsig=[mixedsig1;mixedsig2;mixedsig3;mixedsig4];

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%主成分分析%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

[ED]=pcamat（mixedsig）                          %调用功能函数pcamat得到了特征值和特征向量
warning（‘特征值并没有按照降序排列‘）;
eigenvalues=flipud（sort（diag（D）））               %降序排列特征值

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%计算白化矩阵%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

[nvwmdwm]=whitenv（mixedsigED）               %调用功能函数whitenv得到白化后的信号，白化矩阵，白化矩阵的逆矩阵

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%不动点迭代算法%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

[AW]=fpica（mixedsigwmdwm）;                   %调用功能函数fpica，用定点迭代算法的出分离矩阵W
[icasigAW]=fastica（nv）                        %调用功能函数fastica得到了对源信号的估计信号icasig
is1=icasig（1:）;
is2=icasig（2:）;
is3=icasig（3:）;
is4=icasig（4:）;
figure（‘Name‘‘（3）ICAsig‘）                      %画出用FastICA算法得出的信号，结果显示算法性能良好
subplot（411）plot（is1）
subplot（412）plot（is2）
subplot（413）plot（is3）
subplot（414）plot（is4）
xlabel（‘采样点数‘）

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

     文件       2358  2013-10-28 20:29  FastICA_25\aaa.m

     文件       1307  2005-10-19 16:05  FastICA_25\Contents.m

     文件        794  2005-10-19 16:09  FastICA_25\CVS\Entries

     文件          8  2005-10-19 16:09  FastICA_25\CVS\Repository

     文件         21  2005-10-19 16:09  FastICA_25\CVS\Root

     文件        746  2013-10-28 18:55  FastICA_25\demosig.m

     文件        402  2003-04-05 17:23  FastICA_25\dispsig.m

     文件      18627  2013-10-20 20:20  FastICA_25\fastica.asv

     文件      18436  2005-10-19 16:05  FastICA_25\fastica.m

     文件      19214  2005-10-19 16:05  FastICA_25\fasticag.m

     文件      26041  2005-06-16 15:52  FastICA_25\fpica.m

     文件      13126  2004-07-27 16:09  FastICA_25\gui_adv.m

     文件       7411  2003-09-08 14:28  FastICA_25\gui_advc.m

     文件      19416  2003-09-10 13:33  FastICA_25\gui_cb.m

     文件       2704  2003-04-05 17:23  FastICA_25\gui_cg.m

     文件      14536  2005-10-19 16:05  FastICA_25\gui_help.m

     文件       5129  2004-07-27 16:09  FastICA_25\gui_l.m

     文件       3665  2003-09-11 15:01  FastICA_25\gui_lc.m

     文件       5017  2004-07-27 16:09  FastICA_25\gui_s.m

     文件       2402  2003-09-08 14:29  FastICA_25\gui_sc.m

     文件      13259  2003-04-05 17:23  FastICA_25\icaplot.m

     文件      12075  2003-12-15 20:24  FastICA_25\pcamat.m

     文件        461  2003-04-05 17:23  FastICA_25\remmean.m

     文件     160192  2013-09-30 18:18  FastICA_25\Speech4.mat

     文件       2842  2003-10-12 12:04  FastICA_25\whitenv.m

     文件       2320  2013-10-28 20:26  FastICA_25\zhuhanshu.m

     文件       2320  2013-10-28 20:26  FastICA_25\有关FastICA算法\zhuhanshu.m

     文件      18432  2013-10-30 18:47  FastICA_25\有关FastICA算法\关于FastICA大程序的解释.doc

     文件      39424  2013-10-28 20:40  FastICA_25\有关FastICA算法\编写主函数后运用FastICA的运行结果.doc

     目录          0  2005-10-19 16:09  FastICA_25\CVS

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