text1.m

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

function G=text()
%author lucky zhang
%2005.6.29
t=1:3000;
s(1,:)=sqrt(2)*sin(2*pi*t/97);
s(2,:)=0.2*square(2*pi*t/233);
s(3,:)=0.2*sin(2*pi*t/103).*cos(2*pi*t/277);
s(4,:)=sign(sin((2*pi*t/40)+9*cos(2*pi*t/600)));
A=rand(4);
x=A*s;
%W=jadeR(x);
%y=W*x;
[y,W]=slfilter(x);
%G=W'*A;
G=xsxs(s,y);
figure(1);
mplot(s);
figure(2);
mplot(y);