io.m

来自「一种新的时频分析方法对应的程序送给搞EMD或者HHT但没有下载到完整程序的朋友。」· M 代码 · 共 22 行

function ort = io(x,imf)

% ort = IO(x,imf) computes the index of orthogonality
%
% inputs : - x    : analyzed signal
%          - imf  : empirical mode decomposition

lx = size(imf,2);
n = size(imf,1);

s = 0;

for i = 1:n
  for j =1:n
    if i~=j
      s = s + abs(sum(imf(i,:).*imf(j,:))/sum(x.^2));
    end
  end
end

ort = 0.5*s;