clear?all;??
y1=load?（‘20141223.TXT‘）;?
hh=0;
??[row1g1]=size（y1）;
for?g=1:g1
????c1=y1（:g）;
????[rowclu]=size（c1）;
???for?h=1:row/1024
???????x=c1（（1024*（h-1）+1）:1024*h:）;
???hh=hh+1;
??fs=100;?
??N=1024;?
??t=0:1/fs:（N-1）/fs;?
??z=x;??
??c=emd（z）;??
%計算每個IMF分量及最后一個剩余分量residual與原始信號的相關(guān)性?
??[mn]=size（c）;?
????for?i=1:m;??
??????a=corrcoef（c（i:）z）;?
??????xg（i）=a（12）;?%xg是與原始信號的相關(guān)性?
????end?
???xg;??
???for?i=1:m-1??
%--------------------------------------------------------------------?
%計算各IMF的方差貢獻(xiàn)率
%定義：方差為平方的均值減去均值的平方
%均值的平方
%?bghhvchg.imfp2=mean（c（i:）2）.^