HHT健康诊断使用相关问题

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图候…想得瞬时频率图、瞬时幅值图》
感谢

junzifei

看不懂你在说什么
怎么那么多错字呀？
很费神的