function  qq=eh4ph(x)
%對EH4高頻進行EMD分解并對各模態分量進行FFT變化
t=0.02/4096:0.02/4096:0.02;
f =semd(x);
imf=f';
tn=length(t);

t1=min(t(1),t(tn));
t2=max(t(1),t(tn));

%%%%%%%%%%%%%%%%%%% Plot the EMD %%%%%%%%%%%%%%%%%%%
figure;
clf;

imfStart = 1;
imfStop = size(imf,2);
imfLast = imfStop+1;

for j=imfStart:imfLast;
    subplot(imfLast-imfStart+1,1,j+1-imfStart);
    if j==imfStart
        k=fftph(x,1);
        yMax = max(k);
        yMin = min(k);
            else
       k=fftph(imf(:,j-1),1);
       yMax = max(k);
        yMin = min(k);
    end
      axis([0 96000 yMin yMax]);   
      s1='imf';
   if j==1
       ylabel('signal')
   else
   s2=[s1 num2str(j-1)];
   ylabel(s2);
   end
      if( j == 1 )
       title( 'Spectrum of Signal and IMF Components' );
   end
   if( j < imfLast );
      set(gca,'xticklabel','');
   end;
end;
xlabel('f/Hz');

figure;
for j=imfStart:imfLast;
    subplot(imfLast-imfStart+1,1,j+1-imfStart);
    if j==imfStart
        z=sum(imf(:,imfStart:imfStop),2);
        plot(t,z);
        yMax = max(z);
        yMin = min(z);
    else
        plot(t,imf(:,j-1));
        yMax = max(imf(:,j-1));
        yMin = min(imf(:,j-1));
    end
    if( abs(yMax-yMin) > 1e-10 )
        axis([t1 t2 yMin yMax]);
    else
        axis([t1 t2 -1 1]);
    end
   s1='imf';
   if j==1
       ylabel('signal')
   else
   s2=[s1 num2str(j-1)];
   ylabel(s2);
   end
   if( j == 1 )
       title( 'Time Series and IMF Components' );
   end
   if( j < imfLast );
      set(gca,'xticklabel','');
   end;
end;
xlabel('Time');