eugene1111 发表于 2007-5-20 18:54

缩尺模型加速度记录经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图候…想得瞬时频率图、瞬时幅值图》
感谢

zhlong 发表于 2007-5-20 19:09

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

zhangnan3509 发表于 2007-5-20 19:14

回复 #1 eugene1111 的帖子

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

eugene1111 发表于 2007-5-20 22:11

回复 #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:

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

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

zhlong 发表于 2007-5-21 09:52

回复 #4 eugene1111 的帖子

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

sharc 发表于 2007-5-24 22:26

新建EMD的qq群 主要是初学者的 希望大家参加 40652876

页: [1]
查看完整版本: 缩尺模型加速度记录经emd这,一些疑问