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