io.m

来自「新的信号处理方法－希尔伯特黄变换的EMD分解MATLAB程序 把程序看懂 对希尔」· M 代码 · 共 22 行

function ort = io(x,imf)

% 调用形式为：ort = IO(x,imf) ，用来计算正交性
%
% 输入变量有 : - 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;