希尔伯特变换简介

本文深入探讨了Hilbert变换的定义与应用,包括其实现原理、物理意义及在信号处理中的作用。介绍了Hilbert变换如何将实值信号转换为解析信号,实现信号的包络和瞬时相位计算。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

Hilbert变换简介

在数学与信号处理的领域中,一个实值函数x(t)x(t)x(t),其Hilbert变换记作x^(t)\hat{x}(t)x^(t)定义为:
x^(t)=H[x(t)]=x(t)∗1πt=1π∫−∞∞x(τ)t−τdτ \begin{array}{ll} \hat{x}(t)& =H[x(t)]\\ &=x(t)*\frac{1}{\pi t}\\ &=\frac{1}{\pi }\int^{\infty}_{-\infty}\frac{x(\tau)}{t-\tau}d\tau \end{array} x^(t)=H[x(t)]=x(t)πt1=π1tτx(τ)dτ
在此H标示为希尔伯特变换。
希尔伯特反变换:
x(t)=H−1[x(t)]=−1π∫−∞∞x^(τ)t−τdτ \begin{array}{ll} x(t)&=H^{-1}[x(t)]\\ &=-\frac{1}{\pi }\int^{\infty}_{-\infty}\frac{\hat{x}(\tau)}{t-\tau}d\tau \end{array}x(t)=H1[x(t)]=π1tτx^(τ)dτ
因此,希尔伯特变换结果x′(t)x'(t)x(t)可以被解读为输入是x(t)x(t)x(t)的线性时不变系统的输出,而此系统的脉冲响应为
h(t)=1πth(t)=\frac{1}{\pi t}h(t)=πt1
因此,Hilbert变换可以看成是将原始信号通过一个滤波器,或者一个转向器,这个系统的冲击响应为h(t)。
在这里插入图片描述
h(t)h(t)h(t)做傅里叶变换,得到
H(ω)=−jsgn(ω)=∣H(ω)∣ejϕ(ω)H(\omega)=-jsgn(\omega)\\ =\lvert H(\omega)\rvert e^{j\phi(\omega)}H(ω)=jsgn(ω)=H(ω)ejϕ(ω)
或者写成
H(ω)={−jω&gt;00ω=0+jω&lt;0H(\omega)= \left\{ \begin{array}{lcl} {-j} &amp;\omega&gt;0\\ 0&amp;\omega=0\\ {+j} &amp;\omega&lt;0 \end{array} \right.H(ω)=j0+jω>0ω=0ω<0
其中sgn()sgn()sgn()为符号函数。
sgn(ω)={1ω&gt;00ω=0−1ω&lt;0sgn(\omega)= \left\{ \begin{array}{ll} {1} &amp;\omega&gt;0\\ 0&amp;\omega=0\\ {-1} &amp;\omega&lt;0 \end{array} \right. sgn(ω)=101ω>0ω=0ω<0
因此x(ω)x(\omega)x(ω)的虚部x^(ω)\hat{x}(\omega)x^(ω)的频率部分可以写成下式:
X^(ω)=X(ω)∗H(ω)=X(ω)∗(−jsgn(ω))=−jX(ω)∗sgn(ω) \begin{array}{ll} \hat{X}(\omega)&amp;=X(\omega)*H(\omega)\\ &amp;=X(\omega)*(-jsgn(\omega))\\ &amp;=-jX(\omega)*sgn(\omega)\\ \end{array} X^(ω)=X(ω)H(ω)=X(ω)(jsgn(ω))=jX(ω)sgn(ω)
即:
X^(ω)={−jX(ω)ω&gt;00ω=0jX(ω)ω&lt;0 \hat{X}(\omega)= \left\{ \begin{array}{ll} -jX(\omega) &amp;\omega&gt;0\\ 0&amp;\omega=0\\ jX(\omega) &amp;\omega&lt;0 \end{array} \right. X^(ω)=jX(ω)0jX(ω)ω>0ω=0ω<0
也可以换成另一个写法()
X^(ω)={−jω∣ω∣X(ω)ω≠00ω=0 \hat{X}(\omega)= \left\{ \begin{array}{ll} -j\frac{\omega}{|\omega|}X(\omega) &amp;\omega\neq0\\ 0&amp;\omega=0\\ \end{array} \right. X^(ω)={jωωX(ω)0ω̸=0ω=0
因为X(t)X(t)X(t)的频谱X(ω)X(\omega)X(ω)的导数为
X′(ω)=jωX(ω)X&#x27;(\omega)=j\omega X(\omega)X(ω)=jωX(ω)
所以一个信号的HIlbert变换可以理解为一个信号的平滑梯度(smoothed gradient)。平滑核函数为频域中的∣ω∣−1|\omega|^{-1}ω1

在这里插入图片描述
从频谱上来看,这个滤波器将我们的原始信号的正频率部分乘以−j-jj,也就是说,保持幅度不变的条件下,将相位移动了−pi/2-pi/2pi/2,而对于负频率成分,移动了pi/2pi/2pi/2
下面这个示意图很直观地表示了Hilbert变换,在这里我画出了对原始信号做1到4次Hilbert变换的频谱示意图,是为了说明Hilbert变换的几个性质:
在这里插入图片描述

首先,可以看到,两次希尔伯特变换后,原信号相位翻转了180°,所以,Hilbert逆变换的公式显而易见,就是将正变换加一个符号即可。另外,还可以看到,Hilbert变换四次后就变回本身了。
Hilbert变换的性质有:

  • It is anti-symmetric: H(−x)=−H(x)H(-x)=-H(x)H(x)=H(x)
  • It suppresses the DC component: H(0)=0H(0)=0H(0)=0
  • Its energy is equal to one for all nonzero frequencies: ∣H(x)∣=1∀x≠0\begin{array}{ll} \lvert H(x)\rvert=1 &amp; \forall x\ne0\end{array}H(x)=1x̸=0
  • x(t)x(t)x(t)的Hilbert变换信号x^(t)\hat{x}(t)x^(t)x(t)x(t)x(t)正交:∫ωx^(t)x(t)dt=0\int_{\omega}\hat{x}(t)x(t)dt=0ωx^(t)x(t)dt=0

Hilbert变换物理意义

解析过程概念

一个实值函数x(t)x(t)x(t),其Hilbert变换记作x^(t)\hat{x}(t)x^(t),则x(t)x(t)x(t)的解析信号为:
x~=x(t)+jx^(t)\widetilde{x}=x(t)+j\hat{x}(t)x=x(t)+jx^(t)
这个过程有如下特点,首先,实部和虚部功率谱相同,自相关函数相同;另外,实部和虚部的互相关函数是一个奇函数。
其他的还有:
RX^X(τ)=−R^X(τ)RXX^(τ)=R^X(τ)RX(τ)=2[RX(τ)+jR^X(τ)] R_{\hat{X}X}(\tau)=-\hat{R}_{X}(\tau)\\ R_{X\hat{X}}(\tau)=\hat{R}_{X}(\tau)\\ R_{X}(\tau)=2[R_{X}(\tau)+j\hat{R}_{X}(\tau)] RX^X(τ)=R^X(τ)RXX^(τ)=R^X(τ)RX(τ)=2[RX(τ)+jR^X(τ)]
解析信号的频谱只有正频段且幅度值为原来的两倍,实现了信号由双边谱转换成单边谱。
解析信号的功率谱也只有正频段,强度为原来的四倍。
X~(ω)=X(ω)2u(ω)={2X(ω)ω&gt;0X(0)ω=00ω&lt;0\widetilde{X}(\omega)=X(\omega)2u(\omega)=\\ \{ \begin{array}{lcl} {2X(\omega)} &amp;\omega&gt;0\\ {X(0)} &amp;\omega=0\\ {0} &amp;\omega&lt;0 \end{array} X(ω)=X(ω)2u(ω)={2X(ω)X(0)0ω>0ω=0ω<0

欧拉公式(Euler‘s formula)的启发

eix=cos(x)+isin(x)e^{ix} = cos(x)+isin(x)eix=cos(x)+isin(x)
这个公式说明,用复指数信号可以表示成一个实数信号和一个虚数信号的和的形式。而且,这个实部和虚部是有关系的,一个是cos,一个是sin,两者相差pi/2,看sin和cos的傅里叶变换:

cos(ω0t)cos(\omega_0t)cos(ω0t)π[δ(ω+ω0)+δ(ω−ω0)]\pi[\delta(\omega+\omega_0)+\delta(\omega-\omega_0)]π[δ(ω+ω0)+δ(ωω0)]
sin(ω0t)sin(\omega_0t)sin(ω0t)jπ[δ(ω+ω0)−δ(ω−ω0)]j\pi[\delta(\omega+\omega_0)-\delta(\omega-\omega_0)]jπ[δ(ω+ω0)δ(ωω0)]

可以看出,在正频率上和负频率上两者的相位上的先后顺序刚好相反,但是都是保持90°的差值。
看到这里,大概可以理解Hilbert变换的用意了吧。欧拉公式实际上是一种特殊的,或者说,最简单的Hilbert变换。
复指数信号,就是等号左边的那个,频谱就是一个脉冲,而且是2πδ(ω−ω0)2\pi\delta(\omega-\omega_0)2πδ(ωω0)。只有正频率,且是两倍。虽然时域上是复数,但是在频域只有正分量,实际上是一种简化。

Hilbert解调原理

设窄带信号为 x(t)=a(t)cos(2πfst+ϕ(t))x(t)=a(t)cos(2\pi f_st+\phi(t))x(t)=a(t)cos(2πfst+ϕ(t))
其中fsf_sfs是载波频率,a(t)a(t)a(t)x(t)x(t)x(t)的包络,ϕ(t)\phi(t)ϕ(t)x(t)x(t)x(t)的相位调制信号。
由于x(t)x(t)x(t)是窄带信息,可设a(t)=[1+∑m=1Mxmcos(2πfmt+γm)]a(t)=[1+\sum^{M}_{m=1}x_mcos(2\pi f_mt+\gamma_m)]a(t)=[1+m=1Mxmcos(2πfmt+γm)]fmf_mfm为调幅信号a(t)a(t)a(t)的频率分量,γm\gamma_mγmfmf_mfm的各初相角。因此a(t)a(t)a(t)的各成分频率fm&lt;&lt;fsf_m&lt;&lt;f_sfm<<fs
x(t)x(t)x(t)进行Hilbert变换,并求解析信号,并使用欧拉公式表示:
z(t)=x(t)+x^(t)=A(t)ejΦ(t)z(t)=x(t)+\hat{x}(t)=A(t)e^{j\Phi(t)}z(t)=x(t)+x^(t)=A(t)ejΦ(t)
其中A(t)=a(t)=[1+∑m=1Mxmcos(2πfmt+γm)]A(t)=a(t)=[1+\sum^{M}_{m=1}x_mcos(2\pi f_mt+\gamma_m)]A(t)=a(t)=[1+m=1Mxmcos(2πfmt+γm)]Φ(t)=2πfst+ϕ(t)\Phi(t)=2\pi f_st+\phi(t)Φ(t)=2πfst+ϕ(t)
容易发现
a(t)=A(t)=x2(t)+x^2(t)a(t)=A(t)=\sqrt{x^2(t)+\hat{x}^2(t)}a(t)=A(t)=x2(t)+x^2(t)
ϕ(t)=Φ(t)−2πfst=actan(x(t)x^(t))−2πfst\phi(t)=\Phi(t)-2\pi f_st=actan(\frac{x(t)}{\hat{x}(t)})-2\pi f_stϕ(t)=Φ(t)2πfst=actan(x^(t)x(t))2πfst
并可以利用相位求解出解调频率:f(t)=12πdϕ(t)dt=12πdΦ(t)dt−fsf(t)=\frac{1}{2\pi}\frac{d\phi(t)}{dt}=\frac{1}{2\pi}\frac{d\Phi(t)}{dt}-f_sf(t)=2π1dtdϕ(t)=2π1dtdΦ(t)fs

希尔伯特变换的意义

首先,将实数信号变换成解析信号的结果就是,把一个一维的信号变成了二维复平面上的信号,复数的模和幅角代表了信号的幅度和相位,如图所示
在这里插入图片描述
这样看来,似乎复数信号才是完整的,而实信号只是在复平面的实轴上的一个投影。我们知道,解析信号可以计算包络(瞬时振幅)和瞬时相位。在上图中可以看到,实际上我们计算的包络就是黑色的线围成的立体图形的边界在实部的投影,而计算这个边的投影也很简单,就是在复平面上的螺旋线中的每一个点的模值,也就是A(t) = sqrt(x^2(t) + Hilbert(x(t))^2),而瞬时相位就是虚部(Hilbert变换后的)和实部(原始信号)在某一时间点的比值的arctan,瞬时频率就是它的导数。

评论 9
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值