本帖最后由 牛小贱 于 2014-5-22 12:47 编辑
- clear;clc
- close all;
- t=1:1000;
- t1=t/100*2*pi;
- a1=sin(t1);
- t2=t/10*2*pi;
- b1=linspace(0,0,1000);
- for i=250:350
- b1(i)=0.2*sin(t2(i));
- end
- for i=750:850
- b1(i)=0.2*sin(t2(i));
- end
- x=a1+b1; %x是原信号
- subplot(311);plot(t,a1);
- subplot(312);plot(t,b1);
- subplot(313);plot(t,x);
- Nstd=0.004;
- NE=100; % # of ensemble
- numImf =2; % # of imfs
- [imf]= rcada_eemd(x,Nstd,NE,numImf); % run EEMD [3]
- figure;
- nimf = size(imf,1);
- for m=1:nimf
- subplot(nimf,1,m); plot(t,imf(m,:));
- end
- subplot(nimf,1,1); title('IMFs');
复制代码 国立大学新出的FEEMD很好用,可以限定IMF个数(numImf =2; ),结果如下
|