hht时频图三维_Matlab怎样画HHT三维时频图？

关于三维图我看有人用nspab来画的，具体程序：

function [ output_args ] = Untitled1( input_args )

%UNTITLED1 Summary of this function goes here

%  Detailed explanation goes here

%nspab:

function h1=nspab(data,nyy,minw,maxw,dt)%

% The function NSPAB generates a smoothed HHT spectrum of data(n,k)

% % in time-frequency space, where

% n specifies the length of time series, and

% k is the number of IMF components.

% The frequency-axis range is prefixed.

% Negative frequency sign is reversed.

%

% MATLAB Library function HILBERT is used to calculate the Hilbert transform.

%

% Example, [h,xs,w] = nspab(lod78_p',200,0,0.12,1,3224).

%

% Functions CONTOUR or IMG can be used to view the spectrum,

%    for example contour(xs,w,h) or img(xs,w,h).

%

% Calling sequence-

% [