matlab经验模态分解代码_经验模态分解法

注：以黄河流域花园口水文站年径流序列为例

经验模态分解(EMD)

clc;clear;close alldata=xlsread('testdata.xlsx');t=data(:,1);x=data(:,3);imf=emd(x)';imf2=cell2mat(imf)';%https://blog.youkuaiyun.com/cxj540947672/article/details/82663189k=[x,imf2];%合并矩阵x=k;[M,N]=size(x);%M=743(行)；N=10(列) for i=1:N %i从1到5  subplot(N,1,i);%5行1列，第i个区    %subplot(m，n，p)是将多个图画到一个平面上的工具。    %m和n代表在一个图像窗口中显示m行n列个图像，后面的p代表现在选定第p个图像区域，即在第p个区域作图。  plot(t,x(:,i));%x轴为时间，y轴为第二列的数据    % axis([0 800 0 1.1]); %axis一般用来设置axes的样式，包括坐标轴范围，可读比例等.     %axis([tmin tmax ymin ymax])     axis([1919 2018,-inf,inf])    %xlabel('year');    if i==N    xlabel('year');%第一列为时间，第二列为原始数据。    end  if i==1    ylabel('Y');%第一列为时间，第二列为原始数据。  elseif i==N    ylabel('R');  else    ylabel(strcat('Imf',num2str(i-1)));  endend

集合经验模态分解法(EEMD)

% data为你的一列数据clc;clear;close allData=xlsread('testdata.xlsx');t=Data(:,1);data=Data(:,3);rslt=eemd(data,0.2,100); % 0.2可以改为0.2-0.3，表示信噪比；100可以表示为100-1000；若0.2和100变为0,1；那么EEMD简化为EMD[m,n]=size(rslt);imfs=rslt(:,2:n-1);[sigline50,logep]=significance(imfs,0.5);[sigline80,logep]=significance(imfs,0.2);[sigline90,logep]=significance(imfs,0.1);[sigline95,logep]=significance(imfs,0.05);figure(1)plot(sigline50(:,1),sigline50(:,2),'b');  %  50percenta linehold onplot(sigline80(:,1),sigline80(:,2),'g');  %  80percenta lineplot(sigline90(:,1),sigline90(:,2),'y');  %  90percenta lineplot(sigline95(:,1),sigline95(:,2),'m-');  % 95 percenta lineplot(logep(:,1),logep(:,2),'r*');  plot(logep(1,1),logep(1,2),'b*');grid;xlabel('LOG2 ( Mean Period )');ylabel('LOG2 ( Mean Normalized Energy )');title('Significance test of IMFs of white noise');% axis([16-70]) %1,6表示横坐标范围；-7,0表示纵坐标范围；横纵坐标都可以改figure(2)x=rslt;[M,N]=size(x);%M=743(行)；N=10(列) for i=1:N %i从1到5  subplot(N,1,i);%5行1列，第i个区    %subplot(m，n，p)是将多个图画到一个平面上的工具。    %m和n代表在一个图像窗口中显示m行n列个图像，后面的p代表现在选定第p个图像区域，即在第p个区域作图。  plot(t,x(:,i));%x轴为时间，y轴为第二列的数据    % axis([0 800 0 1.1]); %axis一般用来设置axes的样式，包括坐标轴范围，可读比例等. axis([xmin xmax ymin ymax])        axis([1919 2018,-inf,inf])    if i==N    xlabel('year');  end  if i==1    ylabel('Y');%第一列为时间，第二列为原始数据。  elseif i==N    ylabel('R');  else    ylabel(strcat('Imf',num2str(i-1)));  endend

完备经验模态分解方法(CEEMD)

clc;clear;close alldata=xlsread('testdata.xlsx');Y=data(:,3);t=data(:,1);% 输入时间序列Nstd=0.2;NE=20;            % 迭代次数TNM=8;            % 分量，开头一列原始数据，结尾一列趋势项A=ceemd(Y,Nstd,NE,TNM);x=A;[M,N]=size(x);%M=743(行)；N=10(列) for i=1:N %i从1到5  subplot(N,1,i);%5行1列，第i个区    %subplot(m，n，p)是将多个图画到一个平面上的工具。    %m和n代表在一个图像窗口中显示m行n列个图像，后面的p代表现在选定第p个图像区域，即在第p个区域作图。  plot(t,x(:,i));%x轴为时间，y轴为第二列的数据    % axis([0 800 0 1.1]); %axis一般用来设置axes的样式，包括坐标轴范围，可读比例等. axis([xmin xmax ymin ymax])        axis([1919 2018,-inf,inf])    if i==N    xlabel('year');  end  if i==1    ylabel('Y');%第一列为时间，第二列为原始数据。  elseif i==N    ylabel('R');  else    ylabel(strcat('Imf',num2str(i-1)));  endend

极点对称模态分解(ESMD)