%%?EMD-WP方法：

%?作者：************
%?時(shí)間：2012年************



function?Emd_Wp（dataNNk）

X=data（1:NN）;

%%?參數(shù)設(shè)置
fs=12000;?%采樣頻率
N?=?length（X）;
y_imf=emd（X）;?%Emd處理
[y_IMFIMF_lastrest]=IMF_Choose（Xy_imfk）;%本征模分量的選擇
[mn]=size（y_IMF）;
f=（0:N-1）*fs/N;?%頻率的變化范圍。
t=（0:N-1）/N*（N/fs）;?%時(shí)間的變化范圍。
Fs=1/（t（2）-t（1））;


%%?EMD處理后的IMF圖
if?（m-1）>=4???
????figure;?
????for?i=1:4?%1~4
????????subplot（41i）;
????????plot（ty_IMF（i:））;
????????title（[‘IMF‘num2str（i）]）;%num2str：數(shù)字轉(zhuǎn)化為字符串
????end????
????if?floor（（m-1）/4）>=2?%5~8。floot：取商???
????????figure;??
????????for?i=5:8????????
????????????subplot（41i-4）;????????
????????????plot（ty_IMF（i:））;????
????????????title（[‘IMF‘num2str（i）]）;
????????end???????
????????if?floor（（m-1）/4）>=3%9~12
????????????figure;?
????????????for?i=9:12????????
????????????????subplot（41i-8）;????????
????????????????plot（ty_IMF（i:））;?
????????????????title（[‘IMF‘num2str（i）]）;
????????????end
????????????if?floor（（m-1）/4）==3%13~...
????????????????figure;
????????????????for?i=13:（mod（（m-1）4）+12）??????
????????????????????subplot（41i-12）;????????
????????????????????plot（ty_IMF（i:））;
????????????????????title（[‘IMF‘num2str（i）]）;
????????????????end
????????????end
????????else
????????????if?m>9
????????????figure;?
????????????for?i=9:（rem（（m-1）4）+8）%9~...。?取余
????????????????subplot（41i-8）;????????
????????????????plot（ty_IMF（i:））;????
????????????????title（[‘IMF‘num2str（i）]）;
????????????end
????????????end
????????end?????
????else
????????if?m-1~=4
????????????figure;
????????????for?i=5:（rem（（m-1）4）+4）%5~...。
????????????????subplot（41i-4）;????????
????????????????plot（ty_IMF（i:））;????
????????????????title（[‘IMF‘num2str（i）]）;
????????????end
????????end
????end
else
????figure;
????for?i=1:m-1?
????????subplot（31i）;
????????plot（ty_IMF（i:））;
????????title（[‘IMF‘num2str（i）]）;%num2str：數(shù)字轉(zhuǎn)化為字符串
????end
end
%最后一個(gè)剩余分量residual的圖形
figure;
plot（ty_IMF（m:））;
title（‘殘余量‘）;


%%?對(duì)IMF進(jìn)行一維小波包的分解，返回小波包結(jié)構(gòu)樹T
%?T=wpdec（XN‘wname‘EP）
%?wname=input（‘請(qǐng)輸入進(jìn)行小波包分解的小波函數(shù)（僅限db小波）：‘‘s‘）;
%?num=input（‘請(qǐng)輸入進(jìn)行小波包分解的層數(shù)：‘）;
%?for?i=1:m-1
%?????T=wpdec（y_IMF（i:）numwname‘shannon‘）;
%?????plot（T）;
%?end


%%?選擇感興趣的IMF。
int_IMF=input（‘請(qǐng)輸入對(duì)第幾個(gè)IMF感興趣（數(shù)字）：‘）;


%%?用W