xbp.m

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

function [omega,amp,mag,phase,w]=xbp(x,t,n)
%先得到信号X的小波分解系数，再对各系数序列HIBERT变换，得到三维H-xiaobo谱图
wname='db5';
[c,l] = wavedec(x,n,wname);
for i=1:n;
    d(i,:)=wrcoef('d',c,l,wname,i);
    if i==n;
        d(i+1,:)=wrcoef('a',c,l,wname,i);
    end
end
 
  imf=d';
  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