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