缩尺模型加速度记录经emd这，一些疑问

eugene1111

小弟将加速度记录经置顶文章内下载emd相关matlab码候，
小弟下载置顶emd matlab码候…求得付图的imf

fn=strcat('D:\Studay\ID_damge\k008_nnn_acc.txt') % Hz=200; T=0.005; 笔数5100
s=load(fn);
acc4=s(:,4)'
[imf,ort,nbits] = emd(acc4);

再经 [A,fa,tt]=hhspectrum(imf);

问题产生了…小弟不懂…求大大帮解…

1_欲求得每个imf分量的瞬时频率得付图…感觉很不合理…指令 plot(tt,fa)
    小弟之前用类神经求得结构自然频率约 (1.52 ，4.64，7.33，9.23)，在此图看不出来…
    在图上好乱，错掉了…不知何处有问题…
    理因会在y座标 (1.52 ，4.64，7.33，9.23)处，会很明显…但小得不到
    另外y座标好奇怪…范围0~0.5…y座标不是频率吗…范围0~10…

2_也想得imf分量瞬时幅值…如附图…指令 plot(tt,A)
    这里因该合理…有四个明显趋势…

3_另外小弟想画出 边际谱…翻过讨论区范例…但一直搞不懂…

加速度记录上传於《http://baggio.myweb.hinet.net/k008_nnn_acc.rar》
格式 time floor1 floor2 floor3 floor4 base


不知小弟那里出问题…希望大大帮帮我…
《小弟得到imf图候…想得瞬时频率图、瞬时幅值图》
感谢

zhlong

感觉你的信号没怎么分解。
从你给的瞬时频率图看，分量的频率范围都很宽，说明EMD分解不怎么有效。
可能要结合其它方法同时使用。

zhangnan3509

另外因为瞬时频率是归一化频率因此在0-0.5之间也不足为奇，要乘以采样频率才能行:lol

eugene1111

请问zhlong …我直接将加速度记录进行 [imf,ort,nbits] = emd(acc4); 这样的程序有误吗
有需要调整参数吗?  我用了 新版将舊版产生imf…画出的图都一样耶…

照了破凤文章
fs=1000;    %fs为采样频率;
N=1000;     %采样点数
t=1/fs:1/fs:1;
y1=2*sin(60*pi*t);
y2=5*sin(90*pi*t);
y=[y1;y2;zeros(size(y1))]; %IMF集
%%%%%%%%%%%%%求边际谱
[A,fa,tt]=hhspectrum(y);    plot(tt,fa*1000);   plot(tt,A);

画出的图可以很清楚看出 频率在30、45  而振富在2、5…:victory:

所以小弟仿照「破凤」大大文章的做法…但画出来好怪耶…
小弟另外试其它振动台资料…问题也一样><~~:'(
试过 数值模拟方式…想说用模拟方式可准确…问题也是:'(

不知大大们能否和我找出问题呢…
感恩

zhlong

仿真有效不能说明用于实际信号就有效啊
也有不少文章是用EMD研究你的这种结构振动的，你可以参考一下别人的做法。
我也不是做这个的

sharc