每日精学 126｜三阶Lagrange分数时延插值

2026-09-25

← 技术专题 · 每日精学目录

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

三阶Lagrange插值利用四个相邻样本构造三次多项式，在低频对理想分数时延实现四项矩匹配。

01
模型与符号

使用节点0、1、2、3，所需延迟为 D；固定系数设计以 0≤ D≤3 为插值区间，重点验证 D=1.5。

D：分数时延（sample）

h_k(D)：插值权重（无量纲）

P(t)：插值多项式（输入幅度）

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

02
从模型到公式

基函数构造

节点基函数满足

$\ell_k(m)=\delta_{km}$

其显式形式为

$\ell_k(t)=\prod_{m=0,m\ne k}^{3}\frac{t-m}{k-m}.$

因此

$P(t)=\sum_{k=0}^{3}x[n-k]\ell_k(t)$

延迟输出取 y[n]=P(D)，系数为 h_k(D)=ℓ_k(D)。

多项式再现

任何次数不超过3的多项式都能精确插值，故

$\sum_{k=0}^{3}h_k(D)k^p=D^p,\qquad p=0,1,2,3.$

将指数展开

$e^{-\mathrm j\omega k}=\sum_p(-\mathrm j\omega k)^p/p!$

前四项与理想延迟一致，得到

$H(\omega)=e^{-\mathrm j\omega D}+O(\omega^4).$

半采样中心

当 D=1.5 时，权重为

系数和为1，且关于1.5对称，因而相位延迟精确为1.5样本，但幅度只在低频接近1。高频误差来自有限长度插值，而非索引错误。

低频误差阶数与多项式余项

四点插值在节点零至三上构造三次多项式，目标位置为 d=3/2。对具有四阶连续导数的函数，插值余项为

$f(d)-p_3(d)=\frac{f^{(4)}(\xi)}{4!}\prod_{k=0}^3(d-k) =\frac{3}{128}f^{(4)}(\xi).$

取复指数并作局部级数分析，误差从四阶频率项开始；

这与系数满足零至三阶矩约束一致：

$\sum_{k=0}^3h_k k^m=d^m,\qquad m=0,1,2,3.$

精确再现三次多项式不表示对所有带限信号都精确。

高频复指数的高阶导数随频率增长，Nyquist 附近误差可显著增加。余项界说明低频精度来源于局部平滑性，而不是分数时延参数本身拥有无限分辨率。

03
物理含义

矩匹配保证低阶时间多项式精确再现，也等价于理想频率响应在直流处若干导数相等。

04
固定参数算例

验证四个矩约束及三次多项式再现；同时比较分数时延频响与理想纯延迟，显示低频准确和高频幅度误差。

固定参数数值验证结果

量	数值
时延 / 样本	1.5
权重和	1
max moment residual	0
cubic reproduction error	0
lowband max error	0.00022511393567
nyquist magnitude	9.18485099361e-17

图 1　三阶Lagrange分数时延插值的理论关系与固定参数验证

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

05
数字实现的边界

RTL实现整数四抽头 [−1,9,9,−1]，输出保留分母16之外的精确分子；外部统一缩放避免各抽头分别截断。 RTL 采用同步有效信号和异步低有效复位；

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

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

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

06
适用范围

Lagrange分数时延可用于采样时刻对齐、异步数据插值和有限带宽的时延扫描。四抽头结构提供固定且较短的数据历史，适合中低带宽近似。若工作带宽接近Nyquist频率，应根据幅度和相位误差重新选择阶数或滤波器设计方法。

07
工程上需要注意

有限阶多项式插值不是全带宽理想sinc插值。

远离节点区间的外推会显著增大权重和误差。

快速变化的延迟参数会形成时变调制。

三阶Lagrange插值利用四个相邻样本构造三次多项式，在低频对理想分数时延实现四项矩匹配。

内容依据：每日精学第 126 课《三阶Lagrange分数时延插值》。图表沿用原课固定参数数值结果，属于数值验证，不代表设备实测。

阅读7

写评论...

每日精学 126｜三阶Lagrange分数时延插值

01模型与符号

02从模型到公式

基函数构造

多项式再现

半采样中心

低频误差阶数与多项式余项

03物理含义

04固定参数算例

05数字实现的边界

06适用范围