你贴上的程序我没有细看。我回答你文字写的问题。
首先是引用Hilbert变换将信号变换成分析信号的形式，然后求绝对值这样就能提取出瞬时幅值即IA，那么余下来的部分就用来求取角度，再调用instfreq.m求瞬时频率（默认的情况下是只利用两点的数据，导数）。做归一化是为了去除不同采样频率，程序需要做相应的调整（将时间变量全变成1,2,3,4，...).要想不做归一化，那你自己编写instfreq.m子程序，但是这不是问题，可以还原回来，做出来的瞬时频率乘以采样频率还原回来了。 回复 16 # yghit08 的帖子
多谢！！
对于Grill的程序包进行Hilbert变换的过程还有一点儿疑问，看了好多帖子都没有找到答案~~
在利用IMFs求Hilbert Spectrum时用到的toimage.m，调用语法是=toimage(A,f)，其中f是hhspectrum.m得到的瞬时频率。看toimage.m源程序是发现程序将f的范围限定在之间（大概在94和95行左右，可能版本不同），不明白为什么要做这种限定？求解答~~~
谢谢~~~O(∩_∩)O 本帖最后由 yghit08 于 2012-8-16 14:29 编辑
回复 17 # vera吧噗 的帖子
我不清楚你对Grill中涉及到Hilbert变换的问题，具体可以贴出来。我只说我知道的：在Matlab中Hilbert变换的实现是基于傅里叶变换的（将时域信号变换到频域，然后将负频率部分强制为0再傅里叶逆变换，论坛帖子里看到的），所以Matlab中的Hilbert变换是有边界效应的，所以有一些算法在处理EMD的边界效应时也会应用到Hilbert变换中用来处理边界效应。
具体不清楚为什么限定在之间，可能和采样定理有关。补充的一句是：这就是做归一化的结果：都分布在0-0.5之间。 回复 18 # yghit08 的帖子
明白了~~多谢！！{:{23}:}
又结合论坛里一些帖子说，在MATLAB中频率有的函数是用Fs/2进行归一的，有的函数是用Fs进行归一的。例如数字滤波器设计中，wp和ws都要用Fs/2进行归一，用freqz画出频率响应曲线图中频率域便是用Fs/2进行归一的；而在计算功率谱时，对应的频率却是用Fs进行归一的。所以在调用MATLAB相应函数时要了解它的特性，这样把归一化的频率还原为原始频率才不会发生错误。 DING DING! 怎么把繁体字改成看的懂的中文 ？求详细指导 lcy1989 发表于 2012-4-19 00:23
怎么不能上传txt附件啊？
下不了哦楼主 好想下载看看 我现在也是同样的问题新手刚到IMF这步而且还想请教一下各位高手们 我做的铁路的路肩振动监测 手里的数据是路基各测点的振动加速度数据算不算非平稳信号 做HHT分析会不会比FFT靠谱一些呢 多谢 梦泉 发表于 2012-4-19 10:21
第一张图是用论坛下载的EMD、hhspectrum、toimage、disp_hhs来生成Hilbert谱，第二张图是用诸福磊等编写《 ...
好强大 HIT_LIPENG 发表于 2015-11-10 10:15
我现在也是同样的问题新手刚到IMF这步而且还想请教一下各位高手们 我做的铁路的路肩振动监测 手里 ...
路基振动信号具有信源时变、长度时变、受车速调制的变频非平稳等特征
所以传统的时频域信号分析方法只能分析合成振动的总体强度分布与频谱结构
可以尝试采用HHT分析 Generation 发表于 2015-11-20 13:17
路基振动信号具有信源时变、长度时变、受车速调制的变频非平稳等特征
所以传统的时频域信号分析方法只能 ...
多谢解释 谢谢楼主，学习学习
页:
1
[2]