sep_quality.m

来自「这是盲信号的代码 都已经通过编译了 做这方面的同仁可以参考一下 我觉得蛮惯用的」· M 代码 · 共 13 行

function q=qualilty(V)

  [M,N]=size(V);
  q=zeros(1,M);

  for i=1:M,
    
    [maxsq,imax]=max(V(i,:).^2);
    sumsq=sum(V(i,[1:imax-1 imax+1:N]).^2);
    %q(i)=100*sqrt(1-maxsq/sumsq);
    q(i)=sqrt(sumsq/maxsq);

  end