bss_main.m

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

d1=(s1-y1);                                               %计算s与y的绝对误差
d2=(s2-y2);
r1=(d1./y1);                                              %计算相对误差
r2=(d2./y2);
subplot(5,2,1),plot(s1)                                   %画图表示s1,s2
subplot(5,2,2),plot(s2)
subplot(5,2,3),plot(x1)                                   %画图表示x1,x2
subplot(5,2,4),plot(x2)
subplot(5,2,5),plot(y1)                                   %画图表示y1,y2
subplot(5,2,6),plot(y2)
subplot(5,2,7),plot(r1)                                   %画图表示相对误差
subplot(5,2,8),plot(r2)
C=abs(Y*S')/sqrt((1000*mean(Y'.^2))*(1000*mean(S'.^2))')  %计算估计信号Y与源信号S的相似系数矩阵
PI=0.5*(abs(G(1,1))/max(abs(G(1,:)))+abs(G(2,1))/max(abs(G(2,:)))+abs(G(1,2))/max(abs(G(1,:)))+abs(G(2,2))/max(abs(G(2,:)))+abs(G(1,1))/max(abs(G(:,1)))+abs(G(1,2))/max(abs(G(:,2)))+abs(G(2,1))/max(abs(G(:,1)))+abs(G(2,2))/max(abs(G(:,2)))-4)