S变换（S-Transform）的数学理论详解

S变换（S-Transform）的数学理论详解

引言

S变换（S-Transform）是由加拿大学者Stockwell等人于1996年提出的一种时频分析方法，它巧妙地结合了短时傅里叶变换（STFT）的相位特性和连续小波变换（CWT）的多分辨率分析能力。S变换既保留了与傅里叶变换直接联系的绝对相位信息，又具备小波变换随频率变化的分辨率特性，因此得到了广泛应用。

基本原理

S变换可以看作是短时傅里叶变换的一种扩展，其核心思想是使用高斯窗函数作为基本窗口，但窗口宽度会随着频率自适应调整。低频成分使用宽窗口以获得更好的频率分辨率，高频成分则使用窄窗口以获得更好的时间分辨率。这种自适应性使S变换非常适合分析非平稳信号，尤其是那些同时包含高频和低频成分的复杂信号。

从概念上讲，S变换融合了两个重要的信号分析工具：短时傅里叶变换提供了与时间有关的频率信息，小波变换则提供了多尺度分析能力。S变换可以看作是"相位修正的连续小波变换"，也可以看作是"频率依赖的短时傅里叶变换"。

一、基本定义与数学表达

S变换作为一种时频分析工具，其严格的数学定义源于连续小波变换与短时傅里叶变换的理论基础。

连续S变换的定义

对于一维时间信号 $h(t)$，其连续S变换定义为：

$$S(\tau ,f)={\int }_{-\infty }^{\infty }h(t)\cdot \frac{|f|}{\sqrt{2\pi }}{e}^{-\frac{(t-\tau {)}^2{f}^2}{2}}\cdot e^{-i2\pi ft}dt$$

这里 $\tau$ 表示时间定位参数，$f$ 表示频率参数，指数函数 $e^{-\frac{(t-\tau {)}^2{f}^2}{2}}$ 是频率依赖的高斯窗函数，系数 $\frac{|f|}{\sqrt{2\pi }}$ 确保窗函数的能量为常数。更紧凑的形式可表示为：

$$S(\tau ,f)={\int }_{-\infty }^{\infty }h(t)\cdot w(t-\tau ,f)\cdot e^{-i2\pi ft}dt$$

其中窗函数为：

$$w(t-\tau ,f)=\frac{|f|}{\sqrt{2\pi }}{e}^{-\frac{(t-\tau {)}^2{f}^2}{2}}$$

离散S变换的推导

对于离散时间信号 $h[kT]$ (其中 $k=0,1,...,N-1$，$T$ 是采样间隔)，首先计算其离散傅里叶变换：

$$H[n]=\frac{1}{N}\sum _{k=0}^{N-1}h[kT]\cdot e^{-i2\pi nk/N}$$

对于 $n\ne 0$（非零频率），S变换可以表示为：

$$S[jT,\frac{n}{NT}]=\sum _{m=0}^{N-1}H\left[\frac{(m+n)\mathrm{mod}N}{N}\right]\cdot e^{-2{\pi }^2{m}^2/n^2}\cdot e^{i2\pi mj/N}$$

对于 $n=0$（零频率），S变换退化为信号的直流分量：

$$S[jT,0]=\frac{1}{N}\sum _{k=0}^{N-1}h[kT]=H[0]$$

二、S变换的数学性质

2.1 傅里叶变换性质

S变换与傅里叶变换之间存在严格的数学联系。对S变换结果在时间上积分，可得：

$${\int }_{-\infty }^{\infty }S(\tau ,f)d\tau =H(f)={\int }_{-\infty }^{\infty }h(t)e^{-i2\pi ft}dt$$

证明过程涉及窗函数的性质：

$${\int }_{-\infty }^{\infty }\frac{|f|}{\sqrt{2\pi }}{e}^{-\frac{(t-\tau {)}^2{f}^2}{2}}d\tau =1$$

将此代入S变换定义：

$$\begin{array}{rl}{\int }_{-\infty }^{\infty }S(\tau ,f)d\tau & ={\int }_{-\infty }^{\infty }{\int }_{-\infty }^{\infty }h(t)\cdot \frac{|f|}{\sqrt{2\pi }}{e}^{-\frac{(t-\tau {)}^2{f}^2}{2}}\cdot e^{-i2\pi ft}dtd\tau \\ & ={\int }_{-\infty }^{\infty }h(t)\cdot e^{-i2\pi ft}\left({\int }_{-\infty }^{\infty }\frac{|f|}{\sqrt{2\pi }}{e}^{-\frac{(t-\tau {)}^2{f}^2}{2}}d\tau \right)dt\\ & ={\int }_{-\infty }^{\infty }h(t)\cdot e^{-i2\pi ft}dt\end{array}$$