带通采样定理知乎｜深入解析带通信号采样原理与工程实践
从奈奎斯特采样定理的局限性出发，全面讲解带通采样（Bandpass Sampling）的核心机制、数学推导、典型应用场景及避坑指南，助您掌握数字通信系统设计关键技能。
立即探索带通采样什么是带通采样？——突破奈奎斯特的限制
传统奈奎斯特采样定理的困境
在信号处理领域，奈奎斯特采样定理告诉我们：为避免混叠（aliasing），采样频率 fs > 2 × f_max，其中 f_max 是信号最高频率分量。这个准则适用于基带信号（Baseband Signal），即频谱集中在 0Hz 附近的信号。
但现实中的通信系统往往处理的是带通信号（Bandpass Signal）——频谱能量集中在某个中心频率 f_c 附近，带宽为 B，如 f_c = 10, text{MHz},, B = 2, text{MHz}（即 9–11MHz）。此时若强行套用奈奎斯特准则，采样率需大于 22MHz，对 ADC 要求极高，成本与功耗陡增。
此时，带通采样定理（Bandpass Sampling Theorem）便成为破局关键。其核心思想是：只要满足特定的频率关系，即使采样频率 fs < 2 × f_max，也能无混叠地重构原始信号。
带通采样原理：频谱搬移的数学魔法
频谱周期延拓与混叠规避
当以频率 fs 对信号 x(t) 采样时，其频谱 X(f) 会以 fs 为周期进行周期延拓，形成 X_s(f) = frac{1}{T} sum_{k=-infty}^{infty} X(f - k f_s)。
对基带信号，所有延拓谱带以 0Hz 为中心对称，易重叠；而对带通信号，若选择合适的 fs，可使相邻谱带仅在非信号频段接触（甚至不接触），从而避免混叠。
如下图所示：原始带通信号频谱（蓝色）位于 [f_1, f_2]，采样后产生的镜像谱带（红色）若不与原谱带重叠，则可滤波恢复。
理想采样频谱分布示意图（频域）
原信号：9–11 MHz
采样率：fs = 5, text{MHz}
镜像谱带位置：k × 5, text{MHz} + [9,11]，即：
• k=0：9–11 MHz
• k=1：14–16 MHz
• k=-1：-4–-2 MHz → 镜像至 4–6 MHz（与原信号无重叠）
✅ 无混叠！可用带通滤波器提取 4–6 MHz 基带等效信号。
采样率下限公式
对带宽为 B、最低频率为 f_L、最高频率为 f_H 的带通信号，满足无混叠采样的采样率 fs 需满足：
其中 n 为整数，表示第 n 个允许的采样区间。
物理意义：每个 n 对应一组可行的 fs 区间，需满足该区间与 [2B, 2f_H] 有交集。
若 f_L = 9, text{MHz},, f_H = 11, text{MHz}（即 B = 2, text{MHz}），则：
• 当 n=5：3.6, text{MHz} leq fs leq 4.5, text{MHz} → 取 fs = 4, text{MHz}？❌ 不满足下限
• 当 n=4：5.5, text{MHz} leq fs leq 6.0, text{MHz} → 取 fs = 5.5, text{MHz} ✅
• 当 n=3：7.33, text{MHz} leq fs leq 9.0, text{MHz} → 取 fs = 8, text{MHz} ✅
带通采样公式详解与推导
频谱周期延拓的数学表达
设连续带通信号为 x(t)，其傅里叶变换为 X(f)，仅在 [f_L, f_H] 非零。采样后信号为：
其频谱为：
要避免混叠，需保证任意两个相邻谱带不重叠，即：
化简后可得：
结合频带位置约束，最终导出前述区间公式。
混叠发生的临界条件
混叠发生当且仅当存在整数 k_1 ne k_2，使得：
即区间长度 B 大于区间中心距：|k_1 - k_2| f_s - Delta f，其中 Delta f 为频带中心偏移。
典型混叠案例：
- f_L = 8, text{MHz},, f_H = 12, text{MHz},, B = 4, text{MHz}
- 若 fs = 9, text{MHz}：
- 原谱：8–12 MHz
- 镜像（k=1）：17–21 MHz
- 镜像（k=−1）：−1–3 MHz → 镜像至 6–9 MHz（与原谱 8–12 MHz 重叠于 8–9 MHz）❌
实际工程中，因 ADC 非理想性（如带通响应不平坦、时钟抖动），建议在理论区间内预留 10%~15% 安全余量。
信号重构方法：正交下变频 + 低通滤波
带通采样后得到的离散序列 x[n] = x(nT) 并非原始基带信号，而是带通信号经频谱搬移后的等效低通信号。重构需两步：
- 正交下变频：将 x[n] 乘以复指数 e^{j 2pi f_0 n / f_s}（f_0 为希望搬移至的中心频率，通常取 0Hz），实现频谱左移。
- 低通滤波：用 FIR 或 IIR 滤波器提取目标基带带宽 B 内的信号。
实例：9–11MHz 信号以 fs=5, text{MHz} 采样
采样后谱带位于 4–6 MHz（因镜像折叠）
→ 下变频至 0Hz：乘以 e^{-j 2pi cdot 5, text{MHz} cdot n / 5, text{MHz}} = e^{-j 2pi n} = 1（无需额外调制）
→ 低通滤波（截止 2.5MHz）→ 提取 4–6 MHz 的镜像 → 再次镜像得原始 9–11 MHz 信号
注意：若采用实数 ADC，重构后得到的是实信号，需考虑正负频对称性；若用复数 ADC（I/Q 采样），可直接获得解析信号，避免镜像干扰。
实例分析：从人声到无线通信
–200Hz 人声，采样率 800Hz
原始信号：中心频率 f_c = 150, text{Hz}，带宽 B = 100, text{Hz}。
传统奈奎斯特：需 fs > 400, text{Hz}，取 44.1kHz 过度采样。
带通采样方案：
- 取 n = 3：区间 [2×200/3, 2×100/2] = [133.3, 100] → 无解！
- 取 n = 2：区间 [400/2, 200/1] = [200, 200] → 仅临界点 fs=200, text{Hz}，但易混叠。
- 取 n = 1：区间 [400/1, ∞] = [400, ∞] → 回退至奈奎斯特。
? 结论：对极低频窄带信号，带通采样优势不明显，因 f_c 太小，难以满足 f_c gg B。
✅ 改进：用复指数调制将信号搬移到更高频段再采样（如搬至 1–2kHz），再以 fs=4, text{kHz} 采样，可大幅降低数据率。
–15MHz 中频信号，采样率 20.25MHz
调频广播中频：f_c = 10.7, text{MHz}，带宽 B = 200, text{kHz}（含保护带）。
传统：需 fs > 21.4, text{MHz}，取 25MHz。
带通采样方案：
- 取 n = 5：区间 [2×15/5, 2×7/4] = [6, 3.5] → 无效（下限 > 上限）
- 取 n = 4：区间 [30/4, 14/3] = [7.5, 4.67] → 无效
- 取 n = 3：区间 [30/3, 14/2] = [10, 7] → 无效
- 取 n = 2：区间 [30/2, 14/1] = [15, 14] → 无效
- 取 n = 1：区间 [30/1, ∞] = [30, ∞] → 需 >30MHz
⚠️ 发现：该信号带宽相对中心频率过大（B/f_c ≈ 20%），不满足窄带条件，带通采样区间为空！
✅ 实际方案：选用 fs = 20.25, text{MHz}（非整数比），通过数字下变频（DDC）配合多级滤波器组实现精确重构，本质是欠采样（Undersampling）+ 数字校正。
GHz 频段，100MHz 带宽，ADC 采样率 7.68GHz
G NR 常用频段：f_c = 3.5, text{GHz}，带宽 B = 100, text{MHz}。
直接采样：需 fs > 7.1, text{GHz}，现有 ADC 成本极高。
带通采样方案：
- 取 n = 70：计算得 fs in [100.57, 101.43], text{MHz} → 无效（太小）
- 实际采用分频带采样（Sub-Nyquist Sampling）+ 压缩感知（Compressed Sensing）技术，在已知信号稀疏性前提下，以远低于奈奎斯特率的采样率恢复信号。
网友常见误区解析
❌ 误区：带通采样能突破奈奎斯特极限，采样率可趋近于 0
真相：带通采样仍有下限——采样率必须大于信号带宽的两倍（fs > 2B），否则即使频谱不重叠，也无法满足时域采样定理的最小信息率要求（即 Shannon 信息论下限）。
反例：
- 信号：99–101 MHz（B=2, text{MHz}）
- 若取 fs = 3, text{MHz}（>2B=4MHz？❌ 实际 3 < 4）
- 谱带重叠：原谱 99–101，镜像 k=−32：−297–−295 → 镜像至 297–295 MHz（不重叠）
但 k=−31：−286–−284 → 镜像至 286–284 MHz……持续到 k=−50：−499–−497 → 镜像至 499–497 MHz
关键：因 fs < 2B，相邻谱带间距 fs 小于带宽 B，必然重叠！
带通采样不是绕过奈奎斯特，而是在特定条件下“利用”奈奎斯特准则——将带通信号等效为一个带宽为 B 的低通信号，其采样率需 >2B。
❌ 误区：只要满足 fs > 2B 就能无混叠采样
真相：必须同时满足频带位置约束！仅 fs > 2B 是必要不充分条件。
经典反例：
- 信号：10–12 MHz（B=2, text{MHz}）
- 取 fs = 4.5, text{MHz} > 2B=4, text{MHz}
- 谱带位置：
- 原谱：10–12 MHz
- 镜像 k=−2：−9–−7 MHz → 镜像至 9–7 MHz（即 7–9 MHz）❌ 与原谱 10–12 MHz 不重叠？再看 k=−1：−4.5–−2.5 → 4.5–2.5 MHz（2.5–4.5 MHz）
但 k=−3：−13.5–−11.5 → 13.5–11.5 MHz → 与原谱 10–12 MHz 重叠于 11.5–12 MHz！
✅ 正确做法：用前述区间公式筛选可行 fs，如取 fs = 5.5, text{MHz}（n=4 区间内）。
❌ 误区：带通采样后直接用低通滤波即可恢复原信号
真相：带通采样得到的频谱是原信号的镜像，需先通过数字下变频（DDC）将目标频带搬移到基带，再滤波。
流程对比：
| 步骤 | 基带采样 | 带通采样 |
| 1. 采样 | x(t) → x[n] | x(t) → x[n] |
| 2. 频谱 | X(f) 周期延拓 | X(f) 周期延拓 + 镜像 |
| 3. 滤波 | 直接低通滤波 | 先下变频（乘复指数）→ 再低通滤波 |
数学验证：
否则，滤波后得到的是带通信号的镜像，无法还原原始相位与幅度。
高频问答：带通采样知乎热议TOP5
A：二者目标不同但可结合：
- 带通采样：针对已知频带位置的带通信号，利用其频谱稀疏性，在满足无混叠条件下实现确定性采样（仍需 >2B）。
- 压缩感知：针对时域稀疏信号（如稀疏频谱），通过非线性重构算法，以远低于奈奎斯特率的采样率恢复信号。
✅ 实际系统中，常先用带通采样将信号搬移到较低频段，再用压缩感知进一步降采样率，实现“双重欠采样”。
A：取决于是否需要保留正负频信息：
- 实数 ADC：输出为实序列，频谱共轭对称。若目标频带位于 0–fs/2 且无镜像干扰，可直接采样；否则需通过滤波或算法消除镜像（复杂）。
- 复数 ADC（I/Q 采样）：两路正交采样（I/Q），天然获得解析信号，避免镜像问题，是现代软件无线电（SDR）首选方案。
? 工程建议：若频带靠近 fs/2（如 4–6MHz，fs=10MHz），用实数 ADC 易混叠，推荐 I/Q 采样。
A：是！原因有二：
- 带通信号能量集中在高频段，时钟抖动引起的采样时刻误差 Delta t 会导致相位噪声恶化，表现为 Delta phi = 2pi f_c Delta t，f_c 越高，相位误差越大。
- 欠采样时，微小的 fs 偏移可能导致谱带重叠（如 fs=5.5, text{MHz} pm 0.1, text{MHz} 即可能引发混叠）。
• 采用低抖动时钟（如 VCXO + PLL）
• 使用数字时钟校正算法（如基于训练序列的定时恢复）
• 在 ADC 前加抗混叠滤波器（即使带通采样，也需抑制带外干扰）
A：以下为完整示例代码（9–11MHz 信号，fs=5.5MHz）：
t = 0:T:1e-3; % 1ms 时长
x = cos(2pi10e6t); % 10MHz 载波
X = fftshift(fft(x));
f = (-length(x)/2 : length(x)/2-1) fs / length(x);
plot(f/1e6, abs(X)/max(abs(X))); grid on;
xlabel('频率 (MHz)'); ylabel('幅度'); title('带通采样频谱');
运行后可见：频谱主峰在 ±10MHz，镜像峰在 ±(10−5.5)=±4.5MHz，且无重叠。用带通滤波器提取 4–6MHz 区间，再下变频即可还原。
A：前景广阔但挑战巨大：
- 优势：6G 将进入太赫兹频段（0.1–10 THz），带宽超 100GHz，奈奎斯特采样率超 20THz，现有 ADC 无法实现；带通采样+数字预处理是唯一可行路径。
- 挑战：太赫兹信号相位噪声极大，传统带通采样难以稳定工作；需结合模拟预滤波与AI 辅助重构算法。
✅ 最新进展：2024 年 IEEE 提出“可重构带通采样器”（Reconfigurable BPS），通过可调滤波器动态匹配信号频带，已在 300GHz 实验中实现 10Gbps 传输。