amari.m

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

t=1:1000;
s(1,:) = unifrnd(-1,1,1,1000);
s(2,:) = 0.1*sin(400*t).*cos(30*t);
 s(3,:) = 0.01*sign(sin(500*t+9*cos(40*t)));
 A=randn(3);
 x=A*s;
 figure(1)
 subplot(311);
 plot(s(1,:))
subplot(312)
 plot(s(2,:))
 subplot(313)
 plot(s(3,:))
figure(2)
subplot(311);
 plot(x(1,:))
subplot(312)
 plot(x(2,:))
 subplot(313)
 plot(x(3,:))