contrastratio.m

来自「完整的独立分量分析(ICA)程序包,可以用来进行盲源分离」· M 代码 · 共 21 行

function [valueCR,valueENL] = ContrastRatio(I1)

% EQIM --- Return the ContrastRatio of input vector X
% Date: April,20,2004
%***********************X： is a vector 　
%***********************Y：is ContrastRatio


for i=1:100
    for j=1:128
        temp(i,j)=I1(i+64,j+126);
    end
end
vecorX=reshape(temp',1,128*100);


temp1=std(vecorX);
temp2=mean(vecorX);
valueCR=temp1/temp2;
valueENL=(1/valueCR)^2;