function?imf?=?emd（x）
%?Empiricial?Mode?Decomposition?（Hilbert-Huang?Transform）
%?imf?=?emd（x）
%?Func?:?findpeaks

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;

%?FUNCTIONS

function?u?=?ismonotonic（x）

u1?=?length（findpeaks（x））*length（findpeaks（-x））;
if?u1?>?0?u?=?0;
else??????u?=?1;?end

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

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件??????958087??2007-01-25?14:59??HHT.pdf
?????文件?????????965??2008-04-24?13:36??plot_hht.m
?????文件?????????893??2008-03-13?11:49??emd.m
?????文件?????????146??2007-08-26?00:42??findpeaks.m
?????文件??????102892??2007-04-28?15:57??Hum.wav
?????文件??????????12??2017-06-30?17:28??read?me.txt