matlab程序hht

function plot_hht（xTs）
% Plot the HHT.
% plot_hht（xTs）
% 
% :: Syntax
%    The array x is the input signal and Ts is the sampling period.
%    Example on use: [xFs] = wavread（‘Hum.wav‘）;
%                    plot_hht（x（1:6000）1/Fs）;
% Func : emd

% Get HHT.
imf = emd（x）;
for k = 1:length（imf）
   b（k） = sum（imf{k}.*imf{k}）;
   th   = angle（hilbert（imf{k}））;
   d{k} = diff（th）/Ts/（2*pi）;
end
[uv] = sort（-b）;
b     = 1-b/max（b）;

% Set time-frequency plots