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

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

经验模态分解(EMD)

clc;clear;close alldata=xlsread('testdata.xlsx');t=data(:,1);x=data(:,3);imf=emd(x)';imf2=cell2mat(imf)';%https://blog.csdn.net/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')