随机变量的变换的Hilbert变换与Kramers-Kronig关系
1. 从信号分析的一个基本问题引入
在信号处理和系统理论中,我们经常需要分析信号的幅度和相位,或者实部与虚部之间的关系。对于一个实值信号(可视为随机过程的一次实现),其频域表示(如傅里叶变换)的实部和虚部是独立的吗?实际上,对于一大类因果信号(或更一般地,具有解析性质的信号),其实部和虚部通过一种称为希尔伯特变换(Hilbert transform) 的算子紧密联系。这种关系在概率论与统计中,特别是在研究随机过程的解析性质和谱表示时,有重要应用。
2. 希尔伯特变换的定义
设 \(f(t)\) 是一个实值函数(或随机过程的一条样本路径),其希尔伯特变换 \(\mathcal{H}[f](t)\) 定义为:
\[\mathcal{H}[f](t) = \frac{1}{\pi} \, \text{p.v.} \int_{-\infty}^{\infty} \frac{f(\tau)}{t - \tau} \, d\tau \]
其中 p.v. 表示柯西主值积分。直观上,这可以看作 \(f\) 与核函数 \(1/(\pi t)\) 的卷积,该核函数是所有频率的 −90 度相移器。
关键性质:希尔伯特变换是一个线性算子,且对函数进行两次变换得到 \(\mathcal{H}^2[f] = -f\)。
3. 频域解释与解析信号
在频域中,希尔伯特变换对应一个简单的乘法操作。设 \(F(\omega)\) 是 \(f(t)\) 的傅里叶变换,则有:
\[\mathcal{F}\{ \mathcal{H}[f] \}(\omega) = -i \, \text{sgn}(\omega) F(\omega) \]
其中 \(\text{sgn}(\cdot)\) 是符号函数。这说明希尔伯特变换在频域中相当于对正频率乘以 \(-i\)(相位延迟 \(\pi/2\)),对负频率乘以 \(i\)(相位超前 \(\pi/2\))。
利用这个性质,我们可以构造解析信号(analytic signal):
\[f_a(t) = f(t) + i \mathcal{H}[f](t) \]
其傅里叶变换在负频率处为零,正频率处为 \(2F(\omega)\)。解析信号是研究信号包络和瞬时频率的基础工具。
4. 随机过程的希尔伯特变换
当 \(f(t)\) 是一个零均值的宽平稳随机过程 \(X(t)\) 时,其希尔伯特变换 \(Y(t) = \mathcal{H}[X](t)\) 也是一个零均值的宽平稳随机过程,且与原过程联合平稳。
它们的相关函数满足:
\[R_{YY}(\tau) = R_{XX}(\tau), \quad R_{XY}(\tau) = \mathcal{H}[R_{XX}](\tau) \]
其中 \(R_{XY}\) 是互相关函数。
重要结论:希尔伯特变换保持功率谱密度不变,但改变了相位谱,使得 \(X(t)\) 和 \(Y(t)\) 在同一时刻正交(即 \(R_{XY}(0)=0\))。
5. 克拉默斯-克勒尼希(Kramers-Kronig)关系
在物理光学和系统响应理论中,系统的复值响应函数(如介电常数)的实部与虚部之间存在一种积分关系,称为克拉默斯-克勒尼希关系。本质上,它是希尔伯特变换在频域中的体现。
设 \(\chi(\omega) = \chi_1(\omega) + i\chi_2(\omega)\) 是因果线性系统的复响应函数(即其逆傅里叶变换是因果信号),则其实部和虚部满足:
\[\chi_1(\omega) = \frac{1}{\pi} \, \text{p.v.} \int_{-\infty}^{\infty} \frac{\chi_2(\omega')}{\omega' - \omega} d\omega', \quad \chi_2(\omega) = -\frac{1}{\pi} \, \text{p.v.} \int_{-\infty}^{\infty} \frac{\chi_1(\omega')}{\omega' - \omega} d\omega' \]
这正是希尔伯特变换对在频域中的表现。它说明,对于因果系统,实部和虚部不是独立的,知道其中一个在全部频率上的值就能确定另一个。
6. 在概率与统计中的应用
- 解析随机过程:由宽平稳过程 \(X(t)\) 构造的解析过程 \(Z(t) = X(t) + i Y(t)\) 是复值平稳过程,其谱表示仅包含正频率,便于分析包络和相位。
- 包络与瞬时频率:在信号检测和通信中,解析信号的模 \(|Z(t)|\) 称为 \(X(t)\) 的包络,其相位的时间导数为瞬时频率,常用于分析窄带随机过程(如调频噪声)。
- 谱分析:克拉默斯-克勒尼希关系允许从系统响应的损耗部分(虚部)推算其色散部分(实部),在统计物理和计量经济学中用于构建满足因果性的谱模型。
- 希尔伯特-黄变换:在经验模态分解中,希尔伯特变换用于计算瞬时频率,是非平稳随机信号时频分析的工具之一。
7. 数学的深层联系
希尔伯特变换和克拉默斯-克勒尼希关系本质上是柯西积分公式和解析函数边界性质的体现。如果函数在上半复平面解析且衰减适当,其实部与虚部在实轴上的限制就构成一对希尔伯特变换。在概率论中,这联系到随机过程的解析延拓和谱表示定理,是调和分析与随机过程交叉的经典范例。