function?imf?=?emd（x）
%?Empiricial?Mode?Decomposition?（Hilbert-Huang?Transform）
%?EMD分解或HHT變換
%?返回值為cell類型，依次為一次IMF、二次IMF、...、最后殘差
x???=?transpose（x（:））;%TRANSPOSE函數可返回轉置單元格區域，即將行單元格區域轉置成列單元格區域
imf?=?[];

while?~ismonotonic（x）?%當x不是單調函數，分解終止條件
????x1?=?x;
sd?=?Inf;%均值
%直到x1滿足IMF條件，得c1
????while?（sd?>?0.1）?||?~isimf（x1）?%當標準偏差系數sd大于0.1或x1不是固有模態函數時，分量終止條件
????????s1?=?getspline（x1）;?????????%?極大值點樣條曲線
????????s2?=?-getspline（-x1）;???????%?極小值點樣條曲線
????????x2?=?x1-（s1+s2）/2;
???????
????????sd?=?sum（（x1-x2）.^2）/sum（x1.^2）;
????????x1?=?x2;
????end
???
????imf{end+1}?=?x1;
????x??????????=?x-x1;
end
imf{end+1}?=?x;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2018-04-01?17:48??EMD\
?????文件?????????763??2018-03-31?21:06??EMD\emd.m
?????文件?????????206??2018-03-31?21:08??EMD\FFTAnalysis.m
?????文件?????????992??2018-03-31?22:02??EMD\findpeaks.m
?????文件?????????125??2018-03-31?21:07??EMD\getspline.m
?????文件?????????235??2018-03-31?21:08??EMD\HilbertAnalysis.m
?????文件?????????294??2018-03-31?21:09??EMD\isimf.m
?????文件?????????176??2018-03-31?21:07??EMD\ismonotonic.m
?????文件????????1688??2018-03-31?21:09??EMD\plot_hht.m
?????文件????????1784??2018-03-31?21:46??EMD\test.m