每日精学 127｜线性插值的幅频下垂

2026-09-25

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

两点线性插值虽能精确再现常数和斜坡，但其频率响应存在由分数位置决定的高频幅度下垂。

01
模型与符号

插值权重 0≤μ≤1 固定，输出采用 (1−μ)x[n]+μ x[n−1]。

μ：分数位置（sample）

H(ω)：插值响应（无量纲）

D0：低频延迟（sample）

Aπ：Nyquist幅度（无量纲）

02
从模型到公式

两点模型

直线插值对应

H(\omega)=1-\mu+\mu e^{-\mathrm j\omega}.

在直流处 H(0)=1，一阶展开为

H(\omega)=1-\mathrm j\mu\omega+O(\omega^2)

因此低频延迟为 D0=μ。


幅度平方

乘以复共轭并收集余弦项，得到

|H(\omega)|^2=(1-\mu)^2+\mu^2+2\mu(1-\mu)\cos\omega
=1-4\mu(1-\mu)\sin^2(\omega/2).

最严重的下垂发生于 μ=1/2。

边界检查

Nyquist频率幅度为

A_\pi=|1-2\mu|.

当 μ=0 或1时退化为整数时延，不再下垂；当 μ=1/2 时

H=e^{-\mathrm j\omega/2}\cos(\omega/2)

在 ω=π 处为零。


幅度下垂的二阶展开

两点线性插值的响应为

H=(1-\mu)+\mu e^{-\mathrm j\omega}

将模平方展开并整理，可得到

|H|^2=1-4\mu(1-\mu)\sin^2(\omega/2).

在低频利用正弦和平方根展开，得到

|H|=1-\frac12\mu(1-\mu)\omega^2+O(\omega^4).

因此直流增益准确不意味着整个通带幅度准确；

误差的二阶系数在半样本插值时最大，在整数端点消失。特别地，半样本时 Nyquist 响应为零，这不是数值舍入而是两点加权产生的相消。用更多输出小数位只能减小算术误差，不能消除这一模型幅频下垂。


03
物理含义

两个邻近样本在高频可能具有相反相位，加权平均会抵消高频分量；这不是提高采样率后自动消失的误差。

04
固定参数算例

对分数位置0、0.25、0.5、1计算解析与直接频率响应，验证Nyquist端点和线性序列的精确插值。

固定参数数值验证结果
量数值
分数位置0.25
magnitude identity error3.33066907388e-16
linear reproduction error0
nyquist mu01
nyquist mu0250.5
nyquist mu056.12323399574e-17
nyquist mu11
图 1：线性插值的幅频下垂的理论关系与固定参数验证图 1：线性插值的幅频下垂的理论关系与固定参数验证

图 1　线性插值的幅频下垂的理论关系与固定参数验证

表中数值和图中曲线由 MATLAB 实际运行得到，随机种子为 20261046；运行版本与完整结果另存为独立文件。

05
数字实现的边界

RTL为固定 μ=1/4 的两点插值分子 3a+b，a、b为当前和前一有效样本，由外部统一除以4。 RTL 采用同步有效信号和异步低有效复位；

空拍保持状态且输出无效。除另有说明外，输入为16位有符号整数，结果为64位有符号整数，测试范围内不截断或饱和。输出在采样时钟上升沿更新。自检覆盖连续样本、空拍和复位；

未进行器件布局布线、时序收敛或功耗验证。

独立自检程序已在 Icarus Verilog 13.0 中实际运行并通过。该结果不包含器件综合、静态时序分析或板级测试。

06
适用范围

线性插值可用于低带宽波形的细时间对齐和查表值之间的连续过渡。其优点是权重简单且状态少，但幅度下垂随输入带宽和分数位置变化。对高频精度有要求的链路，应把解析下垂计入误差预算，而不是只检查插值时刻。

07
工程上需要注意

除半采样等特殊位置外，相位并非全带线性。

补偿下垂需要额外滤波且不能恢复精确零点处的信息。

时变分数位置还可能引入幅度调制。

两点线性插值虽能精确再现常数和斜坡，但其频率响应存在由分数位置决定的高频幅度下垂。

内容依据：每日精学第 127 课《线性插值的幅频下垂》。图表沿用原课固定参数数值结果，属于数值验证，不代表设备实测。

阅读8
分享
写评论...