Matlab实现HHT程序(源码-非常珍贵)(共5页).doc

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);endxg;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-1mse(i)=mean(c(i,:).