hht时频谱 matlab 乱序_Matlab实现HHT程序(源码,非常珍贵)

clear all;

x=load ('06514135360001170106.TXT');

fs=1000000;

N=length(x);

t=0:1/fs:(N-1)/fs;

z=x;

c=emd(z);

%

计算每个

IMF

分量及最后一个剩余分量

residual

与原始信号的相关性

[m,n]=size(c);

for i=1:m;

a=corrcoef(c(i,:),z);

xg(i)=a(1,2);

end

xg;

for i=1:m-1

%--------------------------------------------------------------------

%

计算各

IMF

的方差贡献率

%

定义：方差为平方的均值减去均值的平方

%

均值的平方

%imfp2=mean(c(i,:),2).^2

%

平方的均值

%imf2p=mean(c(i,:).^2,2)

%

各个

IMF

的方差

mse(i)=mean(c(i,:).^2,2)-mean(c(i,:),2).^2;

end;

mmse=sum(mse);

for i=1:m-1

mse(i)=mean(c(i,:).