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??2018-07-20?15:02??經驗模態分解\
?????文件????????1089??2018-07-20?15:01??經驗模態分解\emd.m
?????文件?????????208??2018-07-20?15:01??經驗模態分解\FFTAnalysis.m
?????文件?????????911??2018-07-20?15:00??經驗模態分解\findpeaks.m
?????文件?????????235??2018-07-20?15:02??經驗模態分解\HilbertAnalysis.m
?????文件????????1698??2018-07-20?15:00??經驗模態分解\plot_hht.m
?????文件????????1131??2018-07-20?14:59??經驗模態分解\test1.m