帕塞瓦尔定理

在数学中，帕塞瓦尔定理（或称帕塞瓦尔等式），经常指“傅里叶转换是幺正算符”这一结论；简而言之，就是说函数平方的和（或积分）等于其傅里叶转换式平方之和（或者积分）。这个定理产生于法国数学家马克-安托万·帕塞瓦尔（Marc-Antoine Parseval）在1799年所得到的一个有关级数的定理，该定理随后被应用于傅里叶级数。它也被称为瑞利能量定理或瑞利恒等式，以物理学家瑞利命名。

虽说帕塞瓦尔定理这一术语常用来描述任何傅里叶转换的幺正性，尤其是在物理学和工程学上，但这种属性最一般的形式还是称为普朗歇尔定理而不是帕塞瓦尔定理才更合适。

该定理是勾股定理在希尔伯特空间或更广泛的内积空间中的推广，或者说勾股定理是帕塞瓦尔定理在定义了内积的二维欧氏空间中的特例。

## 帕塞瓦尔定理的陈述

在一般的欧氏平面几何中，勾股定理说明直角三角形的两个直角边之长度的平方加起来等于斜边的平方。从另一种角度来看，若在平面上定义了一个直角坐标系xOy（单位向量分别是$(e_x, e_y)$），那么一个向量和它在这两个坐标轴方向上的投影构成一个直角三角形，因此，向量的长度的平方等于它在两个坐标轴方向上的投影的长度的平方之和。

对于一个有限维的欧几里得空间$\mathbb{R}^n$ 以及其中的标准规范正交基$(e_1, e_2, \cdots , e_n)$，空间中的一个向量$v = (v_1, v_2, \cdots , v_n)$ 的长度的平方等于它在各个基向量上的投影的长度的平方之和：

$$\left\| v \right\|^2 = v_1^2 + v_2^2 + \cdots + v_n^2$$

在一般的希尔伯特空间之中，也有类似的等式。设$\mathcal{H}$ 是一个装备了内积：$\left \langle \cdot , \cdot \right \rangle$ 的希尔伯特空间。考虑$\mathcal{H}$ 中的一组规范正交基：$(e_1, e_2, \cdots , e_n, \cdots)$，那么$\mathcal{H}$ 中的每一个向量的范数的平方都等于它在各个基向量上的投影的平方之和。

$\sum_{k} \left| \left \langle x , e_k \right \rangle \right|^2 = \left \| x \right \|^2$

假定A(x)和B(x)都是平方可积的（参照勒贝格测度）复变函数，且定义在R上周期为2π的区间上，分别写成傅里叶级数的形式:

$$A(x)=\sum_{n=-\infty}^\infty a_ne^{inx}$$

和

$$B(x)=\sum_{n=-\infty}^\infty b_ne^{inx}.$$

然后

$$\sum_{n=-\infty}^\infty a_n\overline{b_n} = \frac{1}{2\pi} \int_{-\pi}^\pi A(x)\overline{B(x)} \, dx,$$

这里的i是虚数单位而上划线（horizontal bars）表示复共轭运算。

一般地，给定一个交换的拓扑群 G 和它的Pontryagin对偶 G^, 帕塞瓦尔定理 says the Pontryagin–Fourier transform is a unitary operator between Hilbert spaces L2(G) and L2(G^) (with integration being against the appropriately scaled Haar measures on the two groups.) When G is the unit circle T, G^ is the integers and this is the case discussed above. When G is the real line R, G^ is also R and the unitary transform is the Fourier transform on the real line. When G is the cyclic group Zn, again it is self-dual and the Pontryagin–Fourier transform is what is called discrete-time Fourier transform in applied contexts.

## 物理学和工程学上使用的记号

在 物理学 和 工程学 中, 帕塞瓦尔定理通常描述如下:

$$\int_{-\infty}^\infty | x(t) |^2 \, dt = \int_{-\infty}^\infty | X(f) |^2 \, df$$

其中$X(f) = \mathcal{F} \{ x(t) \}$ 为 x(t) 的连续傅立叶变换（以归一化酉形式），而f代表x的频率分量（非角频率）

帕塞瓦尔定理的此表达形式解释了波形x(t)依时间域t累积的总能量与该波形的傅立叶变换X(f)在频率域f累积的总能量相等。

对于离散时间信号，该理论表达式变换为：

$$\sum_{n=-\infty}^\infty | x[n] |^2 = \frac{1}{2\pi} \int_{-\pi}^\pi | X(e^{i\omega}) |^2 d\omega$$

其中，X为x的离散时间傅立叶变换(DTFT)，而$\omega$为x的角频率（度每样本）。

此外，对于离散傅立叶变换 (DFT)，表达式变换为：

$$\sum_{n=0}^{N-1} | x[n] |^2 = \frac{1}{N} \sum_{k=0}^{N-1} | X[k] |^2$$

其中，Xk为xn的DFT变换，变换前后样本长度皆为N。

## 证明

### 连续傅立叶变换(CTFT)的帕塞瓦尔定理

$$\int_{-\infty}^{\infty}|x(t)|^2dt$$

$$=\int_{-\infty}^{\infty}x(t)x^*(t)dt$$

$$=\int_{-\infty}^{\infty}x(t)[\int_{-\infty}^{\infty}X^*(f)e^{-j2{\pi}ft}df]dt$$

$$=\int_{-\infty}^{\infty}X^*(f)[\int_{-\infty}^{\infty}x(t)e^{-j2{\pi}ft}dt]df$$

$$=\int_{-\infty}^{\infty}X^*(f)X(f)df$$

$$=\int_{-\infty}^{\infty}|X(f)|^2df$$

其中，$x^*(t)$是$x(t)$的共轭复数。

### 离散时间傅立叶变换(DTFT)的帕塞瓦尔定理

$$\sum_{n=-\infty}^{\infty}|x[n]|^2$$

$$=\sum_{n=-\infty}^{\infty}x[n]x^*[n]$$

$$=\sum_{n=-\infty}^{\infty}x[n][\frac{1}{2{\pi}}\int_{0}^{2{\pi}}X^*(e^{j\omega})e^{-j{\omega}n}d\omega]$$

$$=\frac{1}{2{\pi}}\int_{0}^{2{\pi}}[\sum_{n=-\infty}^{\infty}x[n]e^{-j{\omega}n}]X^*(e^{j\omega})d\omega$$

$$=\frac{1}{2{\pi}}\int_{0}^{2{\pi}}X(e^{j\omega})X^*(e^{j\omega})d\omega$$

$$=\frac{1}{2{\pi}}\int_{0}^{2{\pi}}|X(e^{j\omega})|^2d\omega$$

其中，$x^*[n]$是$x[n]$的共轭复数。

### 连续时间傅立叶级数(CTFS)的帕塞瓦尔定理

令x(t)是周期为$T_0=\frac{1}{f_0}$的连续时间函数。

$c_n$是其连续时间傅立叶级数。$c_n=\frac{1}{T_0}\int_{0}^{T_0}x(t)e^{-j2\pi{n}{f_0}t}dt$

$$\sum_{n=-\infty}^{\infty}|c_n|^2$$

$$=\sum_{n=-\infty}^{\infty}c_n{c_n}^*$$

$$=\sum_{n=-\infty}^{\infty}c_n[\frac{1}{T_0}\int_{0}^{T_0}x^*(t)e^{j2\pi{n}{f_0}t}dt]$$

$$=\frac{1}{T_0}\int_{0}^{T_0}x^*(t)[\sum_{n=-\infty}^{\infty}c_ne^{j2\pi{n}{f_0}t}]dt$$

$$=\frac{1}{T_0}\int_{0}^{T_0}x^*(t)x(t)dt$$

$$=\frac{1}{T_0}\int_{0}^{T_0}|x(t)|^2dt$$

### 离散时间傅里叶级数(DTFS)的帕塞瓦尔定理

xn是长度为N的离散时间信号，$a_k$为其离散时间傅立叶级数，亦即$a_k=\frac{1}{N}\sum_{n=0}^{N-1}x[n]e^{-j\omega_0kn}$。

其中$\omega_0$是角基频，$\omega_0=\frac{2\pi}{N}$。

$$\sum_{k=0}^{N-1}|a_k|^2$$

$$=\sum_{k=0}^{N-1}a_k{a_k}^*$$

$$=\sum_{k=0}^{N-1}a_k[\frac{1}{N}\sum_{n=0}^{N-1}x^*[n]e^{j\omega_0kn}]$$

$$=\frac{1}{N}\sum_{n=0}^{N-1}x^*[n][\sum_{k=0}^{N-1}a_ke^{j\omega_0kn}]$$

$$=\frac{1}{N}\sum_{n=0}^{N-1}x^*[n]x[n]$$

$$=\frac{1}{N}\sum_{n=0}^{N-1}|x[n]|^2$$

### 离散傅立叶变换(DFT)的帕塞瓦尔定理

令$x[n]$为一长度是N点的离散时间信号，仅在0≤n≤N-1有值，$x[n]=0$ for $n<0$ or $n>N-1$。

其DFT为$X[k]$，亦为一长度是N点的离散时间信号，仅在0≤k≤N-1有值，$X[k]=0$ for $k<0$ or $k>N-1$。

设$W_N=e^{j\frac{2\pi}{N}}$。

$$\sum_{n=0}^{N-1}|x[n]|^2$$

$$=\sum_{n=0}^{N-1}x[n]x^*[n]$$

$$=\sum_{n=0}^{N-1}x[n][\frac{1}{N}\sum_{k=0}^{N-1}X^*[k]{W_N}^{-kn}]$$

$$=\frac{1}{N}\sum_{k=0}^{N-1}X^*[k][\sum_{n=0}^{N-1}x[n]{W_N}^{-kn}]$$

$$=\frac{1}{N}\sum_{k=0}^{N-1}X^*[k]X[k]$$

$$=\frac{1}{N}\sum_{k=0}^{N-1}|X[k]|^2$$