function?[?y_quzao?]?=?EMD_quzao（?y?）
%%%利用EMD去噪，y為行向量
%t=1:4000;
%y=sin（t.*（2*pi）/1200）+sin（t.*（2*pi）/600）+sin（t.*（2*pi）/300）+normrnd（00.5^0.514000）;

imf=emd（y）;
[ab]=size（imf）;
for?i=1:a
????En（i）=sum（imf（i:）.^2）/b;
????
????[indmin?indmax?indzer]?=?extr（imf（i:））;
????if?（size（indmin2）+size（indmax2））~=0
????Tn（i）=2*b/（size（indmin2）+size（indmax2））;
????else
?????Tn（i）=100000;
????end
????
end
ET=En.*Tn;
for?j=2:a
?Rn（j）=abs（（ET（j）-ET（j-1））/（sum（ET（1:j-1））/（j-1）））;
??if???Rn（j）>3%一般去2或者3
??y_quzao=sum（imf（j:end:）1）;%%%%去除噪聲的真實信號
??????break
??end
??
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????????597??2016-11-13?22:56??EMDLMD_Denoising\EMD_quzao.m

?????文件????????268??2017-04-25?22:44??EMDLMD_Denoising\Example.m

?????文件???????3658??2016-11-13?22:56??EMDLMD_Denoising\lmd.M

?????文件????????618??2016-11-13?22:56??EMDLMD_Denoising\LMD_quzao.m

?????目錄??????????0??2017-04-25?22:43??EMDLMD_Denoising

-----------?---------??----------?-----??----

?????????????????5141????????????????????5