ratioa.m

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

function r = ratioa(h)
 
% The function RATIOA calculates the local non-orthogonal leakage for pair of adjoining components.
%
% Calling sequence-
% r=ratioa(h)
%
% Input-
%	h	- 2-D matrix h(n,k) of IMF components
%		 excluding the trend component.
% Output-
%	r	- vector of the ratio for adjoining components
%		(h(1,i)*h(1,j)+...+h(n,i)*h(n,j))/
%		(h(1,i)^2+h(1,j)^2)+...+(h(n,i)^2+h(n,j)^2)
%

%----- Get dimensions and initialize
[npt,knb] = size(h);
r1=zeros(knb-1,1);

%----- Calculate the ratio for each pair
for i=1:knb-1,
   j=i+1;
   a1=abs(sum(h(:,i).*h(:,j)));
   b1=sum(h(:,i).*h(:,i))+sum(h(:,j).*h(:,j));
   r1(i)=a1/b1;
end
%r=mean(r1);
%s=std(r1);
%q=max(r1-r);
r=r1;