xnzb.m

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

function e=xnzb(y)
%性能指标计算函数
n=length(y);
a=0;
b=0;
c=0;
for i=1:n
   a=sum(abs(y(i,:)))/max(abs(y(i,:)))-1;
   b=sum(abs(y(:,i)))/max(abs(y(:,i)))-1;
   c=c+a+b;
end
e=c/(n*(n-1));