计算机控制系统」4. 从Z平面到性能指标:三大分析的工程实践指南

1. Z平面映射:离散系统的坐标系转换

第一次接触Z平面时,我盯着那些在单位圆内外跳动的极点发懵——它们就像电子游戏里的角色,位置不同居然能决定整个系统的"性格"。后来才明白,Z平面其实就是离散控制系统专属的GPS坐标系,而理解S平面到Z平面的映射关系,就像掌握了一套加密地图的翻译规则。

最关键的映射发生在S平面虚轴(s=jω)到Z平面单位圆(z=e^(jωT))的转换。这里有个工程师容易踩的坑:采样频率ωs=2π/T会直接影响映射精度。比如当ω=ωs/2时,对应的Z平面角度θ=π,这意味着高频信号会被"折叠"成低频表现,就像高速旋转的轮子看起来在倒转一样。我在调试电机控制系统时就遇到过,采样周期T选得太大,导致实际50Hz的振动在系统中显示成10Hz的假信号。

主带映射规律可以总结为三句口诀:

  • 左半平面居民(稳定极点)→单位圆内购房
  • 右半平面流浪汉(不稳定极点)→单位圆外游荡
  • 虚轴上的临界分子单位圆上跳广场舞

等频率线的映射特别有意思。在S平面里,ω=2rad/s和ω=2+ωs rad/s的两个信号,到了Z平面会完全重合。这就解释了为什么离散系统总有混叠现象,就像不同楼层的住户挤进同一部电梯,系统根本分不清谁是谁。

2. 稳定性判据:离散世界的交通规则

记得刚入行时,我拿着连续系统的劳斯判据直接套用在离散系统上,结果被导师骂得狗血淋头。后来才搞懂,离散系统需要特别的"安检设备"——要么用W变换把单位圆"压扁"成左半平面,要么用Jury判据直接检测。

W变换就像给Z平面戴了副哈哈镜,把单位圆内部映射到W平面的左侧。具体操作时有个易错点:变换公式z=(1+w)/(1-w)会把z=1这个点映射到w=∞,所以计算时要特别注意极限情况。有次我分析一个二阶系统,没注意这个细节,导致稳定性判断完全错误,现场设备直接啸叫抗议。

Jury判据更适合手算,它的判定表格就像填数独:

  1. 第一行放特征方程系数
  2. 奇数行正序填写
  3. 偶数行倒序填写
  4. 计算各阶子行列式

但要注意,当遇到某行全零时,需要像处理多项式重根那样构造辅助方程。我在教新人时总强调:判据只是工具,物理直觉才是核心。比如看到极点聚集在z=1附近,就该警惕系统可能有临界振荡倾向。

3. 稳态误差:系统精度的温度计

曾有个项目让我印象深刻:设计的控制器在仿真时误差为零,实际运行却总有0.5mm的恒定偏差。这就是典型没吃透静态误差系数的教训——离散系统的稳态误差计算比连续系统多几个"陷阱"。

首先得记住铁律:先判稳再算误差!就像不能给危房做精装修,不稳定系统的误差计算毫无意义。三种典型输入对应的误差系数:

  • 位置误差系数Kp:对付阶跃输入的"镇定剂"
  • 速度误差系数Kv:追踪斜坡信号的"滑翔伞"
  • 加速度误差系数Ka:应对抛物线输入的"减震器"

计算时最容易栽在z→1的极限运算上。比如处理型别Ⅰ系统时,(1-z^-1)这个因子就像放大镜,能把隐藏在传递函数里的积分特性找出来。有个实用技巧:把z表示为1+Δz,用泰勒展开处理,比直接洛必达更稳妥。

采样周期T在这里扮演双面角色:

  • 减小T能提升Kv和Ka,改善动态跟踪
  • 但过小的T会引发数值计算问题 经验值是选择T≤1/(10ωc),其中ωc是截止频率。

4. 动态特性:极点的舞蹈编排

闭环极点在Z平面的分布就像编舞队形,决定了系统响应的"舞步风格"。我习惯把单位圆分成几个关键区域:

实轴极点区(0<z<1):

  • 单调衰减的"老干部舞步"
  • 越靠近原点,衰减速度越快
  • 负数极点会产生震荡,像机械舞的顿挫感

共轭复数极点区

  • 震荡衰减的"华尔兹"
  • 幅角θ决定频率,模ρ决定衰减速度
  • 最佳阻尼比区域在β=45°射线附近

特别要注意多采样周期系统的响应特点。有次调试一个三阶系统,明明所有极点都在单位圆内,响应却出现周期性波动。后来发现是采样周期与系统自然频率产生谐振,就像跳舞时踩到自己脚。解决方法要么调整T,要么在z=-1附近加零点抵消。

最小拍系统的设计就像编排出道即巅峰的舞蹈:

  1. 所有极点严格控制在原点
  2. 响应在有限拍数内结束
  3. 但对模型误差零容忍 实际工程中我更推荐准最小拍设计,留出10%的调节余量。

5. 频率特性:离散世界的彩色眼镜

虽然离散频率特性不如连续的Bode图直观,但掌握后就像获得夜视仪。关键点在于:

  • 有效频率范围只在0到ωs/2(奈奎斯特频率)
  • 频率响应具有周期性,像万花筒图案
  • 对数坐标依然适用,但刻度要换算

有个实用技巧:分析抗混叠时,可以观察ωs/2处的增益。某次验收滤波器时,就是发现这个频点衰减不足,才避免了一场事故。对于高阶系统,建议先用MATLAB的freqz函数生成幅相曲线,再结合理论分析。

调试数字控制器时,我常关注两个特征频率点:

  1. 相位穿越频率ωc:增益余量的温度计
  2. 谐振峰值频率ωr:系统灵敏度的警报器 经验表明,保持ωc在ωs/6到ωs/4之间,通常能兼顾响应速度与稳定性。