每日精学 142|实序列DFT的共轭对称

2026-09-25

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

实序列DFT的负频率是正频率的复共轭,偶数长度数据只需保存直流至Nyquist的N/2+1个复数频点即可无损恢复。

01
模型与符号

时域输入严格为实数,N为偶数,频率索引按模N解释。

X[k]:DFT系数(幅度谱)

N:偶数长度(sample)

k:频率索引(bin)

P1[k]:单边能量贡献(功率)

02
从模型到公式

共轭关系

因 x[n]=x*[n],取DFT共轭得到

X^*[k]=\sum_nx[n]e^{+\mathrm j2\pi kn/N}=X[(-k)\bmod N].

所以 X[N−k]=X*[k]。

特殊频点

k=0与k=N/2各自映射到自身,因此

\operatorname{Im}X[0]=0,\qquad\operatorname{Im}X[N/2]=0.

其余频点按共轭成对,所需独立存储数量为 N/2+1。

单边能量

利用配对与Parseval,时域能量为

E_t=\frac1N\left(|X[0]|^2+|X[N/2]|^2+2\sum_{k=1}^{N/2-1}|X[k]|^2\right).

只对内部频点乘2,不能把直流和Nyquist重复计数。

独立自由度与自共轭端点

实序列的负频率由正频率共轭确定。偶数 N 时,直流及 Nyquist 两个点必须为实数,其余 N/2−1 个正频率点各含两个实自由度,因此

2+2(N/2-1)=N.

存储半谱不是丢弃一半信息,而是利用已知实值约束消除重复。

若独立修改某个正频率点,负频率对应点也必须作共轭修改:

\widetilde X[N-k]=\widetilde X[k]^*.

否则逆变换一般产生非零虚部。

奇数长度没有独立的 Nyquist 点,半谱端点处理需相应改变。对原本为复数的输入,上述压缩不成立;仅观察某次记录的虚部很小也不足以作为普遍的实信号假设。


03
物理含义

共轭对称来自实数序列的代数约束,而不是频谱幅度碰巧对称。虚部符号是恢复时域实数的必要信息。

04
固定参数算例

对N=128的实多音加噪序列,检查全频共轭残差,用单边数据重建时域并核对单边能量。

固定参数数值验证结果
量数值
length128
independent bins65
conjugate max error0
reconstruction max error6.66133814775e-16
endpoint imag max0
single sided energy error0
时域能量107.709986993
图 1:实序列DFT的共轭对称的理论关系与固定参数验证图 1:实序列DFT的共轭对称的理论关系与固定参数验证

图 1 实序列DFT的共轭对称的理论关系与固定参数验证

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

05
数字实现的边界

RTL实现复共轭打包子核,输入a、b为实虚部,输出a与-b;负满量程取反在64位输出中精确保留,不回绕。 RTL 采用同步有效信号和异步低有效复位;空拍保持状态且输出无效。

除另有说明外,输入为16位有符号整数,结果为64位有符号整数,测试范围内不截断或饱和。输出在采样时钟上升沿更新。自检覆盖连续样本、空拍和复位;未进行器件布局布线、时序收敛或功耗验证。


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

06
适用范围

实序列共轭对称可减少实ADC记录的频谱存储,并用于单边谱显示与数据完整性检查。重建负频率时必须恢复虚部相反符号,同时保留直流和Nyquist的独立实值。对复I/Q记录不应使用这一压缩假设。

07
工程上需要注意

复基带I/Q数据通常不满足实序列共轭对称。

奇数N不存在独立Nyquist点,计数公式需要改变。

仅保存幅度而丢弃相位无法重建原序列。

实序列DFT的负频率是正频率的复共轭,偶数长度数据只需保存直流至Nyquist的N/2+1个复数频点即可无损恢复。

内容依据:每日精学第 142 课《实序列DFT的共轭对称》。图表沿用原课固定参数数值结果,属于数值验证,不代表设备实测。

阅读8
分享
写评论...