hhtmag.m

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

function [h,w]=hhtmag(x,t)
%得到瞬时频率和其对应的频率的复信号
f=emd(x);
imf=f';
starttime=min(t);
endtime=max(t);
[omega,amp,m,phase,w]=myfhilbert(imf,starttime,endtime,[]);
mag=full(m);

[l,n]=size(mag);
h=full(sum(m,2)); %求得个频率点得自HIBERT复信号，累加各频点得瞬时HIBERT复信号
end