劳斯判据在离散系统中的妙用：一个案例讲透双线性变换

劳斯判据在离散系统中的妙用一个案例讲透双线性变换在数字控制系统的设计与分析中稳定性始终是工程师首要关注的焦点。与连续系统不同离散系统的稳定域从s平面的左半平面转变为z平面的单位圆内这一根本差异使得传统的劳斯判据无法直接应用。本文将聚焦于双线性变换这一关键桥梁通过一个完整的案例演示如何将z域问题转化为w域问题最终利用劳斯判据完成稳定性判定。1. 离散系统稳定性的数学本质离散系统的稳定性分析建立在z变换理论基础上。当采样周期为T时s平面与z平面的映射关系由z e^(sT)定义。这个指数映射将s左半平面稳定区域转换为z平面单位圆内部区域而s平面的虚轴则对应z平面单位圆边界。稳定性充要条件离散系统稳定的充要条件是系统特征方程的所有根都位于z平面单位圆内。数学表达为D(z) 0 的所有根满足 |z_i| 1与连续系统相比离散系统的稳定边界从直线虚轴变为曲线单位圆这使得直接应用劳斯判据变得困难。下表对比了两种系统的稳定性特征特征连续系统离散系统稳定域s左半平面z单位圆内边界虚轴单位圆判据适用性直接适用劳斯判据需双线性变换后适用2. 双线性变换的核心原理双线性变换又称w变换的精妙之处在于它建立了一个从z平面到w平面的可逆映射w (z 1)/(z - 1)这个分式线性变换具有以下关键性质单位圆内部映射z平面单位圆内的区域对应w平面的左半平面单位圆边界映射z平面单位圆上的点对应w平面的虚轴单位圆外部映射z平面单位圆外的区域对应w平面的右半平面变换推导过程从映射关系w (z 1)/(z - 1)出发解出z得到逆变换z (w 1)/(w - 1)将z e^(jθ)代入可验证单位圆边界映射到虚轴注意不同文献可能采用w (z - 1)/(z 1)的变体形式本质相同但左右半平面对应关系相反3. 完整案例分析离散系统稳定性判定考虑一个离散系统的特征方程为D(z) z^3 - 1.2z^2 0.45z - 0.05 03.1 执行双线性变换将z (w 1)/(w - 1)代入特征方程计算各次幂z^3 (w 1)^3 / (w - 1)^3 z^2 (w 1)^2 / (w - 1)^2 z^1 (w 1) / (w - 1)通分后得到(w 1)^3 - 1.2(w 1)^2(w - 1) 0.45(w 1)(w - 1)^2 - 0.05(w - 1)^3 0展开并合并同类项w^3 3w^2 3w 1 -1.2(w^3 w^2 - w - 1) 0.45(w^3 - w^2 - w 1) -0.05(w^3 - 3w^2 3w - 1) 0最终得到w域特征方程0.2w^3 2.1w^2 1.5w 0.2 03.2 构建劳斯表基于变换后的特征方程构建劳斯表0.2w^3 2.1w^2 1.5w 0.2 0劳斯表结构如下w^30.21.5w^22.10.2w^1(2.11.5 - 0.20.2)/2.1 ≈ 1.490w^00.2-计算过程第一行奇次幂系数0.2 1.5第二行偶次幂系数2.1 0.2第三行b1 (a1*a2 - a0*a3)/a1 (2.1*1.5 - 0.2*0.2)/2.1 ≈ 1.49第四行c1 (b1*a3 - a1*0)/b1 a3 0.23.3 稳定性判定观察劳斯表第一列0.2 (正)2.1 (正)1.49 (正)0.2 (正)结论劳斯表第一列无符号变化说明w域系统稳定即原z域系统稳定。4. 实践中的常见问题与技巧4.1 数值稳定性问题当特征方程阶数较高时双线性变换可能导致系数剧烈变化。例如# Python示例高阶多项式变换后的系数对比 import numpy as np # 原始z域多项式系数 z_coeff [1, -2.5, 2.35, -0.925, 0.15] # 执行双线性变换需实现变换函数 w_coeff bilinear_transform(z_coeff) print(f变换后w域系数{w_coeff})处理建议使用符号计算工具如SymPy保持精度对高阶系统可分阶段变换检查变换前后多项式阶数是否一致4.2 临界稳定情况当系统处于临界稳定状态时有根在单位圆上劳斯表会出现特殊情形全零行表明存在大小相等、位置相反的根首列零元素可用无穷小量ε法处理示例处理流程出现全零行时用上一行构建辅助方程对辅助方程求导得到系数继续完成劳斯表解辅助方程可得单位圆上的根4.3 变换验证技巧为确保双线性变换正确执行推荐以下验证步骤极点位置验证% MATLAB示例直接求根验证 z_roots roots([1 -1.2 0.45 -0.05]); abs_z abs(z_roots) % 所有模应小于1边界案例测试取z1对应w→∞验证方程最高次项取z-1对应w0验证常数项系数和检查变换前后多项式在z1和w0的值应保持比例关系5. 扩展应用采样周期对稳定性的影响双线性变换不仅用于稳定性分析还可研究采样周期T的影响。考虑系统G(s) 1/(s(s1))采用零阶保持器离散化后不同T值下的z域传递函数T0.1s: G(z) (0.0048z 0.0047)/(z^2 - 1.905z 0.9048) T1s: G(z) (0.3679z 0.2642)/(z^2 - 1.3679z 0.3679)构建特征方程并应用双线性变换T0.1s时劳斯判据显示稳定T1s时可能进入临界状态工程启示过大的采样周期可能导致原本稳定的连续系统离散化后不稳定双线性变换结合劳斯判据可快速评估最大允许采样周期