bac.m

这是盲信号的代码 都已经通过编译了 做这方面的同仁可以参考一下 我觉得蛮惯用的

function bac(s)
A =[0.1509 0.0784;0.6979 0.8600];
 x=A*s;
 subplot(521);
 plot(s(1,:))
 subplot(522)
 plot(s(2,:))
%********************
 ys=slfilter(x);
%********************
 subplot(523)
 plot(ys(1,:))
 subplot(524)
 subplot(523)
 plot(x(1,:))
 subplot(524)
 plot(x(2,:))
 subplot(525)
 plot(ys(1,:))
 subplot(526)
 plot(ys(2,:))
%*******************
 y=fastica(x);
%*******************
 subplot(527)
 plot(y(1,:));
 subplot(528)
 plot(y(2,:));
%**************
 W=jadeR(x);
%**************
 yj=W*x;
 subplot(529)
 plot(yj(1,:));
 subplot(5,2,10)
 plot(yj(2,:));