离散时间信号处理：理论、算法与工程实践
本文面向数字通信、音频、控制和嵌入式DSP领域的学生与工程师，从序列、线性系统、变换和滤波设计出发系统讨论离散时间信号处理。内容包括离散时间信号与序列、卷积与相关、DTFT/DFT/FFT、FIR与IIR滤波器、加窗与谱分析、多采样率处理、量化与定点实现等。文档篇幅较长，适合作为课程参考或工程设计背景资料。
图 1 ：离散时间信号处理中常见的示例序列。
图 2 ：基于窗函数的 FIR 低通滤波器幅频特性简化示意。
图 3 ：矩形窗、 Hann 窗和 Blackman 窗的时域形状及幅度谱示意。
|-----------|---------------------|------------|---------------------------|--------------------|
| 序列类型 | 定义 | 支撑集 | 示例 | 备注 |
| 有限长序列 | 仅在有限范围内非零。 | 常为0≤n≤N-1。 | 滤波器冲激响应、数据块。 | 适合采用DFT进行处理。 |
| 无限长序列 | 在无限多个n上非零。 | 支撑集无限。 | xn=a^n un、噪声过程等。 | 通常需要DTFT/Z变换或截断近似。 |
| 因果序列 | 对n<0为零。 | 支撑在n≥0。 | 大多数物理信号和系统响应。 | 与离散时间系统因果性相关。 |
| 反因果序列 | 对n>0为零。 | 支撑在n≤0。 | 多用于理论分析。 | 不直接对应物理可实现输出。 |
| 周期序列 | 满足xn+N0=xn。 | 周期为N0。 | 离散时间正弦、DFT基序列。 | DFT分析中假设序列周期延拓。 |
表1：离散时间信号处理中的常见序列类型及特征。
|---------------|-----------------|-------------------|-------------------|
| 滤波器类型 | 描述 | 典型用途 | 说明 |
| 低通滤波器 | 通过低频成分，抑制高频成分。 | 抗混叠、降噪、平滑。 | 由截止频率、过渡带和波纹指标刻画。 |
| 高通滤波器 | 通过高频成分，抑制低频成分。 | 边缘检测、强调快速变化。 | 可由低通滤波器经谱反转构造。 |
| 带通滤波器 | 通过某一频带内的成分。 | 信道选择、音调提取。 | 可由低通和高通滤波器组合实现。 |
| 带阻（陷波）滤波器 | 抑制窄频带信号，保留其他频率。 | 干扰抑制、电源噪声（如工频）去除。 | 广泛用于窄带噪声抑制。 |
| 全通滤波器 | 幅度响应为1，相位响应非平凡。 | 群延迟整形、相位均衡。 | 在IIR设计和控制应用中较常见。 |
表2：典型滤波器类型及应用场景。
|-----------------|-----------------|-------------------|-----------------------------|
| 指标 | 含义 | 关注点 | 说明 |
| 幅度响应 | 滤波器/系统对频率的幅度特性。 | 需要满足通带/阻带波纹与衰减指标。 | 决定频率选择性行为。 |
| 相位/ 群延迟 | 相位随频率的变化及其斜率。 | 与时间对齐和波形保持性相关。 | 线性相位滤波器可较好保持信号形状。 |
| 计算复杂度 | 单位样本或数据块所需运算量。 | 影响实时性和功耗。 | FIR vs IIR、FFT vs 直接DFT等权衡。 |
| 量化噪声 | 有限精度引入的噪声。 | 在定点实现中尤为关键。 | 通过字长选择和缩放策略进行控制。 |
表3：评估离散时间信号处理系统和滤波器的关键指标。
1. 离散时间信号与序列
离散时间信号用序列xn表示，仅在整数时间索引上定义。它们来源于连续时间信号的采样，或源自本身离散的系统（如数字通信比特流）。
理解序列的有限/无限长度、因果性和周期性等特征，是选择分析工具和判断实际可实现性的基础。
2. 线性时不变系统与差分方程
离散时间LTI系统常通过线性常系数差分方程刻画，将当前输出样本与当前及过去的输入/输出样本联系起来。这类系统构成数字滤波器的基础。
冲激响应hn可以完整描述LTI系统。系统分析包括求解差分方程、得到hn并利用卷积计算对任意输入的输出。
3. 卷积、相关与系统性质
卷积表达LTI系统输出yn=xn*hn，是离散时间信号处理中的核心运算，广泛应用于滤波和均衡等任务。
相关用于衡量序列之间的相似性，是检测、估计和同步的关键工具。自相关刻画信号功率和周期性；互相关则度量不同序列之间的对齐程度。
4. DTFT、DFT与FFT：频域分析
离散时间傅里叶变换（DTFT）提供序列的连续频率表示，用于观察信号谱和系统频率响应。冲激响应的DTFT给出H(e^{jω})，即LTI系统的频率响应。
离散傅里叶变换（DFT）及其快速实现（FFT）是计算频谱和通过频域乘法实现卷积的实用工具。合理选择窗函数和零填充方式是正确理解DFT结果的前提。
5. FIR与IIR滤波器设计
有限冲激响应（FIR）滤波器具有有限长度的冲激响应，天生稳定，易于获得线性相位，适用于对相位敏感的应用。
无限冲激响应（IIR）滤波器通过反馈形成有效无限长的冲激响应，用较少系数即可实现陡峭的频率特性，但需要在稳定性和定点实现方面进行仔细分析。经典模拟原型（Butterworth、Chebyshev、Elliptic等）经双线性变换等方法用于IIR设计。
6. 加窗、谱泄漏与实用谱分析
在有限长度数据上进行DFT分析会导致谱泄漏，即能量向邻近频点扩散。窗函数（如Hann、Hamming、Blackman等）通过修改时域数据形状来降低泄漏，但会加宽主瓣。
选择窗函数时需要在分辨率和泄漏之间权衡。工程实践中必须了解窗函数特性，避免对实验和仿真中的频谱结果做出错误解读。
7. 多采样率处理与采样率转换
多采样率处理通过抽取（降采样）和插值（升采样）改变序列的采样率，支持采样率转换、高效实现和多带处理等应用。
合理的滤波器设计对于防止多采样率系统中的混叠和镜像至关重要。多相结构和"高尚恒等式"（noble identities）可以显著提高实现效率。
8. 量化、定点DSP与实现问题
实际实现中使用有限精度，导致量化噪声和非线性效应。字长选择、缩放策略和溢出处理对性能和稳定性有重要影响。
工程师需要考虑定点效应、IIR滤波器中的极限环现象以及精度、存储、速度和功耗之间的权衡。在部署前应采用包含量化模型的仿真对算法进行评估。
9. 离散时间信号处理工程实践扩展笔记
扩展笔记 1：关于验证策略：在真实输入场景下同时验证DSP模块的时域和频域行为。
扩展笔记 2：关于窗函数选择：根据谱泄漏和分辨率需求为每类测量任务选择合适窗函数。
扩展笔记 3：关于多采样率设计：清晰记录采样率变化，并确保时钟和数据路径的一致性。
扩展笔记 4：关于定点分析：在实现早期进行字长优化和溢出检查，而不只在编码阶段处理。
扩展笔记 5：关于跨平台迁移：在不同DSP平台间迁移算法时考虑指令集和数值格式差异。
扩展笔记 6：关于维护与复现：将滤波器系数、设计脚本和测试向量纳入版本管理系统。
扩展笔记 7：关于教学与直觉：通过观察简单滤波器和变换对典型信号的作用来增强直觉理解。
扩展笔记 8：关于规格定义：从系统级SNR、带宽和时延需求出发制定滤波和DSP指标，而非仅依赖通用模板。