希尔伯特变换（Hilbert Transform）简介及其物理意义-程序员宅基地

技术标签：信号处理

希尔伯特变换（Hilbert Transform）简介及其物理意义

Hilbert变换简介

希尔伯特变换是信号处理中的一种常用手段，数学定义如下：

这里写图片描述

与卷积的概念进行对比，可以发现，上面的Hilbert变换的表达式实际上就是将原始信号和一个信号做卷积的结果。这个用来卷积的信号就是

$\frac{1}{\pi t}$

因此，Hilbert变换可以看成是将原始信号通过一个滤波器，或者一个系统，这个系统的冲击响应为h(t)。

这里写图片描述

对h(t)做傅里叶变换，可以得到：

$H(j\omega) = -j sgn(\omega)$

或者写成：

这里写图片描述

sgn()是符号函数。从频谱上来看，这个滤波器将我们的原始信号的正频率部分乘以-j，也就是说，保持幅度不变的条件下，将相位移动了-pi/2，而对于负频率成分，移动了pi/2。

下面这个示意图很直观地表示了Hilbert变换，在这里我画出了对原始信号做1到4次Hilbert变换的频谱示意图，是为了说明Hilbert变换的几个性质：

这里写图片描述

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

如果一个信号是两个信号的卷积，即 y = conv(v,x) ，那么Hilbert(y) = conv(Hilbert(v),x) = conv(v,Hilbert(x))

这个性质，只要意识到Hilbert变换本质上是卷积就可以明白。

x(t) 和 Hilbert(x(t))的能量以及平均功率相等，相关函数和功率谱相同。

Hilbert变换物理意义

解析过程概念

对于一个实随机过程 $x (t)$ ， $\hat{x}(t) = Hilbert(x(t))$ 为它的希尔伯特变换，那么定义复随机过程：

$\tilde{x}(t) = x(t) + \hat{x}(t)$

为x(t)的解析过程。

这个过程有如下特点，首先，实部和虚部功率谱相同，自相关函数相同；另外，实部和虚部的互相关函数是一个奇函数。其他的还有：

$R_{\hat X X}(\tau) = -\hat{R}_X(\tau)$

$R_{X\hat X}(\tau) = \hat{R}_X(\tau)$

$R_{\tilde{X}}(\tau) = 2[R_X(\tau) + j \hat{R}_X(\tau)]$

以及一个最重要的特点，就是解析信号的功率谱只有正频段，强度为原来的四倍。或者说是只有正频段且幅度值为原来的两倍：

这里写图片描述

欧拉公式（Euler‘s formula）的启发

欧拉一生中有公式和定理，这里说的欧拉公式指的是这个：

这里写图片描述

这个公式说明，用复指数信号可以表示成一个实数信号和一个虚数信号的和的形式。而且，这个实部和虚部是有关系的，一个是cos，一个是sin，两者相差pi/2，看sin和cos的傅里叶变换：

![这里写图片描述](https://imgconvert.csdnimg.cn/aHR0cDovL2ltZy5ibG9nLmNzZG4ubmV0LzIwMTgwMjI1MDEyOTAxMzY?x-oss-process=image/format,png)

可以看出，在正频率上和负频率上两者的相位上的先后顺序刚好相反，但是都是保持90°的差值。

看到这里，大概可以理解Hilbert变换的用意了吧。欧拉公式实际上是一种特殊的，或者说，最简单的Hilbert变换。

复指数信号，就是等号左边的那个，频谱就是一个脉冲，而且是 $2\pi\delta(\omega - \omega_0)$ 。只有正频率，且是两倍。虽然时域上是复数，但是在频域只有正分量，实际上是一种简化。

希尔伯特变换的意义

首先，将实数信号变换成解析信号的结果就是，把一个一维的信号变成了二维复平面上的信号，复数的模和幅角代表了信号的幅度和相位，如图所示：（图片来源：https://zhuanlan.zhihu.com/p/25250010）

这里写图片描述

这样看来，似乎复数信号才是完整的，而实信号只是在复平面的实轴上的一个投影。我们知道，解析信号可以计算包络（瞬时振幅）和瞬时相位。在上图中可以看到，实际上我们计算的包络就是黑色的线围成的立体图形的边界在实部的投影，而计算这个边的投影也很简单，就是在复平面上的螺旋线中的每一个点的模值，也就是A(t) = sqrt(x^2(t) + Hilbert(x(t))^2)，而瞬时相位就是虚部（Hilbert变换后的）和实部（原始信号）在某一时间点的比值的arctan，瞬时频率就是它的导数。

2018年02月25日01:29:30

本文链接：https://blog.csdn.net/edogawachia/article/details/79366444

原作者删帖不实内容删帖广告或垃圾文章投诉

智能推荐

FPGA设计经验之图像处理_基于fpga的数字图像处理-程序员宅基地

文章浏览阅读8.1k次，点赞37次，收藏278次。今天和大侠简单聊一聊基于FPGA的图像处理，之前也和各位大侠聊过相关的图像处理，这里面也超链接了几篇，具体如下：图像边缘检测算法体验步骤（Photoshop，Matlab）算法系列：基于 FPGA 的图像边缘检测系统设计（sobel算法） FPGA设计中 Verilog HDL实现基本的图像滤波处理仿真基于FPGA的类脑计算平台 —PYNQ 集群的无监督图像识别类脑计算..._基于fpga的数字图像处理

软件质量保证的过程_软件质量保证流程是怎么样的?-程序员宅基地

文章浏览阅读1.8k次。1、软件测试的定义的两面性：2、软件质量的定义：是人们实践产物的属性和行为，是可以认识以及科学描述的，并且通过一些方法和人类的活动来改进的。3、软件质量的模型：McCall模型、Boehm模型、ISO9126模型4、软件缺陷的产生：（1）技术问题：算法错误，语法错误，接口问题（2）团队问题：沟通不充分，误解（3）软件本身：文档错误，用户使用场合不正确；时间上不协调、或不一致；系统自我恢复或数据的异地备份、灾难性恢复等问题。缺陷构成：规格说明书>设计>代码在真正的程序测试之前，_软件质量保证流程是怎么样的?

Debian下配置IPV6和静态路由_debian 配置ipv6-程序员宅基地

文章浏览阅读8.2k次，点赞2次，收藏7次。一、编辑网卡文件vi /etc/network/interfaces二、修改网卡文件# This file describes the network interfaces available on your system# and how to activate them. For more information, see interfaces(5).# The loop..._debian 配置ipv6

hsql 初试之 exception-程序员宅基地

文章浏览阅读130次。hsql 是一个非常轻量的数据库，大小只有六百多K。我看中的是他哪 In-Process 模式，可能把数据内嵌在 java 应用里边，不用看成一个服务器。下载了最新版本 hsqldb_1_8_0_10, 里边的 hsql.jar 可以直接用的了..但我写了一个test，发觉运行时出下边的 exception :[quote]java.lang.NoClassDefFoundErr..._hsql exception

联合概率条件概率边缘概率之间关系_条件=边缘/联合-程序员宅基地

文章浏览阅读367次。联合概率、边缘概率与条件概率之间的关系P(X=a|Y=b)=P(X=a,Y=b)P(Y=b)https://blog.csdn.net/tick_tock97/article/details/79885868_条件=边缘/联合

快递鸟-各种电子面单_Api接口_快递鸟api面单-程序员宅基地

文章浏览阅读3.6k次。电子面单是一种通过热敏纸打印输出纸质物流面单的物流服务。通过热感应显示文字，打印速度比传统针式打印速度提升4~6倍。电子面单以接口形式嵌入到自己的系统、网站上，可以在自己的平台操作打印电子面单。一.电子面单接口类型及定义 1.快递电子面单接口：快递公司自己开发的电子面单服务，商家使用必须快递公司上门做系统对接，使用一家快递则需要对接一次。 2.菜鸟电子面单接口：可一次对接1_快递鸟api面单

随便推点

如何修改antd里Select的样式？？？？？_antd ant-select-selector-程序员宅基地

文章浏览阅读1.9w次，点赞15次，收藏33次。如何修改antd里Select的样式？？？？？antd的Selectantd的Select附一个官网链接 antd Select官方给出的框架用起来是很方便，但是美观程度差强人意，给出的API里虽有className属性，但是这是进行修改Option的样式那如何修改Select的样式呢？废话不多说，直接上代码 :global{ .ant-select-selection{ // antd 版本3.* // 你想要的样式 } .ant-select-selector{ /_antd ant-select-selector