ratio1.m

这是一个关于hht变换很有用的工具箱

function r = ratio1(h)
  
% The function RATIO1 calculates the non-orthogonal leakage of components.
%
% Calling sequence-
% r=ratio1(h)
%
% Input-
%	h	- 2-D matrix h(n,k) of IMF components
%		 excluding the trend component.
% Output-
%	r	- value of the ratio:
%		(|h1*h2|+...+|hk-1*hk|)/(h1+...hk)^2.
%

%----- Get dimensions and initialize
[npt,knb] = size(h);
a1=0;

%----- Calculate the ratio
for i=1:knb-1,
   for j=i+1:knb,
      a1=a1+sum(h(:,i).*h(:,j));
   end
end
h=sum(rot90(h));
h=h.*h;
a2=sum(h);
r=a1/a2;