每日精学 129｜整数抽取的频谱混叠求和
雷达 · 通信 · 电子战 / 基础精学 129
整数抽取把多个频率区间叠加到同一数字频率，抽取后仅靠数字处理无法分离已经相加的频谱分量。
01
模型与符号
抽取倍数M为正整数，取相位0即 y[n]=x[Mn]；公式适用于DTFT存在或广义频谱可定义的序列。
M：抽取倍数（无量纲）
Y(ω)：输出DTFT（幅度谱）
X(ω)：输入DTFT（幅度谱）
r：混叠支路编号（整数）
02
从模型到公式
周期选择器
利用单位根之和
可从原序列中选择每M个样本。
混叠公式
将选择器代入变换并把索引 n=Mm 代回，得到
因DTFT周期性，使用负号频移也得到同一支路集合。
无混叠条件
在主频率区间 −π≤Ω<π 内，若预滤波后 X(Ω)=0 对所有 |Ω|≥π/M 成立，则保留带内只有一支路贡献，
幅度谱中的 1/M 不能误解成单音时域幅度必然衰减；
谱线的频率尺度变换同时改变冲激权重。
混叠的相干相加与可逆性
两条输入复指数若频率相差 2π/M，抽取后在每个保留样本上完全相同。令
输出幅度由复系数之和决定，因此同相分量增强、反相分量可完全抵消。
对于 B=−A，非零输入在所有抽取位置均为零，说明抽取映射在不限制带宽的空间上存在非平凡零空间：
因此事后提高运算精度或增加输出字长不能恢复已经折叠的信息。
抗混叠滤波的作用是预先限制允许输入集合，使不同谱副本不在输出频带重叠；这是一项输入模型约束，而非地址选择逻辑自身的能力。
03
物理含义
抽取后样本无法区分相差新采样率整数倍的输入频率。预滤波的作用是阻止这些成分进入同一输出频率。
04
固定参数算例
以192点确定性记录、M=3验证抽取DFT等于三个输入频段求和除以3，并用两条混叠同频的复指数验证输出相干叠加。
| 量 | 数值 |
|---|---|
| 输入长度 | 192 |
| 抽取倍数 | 3 |
| 输出长度 | 64 |
| alias sum max error | 5.68711676635e-15 |
| coherent alias error | 2.76917754872e-14 |
| aliased amplitude | 1.5 |
图 1 整数抽取的频谱混叠求和的理论关系与固定参数验证
表中数值和图中曲线由 MATLAB 实际运行得到，随机种子为 20261048；运行版本与完整结果另存为独立文件。
05
数字实现的边界
RTL输出当前整数样本和外部样本索引b的相位0选择标志。实际抽取接口应仅把标志为1的样本计为输出；该核不包含抗混叠低通。 RTL 采用同步有效信号和异步低有效复位；
空拍保持状态且输出无效。除另有说明外，输入为16位有符号整数，结果为64位有符号整数，测试范围内不截断或饱和。输出在采样时钟上升沿更新。自检覆盖连续样本、空拍和复位；
未进行器件布局布线、时序收敛或功耗验证。
独立自检程序已在 Icarus Verilog 13.0 中实际运行并通过。该结果不包含器件综合、静态时序分析或板级测试。
06
适用范围
整数抽取用于把宽带采样流转换为低带宽测量通道。混叠支路求和公式可把各阻带噪声贡献映射到输出带内，从而确定抗混叠滤波需求。抽取后的异常谱线也可依据该映射检查是否来自上游未充分衰减的频率分量。
07
工程上需要注意
非理想抗混叠滤波使阻带噪声折回。
抽取相位不同会给各支路附加相位因子。
有限块长度必须能与所用DFT支路划分相容。
整数抽取把多个频率区间叠加到同一数字频率，抽取后仅靠数字处理无法分离已经相加的频谱分量。
内容依据：每日精学第 129 课《整数抽取的频谱混叠求和》。图表沿用原课固定参数数值结果，属于数值验证，不代表设备实测。