hht变换中emd批量处理数据

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与原始信号的相关性 
  [mn]=size（c）; 
    for i=1:m;  
      a=corrcoef（c（i:）z）; 
      xg（i）=a（12）; %xg是与原始信号的相关性 
    end 
   xg;  
   for i=1:m-1  
%-------------------------------------------------------------------- 
%计算各IMF的方差贡献率
%定义：方差为平方的均值减去均值的平方
%均值的平方
% bghhvchg.imfp2=mean（c（i:）2）.^