Matlab实现HHT程序源码-非常珍贵

1、Mat lab实现HHT程序（源码非常珍贵）clear all;x=load ();fs=1000000;N=length(x);t=O:l/fs:(N-l)/fs;z=x;c=emd(z);%计算每个IMF分量及最后一个剩余分量residual与原始信号的相尖性m,n=size(c);for i=l:m;a=corrcoef(c(i,:)/z);xg(i)=a(l,2)；endxg;for i=l:m-l%计算各IMF的方差贡献率%定义：方差为平方的均值减去均值的平方%均值的平方%imfp2=mean(c(i,:)/2).A2%平方的均值%imf2p=mean(c(i,:).A2/2)%各个

2、IMF的方差mse(i)=mean(c(i/：).A2/2)-mean(c(i,:)z2).A2;end;mmse=sum(mse);for i=l:m-lmse(i)=mea n(c(i/：).A2/2)-mean(c(i/：)z2).A2;%方差百分比，也就是方差贡献率mseb(i)=mse(i)/mmse*100;%显示各个IMF的方差和贡献率end;%画出每个IMF分量及最后一个剩余分量residual的图形figure(l)for i=l:m-ldisp(imf/int2str(i) ;disp(mse(i) mseb(i);end;subplot(m+IU)plot(tz)set(

3、gca/fontname7times New RomarV)set(gca/f on tsize1,ylabelfsign ar/Amplitude)for i=l:m-lsubplot(m+lzl,i+l);setfgcf/colorVw1)plot(t,c(ijk)set(gca/fontname叮times New RomarV)set(gca/f on tsize1,ylabelimf/int2str(i)endsubplot(m+lzlzm+l);setlgcf/color/w)plot(tjC(m/:)；k/)set(gca/fontname叮times New Roman)set

4、(gca/f on tsize1,ylabel(rint2str(m-l)%画出每个IMF分量及剩余分量residual的幅频曲线figure(2)subplot(m+blj)setfgcf/color/w1)f,z=fft(t,z);plot(f,z； k)set(gca/fontnametimes New Roman)set(gca/f on tsizeylabelinitial signal1,int2str(m-l),Amplitude)for i=l:m-lsubplot(m+lzl,i+l);settgcf/color/w1)f,z=fft(t/c(i/：);plot(fzk)se

5、t(gca/fontname7times New RomarV)set(gca/f on tsize1,ylabel(imf/int2str(i)/Amplitude= 0errorfinf doit etre 01)endM=max(ma x(im);im = Iogl0(im/M+Ie-300);inf=inf/10;imagesc(t,fliplr(l:size(im/l)/(2*size(im/l)/im,inf/0);set(gca/YDir/norma I1)xlabeKCtime)ylabel(normalized frequency”)titlefAilbert-Huang spectrum)function f,z=fftfenxi(t,y)L=length(t);N=2Anextpow2(L);%fft默认计算的信号是从0开始的t=linspace(t(l)/t(L),N);deta=t(2)-t(l);m=0:N-l;f=l./(N*deta)*m;%下面计算的Y就是x(t)的傅里叶变换数值%Y=ex