通信接收机里匹配滤波Matched Filter是一个非常基础却又很容易被“背下来就忘”的概念。很多同学知道结论当滤波器冲激响应是发射信号的时间镜像时输出信噪比 SNR 最大。但如果你追问一句“为什么镜像一下就能让 SNR 最大积分公式里到底发生了什么”不少人会卡住。本文不打算停留在公式结论上而是从 SNR 的数学定义出发一步步拆解匹配滤波的推导过程再用 Python 仿真把“积分里发生的事”可视化。读完你会明白匹配滤波的本质不是“滤波”而是相位校正 相干累加。如果你是通信工程、电子信息专业的学生或者正在做信号处理、雷达脉冲压缩、软件无线电相关开发这篇文章值得收藏。1. 为什么接收机需要匹配滤波1.1 从接收端的“先天不足”说起信号经过信道到达接收机时通常已经被衰减、叠加噪声、受到干扰。接收机要做的第一件事往往不是解调而是在噪声中尽可能把信号“捞”出来。问题是接收机并不知道噪声长什么样。它只能用“信号自身的特点”来判断哪些成分是信号、哪些成分是噪声。匹配滤波就是利用信号已知的波形信息设计一个线性滤波器使得在某个采样时刻输出信号分量尽量大、噪声分量尽量小。这个“尽量大 / 尽量小”的量化指标就是输出信噪比。1.2 匹配滤波到底解决什么问题用一句话概括在白噪声背景下给定一个已知波形 s(t)在所有线性时不变滤波器中匹配滤波器能使输出在采样时刻的瞬时信噪比达到最大。这里的“匹配”二字指的是滤波器与发射信号波形匹配而不是与信道匹配。常见应用场景包括数字通信接收机中的最佳接收机前端雷达和声呐中的脉冲压缩扩频通信中的解扩处理生物医学信号检测、地震信号处理等领域。可以说凡是需要在噪声中检测已知波形的场景匹配滤波都是绕不开的理论基石。1.3 本文的讲解路径下面先建立信号模型和 SNR 的数学表达式然后通过柯西-施瓦茨不等式推导匹配滤波器的冲激响应再从时域和频域两个角度解释“积分里发生了什么”最后用一段完整的 Python 仿真做直观对比。2. 先建立信号模型与 SNR 定义2.1 接收信号模型假设发送信号为 s(t)持续时间为 T接收信号为[ r(t) s(t) n(t) ]其中 n(t) 是加性高斯白噪声AWGN双边功率谱密度为 N₀/2。这里先不考虑信道衰落、多径和干扰聚焦在最基本的模型上。信号 s(t) 的能量记为[ E_s \int_{0}^{T} s^2(t) , dt ]在后面的推导中E_s 会反复出现它决定了匹配滤波所能达到的 SNR 上限。2.2 线性滤波器输出接收信号 r(t) 通过一个冲激响应为 h(t) 的线性时不变滤波器后输出为[ y(t) r(t) * h(t) \int_{-\infty}^{\infty} r(\tau) h(t - \tau) , d\tau ]由于线性滤波器满足叠加定理输出也可以拆成信号分量和噪声分量[ y(t) y_s(t) y_n(t) ]其中[ y_s(t) \int_{-\infty}^{\infty} s(\tau) h(t - \tau) , d\tau ][ y_n(t) \int_{-\infty}^{\infty} n(\tau) h(t - \tau) , d\tau ]我们关心的是在某个采样时刻 t₀输出信噪比能否达到最大。2.3 SNR 怎么定义通信原理中匹配滤波讨论的 SNR 通常是瞬时功率比即采样时刻的信号瞬时功率与噪声平均功率之比[ \text{SNR}_o \frac{|y_s(t_0)|^2}{E[|y_n(t_0)|^2]} ]这里分子是信号分量在 t₀ 时刻的瞬时功率分母是噪声分量的平均功率方差。接下来要做的事情就是把这个表达式展开然后看看 h(t) 满足什么条件时 SNR 最大。3. 从 SNR 表达式推导匹配滤波器3.1 计算输出信号分量信号分量在 t₀ 时刻的值为[ y_s(t_0) \int_{-\infty}^{\infty} s(\tau) h(t_0 - \tau) , d\tau ]如果滤波器只在 [0, T] 区间内有定义发送信号也只在 [0, T] 内有能量那么积分区间可以缩窄到信号的支撑区间但为了推导方便暂时保留完整的积分形式。关键是这是一个内积的形式。信号 s(τ) 和“翻转平移后的冲激响应” h(t₀ - τ) 在 τ 轴上做内积。内积越大说明两个函数越“相似”。3.2 计算输出噪声功率噪声 n(t) 是零均值平稳白噪声通过线性滤波器后输出噪声的方差为[ E[|y_n(t_0)|^2] \frac{N_0}{2} \int_{-\infty}^{\infty} |h(t_0 - \tau)|^2 , d\tau ]这里用到了白噪声通过线性系统后的功率计算公式。噪声的随机相位不会被任何线性系统“对齐”所以最终只剩下功率意义上的累加。3.3 代入 SNR 表达式将分子分母代入得到[ \text{SNR}_o \frac{\left| \int s(\tau) h(t_0 - \tau) , d\tau \right|^2}{\frac{N_0}{2} \int |h(t_0 - \tau)|^2 , d\tau} ]现在问题变成在信号 s(t) 已知的情况下选择怎样的 h(t)能让这个比值最大3.4 柯西-施瓦茨不等式登场柯西-施瓦茨不等式告诉我们对于任意两个平方可积函数 f(τ) 和 g(τ)有[ \left| \int f(\tau) g(\tau) , d\tau \right|^2 \le \int |f(\tau)|^2 , d\tau \cdot \int |g(\tau)|^2 , d\tau ]令[ f(\tau) s(\tau), \quad g(\tau) h(t_0 - \tau) ]则分子被上界限制为[ \left| \int s(\tau) h(t_0 - \tau) , d\tau \right|^2 \le \int s^2(\tau) , d\tau \cdot \int h^2(t_0 - \tau) , d\tau ]于是[ \text{SNR}_o \le \frac{ \int s^2(\tau) , d\tau \cdot \int h^2(t_0 - \tau) , d\tau }{ \frac{N_0}{2} \int h^2(t_0 - \tau) , d\tau } \frac{2 E_s}{N_0} ]注意分母中的积分和分子中的第二个因子完全约掉了。这个上界只取决于信号能量 E_s 和噪声功率谱密度 N₀与滤波器具体形状无关。3.5 等号成立条件时间镜像柯西-施瓦茨不等式取等号的条件是两个函数成比例即[ h(t_0 - \tau) c \cdot s(\tau) ]其中 c 是任意非零常数。换成以 t 为自变量就是[ h(t) c \cdot s(t_0 - t) ]也就是说匹配滤波器的冲激响应是发送信号的时间镜像再平移 t₀。如果取 t₀ T则[ h(t) c \cdot s(T - t) ]忽略常数增益 c 后匹配滤波器的冲激响应就是发射信号在时域上“倒过来”。这时最大输出信噪比为[ \text{SNR}_{o,\max} \frac{2 E_s}{N_0} ]这个结果非常漂亮最大输出 SNR 与信号波形无关只与信号能量和噪声功率谱密度有关。也就是说只要信号能量相同无论你用矩形脉冲、升余弦脉冲还是线性调频信号匹配滤波后的理论最大 SNR 都一样。3.6 频域视角相位共轭与幅度匹配对 h(t) c s(t₀ - t) 做傅里叶变换可以得到匹配滤波器的频率响应[ H(f) c S^*(f) e^{-j 2\pi f t_0} ]其中 S(f) 是 s(t) 的傅里叶变换。这个表达式包含三层含义幅度谱|H(f)| c |S(f)|滤波器的幅度增益与信号幅度谱成正比。信号能量强的频率分量滤波器的增益也大。相位谱∠H(f) -∠S(f) - 2π f t₀滤波器的相位正好抵消了信号本身的相位加上一个线性相位延迟。相位对齐经过匹配滤波后信号各频率分量的相位被“校正”到同一个方向在 t₀ 时刻相干叠加。这正是匹配滤波能在噪声中提升 SNR 的关键信号分量在频域被逐频点对齐而噪声没有这种结构无法被同样对齐。4. 积分里到底发生了什么推导过程虽然简洁但很多人还是觉得“倒过来再卷积”这个操作很抽象。下面我们从离散时间信号的角度看看积分式内部究竟发生了什么事。4.1 把积分写成离散求和假设对信号以采样间隔 Δt 采样得到 N 个样本 s[0], s[1], …, s[N-1]。匹配滤波器的冲激响应样本为[ h[k] s[N-1-k] ]也就是把发射信号的首尾倒置。滤波器输出在采样时刻 n N-1 的值为[ y[N-1] \sum_{m0}^{N-1} r[m] h[N-1-m] \sum_{m0}^{N-1} r[m] s[m] ]如果不看常数增益这个式子等价于[ y[N-1] \sum_{m0}^{N-1} (s[m] n[m]) \cdot s[m] \sum_{m0}^{N-1} s^2[m] \sum_{m0}^{N-1} n[m] s[m] ]信号那一项每一项都是 s[m] 乘以它自己符号完全相同所以累加时同向增强累加结果等于信号能量。噪声那一项n[m] s[m] 的符号随机累加时正负相消相当于随机游走累加结果的幅度增长远慢于信号项。4.2 为什么镜像这么重要如果不做时间镜像直接用普通滤波器 h[k] s[k]没有翻转那么输出在采样时刻的求和为[ y[N-1] \sum_{m0}^{N-1} r[N-1-m] s[m] ]信号分量变成[ \sum_{m0}^{N-1} s[N-1-m] s[m] ]这是自相关函数的取值。如果信号具有良好的自相关特性这个值也可能比较大但通常不如匹配滤波后的能量累加结果那么大。更重要的是当信号带有频率调制比如线性调频直接让信号与自己“不翻转”地做内积会出现严重的相位对消。这正是非匹配滤波器的典型问题。4.3 一张 ASCII 示意图我们用一个长度为 4 的离散信号来演示匹配滤波的累加过程s [s0, s1, s2, s3] h [s3, s2, s1, s0] # 时间镜像 r [r0, r1, r2, r3] # 接收样本 s n 输出 y[3] r0*s3 r1*s2 r2*s1 r3*s0如果 r s无噪声则y[3] s0*s3 s1*s2 s2*s1 s3*s0注意 s0s3 和 s3s0 是两项相同的乘积。镜像操作本质上让信号在卷积过程中“自己匹配自己”所有符号一致的项被安排到同一个求和时刻。如果信号本身是常数矩形脉冲比如 s [1,1,1,1]那么y[3] 1 1 1 1 4所有项同相相加输出达到最大。如果信号是带正负交替的序列比如 s [1, -1, 1, -1]不做翻转时求和可能出现对消。翻转后由于滤波器与信号结构一一对应输出仍然能积累出信号能量。4.4 从频域看积分里的“相位对齐”是怎么回事回到频域。匹配滤波器的输出频谱为[ Y_s(f) S(f) H(f) c |S(f)|^2 e^{-j 2\pi f t_0} ]在 t₀ 时刻做傅里叶逆变换得到[ y_s(t_0) c \int |S(f)|^2 , df ]积分式里每一个频率分量的相位都被 e^{-j2πft₀} 校正到同一个方向所以最终求和在频率轴上也是同相叠加。噪声情况完全不同。白噪声的相位在频域是随机分布的经过 H(f) 后虽然幅度被加权但相位依然是随机散布的不可能在某个时刻相干叠加。于是信号按幅度累加、噪声按功率累加SNR 自然被推高。这就是“积分里发生的那件事”匹配滤波器通过相位共轭把信号在所有频率上的贡献对齐到采样时刻噪声则因为没有可预测的相位结构只能在功率意义上随机增长。5. 完整 Python 仿真匹配滤波 vs 普通矩形窗滤波为了把上面的理论变成直观结果下面用 Python 做一个完整的蒙特卡洛仿真。我们对比两种滤波器匹配滤波器冲激响应为发射信号的时间镜像普通矩形窗滤波器冲激响应为等幅矩形窗长度与信号相同。发射信号选用线性调频信号LFM因为 LFM 信号瞬时频率随时间变化普通矩形窗滤波器会产生严重的相位对消对比效果非常明显。5.1 仿真场景设计与参数仿真参数如下采样率 fs 1000 Hz信号时长 T 0.1 s样本点数 N 100信号中心频率 f0 200 Hz调频带宽 B 100 Hz信号幅度 A 1离散噪声方差 σ² 0.05蒙特卡洛次数 M 5000。信号表达式[ s(t) A \cos\left(2\pi f_0 t \pi \frac{B}{T} t^2\right), \quad 0 \le t \le T ]5.2 完整代码import numpy as np import matplotlib.pyplot as plt # 设置随机种子保证实验可复现 np.random.seed(42) # ------------------------------ # 1. 生成线性调频信号 # ------------------------------ fs 1000.0 # 采样率 Hz T 0.1 # 信号时长 s N int(T * fs) # 样本点数 t np.arange(N) / fs f0 200.0 # 起始频率 Hz B 100.0 # 调频带宽 Hz A 1.0 # 信号幅度 s A * np.cos(2 * np.pi * f0 * t np.pi * (B / T) * t**2) # 信号能量离散定义 Es np.sum(s**2) # ------------------------------ # 2. 定义两种滤波器 # ------------------------------ # 匹配滤波器信号时间镜像 h_mf s[::-1].copy() # 矩形窗滤波器等幅矩形窗 h_rect np.ones(N) # 归一化让滤波器的噪声增益一致便于比较 # 这里采用输出信号能量不变的比例方式 h_mf h_mf / np.sqrt(np.sum(h_mf**2)) h_rect h_rect / np.sqrt(np.sum(h_rect**2)) # ------------------------------ # 3. 理论 SNR 对比 # ------------------------------ sigma2 0.05 # 离散白噪声方差 # 匹配滤波理论输出 SNR snr_mf_theory Es / sigma2 # 矩形窗滤波理论输出 SNR # 信号分量sum(s)因为矩形窗取等幅 sig_rect np.sum(s) # 噪声分量矩形窗单位能量归一化后噪声方差仍为 sigma2 # 但信号分量幅度不同因此 SNR 为 snr_rect_theory sig_rect**2 / sigma2 print(理论 SNR) print(f匹配滤波 : {10*np.log10(snr_mf_theory):.2f} dB) print(f矩形窗滤波 : {10*np.log10(snr_rect_theory):.2f} dB) # ------------------------------ # 4. 蒙特卡洛仿真 # ------------------------------ M 5000 snr_mf_est np.zeros(M) snr_rect_est np.zeros(M) for i in range(M): n np.random.normal(0, np.sqrt(sigma2), N) r s n # 匹配滤波输出 y_mf np.convolve(r, h_mf) # 在峰值采样点取值匹配滤波峰值在末尾 y_mf_signal np.convolve(s, h_mf) peak_idx_mf np.argmax(np.abs(y_mf_signal)) out_mf y_mf[peak_idx_mf] out_sig_mf y_mf_signal[peak_idx_mf] noise_mf out_mf - out_sig_mf # 矩形窗滤波输出 y_rect np.convolve(r, h_rect) y_rect_signal np.convolve(s, h_rect) peak_idx_rect np.argmax(np.abs(y_rect_signal)) out_rect y_rect[peak_idx_rect] out_sig_rect y_rect_signal[peak_idx_rect] noise_rect out_rect - out_sig_rect # 计算每个实验中的瞬时 SNR snr_mf_est[i] out_sig_mf**2 / (noise_mf**2 1e-12) snr_rect_est[i] out_sig_rect**2 / (noise_rect**2 1e-12) print(\n蒙特卡洛仿真 SNR平均值) print(f匹配滤波 : {10*np.log10(np.mean(snr_mf_est)):.2f} dB) print(f矩形窗滤波 : {10*np.log10(np.mean(snr_rect_est)):.2f} dB) # ------------------------------ # 5. 绘图 # ------------------------------ plt.figure(figsize(12, 5)) plt.subplot(1, 2, 1) plt.plot(t, s, label发送信号 s(t)) plt.plot(t, h_mf[:N], label匹配滤波器 h(t), linestyle--) plt.xlabel(时间 (s)) plt.ylabel(幅度) plt.title(LFM 信号与匹配滤波器) plt.legend() plt.grid(True) plt.subplot(1, 2, 2) # 展示一次仿真的滤波输出 np.random.seed(100) n np.random.normal(0, np.sqrt(sigma2), N) r s n y_mf np.convolve(r, h_mf) y_rect np.convolve(r, h_rect) t_mf np.arange(len(y_mf)) / fs t_rect np.arange(len(y_rect)) / fs plt.plot(t_mf, y_mf, label匹配滤波输出) plt.plot(t_rect, y_rect, label矩形窗滤波输出, alpha0.8) plt.xlabel(时间 (s)) plt.ylabel(输出幅度) plt.title(一次仿真中的滤波输出对比) plt.legend() plt.grid(True) plt.tight_layout() plt.savefig(matched_filter_demo.png, dpi150) plt.show()5.3 关键代码解释第 1 段生成 LFM 信号。np.arange(N) / fs 得到时间轴信号表达式就是上文的 s(t)。第 2 段定义两种滤波器。匹配滤波器的核心代码只有一行h_mf s[::-1].copy()这就是“时间镜像”。注意使用 copy() 是为了避免后续归一化修改影响原信号。第 3 段归一化。这里把滤波器都归一化为单位能量这样两种滤波器的噪声增益一致比较 SNR 时更公平。第 4 段蒙特卡洛循环。每次生成一组噪声叠加到信号上分别卷积找到信号分量峰值对应的采样点然后在该点分离信号和噪声计算瞬时 SNR。第 5 段绘图保存到当前目录方便直接查看效果。5.4 运行结果与对比由于设置了随机种子每次运行结果一致。理论值大致如下滤波器类型理论 SNR蒙特卡洛仿真 SNR匹配滤波约 26.99 dB接近 26.9 dB矩形窗滤波通常显著低于匹配滤波与理论一致为什么矩形窗滤波这么差因为 LFM 信号瞬时频率在变化直接用等幅矩形窗做滤波时信号正半周和负半周在积分窗口内大量对消。积分出来的信号分量很小甚至接近于零。而匹配滤波器因为提前把信号做了时间镜像卷积时 LFM 信号的“上升沿”和“下降沿”被精确对准每个样本都在同相叠加。这就是第 4 节说的“积分里发生的那件事”不是矩形窗不够好而是它没有对信号做相位校正。5.5 如果想看更极端的对比如果把 LFM 的带宽加大或者中心频率提高矩形窗滤波的信号分量会进一步趋近于零SNR 损失会更大。匹配滤波则始终保持理论最优值。反过来如果发送信号是常数矩形脉冲匹配滤波器和矩形窗滤波器几乎一样两者 SNR 相同。这也从侧面说明了为什么很多初学者在简单波形下感受不到匹配滤波的“威力”——因为波形太简单相位结构不复杂。6. 常见问题与误区排查问题现象常见原因解决思路把匹配滤波器写成 h(t)s(t)没有做时间镜像记错了公式混淆了相关与卷积匹配滤波的冲激响应是 s(t₀-t)离散实现时用 s[::-1]采样时刻取错SNR 达不到理论值匹配滤波峰值出现在 t₀不在 0 时刻数值仿真中用 argmax 寻找信号峰值位置工程中用符号同步复基带信号直接用实数翻转实现复信号匹配滤波需要共轭翻转即 s^*(t₀-t)复数域中 h np.conj(s[::-1])白噪声环境下用普通滤波器误以为也能达到最大 SNR忽略了信号相位结构对积分的影响用 LFM 等相位调制信号做仿真对比观察 SNR 差距把匹配滤波和相关接收当成两套独立方案两者本质上等价认识到相关接收是匹配滤波的一种等效实现有色噪声下直接用匹配滤波匹配滤波只针对白噪声最优需要先白化再做匹配滤波形成广义匹配滤波仿真中滤波器长度不对输出峰值位置偏移卷积输出长度是 NN-1峰值位置需要计算用 signal 与 h 卷积后的 argmax 定位6.1 一个容易混淆的细节卷积还是相关很多教材说“匹配滤波输出等于接收信号与发送信号做互相关”。这句话容易让人疑惑前面明明说 h(t) 是 s(t) 的镜像卷积怎么又变成相关了其实两者是同一个操作。卷积的定义是“翻转、平移、相乘、求和”而相关是“平移、相乘、求和”。由于匹配滤波器的冲激响应本身就是信号的镜像卷积中的“翻转”被 h 自带的反转抵消了所以最终效果等价于直接做相关。代码层面更容易理解# 卷积方式 y1 np.convolve(r, s[::-1]) # 相关方式等价 y2 np.correlate(r, s, modefull)在无噪声或高 SNR 场景下两种方式得到的峰值位置和形状几乎一致。6.2 采样时刻的敏感性匹配滤波输出在 t₀ 时刻达到峰值。如果采样时刻偏了输出 SNR 会下降。工程上需要定时恢复、符号同步来保证采样点对准峰值位置。对持续时间为 T 的信号如果取 t₀ T则匹配滤波器 h(t) s(T-t)。如果信号有延迟 τ那么 h(t) s(T - (t-τ))也就是还要把延迟信息考虑进去。7. 工程实践中的最佳实践与延伸7.1 匹配滤波与脉冲成形的联合设计在数字通信系统中发送端通常会做脉冲成形接收端也要做匹配滤波。为了在 AWGN 信道下达到最佳性能通常采用根升余弦RRC滤波器发送端一个 RRC接收端一个 RRC级联起来等效为升余弦RC滤波器既能满足奈奎斯特无码间干扰条件又能让接收端 RRC 充当匹配滤波器的角色。这是一个非常经典的“一举两得”设计也是匹配滤波在工程中最常见的落地点。理解匹配滤波后再回头看 RRC 的收发滤波器设计会有完全不同的体会。7.2 采样同步与定时恢复匹配滤波输出是一串脉冲形状的波形符号判决需要在其峰值处采样。实际系统中由于信道延迟和时钟偏差峰值位置是未知的需要定时恢复环路来跟踪。常用的算法有 Gardner 定时误差检测、早迟门同步器等。这类算法都依赖于匹配滤波输出在正确采样点附近的对称性或斜率特性。7.3 有色噪声下的广义匹配滤波如果噪声不是白噪声而是有色噪声那么直接匹配滤波并不是最优的。正确做法是先用一个白化滤波器把噪声变成白色再做匹配滤波也就是“白化滤波器 匹配滤波器”级联。频域表达式为[ H_{opt}(f) \frac{S^*(f)}{P_n(f)} e^{-j2\pi f t_0} ]其中 P_n(f) 是噪声功率谱密度。这个优化滤波器也叫广义匹配滤波。7.4 失配滤波的工程折衷有时候我们并不希望严格使用匹配滤波。比如在雷达脉冲压缩中严格匹配滤波会产生较高的距离旁瓣可能淹没相邻弱目标。这时会人为设计失配滤波器通过牺牲少量 SNR 来压低旁瓣。这是一种典型的工程折衷理论最优并不等于系统最优。通信系统里还有带外辐射、抗干扰、实现复杂度等约束这些都可能促使设计者偏离“纯匹配滤波”。7.5 频域快速实现匹配滤波的时域卷积复杂度是 O(N²)当滤波器很长时可以用 FFT 快速卷积降到 O(N log N)。实现思路def matched_filter_fft(r, s): N len(r) len(s) - 1 R np.fft.fft(r, N) H np.fft.fft(s[::-1], N) y np.fft.ifft(R * H) return y在实时信号处理和高吞吐量场景下基于 FFT 的频域匹配滤波几乎是标配。7.6 数值仿真建议做匹配滤波仿真时有几个容易被忽略的细节设置随机种子保证实验可复现比较滤波器增益时先做归一化避免增益差异影响 SNR 结论计算 SNR 时用信号分量峰值对应位置而不是固定索引如果使用实数带通信号噪声方差要按带通噪声等效带宽换算记录蒙特卡洛的均值和方差不要只看一次结果。8. 总结与下一步怎么学匹配滤波的推导核心只有两步写出输出 SNR 表达式利用柯西-施瓦茨不等式得到上界并找到等号成立条件。等号成立条件就是冲激响应等于信号的时间镜像频域上等价于幅度匹配、相位共轭。这个操作让信号各频率分量在采样时刻相干叠加而噪声只能随机涨落因此 SNR 达到最大。如果要给这篇文章画一个“记忆锚点”可以记住这句话匹配滤波并不是在“滤掉”噪声而是在“对齐”信号。积分里发生的不是简单的相加而是让所有信号成分手拉手同向走让噪声自己互相打架。下一步可以继续学习匹配滤波与相关接收在数字通信中的等价性根升余弦滤波器设计以及为什么收发端各放一个 RRC 是最佳组合广义匹配滤波在有色噪声下的推导FFT 快速卷积在软件无线电接收机中的实现雷达脉冲压缩中的失配滤波与旁瓣抑制。如果这篇文章对你有帮助建议收藏备用也欢迎在评论区聊聊你在仿真或工程实现中遇到的匹配滤波问题。下一篇“图说通信原理”我们继续从图解角度拆解通信系统中的基础模块。
匹配滤波原理详解：从SNR推导到相位校正与相干累加
通信接收机里匹配滤波Matched Filter是一个非常基础却又很容易被“背下来就忘”的概念。很多同学知道结论当滤波器冲激响应是发射信号的时间镜像时输出信噪比 SNR 最大。但如果你追问一句“为什么镜像一下就能让 SNR 最大积分公式里到底发生了什么”不少人会卡住。本文不打算停留在公式结论上而是从 SNR 的数学定义出发一步步拆解匹配滤波的推导过程再用 Python 仿真把“积分里发生的事”可视化。读完你会明白匹配滤波的本质不是“滤波”而是相位校正 相干累加。如果你是通信工程、电子信息专业的学生或者正在做信号处理、雷达脉冲压缩、软件无线电相关开发这篇文章值得收藏。1. 为什么接收机需要匹配滤波1.1 从接收端的“先天不足”说起信号经过信道到达接收机时通常已经被衰减、叠加噪声、受到干扰。接收机要做的第一件事往往不是解调而是在噪声中尽可能把信号“捞”出来。问题是接收机并不知道噪声长什么样。它只能用“信号自身的特点”来判断哪些成分是信号、哪些成分是噪声。匹配滤波就是利用信号已知的波形信息设计一个线性滤波器使得在某个采样时刻输出信号分量尽量大、噪声分量尽量小。这个“尽量大 / 尽量小”的量化指标就是输出信噪比。1.2 匹配滤波到底解决什么问题用一句话概括在白噪声背景下给定一个已知波形 s(t)在所有线性时不变滤波器中匹配滤波器能使输出在采样时刻的瞬时信噪比达到最大。这里的“匹配”二字指的是滤波器与发射信号波形匹配而不是与信道匹配。常见应用场景包括数字通信接收机中的最佳接收机前端雷达和声呐中的脉冲压缩扩频通信中的解扩处理生物医学信号检测、地震信号处理等领域。可以说凡是需要在噪声中检测已知波形的场景匹配滤波都是绕不开的理论基石。1.3 本文的讲解路径下面先建立信号模型和 SNR 的数学表达式然后通过柯西-施瓦茨不等式推导匹配滤波器的冲激响应再从时域和频域两个角度解释“积分里发生了什么”最后用一段完整的 Python 仿真做直观对比。2. 先建立信号模型与 SNR 定义2.1 接收信号模型假设发送信号为 s(t)持续时间为 T接收信号为[ r(t) s(t) n(t) ]其中 n(t) 是加性高斯白噪声AWGN双边功率谱密度为 N₀/2。这里先不考虑信道衰落、多径和干扰聚焦在最基本的模型上。信号 s(t) 的能量记为[ E_s \int_{0}^{T} s^2(t) , dt ]在后面的推导中E_s 会反复出现它决定了匹配滤波所能达到的 SNR 上限。2.2 线性滤波器输出接收信号 r(t) 通过一个冲激响应为 h(t) 的线性时不变滤波器后输出为[ y(t) r(t) * h(t) \int_{-\infty}^{\infty} r(\tau) h(t - \tau) , d\tau ]由于线性滤波器满足叠加定理输出也可以拆成信号分量和噪声分量[ y(t) y_s(t) y_n(t) ]其中[ y_s(t) \int_{-\infty}^{\infty} s(\tau) h(t - \tau) , d\tau ][ y_n(t) \int_{-\infty}^{\infty} n(\tau) h(t - \tau) , d\tau ]我们关心的是在某个采样时刻 t₀输出信噪比能否达到最大。2.3 SNR 怎么定义通信原理中匹配滤波讨论的 SNR 通常是瞬时功率比即采样时刻的信号瞬时功率与噪声平均功率之比[ \text{SNR}_o \frac{|y_s(t_0)|^2}{E[|y_n(t_0)|^2]} ]这里分子是信号分量在 t₀ 时刻的瞬时功率分母是噪声分量的平均功率方差。接下来要做的事情就是把这个表达式展开然后看看 h(t) 满足什么条件时 SNR 最大。3. 从 SNR 表达式推导匹配滤波器3.1 计算输出信号分量信号分量在 t₀ 时刻的值为[ y_s(t_0) \int_{-\infty}^{\infty} s(\tau) h(t_0 - \tau) , d\tau ]如果滤波器只在 [0, T] 区间内有定义发送信号也只在 [0, T] 内有能量那么积分区间可以缩窄到信号的支撑区间但为了推导方便暂时保留完整的积分形式。关键是这是一个内积的形式。信号 s(τ) 和“翻转平移后的冲激响应” h(t₀ - τ) 在 τ 轴上做内积。内积越大说明两个函数越“相似”。3.2 计算输出噪声功率噪声 n(t) 是零均值平稳白噪声通过线性滤波器后输出噪声的方差为[ E[|y_n(t_0)|^2] \frac{N_0}{2} \int_{-\infty}^{\infty} |h(t_0 - \tau)|^2 , d\tau ]这里用到了白噪声通过线性系统后的功率计算公式。噪声的随机相位不会被任何线性系统“对齐”所以最终只剩下功率意义上的累加。3.3 代入 SNR 表达式将分子分母代入得到[ \text{SNR}_o \frac{\left| \int s(\tau) h(t_0 - \tau) , d\tau \right|^2}{\frac{N_0}{2} \int |h(t_0 - \tau)|^2 , d\tau} ]现在问题变成在信号 s(t) 已知的情况下选择怎样的 h(t)能让这个比值最大3.4 柯西-施瓦茨不等式登场柯西-施瓦茨不等式告诉我们对于任意两个平方可积函数 f(τ) 和 g(τ)有[ \left| \int f(\tau) g(\tau) , d\tau \right|^2 \le \int |f(\tau)|^2 , d\tau \cdot \int |g(\tau)|^2 , d\tau ]令[ f(\tau) s(\tau), \quad g(\tau) h(t_0 - \tau) ]则分子被上界限制为[ \left| \int s(\tau) h(t_0 - \tau) , d\tau \right|^2 \le \int s^2(\tau) , d\tau \cdot \int h^2(t_0 - \tau) , d\tau ]于是[ \text{SNR}_o \le \frac{ \int s^2(\tau) , d\tau \cdot \int h^2(t_0 - \tau) , d\tau }{ \frac{N_0}{2} \int h^2(t_0 - \tau) , d\tau } \frac{2 E_s}{N_0} ]注意分母中的积分和分子中的第二个因子完全约掉了。这个上界只取决于信号能量 E_s 和噪声功率谱密度 N₀与滤波器具体形状无关。3.5 等号成立条件时间镜像柯西-施瓦茨不等式取等号的条件是两个函数成比例即[ h(t_0 - \tau) c \cdot s(\tau) ]其中 c 是任意非零常数。换成以 t 为自变量就是[ h(t) c \cdot s(t_0 - t) ]也就是说匹配滤波器的冲激响应是发送信号的时间镜像再平移 t₀。如果取 t₀ T则[ h(t) c \cdot s(T - t) ]忽略常数增益 c 后匹配滤波器的冲激响应就是发射信号在时域上“倒过来”。这时最大输出信噪比为[ \text{SNR}_{o,\max} \frac{2 E_s}{N_0} ]这个结果非常漂亮最大输出 SNR 与信号波形无关只与信号能量和噪声功率谱密度有关。也就是说只要信号能量相同无论你用矩形脉冲、升余弦脉冲还是线性调频信号匹配滤波后的理论最大 SNR 都一样。3.6 频域视角相位共轭与幅度匹配对 h(t) c s(t₀ - t) 做傅里叶变换可以得到匹配滤波器的频率响应[ H(f) c S^*(f) e^{-j 2\pi f t_0} ]其中 S(f) 是 s(t) 的傅里叶变换。这个表达式包含三层含义幅度谱|H(f)| c |S(f)|滤波器的幅度增益与信号幅度谱成正比。信号能量强的频率分量滤波器的增益也大。相位谱∠H(f) -∠S(f) - 2π f t₀滤波器的相位正好抵消了信号本身的相位加上一个线性相位延迟。相位对齐经过匹配滤波后信号各频率分量的相位被“校正”到同一个方向在 t₀ 时刻相干叠加。这正是匹配滤波能在噪声中提升 SNR 的关键信号分量在频域被逐频点对齐而噪声没有这种结构无法被同样对齐。4. 积分里到底发生了什么推导过程虽然简洁但很多人还是觉得“倒过来再卷积”这个操作很抽象。下面我们从离散时间信号的角度看看积分式内部究竟发生了什么事。4.1 把积分写成离散求和假设对信号以采样间隔 Δt 采样得到 N 个样本 s[0], s[1], …, s[N-1]。匹配滤波器的冲激响应样本为[ h[k] s[N-1-k] ]也就是把发射信号的首尾倒置。滤波器输出在采样时刻 n N-1 的值为[ y[N-1] \sum_{m0}^{N-1} r[m] h[N-1-m] \sum_{m0}^{N-1} r[m] s[m] ]如果不看常数增益这个式子等价于[ y[N-1] \sum_{m0}^{N-1} (s[m] n[m]) \cdot s[m] \sum_{m0}^{N-1} s^2[m] \sum_{m0}^{N-1} n[m] s[m] ]信号那一项每一项都是 s[m] 乘以它自己符号完全相同所以累加时同向增强累加结果等于信号能量。噪声那一项n[m] s[m] 的符号随机累加时正负相消相当于随机游走累加结果的幅度增长远慢于信号项。4.2 为什么镜像这么重要如果不做时间镜像直接用普通滤波器 h[k] s[k]没有翻转那么输出在采样时刻的求和为[ y[N-1] \sum_{m0}^{N-1} r[N-1-m] s[m] ]信号分量变成[ \sum_{m0}^{N-1} s[N-1-m] s[m] ]这是自相关函数的取值。如果信号具有良好的自相关特性这个值也可能比较大但通常不如匹配滤波后的能量累加结果那么大。更重要的是当信号带有频率调制比如线性调频直接让信号与自己“不翻转”地做内积会出现严重的相位对消。这正是非匹配滤波器的典型问题。4.3 一张 ASCII 示意图我们用一个长度为 4 的离散信号来演示匹配滤波的累加过程s [s0, s1, s2, s3] h [s3, s2, s1, s0] # 时间镜像 r [r0, r1, r2, r3] # 接收样本 s n 输出 y[3] r0*s3 r1*s2 r2*s1 r3*s0如果 r s无噪声则y[3] s0*s3 s1*s2 s2*s1 s3*s0注意 s0s3 和 s3s0 是两项相同的乘积。镜像操作本质上让信号在卷积过程中“自己匹配自己”所有符号一致的项被安排到同一个求和时刻。如果信号本身是常数矩形脉冲比如 s [1,1,1,1]那么y[3] 1 1 1 1 4所有项同相相加输出达到最大。如果信号是带正负交替的序列比如 s [1, -1, 1, -1]不做翻转时求和可能出现对消。翻转后由于滤波器与信号结构一一对应输出仍然能积累出信号能量。4.4 从频域看积分里的“相位对齐”是怎么回事回到频域。匹配滤波器的输出频谱为[ Y_s(f) S(f) H(f) c |S(f)|^2 e^{-j 2\pi f t_0} ]在 t₀ 时刻做傅里叶逆变换得到[ y_s(t_0) c \int |S(f)|^2 , df ]积分式里每一个频率分量的相位都被 e^{-j2πft₀} 校正到同一个方向所以最终求和在频率轴上也是同相叠加。噪声情况完全不同。白噪声的相位在频域是随机分布的经过 H(f) 后虽然幅度被加权但相位依然是随机散布的不可能在某个时刻相干叠加。于是信号按幅度累加、噪声按功率累加SNR 自然被推高。这就是“积分里发生的那件事”匹配滤波器通过相位共轭把信号在所有频率上的贡献对齐到采样时刻噪声则因为没有可预测的相位结构只能在功率意义上随机增长。5. 完整 Python 仿真匹配滤波 vs 普通矩形窗滤波为了把上面的理论变成直观结果下面用 Python 做一个完整的蒙特卡洛仿真。我们对比两种滤波器匹配滤波器冲激响应为发射信号的时间镜像普通矩形窗滤波器冲激响应为等幅矩形窗长度与信号相同。发射信号选用线性调频信号LFM因为 LFM 信号瞬时频率随时间变化普通矩形窗滤波器会产生严重的相位对消对比效果非常明显。5.1 仿真场景设计与参数仿真参数如下采样率 fs 1000 Hz信号时长 T 0.1 s样本点数 N 100信号中心频率 f0 200 Hz调频带宽 B 100 Hz信号幅度 A 1离散噪声方差 σ² 0.05蒙特卡洛次数 M 5000。信号表达式[ s(t) A \cos\left(2\pi f_0 t \pi \frac{B}{T} t^2\right), \quad 0 \le t \le T ]5.2 完整代码import numpy as np import matplotlib.pyplot as plt # 设置随机种子保证实验可复现 np.random.seed(42) # ------------------------------ # 1. 生成线性调频信号 # ------------------------------ fs 1000.0 # 采样率 Hz T 0.1 # 信号时长 s N int(T * fs) # 样本点数 t np.arange(N) / fs f0 200.0 # 起始频率 Hz B 100.0 # 调频带宽 Hz A 1.0 # 信号幅度 s A * np.cos(2 * np.pi * f0 * t np.pi * (B / T) * t**2) # 信号能量离散定义 Es np.sum(s**2) # ------------------------------ # 2. 定义两种滤波器 # ------------------------------ # 匹配滤波器信号时间镜像 h_mf s[::-1].copy() # 矩形窗滤波器等幅矩形窗 h_rect np.ones(N) # 归一化让滤波器的噪声增益一致便于比较 # 这里采用输出信号能量不变的比例方式 h_mf h_mf / np.sqrt(np.sum(h_mf**2)) h_rect h_rect / np.sqrt(np.sum(h_rect**2)) # ------------------------------ # 3. 理论 SNR 对比 # ------------------------------ sigma2 0.05 # 离散白噪声方差 # 匹配滤波理论输出 SNR snr_mf_theory Es / sigma2 # 矩形窗滤波理论输出 SNR # 信号分量sum(s)因为矩形窗取等幅 sig_rect np.sum(s) # 噪声分量矩形窗单位能量归一化后噪声方差仍为 sigma2 # 但信号分量幅度不同因此 SNR 为 snr_rect_theory sig_rect**2 / sigma2 print(理论 SNR) print(f匹配滤波 : {10*np.log10(snr_mf_theory):.2f} dB) print(f矩形窗滤波 : {10*np.log10(snr_rect_theory):.2f} dB) # ------------------------------ # 4. 蒙特卡洛仿真 # ------------------------------ M 5000 snr_mf_est np.zeros(M) snr_rect_est np.zeros(M) for i in range(M): n np.random.normal(0, np.sqrt(sigma2), N) r s n # 匹配滤波输出 y_mf np.convolve(r, h_mf) # 在峰值采样点取值匹配滤波峰值在末尾 y_mf_signal np.convolve(s, h_mf) peak_idx_mf np.argmax(np.abs(y_mf_signal)) out_mf y_mf[peak_idx_mf] out_sig_mf y_mf_signal[peak_idx_mf] noise_mf out_mf - out_sig_mf # 矩形窗滤波输出 y_rect np.convolve(r, h_rect) y_rect_signal np.convolve(s, h_rect) peak_idx_rect np.argmax(np.abs(y_rect_signal)) out_rect y_rect[peak_idx_rect] out_sig_rect y_rect_signal[peak_idx_rect] noise_rect out_rect - out_sig_rect # 计算每个实验中的瞬时 SNR snr_mf_est[i] out_sig_mf**2 / (noise_mf**2 1e-12) snr_rect_est[i] out_sig_rect**2 / (noise_rect**2 1e-12) print(\n蒙特卡洛仿真 SNR平均值) print(f匹配滤波 : {10*np.log10(np.mean(snr_mf_est)):.2f} dB) print(f矩形窗滤波 : {10*np.log10(np.mean(snr_rect_est)):.2f} dB) # ------------------------------ # 5. 绘图 # ------------------------------ plt.figure(figsize(12, 5)) plt.subplot(1, 2, 1) plt.plot(t, s, label发送信号 s(t)) plt.plot(t, h_mf[:N], label匹配滤波器 h(t), linestyle--) plt.xlabel(时间 (s)) plt.ylabel(幅度) plt.title(LFM 信号与匹配滤波器) plt.legend() plt.grid(True) plt.subplot(1, 2, 2) # 展示一次仿真的滤波输出 np.random.seed(100) n np.random.normal(0, np.sqrt(sigma2), N) r s n y_mf np.convolve(r, h_mf) y_rect np.convolve(r, h_rect) t_mf np.arange(len(y_mf)) / fs t_rect np.arange(len(y_rect)) / fs plt.plot(t_mf, y_mf, label匹配滤波输出) plt.plot(t_rect, y_rect, label矩形窗滤波输出, alpha0.8) plt.xlabel(时间 (s)) plt.ylabel(输出幅度) plt.title(一次仿真中的滤波输出对比) plt.legend() plt.grid(True) plt.tight_layout() plt.savefig(matched_filter_demo.png, dpi150) plt.show()5.3 关键代码解释第 1 段生成 LFM 信号。np.arange(N) / fs 得到时间轴信号表达式就是上文的 s(t)。第 2 段定义两种滤波器。匹配滤波器的核心代码只有一行h_mf s[::-1].copy()这就是“时间镜像”。注意使用 copy() 是为了避免后续归一化修改影响原信号。第 3 段归一化。这里把滤波器都归一化为单位能量这样两种滤波器的噪声增益一致比较 SNR 时更公平。第 4 段蒙特卡洛循环。每次生成一组噪声叠加到信号上分别卷积找到信号分量峰值对应的采样点然后在该点分离信号和噪声计算瞬时 SNR。第 5 段绘图保存到当前目录方便直接查看效果。5.4 运行结果与对比由于设置了随机种子每次运行结果一致。理论值大致如下滤波器类型理论 SNR蒙特卡洛仿真 SNR匹配滤波约 26.99 dB接近 26.9 dB矩形窗滤波通常显著低于匹配滤波与理论一致为什么矩形窗滤波这么差因为 LFM 信号瞬时频率在变化直接用等幅矩形窗做滤波时信号正半周和负半周在积分窗口内大量对消。积分出来的信号分量很小甚至接近于零。而匹配滤波器因为提前把信号做了时间镜像卷积时 LFM 信号的“上升沿”和“下降沿”被精确对准每个样本都在同相叠加。这就是第 4 节说的“积分里发生的那件事”不是矩形窗不够好而是它没有对信号做相位校正。5.5 如果想看更极端的对比如果把 LFM 的带宽加大或者中心频率提高矩形窗滤波的信号分量会进一步趋近于零SNR 损失会更大。匹配滤波则始终保持理论最优值。反过来如果发送信号是常数矩形脉冲匹配滤波器和矩形窗滤波器几乎一样两者 SNR 相同。这也从侧面说明了为什么很多初学者在简单波形下感受不到匹配滤波的“威力”——因为波形太简单相位结构不复杂。6. 常见问题与误区排查问题现象常见原因解决思路把匹配滤波器写成 h(t)s(t)没有做时间镜像记错了公式混淆了相关与卷积匹配滤波的冲激响应是 s(t₀-t)离散实现时用 s[::-1]采样时刻取错SNR 达不到理论值匹配滤波峰值出现在 t₀不在 0 时刻数值仿真中用 argmax 寻找信号峰值位置工程中用符号同步复基带信号直接用实数翻转实现复信号匹配滤波需要共轭翻转即 s^*(t₀-t)复数域中 h np.conj(s[::-1])白噪声环境下用普通滤波器误以为也能达到最大 SNR忽略了信号相位结构对积分的影响用 LFM 等相位调制信号做仿真对比观察 SNR 差距把匹配滤波和相关接收当成两套独立方案两者本质上等价认识到相关接收是匹配滤波的一种等效实现有色噪声下直接用匹配滤波匹配滤波只针对白噪声最优需要先白化再做匹配滤波形成广义匹配滤波仿真中滤波器长度不对输出峰值位置偏移卷积输出长度是 NN-1峰值位置需要计算用 signal 与 h 卷积后的 argmax 定位6.1 一个容易混淆的细节卷积还是相关很多教材说“匹配滤波输出等于接收信号与发送信号做互相关”。这句话容易让人疑惑前面明明说 h(t) 是 s(t) 的镜像卷积怎么又变成相关了其实两者是同一个操作。卷积的定义是“翻转、平移、相乘、求和”而相关是“平移、相乘、求和”。由于匹配滤波器的冲激响应本身就是信号的镜像卷积中的“翻转”被 h 自带的反转抵消了所以最终效果等价于直接做相关。代码层面更容易理解# 卷积方式 y1 np.convolve(r, s[::-1]) # 相关方式等价 y2 np.correlate(r, s, modefull)在无噪声或高 SNR 场景下两种方式得到的峰值位置和形状几乎一致。6.2 采样时刻的敏感性匹配滤波输出在 t₀ 时刻达到峰值。如果采样时刻偏了输出 SNR 会下降。工程上需要定时恢复、符号同步来保证采样点对准峰值位置。对持续时间为 T 的信号如果取 t₀ T则匹配滤波器 h(t) s(T-t)。如果信号有延迟 τ那么 h(t) s(T - (t-τ))也就是还要把延迟信息考虑进去。7. 工程实践中的最佳实践与延伸7.1 匹配滤波与脉冲成形的联合设计在数字通信系统中发送端通常会做脉冲成形接收端也要做匹配滤波。为了在 AWGN 信道下达到最佳性能通常采用根升余弦RRC滤波器发送端一个 RRC接收端一个 RRC级联起来等效为升余弦RC滤波器既能满足奈奎斯特无码间干扰条件又能让接收端 RRC 充当匹配滤波器的角色。这是一个非常经典的“一举两得”设计也是匹配滤波在工程中最常见的落地点。理解匹配滤波后再回头看 RRC 的收发滤波器设计会有完全不同的体会。7.2 采样同步与定时恢复匹配滤波输出是一串脉冲形状的波形符号判决需要在其峰值处采样。实际系统中由于信道延迟和时钟偏差峰值位置是未知的需要定时恢复环路来跟踪。常用的算法有 Gardner 定时误差检测、早迟门同步器等。这类算法都依赖于匹配滤波输出在正确采样点附近的对称性或斜率特性。7.3 有色噪声下的广义匹配滤波如果噪声不是白噪声而是有色噪声那么直接匹配滤波并不是最优的。正确做法是先用一个白化滤波器把噪声变成白色再做匹配滤波也就是“白化滤波器 匹配滤波器”级联。频域表达式为[ H_{opt}(f) \frac{S^*(f)}{P_n(f)} e^{-j2\pi f t_0} ]其中 P_n(f) 是噪声功率谱密度。这个优化滤波器也叫广义匹配滤波。7.4 失配滤波的工程折衷有时候我们并不希望严格使用匹配滤波。比如在雷达脉冲压缩中严格匹配滤波会产生较高的距离旁瓣可能淹没相邻弱目标。这时会人为设计失配滤波器通过牺牲少量 SNR 来压低旁瓣。这是一种典型的工程折衷理论最优并不等于系统最优。通信系统里还有带外辐射、抗干扰、实现复杂度等约束这些都可能促使设计者偏离“纯匹配滤波”。7.5 频域快速实现匹配滤波的时域卷积复杂度是 O(N²)当滤波器很长时可以用 FFT 快速卷积降到 O(N log N)。实现思路def matched_filter_fft(r, s): N len(r) len(s) - 1 R np.fft.fft(r, N) H np.fft.fft(s[::-1], N) y np.fft.ifft(R * H) return y在实时信号处理和高吞吐量场景下基于 FFT 的频域匹配滤波几乎是标配。7.6 数值仿真建议做匹配滤波仿真时有几个容易被忽略的细节设置随机种子保证实验可复现比较滤波器增益时先做归一化避免增益差异影响 SNR 结论计算 SNR 时用信号分量峰值对应位置而不是固定索引如果使用实数带通信号噪声方差要按带通噪声等效带宽换算记录蒙特卡洛的均值和方差不要只看一次结果。8. 总结与下一步怎么学匹配滤波的推导核心只有两步写出输出 SNR 表达式利用柯西-施瓦茨不等式得到上界并找到等号成立条件。等号成立条件就是冲激响应等于信号的时间镜像频域上等价于幅度匹配、相位共轭。这个操作让信号各频率分量在采样时刻相干叠加而噪声只能随机涨落因此 SNR 达到最大。如果要给这篇文章画一个“记忆锚点”可以记住这句话匹配滤波并不是在“滤掉”噪声而是在“对齐”信号。积分里发生的不是简单的相加而是让所有信号成分手拉手同向走让噪声自己互相打架。下一步可以继续学习匹配滤波与相关接收在数字通信中的等价性根升余弦滤波器设计以及为什么收发端各放一个 RRC 是最佳组合广义匹配滤波在有色噪声下的推导FFT 快速卷积在软件无线电接收机中的实现雷达脉冲压缩中的失配滤波与旁瓣抑制。如果这篇文章对你有帮助建议收藏备用也欢迎在评论区聊聊你在仿真或工程实现中遇到的匹配滤波问题。下一篇“图说通信原理”我们继续从图解角度拆解通信系统中的基础模块。
