|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
本人按照于德介老师的《机械故障诊断的Hibert-Huang变换方法》一书中P29页的例子:
x(t)=2*sin(30*pi*t)+4*sin(20*pi*t).*sin(2*pi*1/10*t)+sin(2*pi*5*t),t∈[0,1]进行EMD分解,程序如下:
clc;
clear all;
t=0:0.0001:1;
x=2*sin((2*pi*15)*t)+4*sin((2*pi*10)*t).*sin((2*pi*t)/10)+sin((2*pi*5)*t);
figure(1);
subplot(211);
plot(t,x);
xlabel('t/s');ylabel('x(t)');
imf=emd(x);
figure(2);
for i=1:5
subplot(5,1,i);
plot(t,imf(i,:));
end
可是出来的IMF结果与书上的完全对不上,不知何故?EMD用的是Grilling的2007版本。
下面是出来的结果图。
论坛中的相关讨论(http://forum.vibunion.com/thread-61930-1-1.html),我不清楚huangqiang 是怎么设置的。
|
|