每日精学 005｜有限长序列的离散时间傅里叶变换

2026-09-25

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

离散时间 LTI 系统的频率响应描述系统对单个复指数的复增益，但要分析一般有限记录，还需把样本序列分解到连续的离散时间频率轴上。DTFT 对每个归一化角频率计算序列与复指数基函数的相关和，所得复数同时包含该频率处的幅度与相位信息。

有限长序列出现在脉冲响应测量、通信基带突发、雷达回波数据窗和数字中频缓存中。其 DTFT 既可直接描述有限冲激响应滤波器，也可解释有限观测造成的谱展宽，并为离散傅里叶变换提供连续频率参照。

本文只处理有限支撑序列的 DTFT 定义及其由定义直接得到的基本性质，不展开随机过程功率谱、快速傅里叶变换或具体系统检测方案。

01
模型与符号

设复序列 x[n] 仅在整数区间 na≤ n≤ nb 内可能非零。主要结论采用以下假设：

  • 时间索引为整数，归一化角频率 ω 的单位为弧度每样本；

  • 序列具有有限支撑，故 DTFT 是有限项复指数之和，不需要额外收敛条件；

  • 反变换积分取任意长度为 2π 的完整频率区间，本文固定为 [−π,π]；

  • 共轭对称结论要求 x[n] 为实序列，能量关系允许复序列；

  • 数字实现按连续到达的 8 个有效样本划分数据块，输入采用二进制补码有符号整数。

n,m：离散时间索引与反变换索引（ℤ）

na,nb：有限序列支撑的起点与终点（ℤ）

ω：归一化角频率，单位为弧度每样本（实数）

x[n]：有限长离散时间序列（实数或复数）

X(ejω)：x[n] 的 DTFT（复数）

N：矩形序列或数据块的样本数（正整数）

n0：整数时移量，正值表示向右延迟（ℤ）

02
有限长序列的 DTFT 定义

正变换与连续频率含义

有限长序列的 DTFT 定义为

\boxed{ X(e^{\mathrm{j}\omega}) =\sum_{n=-\infty}^{\infty}x[n]e^{-\mathrm{j}\omega n} =\sum_{n=n_a}^{n_b}x[n]e^{-\mathrm{j}\omega n} }.

求和只含 nb−na+1 项，因此 X(ejω) 对全部实数 ω 存在且连续。

变量 ω 连续变化，故 DTFT 本身不是有限个频点组成的向量。每个频点的变换值可写为

X(e^{\mathrm{j}\omega}) =\sum_n x[n]\cos(\omega n) -\mathrm{j}\sum_n x[n]\sin(\omega n),

其实部和虚部分别是样本对余弦与正弦基函数的加权相关。


由于整数 n 满足

e^{-\mathrm{j}2\pi n}=1

有

\begin{aligned} X(e^{\mathrm{j}(\omega+2\pi)}) &=\sum_n x[n]e^{-\mathrm{j}(\omega+2\pi)n}\\ &=\sum_n x[n]e^{-\mathrm{j}\omega n}e^{-\mathrm{j}2\pi n}\\ &=X(e^{\mathrm{j}\omega}). \end{aligned}

因此，DTFT 沿频率轴以 2π 为周期重复，一个长度为 2π 的区间即可包含全部不重复信息。


反变换的正交性推导

将正变换代入下式并交换有限求和与积分次序：

\begin{aligned} \frac{1}{2\pi}\int_{-\pi}^{\pi} X(e^{\mathrm{j}\omega})e^{\mathrm{j}\omega m}\,\mathrm{d}\omega &=\sum_n x[n]\frac{1}{2\pi}\int_{-\pi}^{\pi} e^{\mathrm{j}\omega(m-n)}\,\mathrm{d}\omega. \end{aligned}

对整数 m−n，复指数在完整周期上的积分满足

\frac{1}{2\pi}\int_{-\pi}^{\pi} e^{\mathrm{j}\omega(m-n)}\,\mathrm{d}\omega =\begin{cases} 1, & m=n,\\ 0, & m\ne n. \end{cases}

积分只保留 n=m 的一项，因而得到反变换

\boxed{ x[m]=\frac{1}{2\pi}\int_{-\pi}^{\pi} X(e^{\mathrm{j}\omega})e^{\mathrm{j}\omega m}\,\mathrm{d}\omega }.

该式表明连续且周期的频域函数能够恢复离散样本；

时域离散性与频域周期性相互对应。

03
由定义得到的基本性质

线性与整数时移

对常数 a,b 及有限长序列 x1[n],x2[n]，有限和的分配律直接给出

a x_1[n]+b x_2[n] \quad\longleftrightarrow\quad aX_1(e^{\mathrm{j}\omega})+bX_2(e^{\mathrm{j}\omega}).

若 y[n]=x[n−n0]，令 m=n−n0，则

\begin{aligned} Y(e^{\mathrm{j}\omega}) &=\sum_n x[n-n_0]e^{-\mathrm{j}\omega n}\\ &=\sum_m x[m]e^{-\mathrm{j}\omega(m+n_0)}\\ &=\boxed{e^{-\mathrm{j}\omega n_0}X(e^{\mathrm{j}\omega})}. \end{aligned}

整数时移不改变幅度谱，只增加关于 ω 的线性相位。

正时移 n0>0 对应负相位斜率。

实序列的共轭对称

当 x[n] 为实数时，

\begin{aligned} X^*(e^{\mathrm{j}\omega}) &=\sum_n x[n]e^{\mathrm{j}\omega n} =X(e^{-\mathrm{j}\omega}). \end{aligned}

因此 |X(ejω)| 为偶函数，实部为偶函数，虚部为奇函数。

相位在幅度非零且选取一致分支时呈奇对称；在谱零点处，相位没有定义。

Parseval 能量关系

频域模平方可展开为

|X(e^{\mathrm{j}\omega})|^2 =\sum_n\sum_m x[n]x^*[m]e^{-\mathrm{j}\omega(n-m)}.

对一个完整频率周期积分后，正交关系消除所有 n≠ m 的交叉项，得到

\boxed{ \sum_n|x[n]|^2 =\frac{1}{2\pi}\int_{-\pi}^{\pi} |X(e^{\mathrm{j}\omega})|^2\,\mathrm{d}\omega }.

该关系说明有限序列的总能量在时域求和与频域积分下保持一致。

系数 1/(2π) 来自本文采用的正反变换归一化约定。

04
矩形序列的解析推导

定义长度为 N 的因果矩形序列

r_N[n]=\begin{cases} 1, & 0\le n\le N-1,\\ 0, & \text{其他}. \end{cases}

其 DTFT 是有限几何级数。

当 ω 不是 2π 的整数倍时，

\begin{aligned} R_N(e^{\mathrm{j}\omega}) &=\sum_{n=0}^{N-1}e^{-\mathrm{j}\omega n}\\ &=\frac{1-e^{-\mathrm{j}N\omega}}{1-e^{-\mathrm{j}\omega}}. \end{aligned}

利用恒等式

1-e^{-\mathrm{j}\alpha}=2\mathrm{j}e^{-\mathrm{j}\alpha/2}\sin(\alpha/2)

可化为

\boxed{ R_N(e^{\mathrm{j}\omega}) =e^{-\mathrm{j}\omega(N-1)/2} \frac{\sin(N\omega/2)}{\sin(\omega/2)} }.

在 ω=2πℓ 处，分子与分母同时为零，但原有限和的值为 N；

闭式表达式的连续极限同样为 N。

因子 e−jω(N−1)/2 表示矩形序列的中心位于 (N−1)/2，其相位具有相应的线性斜率。正弦比值为实数，并在变号时引入 π 相位跳变。

除 2π 整数倍外，零点位于

\omega_k=\frac{2\pi k}{N}, \qquad k=1,2,\ldots,N-1 \pmod{2\pi}.

围绕 ω=0 的两个最近零点为 ±2π/N，故主瓣的零点间宽度为 4π/N。

序列越长，归一化幅度谱的主瓣越窄；有限截断仍保留旁瓣。

05
固定参数算例

MATLAB 程序使用 N=9 的矩形序列和 32769 点频率网格。程序分别按有限和与正弦比值闭式计算 DTFT，另行检查 2π 周期性、右移 3 个样本后的线性相位因子、实序列共轭对称、数值反变换和 Parseval 关系。

程序已在 MATLAB R2026a 中实际运行，关键结果列于表 2。

MATLAB 数值验证结果
验证项目结果
有限和与解析闭式的最大复误差3.66×10−15
2π 周期性最大误差1.53×10−14
右移 3 个样本的时移性质最大误差6.83×10−15
实序列共轭对称最大误差0
数值反变换最大重建误差8.88×10−16
时域与频域能量9.000000000000
Parseval 能量绝对误差1.78×10−15
理论零点处的最大剩余幅度6.00×10−15
ω=π/2 单频点结果2+17j
图 1：长度为 9 的矩形序列、归一化 DTFT 幅度、谱零点外的主值相位及 2π 周期重复。（子图 1）图 1：长度为 9 的矩形序列、归一化 DTFT 幅度、谱零点外的主值相位及 2π 周期重复。（子图 2）图 1：长度为 9 的矩形序列、归一化 DTFT 幅度、谱零点外的主值相位及 2π 周期重复。（子图 3）图 1：长度为 9 的矩形序列、归一化 DTFT 幅度、谱零点外的主值相位及 2π 周期重复。（子图 4）

图 1　长度为 9 的矩形序列、归一化 DTFT 幅度、谱零点外的主值相位及 2π 周期重复。

闭式误差和零点残差处于双精度浮点舍入量级。反变换积分在均匀密集网格上采用梯形求积，重建结果与原 9 个单位样本一致；时域能量与频域积分均为 9。周期性误差略高于其他误差，主要来自大角度复指数求值中的有限精度约简，不影响解析恒等式。

06
数字实现的边界

07
适用范围

有限长 DTFT 可直接用于计算 FIR 系数序列的连续频率响应、观察有限数据记录的幅相分布，以及验证离散频率算法是否正确采样连续频谱。若长度为 L 的序列补零或截取到 NDFT≥ L 点，则其 NDFT 点离散傅里叶变换等于 DTFT 在

\omega_k=\frac{2\pi k}{N_{\mathrm{DFT}}}, \qquad k=0,1,\ldots,N_{\mathrm{DFT}}-1

处的均匀样值。

补零增加频率采样密度，但不改变原有限序列的连续 DTFT，也不消除矩形截断形成的旁瓣。

在通信、雷达和电子测量的公开基础处理中，DTFT 可描述短数据块与参考复指数的相关、分析脉冲包络的谱宽、检查数字滤波器系数的幅相响应，并比较不同窗函数对有限记录频谱的影响。使用这些结果时必须同时保留采样率、数据长度和归一化约定，否则归一化角频率无法唯一换算为赫兹。

08
工程上需要注意

  • 连续频率与离散网格：绘图和数值算法只能抽取有限频点，过疏网格可能漏掉窄峰或零点，不能替代解析检查。

  • 有限观测：截取数据等价于乘以时间窗，所得频谱同时包含原信号特征与窗函数影响；矩形窗会产生明显旁瓣。

  • 相位零点：频谱幅度为零或接近零时，相位对舍入和噪声极敏感，主值相位跳变不应直接解释为稳定延迟。

  • 归一化约定：正反变换的 2π 系数可按教材约定重新分配，能量公式必须与所用约定一致。

  • 有限精度：大角度复指数、相近复数相消和数值积分都会积累浮点误差；定点实现还受系数量化、累加器溢出和舍入影响。

  • 数据块边界：硬件单频点累加器依赖明确的 8 个有效样本边界；丢失或重复有效标记会改变频率系数与样本的对应关系。

  • 模型范围：DTFT 描述确定性离散序列，不直接给出随机过程的统计功率谱，也不包含采样时钟抖动、模拟前端失真或非线性效应。

有限长序列的 DTFT 是有限项复指数之和，因而对全部实频率存在、连续并以 2π 为周期。复指数在完整频率周期上的正交性给出反变换，同时导出 Parseval 能量关系；

变量替换和共轭运算则分别给出整数时移的线性相位因子与实序列的共轭对称。长度为 N 的矩形序列具有

R_N(e^{\mathrm{j}\omega}) =e^{-\mathrm{j}\omega(N-1)/2} \frac{\sin(N\omega/2)}{\sin(\omega/2)},

其最近零点确定 4π/N 的主瓣宽度，谱零点同时构成相位解释的边界。

MATLAB 实际运行结果在双精度误差范围内符合定义、闭式、周期性、时移、反变换和能量关系。SystemVerilog 文件把 ω=π/2 的固定复系数化为符号交换与累加，提供了从 DTFT 定义到单频点数字数据通路的直接对应；

动态 HDL 结论仍以可用仿真器上的自检运行为边界。

内容依据：每日精学第 05 课《有限长序列的离散时间傅里叶变换》。图表沿用原课固定参数数值结果，属于数值验证，不代表设备实测。

阅读15
分享
写评论...