clc
clear
load?signal.mat
%%?輸入數據
ts?=?0:0.001:0.3;

fs?=?1000;
x?=?cos（2*pi*20*ts）?+?2*cos（2*pi*100*ts）;
N?=?length（x）;

%%?EMD和HT
figure（）
emd（x）;
[imfresidualinfo]=emd（x‘Interpolation‘‘pchip‘‘Display‘0）;
figure（）
hht（imffs）;
%?橫軸表示時間、縱軸表示頻率，顏色表示能量
[hs?f?t?imfinsf?imfinse]?=?hht（imffs）;
%?hs——信號的希爾伯特譜（Hilbert?Spectrum?）
%?f——信號的頻率向量（Frequency?vector?of?signal）
%?t——信號的時間向量（Time?vector?of?signal）
%?imfinsf——每個imf的瞬時頻率（instantaneous?frequency?of?each?imf）
%?imfinse——每個imf的瞬時能量（instantaneous?energy?of?each?imf）
[mn]?=?size（hs）;
l?=?size（imf2）;

%?畫出所有imf函數。
figure（）
for?i?=?1:l
????plot（1:nimfinsf（:i））;
????hold?on;
end

%%?由于HHT變換，得到的時頻譜hs是錯誤排列，所以需要重建時頻譜
%?重建思路：HHT變換得到每一條imf的瞬時頻率向量，和對應的瞬時能量向量
%?f?=?（0:m-1）/（m-1）*（fs/2），這個是記錄每一個點的記錄頻率，反向操作，

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件?????508741??2020-03-10?16:39??HHT變換思想與malab實現.docx

?????文件???????1776??2020-03-10?16:38??HHT.m

-----------?---------??----------?-----??----

???????????????510517????????????????????2