xsxs.m

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

function rs=xsxs(s,y)
%
%计算相似系数程序
%copyright 2005.5.26
%author lucky zhang
%
[dims,samples]=size(s);
if(dims<samples)
    num_sources=dims;
    s=s';
end
[dimy,sampley]=size(y);
if(dimy<sampley)
    y=y';
end
[dims,samples]=size(s);
[dimy,sampley]=size(y);
if ~((dims==dimy)&&(samples==sampley))
    error('length of s and y must be equal\n');
end

a=[s y]; 
c=corrcoef(a);  
rs=c(1:num_sources,num_sources+1:num_sources*2);
%fprintf('源信号和恢复信号的相似系数 ...\n'); 
rs=abs(rs);