hhtp.m

基于hht的大地电磁信号资料处理MATLAB方法实现

function [omega,amp,mag,phase,w]=hhtp(x,t)
%得到信号X的三维HHT谱图
  f=emd(x);
  imf=f';
  starttime=min(t);
  endtime=max(t);
 [omega,amp,m,phase,w]=fhilbert(imf,starttime,endtime,[]);
 mag=full(m);
 mesh(t,w,mag);
 xlabel('Time/s');ylabel('Frequency/Hz');zlabel('Power');
end