|
楼主 |
发表于 2012-10-18 17:06
|
显示全部楼层
本帖最后由 yghit08 于 2012-10-18 17:16 编辑
danielrichie 发表于 2012-10-18 16:51
首先非常感谢楼主的提点,用EMD中Hilbert谱来画时频图之前有做过,就是因为对LMD原始算法中不用Hilbert谱 ...
我想画时频谱是后处理过程,不清楚楼主是否明白了EMD中利用Hilbert谱显示时频图的各个子函数的功能。如我在回复你的那楼中所示的例子来说,假设分解完成各阶PF后。接着,需要做的就是对MF矩阵做求解瞬时频率的求解。如果这部分不是利用Hilbert谱完成的,那假设是用其他方法完成的,那么对应出来的就是瞬时频率的矩阵。不论是利用Hilbert谱还是其他方法,接下来就是将分解的结果做时频显示了。函数是[im,tt]=toimage(A,F),其中A是瞬时幅值矩阵(对应AF),f是瞬时频率的矩阵(不论你利用什么做的瞬时频率的求解,此处注意建议在利用其他方法做瞬时频率的求解是也做归一化处理,免得后续麻烦。注意每一阶存成一行且和AF对应上),接下来就是显示disp_hhs(im,tt/SampleFrequency,[],SampleFrequency),其中SampleFrequency是采样频率。不知道这样说你是否明白??如果还是不明白,那咱们就用实例说话。你把你的求解瞬时频率的程序贴上来,我做一个例子给你看!
|
|