之前在CSDN上下载这些东西要积分,下载不了,所以搞了个这样的,亲测有效
集成经验模式分解eemd,一种改进的emd的信号分解方法matlab代码实现
function allmode=eemd(Y,Nstd,NE)
% This is an EMD/EEMD program
%
% INPUT:
% Y: Inputted data;1-d data only
% Nstd: ratio of the standard deviation of the added noise and that of
% Y; Nstd = (0.1 ~ 0.4)std(Y).
% NE: Ensemble number for the EEMD, NE = 10-50.
% OUTPUT:
% A matrix of N(m+1) matrix, where N is the length of the input
% data Y, and m=fix(log2(N))-1. Column 1 is the original data, columns 2, 3, …
% m are the IMFs from high to low frequency, and comlumn (m+1) is the
% residual (over all trend).
%
% NOTE:
% It should be noted that when Nstd is set to zero and NE is set to 1, the
% program degenerates to a EMD program.(for EMD Nstd=0,NE=1)
% This code limited sift number=10 ,the stoppage criteria can’t change.
% Re