sep_quality.m

该代码为盲源分离中的相对牛顿法

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