每日精学 132｜多相插值器的等效分解

2026-09-25

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

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

整数插值FIR可按系数索引模L拆成L个低采样率支路，使每个输出相位只计算必要的非零乘积。

01
模型与符号

固定插值倍数L，零初始状态；系数不足整组时以零补齐，不改变传递函数。

Er(z)：第r相子滤波器（无量纲）

L：插值倍数（无量纲）

h[k]：原型系数（无量纲）

y[Lm+r]：第r相输出（幅度）

02
从模型到公式

系数分类

把任意索引写成 k=pL+r，得到

H(z)=\sum_{r=0}^{L-1}z^{-r}E_r(z^L),\quad E_r(z)=\sum_ph[pL+r]z^{-p}.

分解只重新排序有限和，不引入近似。

时域支路

对插零输入，输出第r相满足

y[Lm+r]=\sum_ph[pL+r]x[m-p].

每支路处理原采样率输入，输出交织为高采样率序列。

运算节约

原型长度K且可被L整除时，每输出只需K/L次乘法；每个原始输入对应L个输出，总计

C_{\mathrm{poly}}=K,\qquad C_{\mathrm{dense}}=KL

次乘法每原始输入。该比较不把乘以已知零的运算优化预先计入密集基准。

支路直流增益与输出相位纹波

多相分解保持完整卷积，但各支路系数和不一定相同。对原型 [1,2,3,4]/5，两相直流增益为

G_0=(1+3)/5=4/5,\qquad G_1=(2+4)/5=6/5.

常值输入进入稳态后，输出在两个相位交替取这两个值。

全原型系数和等于二仅保证平均幅度为一，并不能单独保证每个插值相位都平坦：

\frac{G_0+G_1}{2}=1,\qquad G_0\ne G_1.

该原型适合验证分解恒等式，不宜仅凭等价性结论视为高质量重建低通。

实际插值滤波器还需检验镜像抑制和分相直流误差；数学分解正确与原型指标合格是相互独立的判断。

03
物理含义

多相分解把本来乘零的项从数据通路中移除，同时保持完整滤波响应和所有输出时刻。

04
固定参数算例

采用L=2、h=[1,2,3,4]/5，比较两个支路交织和显式插零卷积的完整输出，并核对尾部零延拓。

固定参数数值验证结果
量数值
factor2
原型滤波器长度4
输入长度35
输出长度72
equivalence error0
polyphase mult per input4
dense mult per input8
图 1：多相插值器的等效分解的理论关系与固定参数验证图 1：多相插值器的等效分解的理论关系与固定参数验证

图 1　多相插值器的等效分解的理论关系与固定参数验证

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

05
数字实现的边界

RTL提供两相分子并行核：输入a为当前样本、b为前一样本，输出 a+3b 与 2a+4b，外部除以5并交织两路。 RTL 采用同步有效信号和异步低有效复位；

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

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

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

06
适用范围

多相插值器用于提高采样率并抑制频谱镜像。各分支可串行复用乘法器，也可并行产生多个输出相位；这两种资源选择均应保持相同的系数相位编号与输出交织顺序。原型滤波器的增益决定最终幅度恢复。

07
工程上需要注意

输出交织次序错误会造成周期性失真。

支路长度不同必须补零或独立管理尾部。

固定点各支路的缩放必须一致。

整数插值FIR可按系数索引模L拆成L个低采样率支路，使每个输出相位只计算必要的非零乘积。

内容依据：每日精学第 132 课《多相插值器的等效分解》。图表沿用原课固定参数数值结果，属于数值验证，不代表设备实测。

阅读7
分享
写评论...