eemd分解matlab代码

以下是一个使用MATLAB实现的EEMD分解代码示例：

function [imf,residue]=emd(x)
% x:	信号数据
% imf:	分解出的各个imf成分
% residue:	剩余的趋势项
% n:	当前是第几次分解

if length(x)<3
    imf=x;
    residue=zeros(size(x));
    return;
end

% 计算均值
m=mean(x);
% 计算极值
u=max(x);
d=min(x);

if u==d
    imf=x;
    residue=zeros(size(x));
    return;
end

% 计算极值均值
h=(u+d)/2;

% 构造辅助函数
x1=x(1);
x2=x(2);
if x1>=h && x2>=h
    s=1;
elseif x1<=h && x2<=h
    s=-1;
else
    s=0;
end
p=0;
for i=2:length(x)-1
    x0=x(i);
    if (x1<=h && x0>=h) || (x1>=h && x0<=h)
        p=p+1;
        imf1(p)=(x1+x0)/2;
    end
    x1=x0;
end

if p<=1
    imf=x;
    residue=zeros(size(x));
    return;
end

imf2=emd(imf1);
imf3=emd(x-imf2);
imf=[imf2 imf3];
residue=x-sum(imf,2);

end

这段代码实现了EEMD分解，将信号进行分解得到各个imf成分和剩余的趋势项。