function?imf?=?emd（x）
%?Empiricial?Mode?Decomposition?（Hilbert-Huang?Transform）
%?EMD分解或HHT變換
%?返回值為cell類型，依次為一次IMF、二次IMF、...、最后殘差
x???=?transpose（x（:））;
imf?=?[];
while?~ismonotonic（x）
????x1?=?x;
????sd?=?Inf;
????while?（sd?>?0.1）?||?~isimf（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;
%?是否單調
function?u?=?ismonotonic（x）
u1?=?length（findpeaks（x））*length（findpeaks（-x））;
if?u1?>?0
????u?=?0;
else
????u?=?1;
end
%?是否IMF分量
function?u?=?isimf（x）
N??=?length（x）;
u1?=?sum（x（1:N-1）.*x（2:N）?u2?=?length（findpeaks（x））+length（findpeaks（-x））;????%?極值點的個數
if?abs（u1-u2）?>?1
????u?=?0;
else
????u?=?1;
end
%?據極大值點構造樣條曲線
function?s?=?getspline（x）
N?=?length（x）;
p?=?findpeaks（x）;
s?=?spline（[0?p?N+1][0?x（p）?0]1:N）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2017-11-21?19:51??EMD\
?????文件????????1081??2017-11-21?19:47??EMD\emd.m
?????文件?????????912??2017-11-21?19:50??EMD\findpeaks.m
?????文件?????????235??2017-11-21?19:48??EMD\HilbertAnalysis.m
?????文件????????1688??2017-11-21?19:49??EMD\plot_hht.m
?????文件????????1114??2017-11-21?19:51??EMD\test.m