six_signal.m

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

function [rs,PI] = six_signal()
t=1:512;
s1=2*sin(t/149).*cos(t/8);
s2=sign(sin(12*t+9*cos(2/29)));
s3=(rem(t,79)-17)/23;
s4=((rem(t,23)-11)/9).^5;
s5=5*exp(-t/121).*cos(37*t);
s6=((rand(1,512)<0.5)*2-1).*log(rand(1,512));
subplot(321);
plot(s1);
subplot(322);
plot(s4);
subplot(323);
plot(s2);
subplot(324);
plot(s5);
subplot(325);
plot(s3);
subplot(326);
plot(s6);
figure(2);
s=[s1;s2;s3;s4;s5;s6];
%A=rand(6);
A =[0.1902    0.8586    0.6782    0.5433    0.2421    0.1167;...
    0.6590    0.7318    0.1786    0.4048    0.0777    0.8609;...
    0.1482    0.0916    0.6996    0.8609    0.4826    0.9478;...
    0.6654    0.3480    0.1100    0.8997    0.7154    0.5094;...
    0.7148    0.3194    0.7347    0.3685    0.2073    0.6728;...
    0.7766    0.8016    0.2381    0.9409    0.6049    0.3079];
x=A*s;
%*****************************************
W =  shibbsR(x);
y=W*x;
%*******************************
y1=y(1,:);
y2=y(2,:);
y3=y(3,:);
y4=y(4,:);
y5=y(5,:);
y6=y(6,:);
subplot(321);
plot(y1);
subplot(322);
plot(y4);
subplot(323);
plot(y2);
subplot(324);
plot(y5);
subplot(325);
plot(y3);
subplot(326);
plot(y6);

rs=xsxs(s,y);
PI = xnzb(rs);