test3.m

来自「新的信号处理方法－希尔伯特黄变换的EMD分解MATLAB程序 把程序看懂 对希尔」· M 代码 · 共 16 行

clear
clc
t=0:0.01:10;
a=sin(20*pi*t);
b=5*sin(10*pi*t);
c=5*cos(2*pi*t);
x=a+b+c;
stop = [0.1,0.1,0.1];
[imf,ort,nbits]=myemd(x,t,stop,2);


imfsize=size(imf);
columnNum=imfsize(1)/2+1;
for i=1:columnNum(1)
    subplot(columnNum,2,i);plot(imf(i,:));
end