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）;?????????%?極大值點(diǎn)樣條曲線
????????s2?=?-getspline（-x1）;???????%?極小值點(diǎn)樣條曲線
????????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;

%?是否單調(diào)
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（find

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2017-04-09?15:08??EMD\
?????文件?????????227??2012-10-08?22:32??EMD\FFTAnalysis.m
?????文件?????????250??2012-10-08?22:33??EMD\HilbertAnalysis.m
?????文件????????1108??2012-10-08?22:32??EMD\emd.m
?????文件????????2501??2017-02-02?20:34??EMD\exe_emd_20170201.asv
?????文件????????2129??2017-04-09?15:08??EMD\exe_emd_20170201.m
?????文件????????2595??2017-04-05?09:00??EMD\exe_emd_20170201_moni.asv
?????文件?????????925??2012-10-08?22:30??EMD\findpeaks.asv
?????文件?????????925??2012-10-08?22:30??EMD\findpeaks.m
?????文件????????1445??2007-06-27?22:12??EMD\hhspectrum.m
?????文件????????3427??2005-10-13?23:15??EMD\instfreq.m
?????文件????????1721??2012-10-12?15:48??EMD\plot_hht.asv
?????文件????????1755??2017-02-04?12:04??EMD\plot_hht.m
?????文件?????????595??2012-10-11?22:50??EMD\test_bianjipu.m
?????文件????????1330??2017-01-26?17:01??EMD\test_emd.asv
?????文件?????????442??2012-10-13?09:01??EMD\test_emd_shu.asv
?????文件?????????274??2012-10-11?22:51??EMD\youyige_bianjipu.m