Choi-Williams分布详解

Choi-Williams分布详解

引言

时频分析是信号处理中的重要领域，旨在同时描述信号在时间和频率域的特性。传统的傅里叶变换虽然能将信号从时域转换到频域，但无法同时提供时间和频率的局部信息。而Choi-Williams分布（Choi-Williams Distribution，简称CWD）作为Cohen类时频分布的一种重要形式，提供了一种有效的方法来分析非平稳信号的时频特性，同时抑制交叉项的影响。

时频分析基础

在深入了解Choi-Williams分布之前，我们需要先理解时频分析的基本概念。对于非平稳信号，其频率内容随时间变化，这使得仅使用傅里叶变换不足以完全描述信号的特性。时频分析的目标是创建一个联合时频分布，即时频平面上的二维函数，用于表示信号能量如何在时间和频率上分布。常见的时频分析方法包括短时傅里叶变换(STFT)、Wigner-Ville分布(WVD)以及各种Cohen类分布。在这些方法中，Wigner-Ville分布提供了良好的时频分辨率，但存在严重的交叉项干扰，而Choi-Williams分布则是为了解决这一问题而提出的。

Choi-Williams分布的定义与理论基础

Choi-Williams分布是Cohen类分布的一种特殊形式，通过引入指数核函数来抑制交叉项，同时保留自项的时频特性。从理论上讲，任何满足一定条件的Cohen类二次时频分布可以通过以下一般形式表示：

$C_x(t,\omega )=\frac{1}{4{\pi }^2}{\int }_{-\infty }^{\infty }{\int }_{-\infty }^{\infty }{\int }_{-\infty }^{\infty }{e}^{j\theta (u-t)-j\omega \tau }\cdot \varphi (\theta ,\tau )\cdot x(u+\frac{\tau }{2})\cdot x^{\ast }(u-\frac{\tau }{2})d\theta d\tau du$

其中$\varphi (\theta ,\tau )$是二维核函数，定义了分布的特性。通过引入模糊函数（ambiguity function），可以将上式重写为：

$C_x(t,\omega )=\frac{1}{4{\pi }^2}{\int }_{-\infty }^{\infty }{\int }_{-\infty }^{\infty }\varphi (\theta ,\tau )\cdot e^{j\theta t-j\omega \tau }\cdot A_x(\theta ,\tau )d\theta d\tau$

Choi-Williams分布则特别选择了指数形式的核函数，其完整定义如下：

$CW{D}_x(t,\omega )={\int }_{-\infty }^{\infty }{\int }_{-\infty }^{\infty }\Phi (\theta ,\tau )\cdot e^{j\theta t-j\omega \tau }\cdot A_x(\theta ,\tau )d\theta d\tau$

其中，$A_x(\theta ,\tau )$是信号$x(t)$的模糊函数，定义为：

$A_x(\theta ,\tau )={\int }_{-\infty }^{\infty }x(t+\frac{\tau }{2})\cdot x^{\ast }(t-\frac{\tau }{2})\cdot e^{-j\theta t}dt$

而$\Phi (\theta ,\tau )$是Choi-Williams核函数，定义为：

$\Phi (\theta ,\tau )=e^{-\frac{{\theta }^2{\tau }^2}{\sigma }}$

其中$\sigma$是一个正的参数，用于控制核函数的形状和交叉项抑制的程度。

从另一个角度看，Choi-Williams分布也可以通过信号的Wigner-Ville分布与时频核的二维卷积来表示：

$CW{D}_x(t,\omega )=WV{D}_x(t,\omega )\ast \ast G(t,\omega )$

其中$WV{D}_x(t,\omega )$是信号$x(t)$的Wigner-Ville分布，$G(t,\omega )$是$\Phi (\theta ,\tau )$的二维傅里叶变换，$\ast \ast$表示二维卷积。$G(t,\omega )$可以通过以下公式计算：

$G(t,\omega )=\frac{1}{4{\pi }^2}{\int }_{-\infty }^{\infty }{\int }_{-\infty }^{\infty }{e}^{j\theta t-j\omega \tau }\cdot \Phi (\theta ,\tau )d\theta d\tau =\frac{1}{4{\pi }^2}{\int }_{-\infty }^{\infty }{\int }_{-\infty }^{\infty }$