%當X是一個真正的矢量時，計算X的經驗模態分解，
%每一行包含一個max的經驗模態？最后一個是殘余向量，默認準則在文獻2提出
%本函數用于對原始信號varargin進行分解，imf為分解出的本征模態向量
%函數返回值imf為本征模態向量ort=一個數，nbits為1*4的向量？
function?[imfortnbits]?=?emd（varargin）

[xtsdsd2tolMODE_COMPLEXndirsdisplay_siftingsdtsd2trimfknbitNbItMAXITERATIONSFIXEFIXE_HMAXMODESINTERPmask]?=?init（varargin{:}）;

if?display_sifting
??fig_h?=?figure;
end


%主循環:至少需要3極值才能繼續
while?~stop_EMD（rMODE_COMPLEXndirs）?&&?（k?
??%?當前模式
??m?=?r;
??%?在前一個迭代模式
??mp?=?m;
??%計算平均值和停止準則
??if?FIXE
????[stop_siftmoyenne]?=?stop_sifting_fixe（tmINTERPMODE_COMPLEXndirs）;
??elseif?FIXE_H
????stop_count?=?0;
????[stop_siftmoyenne]?=?stop_sifting_fixe_h（tmINTERPstop_countFIXE_HMODE_COMPLEXndirs）;
??else
????[stop_siftmoyenne]?=?stop_sifting（mtsdsd2tolINTERPMODE_COMPLEXndirs）;
??end

??%?in?case?the?current?mode?is?so?small?that?machine?precision?can?cause
??%如果當前模式太小機器精度的原因
??%?spurious?extrema?to?appear
??%偽極值出現
??if?（max（abs（m）））?????if?~stop_sift
??????warning（‘emd:warning‘‘forced?stop?of?EMD?:?too?small?amplitude‘）
????else
??????disp（‘forced?stop?of?EMD?:?too?small?amplitude‘）
????end
????break
??end


??%?篩選循環
??while?~stop_sift?&&?nbit
????if（~MODE_COMPLEX?&&?nbit>MAXITERATIONS/5?&&?mod（nbitfloor（MAXITERATIONS/10））==0?&&?~FIXE?&&?nbit?>?100）
??????disp（[‘mode?‘int2str（k）‘?iteration?‘int2str（nbit）]）
??????if?exist（‘s‘‘var‘）
????????disp（[‘stop?parameter?mean?value?:?‘num2str（s）]）
??????end
??????[imiM]?=?extr（m）;
??????disp（[int2str（sum（m（im）?>?0））‘?minima?>?0;?‘int2str（sum（m（iM）?????end

????%篩選
????m?=?m?-?moyenne;

????%計算平均值和停止準則
????if?FIXE
??????[stop_siftmoyenne]?=?stop_sifting_fixe（tmINTERPMODE_COMPLEXndirs）;
????elseif?FIXE_H
??????[stop_siftmoyennestop_count]?=?stop_sifting_fixe_h（tmINTERPstop_countFIXE_HMODE_COMPLEXndirs）;
????else
??????[stop_siftmoyennes]?=?stop_sifting（mtsdsd2tolINTERPMODE_COMPLEXndirs）;
????end

????%?display
????if?display_sifting?&&?~MODE_COMPLEX
??????NBSYM?=?2;
??????[indminindmax]?=?extr（mp）;
??????[tmintmaxmminmmax]?=?boundary_conditions（indminindmaxtmpmpNBSYM）;
??????envminp?=?interp1（tminmmintINTERP）;
??????envmaxp?=?interp1（tmaxmmaxtINTERP）;
??????envmoyp?=?（envminp+envmaxp）/2;
??????if?FIXE?||?FIXE_H
????????display_emd_fixe（tmmprenvminpenvmaxpenvmoypnbitkdisplay_sifting）
??????else
????????sxp=2*（abs（envmoyp））./（abs（envmaxp-envminp））;
????????sp?=?mean（sxp）;
????????display_emd（tmmprenvminpenvmaxpenvmoypsspsxpsdtsd2tnbitkdisplay_siftingstop_sift）
??????end
????end

????mp?=?m;
????nbit=nbit+1;
????NbIt=NbIt+1;

????if（nbit==（MAXITERATIONS-1）?&&?~FIXE?&&?nbit?>?100）
??????if?exist（‘s‘‘var‘）
????????warning（‘emd:warning‘[‘forced?stop?of?sifting?:?too?many?iterations...?mode?‘int2str（k）‘.?stop?parameter?mean?value?:?‘num2str（s）]）
??????else
????????warning（‘emd:warning‘[‘forced?stop?of?sifting?:?too?many?ite

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????1701??2015-09-02?21:36??EMD樣本熵特征提取\Emdplus.m
?????文件???????17907??2015-09-02?21:37??EMD樣本熵特征提取\emd.m
?????文件?????????656??2015-09-02?21:35??EMD樣本熵特征提取\samp.m
?????目錄???????????0??2015-09-02?22:36??EMD樣本熵特征提取\