双线性变换 - Ciallo Wiki

1 概述与定义

双线性变换（bilinear transform）是一种把连续时间域的复变量表示（常见为 \(s\) 平面）转换为离散时间域等价表示（常见为 \(z\) 平面）的标准方法。它通过一个线性分式变换，将连续域中的微分算子与离散域中的差分算子联系起来，从而把连续系统的传递函数或状态方程，系统性地离散化到数字实现所用的形式。

1.1 基本思想：连续域到离散域的映射

在连续系统中，导数算子与复变量 \(s\) 对应；在离散系统中，移位与差分操作与 \(z\) 变量相关。双线性变换的核心思路是：用一个“非线性频率标尺关系”的映射，把连续域的 \(s\) 写成 \(z\) 的函数，再把原系统的数学表达式替换到离散域，从而得到离散系统模型。

由于这种替换是系统性的代数操作，工程上通常不需要直接对微分方程进行数值积分，而是从频域/传递函数出发完成离散化。

1.2 典型形式：\(s\)-平面与 \(z\)-平面的关联

最常用的双线性映射是 \[ s=\frac{2}{T}\frac{z-1}{z+1}, \] 其中 \(T\) 为采样周期。该式等价地把单位圆与复平面中的稳定区域、虚轴等几何对象联系起来：当 \(z\) 在单位圆上变化时，\(s\) 落在与频率相关的曲线上；当 \(z\) 在单位圆内部变化时，\(s\) 落在左半平面。

在信号处理语境中，这一变换也常被写成以 \(z^{-1}\) 表示的差分形式，便于得到实现用的差分方程。

1.3 数学要点：线性分式变换与可逆性条件

双线性变换本质上是线性分式变换（Möbius 变换）的一种特例。线性分式形式意味着：

除了少数使分母为零的点外，映射是可逆的；
映射保持某些几何结构（例如圆与直线在复平面中的对应关系）；
系统的极点与零点会按复平面映射规律相应变换。

可逆性通常要求避开映射中的奇异点：例如在上述形式中，\(z=-1\) 对应 \(z+1=0\)，因此在推导或代数替换时需注意不要把该点引入到不适用的表达式里。同时，当连续模型在映射所覆盖的区域内满足常规可实现性条件时，离散化结果才具有良好的物理与数值解释。

2 频域与映射性质

双线性变换之所以广受使用，除了离散化方便外，更关键在于它对稳定性与频率刻度的“结构性继承”。这些性质可从复平面几何映射与虚轴到单位圆的对应关系得到解释。

2.1 稳定性保持：左半平面到单位圆内部

若连续系统在 \(s\) 平面中“稳定”（通常指所有极点位于左半平面），则使用双线性映射后，离散系统的极点会落在 \(z\) 平面单位圆内部。这种单调的几何保持使得离散化不会破坏稳定性，从而避免某些其他离散化方法可能出现的稳定性丢失问题。

直观理解是：双线性变换把左半平面映射到单位圆内部，而右半平面映射到单位圆外部；稳定区域与发散区域之间由单位圆分隔。

2.2 频率映射：从虚轴到单位圆的对应

连续系统的频率响应常由 \(s=j\omega\) 描述。将 \(s=j\omega\) 代入双线性映射后，可得到单位圆上的 \(z\) 与连续频率 \(\omega\) 的关系。于是，数字系统的归一化频率变量与模拟角频率之间存在确定的非线性对应。

这种对应的直接工程意义是：当你把模拟滤波器的频率特性“搬运”到数字域时，每个模拟频率点不会线性对应到数字频率点，而是按双线性变换的频率标尺发生弯曲映射。

2.3 角频率“预畸变”与双线性压缩

由于从 \(j\omega\) 到单位圆的映射是非线性的，数字域的频率轴与模拟域的频率轴之间通常呈现“压缩/预畸变”关系。工程上常采用“预畸变频率”（pre-warping）的技巧：在设计之前先把目标数字频率反解回模拟频率，再在模拟原型上设计，从而让关键通带或阻带边界在离散域对齐。

简单地说：双线性变换在把频响从模拟域移到数字域时，会对频率进行一种可计算的弯曲重标定；通过预畸变可以抵消在关注频段内的偏差。

2.4 映射对特征值/极点零点的影响

若连续传递函数的极点与零点分别为 \(\{p_i\}\) 与 \(\{z_i\}\)（此处 \(z_i\) 作为零点集合记号与 \(z\) 平面变量字母不混淆），则离散系统的极点与零点由同一个映射规则由 \(\{p_i\}\) 与 \(\{z_i\}\) 对应得到。线性分式变换的代数性质使得：

极点位置按映射曲线迁移；
零点也同样映射；
因此频响形态（例如共振峰、衰减斜率的拐点位置等）会以可预测方式改变。

这也是为什么双线性变换常被用于把模拟滤波器结构“迁移”到数字滤波器时：它不会随意重排系统的关键根结构，而是按规则搬运。

3 在系统与控制中的应用

双线性变换不仅用于滤波器，还用于一般控制系统的连续到离散转换。其常见场景包括将连续模型转为离散状态空间或离散传递函数，并用于数字控制器设计与数字仿真。

3.1 连续系统到离散系统：状态空间离散化

在状态空间形式中，连续系统可写为 \[ \dot{x}(t)=Ax(t)+Bu(t),\quad y(t)=Cx(t)+Du(t). \] 双线性变换可通过把导数算子与差分算子之间的关系替换到等价的离散时间形式，从而得到离散系统矩阵（例如将连续传递函数对应到离散传递函数，再回到状态空间实现）。其优势在于输出与输入的离散化保持了稳定性结构，且计算流程通常比某些需要矩阵指数的离散化更直接。

实际工程中，状态空间离散化可能结合“先求传递函数再映射”的路径，或直接对可实现的算子形式进行替换。

3.2 从传递函数进行变换的常用步骤

将连续传递函数 \[ G(s)=\frac{N(s)}{D(s)} \] 映射为离散传递函数的典型步骤包括：

将映射关系 \(s=\frac{2}{T}\frac{z-1}{z+1}\) 代入 \(G(s)\)；
对分子分母做代数化简，使结果表达成关于 \(z\)（或 \(z^{-1}\)）的有理函数；
必要时进行归一化与因式整理，以获得便于实现的结构（例如按幂次排列，保证因子“落在合理的实现方向”）。

对于多项式阶次较高的情况，代数化简需要注意数值规模，避免系数爆炸。

3.3 离散化误差与等效性分析

双线性变换并非在所有频率都保持严格等效。由于频率映射非线性，离散系统对连续系统的响应在频段之间会出现偏差。误差通常随以下因素变化：

采样周期 \(T\) 的大小（采样越密，越接近连续行为）；
关注的频率范围（尤其是接近折返相关的边界区域时）；
系统阶数与极点分布（影响误差在频域的扩散）。

等效性分析常以频率响应或输入输出误差度量进行；当系统或滤波器被设计为在关键频段对齐（如通过预畸变），误差可显著改善。

3.4 与其他离散化方法（如前向/后向差分）的比较

相较于前向差分或后向差分等方法，双线性变换常被认为具有更好的稳定性保持与更平滑的频率特性继承：

前向/后向差分的频率响应映射与极点迁移通常不同，可能在某些情况下导致稳定性或频率畸变更严重；
双线性变换因把左半平面直接对应到单位圆内部，在稳定性方面更“稳健”。

当然，代价在于其映射的频率非线性需要预畸变或接受频率扭曲；若系统只需局部频率近似，选择何种差分格式会影响实现与精度的折中。

4 在数字滤波器设计中的应用

在数字滤波器领域，双线性变换是一种经典的“模拟原型到数字实现”的转换桥梁。工程流程通常围绕频率映射与边界对齐展开。

4.1 模拟低通原型到数字滤波器的流程

常见流程可概括为：

选取模拟低通原型（如巴特沃斯、切比雪夫或椭圆型等）并确定其阶数与通带/阻带指标；
将所需的数字采样率、通带边界与阻带边界换算为模拟域目标频率（必要时使用预畸变）；
通过低通到数字低通的双线性映射，得到离散域低通传递函数；
若需其他类型滤波器（带通、带阻），再进行相应的模拟域变换组合。

这样做的优点是：可复用成熟的模拟原型设计经验，并在数字域得到可实现的 IIR 结构。

4.2 模拟到数字的带通/带阻变换组合

若目标为带通或带阻滤波器，典型思路是先在模拟域完成低通原型的带通/带阻频率变换，再用双线性映射进入数字域。常见的做法包括：

模拟低通 \(G_{LP}(s)\) 经过频率变换得到 \(G_{BP}(s)\) 或 \(G_{BS}(s)\)；
最终将得到的模拟带通/带阻传递函数用双线性映射替换 \(s\) 并化简。

组合中需要特别注意频率边界在映射后的对应关系，否则带通中心频率与带宽可能出现偏移；预畸变在这里通常更重要。

4.3 冲激响应与等效差分方程的构造

双线性变换得到的是离散传递函数。为了实现滤波器，通常要写成差分方程或结构形式。对一般有理函数 \[ H(z)=\frac{B(z)}{A(z)}, \] 可通过乘以 \(A(z)\) 并使用 \(z^{-1}\) 与延时算子的关系，把输出写成输入与历史输出的线性组合，从而得到对应的差分方程。冲激响应则由在离散系统中对单位冲激输入进行迭代计算得到，或由传递函数展开系数求得。

在工程实现中，直接将离散传递函数转为差分方程时，往往选用更数值稳健的实现形式（见后文实现细节）。

4.4 频响特性继承：通带/阻带与幅相表现

滤波器设计目标通常强调幅频特性，如通带纹波、阻带衰减等。双线性变换在一般条件下能够较好继承模拟原型的幅频形状，但由于频率非线性，边界处的实际位置可能发生偏移。通过预畸变让关键频率点对齐后，数字滤波器的通带与阻带特性可以更接近原始指标。

幅相方面，双线性变换对应的相位响应也会随映射改变。工程上若系统对相位或群延迟有特殊要求，设计时需要在目标频段做更精细的对齐与误差评估。

5 实用实现细节

将双线性变换用于具体系统实现时，数值稳定性、结构选择与频率单位处理非常关键。许多问题并非来自理论映射，而来自实现方式与尺度管理。

5.1 数值稳定性与实现形式选择（如直接型/级联型）

得到离散传递函数后，如何实现会显著影响数值表现：

直接型（direct form）可能在系数精度有限时出现数值误差积累；
级联型（cascade of biquads）把高阶系统分解为低阶部分，通常更稳定；
反归一化、对极点零点配对以及归一化处理也能改善数值误差。

在双线性变换后的系数可能跨越较大数量级时，级联结构与适当的缩放通常更可靠。

5.2 参数标定：采样周期对映射的影响

采样周期 \(T\) 决定映射的尺度。改变 \(T\) 会改变同一连续系统在离散域的频率对应关系，也会影响最终离散滤波器的边界频率位置。因此工程设计中必须统一：

连续原型的频率定义方式；
目标数字系统的归一化频率或实际采样率；
单位制（例如“弧度/秒”与“赫兹”的换算）。

标定错误会直接导致滤波器中心频率偏移或通带/阻带错位。

5.3 处理多阶系统与运算复杂度

高阶模型在双线性映射后多项式代数化简与系数计算会变得复杂。常见策略包括：

先求极点与零点，再按映射规则得到离散根结构，最后构造多项式；
使用对多项式因式分解更稳健的数值流程；
在实现端采用分节结构以降低累计误差。

此外，若目标是实时处理，结构的乘加次数与内存需求也会影响实现选择。

5.4 常见陷阱：尺度、单位与频率归一化

实际开发中最常见的错误包括：

忘记把目标频率从赫兹转成弧度/秒，或反过来；
将“归一化数字频率”与“模拟角频率”混用；
预畸变只做了一侧边界或使用了错误的频率反解公式；
在软件工具链中重复应用缩放（例如采样率参数输入方式不同导致重复换算）。

这些陷阱通常会让结果看似“系数正确但指标完全不对”。因此在设计流程的每一步都应核对频率尺度的一致性。

6 误差、局限与扩展

双线性变换具有稳健的稳定性保持与良好的总体频域表现，但其频率非线性映射带来的误差、边界频率限制等仍需理解。相关改进思路也在工程中被广泛讨论。

6.1 近似精度随频率范围的变化

在低频或关注频段内，双线性映射的频率对应关系更接近线性，使得数字系统更能逼近连续系统的响应。但当频率增大、靠近目标边界（尤其接近映射中对频率分辨率不利的区域）时，预期的响应可能出现偏差。

因此，设计时一般不会把“全频带一致性”当作默认前提，而是依据具体任务选择频率范围并配合预畸变。

6.2 对高频与边界频率的局限

由于映射把无限模拟频率压缩到有限的数字频率范围，越高的连续频率在离散域中越容易出现压缩带来的分辨率下降效应。此外，靠近关键边界频率时的误差对最终指标影响更显著，可能导致：

通带边缘更早进入过渡带；
阻带开始衰减的频率偏移；
相位与群延迟偏差加剧。

这也是工程上常强调“关键频率点对齐”的原因。

6.3 加权双线性变换与改进思路（如保通带/保相位取舍）

为改善不同频率段的误差分布，常出现加权或改进型双线性变换：通过引入权重或调整映射参数，让误差更偏向某个目标（例如更准确地保通带纹波形状、或在某些区间保相位特性）。在具体实现上，这类方法会改变映射公式与频率对应关系，从而调整频域误差的形态。

在工程权衡上，“保通带”与“保相位”往往存在取舍：更强调幅度拟合的同时，相位误差可能增加；反之亦然。选择应由应用场景决定。

6.4 面向特定目标的替代方法概览

除了加权双线性变换，离散化与模拟-数字映射领域还有其他替代思路，例如：

更强调时间域或特定响应指标的离散化策略；
在某些结构上保持特定性质（如更精确的幅频匹配或更好的相位保持）的设计流程；
结合数值优化对离散滤波器系数进行直接拟合。

这些方法各自适合不同目标与约束条件，双线性变换通常因其通用性与稳定性而仍然是起点方案。

7 相关概念与术语

理解双线性变换常需借助一些与复平面几何、离散变换以及频率刻度相关的基本概念。

7.1 单位圆、稳定区域与映射图像

单位圆是离散系统分析中最常见的边界：离散系统的稳定性通常与极点是否在单位圆内部有关。双线性变换通过复平面映射把连续系统的左半平面稳定区域映射到单位圆内部，因此常用“映射图像”帮助理解：稳定区域如何转化，以及虚轴对应到单位圆上的哪一段。

在进行可视化或检查结果时，绘制极点迁移到单位圆的示意图能快速判断稳定性是否保持。

7.2 采样、Z变换与差分算子

离散系统中 \(z\) 变量通常由 Z 变换引入，它把序列与复变量对应起来。移位操作对应于乘以 \(z^{-1}\)，从而差分方程可以写成关于 \(z\) 的代数关系。双线性变换所给出的 \(s\)-\(z\) 关联，正是把“导数/微分”这类连续算子与“移位/差分”这类离散算子对应起来。

采样周期 \(T\) 在该对应关系里扮演尺度因子的角色，直接决定频率映射的弯曲程度。

7.3 频率预畸变与模拟-数字对照图

预畸变是把目标数字频率反映回模拟域时使用的补偿步骤。通过对照频率映射曲线（模拟角频率与数字归一化频率之间的函数关系），设计者可计算对应的模拟频率值，使得离散设计后关键边界更准确落在期望位置。

在工程实践中，“对照图 + 关键点反解”是一种效率较高的工作方式。

8 参考与延伸阅读

双线性变换的基础推导、滤波器设计流程与频域性质在许多经典教材与工程手册中都有系统叙述。延伸阅读通常围绕“复变函数映射”“离散滤波器设计”“误差与频率预畸变”等主题展开。

8.1 经典教材与公式来源（按主题归类）

可按主题归类查找资料：

控制与离散化：双线性变换的稳定性保持与系统离散化流程；
数字信号处理：模拟原型到数字滤波器的映射、频率预畸变与设计步骤；
复变函数：线性分式变换的几何性质与映射定理。

不同教材对符号体系与归一化频率约定可能不同，阅读时需对照采样频率与频率变量的定义。

8.2 工程手册中的设计流程示例

工程手册常给出从规格到系数的实际步骤示例，例如：选择原型、确定阶数、计算预畸变频率、进行双线性映射并最终得到可实现的滤波器结构。通过这些示例，读者可以快速建立“从需求到实现”的工程直觉。

8.3 进一步的数学推导与证明资源

进一步推导通常包括：映射的可逆性与奇异点讨论、稳定性映射的严格证明、极点零点迁移的一般形式，以及频率映射中角频率关系的推导。若关注更精细的误差分析，还可查找关于频率响应近似与加权改进方法的数学说明。