缩尺模型加速度记录经emd这,一些疑问
小弟将加速度记录经置顶文章内下载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)'
= emd(acc4);
再经 =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图候…想得瞬时频率图、瞬时幅值图》
感谢 感觉你的信号没怎么分解。
从你给的瞬时频率图看,分量的频率范围都很宽,说明EMD分解不怎么有效。
可能要结合其它方法同时使用。
回复 #1 eugene1111 的帖子
另外因为瞬时频率是归一化频率因此在0-0.5之间也不足为奇,要乘以采样频率才能行:lol回复 #2 zhlong 的帖子
请问zhlong …我直接将加速度记录进行 = 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=; %IMF集
%%%%%%%%%%%%%求边际谱
=hhspectrum(y); plot(tt,fa*1000); plot(tt,A);
画出的图可以很清楚看出 频率在30、45而振富在2、5…:victory:
所以小弟仿照「破凤」大大文章的做法…但画出来好怪耶…
小弟另外试其它振动台资料…问题也一样><~~:'(
试过 数值模拟方式…想说用模拟方式可准确…问题也是:'(
不知大大们能否和我找出问题呢…
感恩
回复 #4 eugene1111 的帖子
仿真有效不能说明用于实际信号就有效啊也有不少文章是用EMD研究你的这种结构振动的,你可以参考一下别人的做法。
我也不是做这个的
新建EMD的qq群 主要是初学者的 希望大家参加 40652876
页:
[1]