调用eemd是三个参数,信号,噪声参数和分解次数。
文件夹包括了emd EEMD 和cemd函数:https://blog.youkuaiyun.com/yyyyllllxxxx/article/details/102760130
(感谢楼主分享!)
下图是对于上篇文章中代码的补充:上篇文章中代码需要调用作者封装的工具箱,这个是直接调用的开源的EMD和EEMD工具箱即可。
fs = 400; %采样频率
t = 0:1/fs:0.75; %时间轴
x = sin(2*pi*4*t); %低频正弦信号
y = 0.5*sin(2*pi*120*t); %高频正弦信号
for i = 1:length(t) %将高频信号处理成间断性
if mod(t(i),0.25)>0.11&&mod(t(i),0.25)<0.12
else
y(i) = 0;
end
end
sig = x+y; %信号叠加
%emd分解画图,我给注释掉了,作图的时候解除注释即可
% imf=emd(sig);
% emd_visu(sig,t,imf) % EMD画图调用
%eemd分解及作图
imf2=eemd(sig,0.2,100);
m = size(imf2,2);
for i = 1:m
subplot(m,1,i)
plot(imf2(:,i))
end
emd分解结果:(第一个图是原信号)
eemd分解结果: