online.m

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

function online()
%2005.6.24
%author lucky zhang
t=1:10000;
t=t/20000;
n=length(t);
s(1,:)=unifrnd(-1,1,1,10000);
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;
w=eye(3);
y=w*x;
for i=1:n
  %   eta=250*exp(-5*t(i));
   eta=0.001;
    fy(:,i)=3/4*y(:,i).^11+25/4*y(:,i).^9-14/3*y(:,i).^7-47/4*y(:,i).^5+29/4*y(:,i).^3;
 %   D=diag(diag(fy(:,i)*y(:,i)'));
    dw=eta*(eye(3)-fy(:,i)*y(:,i)')*w;
    w=w+dw;
    y=w*x;
 end
rs=xsxs(s,y)
figure(1);
mplot(s);
xlabel('原始信号');
figure(2);
mplot(x);
xlabel('混合信号')
figure(3);
mplot(y);
xlabel('分离信号');