%?EMD:??Emprical?mode?decomposition
%
%?imf?=?emd（xn）
%
%?x???-?input?signal?（must?be?a?column?vector）
%?n???-?number?of?intrinsic?mode?functions
%
%?imf?-?Matrix?of?intrinsic?mode?functions?（each?as?a?row）
%
%?See:??Huang?et?al?Royal?Society?Proceedings?on?Math?Physical?
%???????and?Engineering?Sciences?vol.?454?no.?1971?pp.?903-995?
%???????8?March?1998
%
%?Author:?Ivan?Magrin-Chagnolleau??
%?
function?imf?=?emd（xn）;%%最好把函數名改為emd1之類的，以免和Grilling的emd沖突????????????????????????
??????????????????????????????????%%n為你想得到的IMF的個數
c?=?x（:）‘;?%?copy?of?the?input?signal?（as?a?row?vector）
N?=?length（x）;-
%?loop?to?decompose?the?input?signal?into?n?successive?IMFs
imf?=?[];?%?Matrix?which?will?contain?the?successive?IMF?and?the?residuefor?t=1:n?

while