每日精学 141｜DFT的Parseval能量恒等式

2026-09-25

雷达 · 通信 · 电子战　/　基础精学 141

非归一化正DFT的频域能量是时域能量的N倍，因而功率、窗增益和FFT缩放的计量必须遵循一致的变换约定。

01
模型与符号

有限长度N的任意复序列，正DFT不缩放，逆DFT使用1/N；不引入随机过程极限。

E_t：时域能量（幅度平方）

E_f：频域能量（幅度平方）

N：变换长度（sample）

X[k]：非归一化DFT（幅度谱）

02
从模型到公式

展开平方和

把DFT定义代入可得

$\sum_k|X[k]|^2=\sum_n\sum_mx[n]x^*[m]\sum_ke^{-\mathrm j2\pi k(n-m)/N}.$

交换求和合法，因为所有和均有限。

正交消交叉项

单位根几何级数满足

$\sum_{k=0}^{N-1}e^{-\mathrm j2\pi k(n-m)/N}=N\delta_{nm},$

其中 0≤ n,m<N。故

$\sum_n|x[n]|^2=\frac1N\sum_k|X[k]|^2.$

缩放约定

平均样本功率为

$P=\frac{E_t}{N}=\frac{1}{N^2}\sum_k|X[k]|^2.$

若采用单位酉变换

$X_u=X/\sqrt N$

则

$\sum|X_u|^2=E_t$

窗化数据必须以实际 w[n]x[n] 的能量比较，不能用原始未窗化能量替代。

一般DFT归一化与单边能量

若正变换额外乘系数 a，频域能量因子随之改变。由复指数正交性得到

$X_a[k]=a\sum_{n=0}^{N-1}x[n]e^{-2\pi\mathrm jkn/N},\qquad \sum_k|X_a[k]|^2=N|a|^2\sum_n|x[n]|^2.$

因此未归一化 DFT 使用 1/N 补偿，单位酉形式使用

$a=1/\sqrt N$

对于偶数长度实序列，单边能量应写为

$N\sum_n|x[n]|^2=|X[0]|^2+|X[N/2]|^2+2\sum_{k=1}^{N/2-1}|X[k]|^2.$

直流与 Nyquist 是自共轭频点，不能重复加倍。

该恒等式针对有限记录本身；窗函数改变了记录能量，不能未经窗增益校正就把窗后结果当作原记录的物理功率。

03
物理含义

DFT是在正交坐标系间改写有限向量，交叉项抵消后只剩坐标能量；额外N因子完全来自基向量未归一化。

04
固定参数算例

对复随机序列、单音和窗化数据分别验证Parseval以及单位酉形式；另展示错漏1/N时出现的固定倍率误差。

固定参数数值验证结果

量	数值
length	128
complex record energy	247.734305127
unit tone energy	128
windowed energy	88.2858451001
parseval max error	5.68434188608e-14
wrong scaling ratio	128
样本功率	1.93542425881

图 1　DFT的Parseval能量恒等式的理论关系与固定参数验证

表中数值和图中曲线由 MATLAB 实际运行得到，随机种子为 20261060；运行版本与完整结果另存为独立文件。

05
数字实现的边界

RTL实现复样本能量 a²+b²，a、b为整数实虚部，64位输出保留平方和；跨N点积累及频域归一化由外部完成。 RTL 采用同步有效信号和异步低有效复位；

空拍保持状态且输出无效。除另有说明外，输入为16位有符号整数，结果为64位有符号整数，测试范围内不截断或饱和。输出在采样时钟上升沿更新。自检覆盖连续样本、空拍和复位；

未进行器件布局布线、时序收敛或功耗验证。

独立自检程序已在 Icarus Verilog 13.0 中实际运行并通过。该结果不包含器件综合、静态时序分析或板级测试。

06
适用范围

Parseval恒等式可用于检查FFT功率测量、窗化能量和频域滤波器的缩放约定。时域与频域两条独立计算路径能够暴露漏乘N、重复单边加倍或错误窗归一化等问题。该恒等式检验能量一致性，但不单独证明频点索引正确。

07
工程上需要注意

单边谱需要对非直流非Nyquist频点加倍。

功率谱密度还涉及采样率和窗能量，不仅是FFT长度。

有限精度FFT和累加器会产生小残差。

非归一化正DFT的频域能量是时域能量的N倍，因而功率、窗增益和FFT缩放的计量必须遵循一致的变换约定。

内容依据：每日精学第 141 课《DFT的Parseval能量恒等式》。图表沿用原课固定参数数值结果，属于数值验证，不代表设备实测。

阅读8

写评论...

每日精学 141｜DFT的Parseval能量恒等式

01模型与符号

02从模型到公式

展开平方和

正交消交叉项

缩放约定

一般DFT归一化与单边能量

03物理含义

04固定参数算例

05数字实现的边界

06适用范围