经验模态分解Matlab代码

clearclcclose all;

ss=importdata（‘6.mat‘）;
s=awgn（ss60‘measured‘）;

% t=0:1e-3:1;
% s=2*sawtooth（2*pi*15*t）;%+4*sin（2*pi*10*t）.*sin（2*pi*0.1*t）+sin（2*pi*5*t）;

% t=0:1e-3:1;
% s1=2*sin（2*pi*5*t）;
% s2=0.3*sin（2*pi*100*t）;
% s2（1:450）=zeros（1450）;
% s2（550:end）=zeros（1452）;
% s=s1+s2;

sbk=s;
n=length（s）;x=1:n;
q=0;
while（1）
q=q+1;
qq=0;
while（1）
qq=qq+1;
j=0;k=0;p=0;
for i=2:（n-1）
    if（s（i）>s（i+1） && s（i）>s（i-1））
        j=j+1;
        jd（j）=s（i）;xd（j）=i;
    end
    if（s（i）        k=k+1;
        jx（k）=s（i）;xx（k）=i;
    end
    if（s（i-1）*s（i）<=0）
        p=p+1;
    end
end
% if（abs（j+k-p）<=1 || qq>10）
%     clear jd jx xd xx ub db mid;
%     break;
% end

if s（1）>s（2）
    dlp=s（1）;xlp=-s（1）;