function [m,p] = hierxiangwei(imf,N); L = size(imf,1); % Number of components in the decomposition S = []; % Matrix which will contain the time-frequency-amplitude representation clear x z m p freq x = imf'; % now each column is a component z = hilbert(x); % analytic signal m = abs(z); % module of z p = angle(z); % phase of z clear all pp=load( 'g:\matlab程序\2.txt'); n=zeros(6); m=zeros(6); for j=1:6 a=7-j; m=zeros(6); for c=1:a x0=pp(1:10,j); [z0,a0]=hierxiangwei(x0,10) x1=pp(1:10,j+c); [z1,a1]=hierxiangwei(x1,10) r0=abs(1/10*sum(exp(i*(a0-a1)))) m(c,1)=r0; end n(:,j)=m(:,1); end A=n; n=6; for b=1:n for c=1:n if A(b,c)>=0.9 d(b,c)=1; else d(b,c)=0; end end end