huangyong87 发表于 2007-5-4 18:21

若已知瞬时频率与幅值,如何重构该信号?

原来认为这个问题很简单,若已知某一信号成分的瞬时频率F(t)与瞬时幅值A(t),重构信号X(t)=A(t)*cos(F(t)*t*2*pi);即可。不想数值试验时发现,这只对瞬时频率F(t)为常量的情况适用,当瞬时频率F(t)为时变的情况,比如x=(1+0.6*cos(pi*t/100)).* cos(pi*t/5+4*sin(3*pi*t/200));有调频成分的情况,误差就十分大。
已知瞬时频率与幅值,如何重构该信号?请各位高人指点!

eight 发表于 2007-5-4 19:03

本帖最后由 VibInfo 于 2016-10-11 16:03 编辑

原帖由 huangyong87 于 2007-5-4 18:21 发表
原来认为这个问题很简单,若已知某一信号成分的瞬时频率F(t)与瞬时幅值A(t),重构信号X(t)=A(t)*cos(F(t)*t*2*pi);即可。不想数值试验时发现,这只对瞬时频率F(t)为常量的情况适用,当瞬时频率F(t)为时变的情况 ...
通过上述方法对原信号重构,有一个前提,就是调频不能受到调幅的影响,即原信号必须满足 Bedrosian 定理:解析信号的低频成分在幅度调制上,高频成分在相位调制上。

ps:对这个问题没有仔细研究过,所以上述观点不一定正确,仅供参考,呵呵。

playfish 发表于 2007-5-4 21:09

我想过这个问题,其实就是HHT的逆变换问题。其步骤是:由于HHT中的瞬时频率是由瞬时相位微分而来,因此,要重构信号,首先要重建瞬时相位函数,即对瞬时频率积分,然后与幅度函数一起求得各IMF的解析信号,再把这些信号相加得到原信号的解析信号,取实部得到原信号。这时问题就出来了,先不管函数本身是否可积,积分后会带有常数项,这个常数项怎么确定是个令人头疼的问题。此外,采用什么积分方法也是个不好确定的问题。若瞬时频率函数不可积,这时要求瞬时相位就很麻烦。
这是我的思路,请eight兄指点。

eight 发表于 2007-5-4 21:20

本帖最后由 VibInfo 于 2016-10-11 16:03 编辑

原帖由 playfish 于 2007-5-4 21:09 发表
我想过这个问题,其实就是HHT的逆变换问题。其步骤是:由于HHT中的瞬时频率是由瞬时相位微分而来,因此,要重构信号,首先要重建瞬时相位函数,即对瞬时频率积分,然后与幅度函数一起求得各IMF的解析信号,再把 ...

你的问题与楼主的不同,是另外一个问题。积分常数的问题,不会影响最终结果的吧?因为是出现在复部,而重构需要的是实部

playfish 发表于 2007-5-4 22:21

假设f(t)是频率函数,a(t)是幅度函数,th(t)+c是f(t)的积分,c是常数。则imf的解析表达式为a(t)*exp(j*(th(t)+c)),由欧拉公式,此式等于a(t)(cos(th(t)+c)+j*sin(th(t)+c)),其实部是a(t)cos(th(t)+c)。从这里看常数项对实部是有影响的。

eight 发表于 2007-5-4 22:35

本帖最后由 VibInfo 于 2016-10-11 16:04 编辑

原帖由 playfish 于 2007-5-4 22:21 发表
假设f(t)是频率函数,a(t)是幅度函数,th(t)+c是f(t)的积分,c是常数。则imf的解析表达式为a(t)*exp(j*(th(t)+c)),由欧拉公式,此式等于a(t)(cos(th(t)+c)+j*sin(th(t)+c)),其实部是a(t)cos(th(t)+c)。从这里看 ...
哦,对,我想错了。那能否不用求积分的方式恢复相位呢?

playfish 发表于 2007-5-5 09:22

因为瞬时频率就是相位微分而来,在只知道频率函数的情况下,好像只能积分来求相位函数了。从这个角度说,HHT是没有逆变换的。

huangyong87 发表于 2007-5-5 11:05

仍以信号x=(1+0.6*cos(pi*t/100)).* cos(pi*t/5+4*sin(3*pi*t/200));为例,其瞬时频率显然可以通过微分得到,为w(t)=0.1*2*pi+3*2*pi/100*cos(3*2*pi*t/400);如何通过积分得到瞬时相位P(t)=pi*t/5+4*sin(3*pi*t/200)是一个关键。若用数值积分计算w(t)的定积分将得到一个常数,而分微段求面积虽然能得到一个时变值,但结果似乎不对。这个积分该如何做呢?

eight 发表于 2007-5-5 11:18

本帖最后由 VibInfo 于 2016-10-11 16:04 编辑

原帖由 playfish 于 2007-5-5 09:22 发表
因为瞬时频率就是相位微分而来,在只知道频率函数的情况下,好像只能积分来求相位函数了。从这个角度说,HHT是没有逆变换的。

按照这个说法,对于一种时频方法,只要频率是基于相位的求导得到的(事实上这是目前比较公认的定义),那么,它就没有逆变换了。从瞬时频率和赋值上重构信号,其可行性值得探讨

zhlong 发表于 2007-8-10 13:03

回复 #8 huangyong87 的帖子

可以参考由加速度信号积分得到速度信号,由速度信号积分得到位移信号的方法。
http://forum.vibunion.com/forum/thread-49413-1-1.html 里面有这方面帖子的链接!

破凰 发表于 2007-8-12 22:57

要想进行重构,关键是要求出信号的初相位。

glwh 发表于 2010-10-6 19:07

像小波系数一样逆变换?

ohaiou 发表于 2011-3-29 10:29

hilbert变换揭示了实信号的Fouier变换的实部和虚部、幅频响应和相频响应之间的相互依性,即它们彼此不是独立的,存在着制约关系。

潘文杰. 傅里叶分析及其应用. 北京:北京大学出版社,2000
程乾生. 数字信号处理. 北京:北京大学出版社,2002
页: [1]
查看完整版本: 若已知瞬时频率与幅值,如何重构该信号?