function?imf=emd（x）
%EMD分解或HHT變換
%返回值為cell類型，依次為一次IMF、二次IMF、...、最后殘差
x=transpose（x（:））;?????%對(duì)x進(jìn)行轉(zhuǎn)置
imf=[];
while~ismonotonic（x）?%當(dāng)x不是單調(diào)函數(shù)時(shí)
????x1=x;
????sd=Inf;%inf是無(wú)窮大
????while?（sd?>?0.1）?||?~isimf（x1）??%?當(dāng)sd>0.1并且不滿足isimf（x1）這個(gè)條件，循環(huán)
????????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（findpeaks（x））+length（findpeaks（

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2016-09-21?16:59??EMD程序\
?????目錄???????????0??2016-09-09?14:25??EMD程序\EMD\
?????文件????????1090??2016-07-29?17:05??EMD程序\emd.m
?????文件???????18487??2016-09-14?15:42??EMD程序\emd1.m
?????目錄???????????0??2016-09-12?11:31??EMD程序\EMD\復(fù)雜中斷\
?????文件???????12200??2010-06-06?19:07??EMD程序\EMD\復(fù)雜中斷\3435567788.fig
?????文件???????19815??2016-09-12?11:34??EMD程序\EMD\復(fù)雜中斷\emd.m
?????文件????????1133??2010-04-26?23:58??EMD程序\EMD\復(fù)雜中斷\emdyy.m
?????文件????????2850??2010-04-26?23:06??EMD程序\EMD\復(fù)雜中斷\hhspectrum.m
?????文件????????2299??2010-06-13?00:03??EMD程序\EMD\復(fù)雜中斷\IMFeg.asv
?????文件????????2378??2016-09-12?11:31??EMD程序\EMD\復(fù)雜中斷\IMFeg.m
?????文件????????3003??2007-06-27?22:12??EMD程序\EMD\復(fù)雜中斷\toimage.m
?????目錄???????????0??2010-06-23?12:28??EMD程序\EMD\復(fù)雜暫降（幅值）\
?????文件???????17417??2010-04-26?23:21??EMD程序\EMD\復(fù)雜暫降（幅值）\emd.m
?????文件????????1133??2010-04-26?23:58??EMD程序\EMD\復(fù)雜暫降（幅值）\emdyy.m
?????文件????????1915??2010-06-16?20:09??EMD程序\EMD\復(fù)雜暫降（幅值）\hhspectrum.m
?????文件????????2292??2010-06-12?18:13??EMD程序\EMD\復(fù)雜暫降（幅值）\IMFeg.asv
?????文件????????1692??2010-06-18?12:25??EMD程序\EMD\復(fù)雜暫降（幅值）\IMFeg.m
?????文件????????1942??2010-06-16?20:09??EMD程序\EMD\復(fù)雜暫降（幅值）\toimage.m
?????目錄???????????0??2010-06-23?12:28??EMD程序\EMD\復(fù)雜暫降（時(shí)刻）\
?????文件???????17417??2010-04-26?23:21??EMD程序\EMD\復(fù)雜暫降（時(shí)刻）\emd.m
?????文件????????1133??2010-04-26?23:58??EMD程序\EMD\復(fù)雜暫降（時(shí)刻）\emdyy.m
?????文件????????1915??2010-06-16?20:09??EMD程序\EMD\復(fù)雜暫降（時(shí)刻）\hhspectrum.m
?????文件????????2292??2010-06-12?18:13??EMD程序\EMD\復(fù)雜暫降（時(shí)刻）\IMFeg.asv
?????文件????????1632??2010-06-17?21:09??EMD程序\EMD\復(fù)雜暫降（時(shí)刻）\IMFeg.m
?????文件????????1942??2010-06-16?20:09??EMD程序\EMD\復(fù)雜暫降（時(shí)刻）\toimage.m
?????目錄???????????0??2016-09-19?11:27??EMD程序\EMD\復(fù)雜諧波\
?????文件???????17417??2010-04-26?23:21??EMD程序\EMD\復(fù)雜諧波\emd.m
?????文件????????1133??2010-04-26?23:58??EMD程序\EMD\復(fù)雜諧波\emdyy.m
?????文件????????1916??2010-06-16?20:10??EMD程序\EMD\復(fù)雜諧波\hhspectrum.m
?????文件????????1626??2010-06-16?20:10??EMD程序\EMD\復(fù)雜諧波\IMFeg.m
............此處省略22個(gè)文件信息