每日精学 125｜一阶全通滤波器的群时延

2026-09-25

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

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

一阶全通滤波器通过互为倒数的零极点保持单位幅度，同时产生随频率变化的群时延。

01
模型与符号

实系数 |a|<1，传递函数为 (a+z⁻¹)/(1+az⁻¹)，群时延以样本为单位。

a：全通系数（无量纲）

ϕ：展开相位（rad）

τ_g：群时延（sample）

ω：数字频率（rad/sample）

02
从模型到公式

单位幅度

在单位圆上，分子和分母的模平方均为

$1+a^2+2a\cos\omega.$

所以

$|H(e^{\mathrm j\omega})|^2=1.$

极点为 −a，|a|<1 保证因果稳定。

相位导数

将分子写成

$e^{-\mathrm j\omega}(1+ae^{\mathrm j\omega})$

得到

$\phi(\omega)=-\omega-2\arg(1+ae^{-\mathrm j\omega}).$

对角度求导并化简得到

$\tau_g(\omega)=\frac{1-a^2}{1+a^2+2a\cos\omega}.$

低频延迟与实现

令 ω=0 得 D₀=(1−a)/(1+a)，所以指定低频延迟对应

$a=\frac{1-D_0}{1+D_0}.$

交叉相乘给出

此式只在低频近似恒定延迟，不能把 D₀ 当成全带宽时延。

群时延的频率积分约束

对实参数 |a|<1 的一阶全通，群时延可展开为收敛的余弦级数：

$\tau_g(\omega)=\frac{1-a^2}{1+a^2+2a\cos\omega} =1+2\sum_{m=1}^{\infty}(-a)^m\cos(m\omega).$

全频积分时所有非零谐波消失，故

$\frac1{2\pi}\int_{-\pi}^{\pi}\tau_g(\omega)d\omega=1.$

这说明一阶全通可以重新分配不同频率的延迟，却不能使整个频带的平均群时延任意变小。

对 a=1/2，直流处为三分之一，Nyquist 处为三，二者并不矛盾。接近单位圆的极点使群时延在某一频段集中，对频率网格和相位解缠提出更高要求。

03
物理含义

幅度平坦并不意味着波形只发生刚性平移。不同频率具有不同延迟，宽带脉冲仍会展宽或改变形状。

04
固定参数算例

采用 a=0.5，比较相位数值导数与解析群时延，验证模平方为1及直流群时延1/3样本。

固定参数数值验证结果

量	数值
coefficient	0.5
amplitude max error	4.4408920985e-16
delay derivative error	1.47068502709e-07
dc delay	0.333333333333
nyquist delay	3
冲激响应能量	1

图 1　一阶全通滤波器的群时延的理论关系与固定参数验证

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

05
数字实现的边界

RTL实现 a=1/2 的一阶全通递推，乘除二采用算术移位；有限字长截断使幅度严格为1的数学性质只对精确算术成立。 RTL 采用同步有效信号和异步低有效复位；

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

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

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

06
适用范围

一阶全通可用作低频分数延迟和相位均衡的基本单元。幅度不变便于把相位校正与增益校正分开设计，但必须限定工作频带内的群时延起伏。多个全通级联时，群时延相加而瞬态长度和定点状态风险也随之增加。

07
工程上需要注意

极点接近单位圆时群时延峰值和瞬态都增大。

时变系数不再满足固定LTI全通推导。

反馈舍入误差可能造成极限环。

一阶全通滤波器通过互为倒数的零极点保持单位幅度，同时产生随频率变化的群时延。

内容依据：每日精学第 125 课《一阶全通滤波器的群时延》。图表沿用原课固定参数数值结果，属于数值验证，不代表设备实测。

阅读7

写评论...

每日精学 125｜一阶全通滤波器的群时延

01模型与符号

02从模型到公式

单位幅度

相位导数

低频延迟与实现

群时延的频率积分约束

03物理含义

04固定参数算例

05数字实现的边界

06适用范围