4、采样系统稳定性分析：S平面到Z平面的映射，稳定性的充要条件
好，咱们今天聊点硬核的——采样系统的稳定性。
做控制这么多年，我见过太多人把连续系统的稳定性分析背得滚瓜烂熟，一碰到离散系统就懵了。说白了，就是没搞明白S平面和Z平面之间那点「映射关系」。今天我就把这块掰开了揉碎了讲清楚。
4.1 从S平面到Z平面：映射的本质
先问一个问题：为什么连续系统用S平面，离散系统就得用Z平面？
答案其实很简单——采样让系统变了「脾气」。
连续系统里，我们看极点落在S平面的左半平面还是右半平面。左半平面稳定，右半平面不稳定，这个大家应该都熟。但采样之后呢？系统变成了离散时间系统，我们得用Z变换来分析。
Z变换和拉普拉斯变换的关系，核心就是这一条：
z = e^(sT)其中T是采样周期。这个公式看着简单，但内涵很深。我当年刚学的时候，觉得这不就是个指数映射嘛，有啥了不起？后来做项目才发现，这个映射把S平面「卷」成了Z平面。
具体怎么卷的？
- S平面的虚轴（s = jω）映射到Z平面的单位圆上
- S平面的左半平面（σ < 0）映射到Z平面的单位圆内部
- S平面的右半平面（σ > 0）映射到Z平面的单位圆外部
嗯，这里要注意：这个映射不是一对一的。S平面上每隔2π/T的频率，都会映射到Z平面的同一个点。这就是所谓的「频率混叠」——我在项目中吃过这个亏，后面会细说。
4.2 稳定性的充要条件：一句话的事
采样系统稳定的充要条件，其实就一句话：
所有闭环极点必须位于Z平面的单位圆内部。
换句话说，如果系统的特征方程是：
1 + G(z)H(z) = 0那么所有特征根的模必须小于1，即 |z_i| < 1。
为什么？你想想看，连续系统里极点必须在左半平面，对应的是指数衰减。而离散系统里，极点模小于1对应的是几何级数衰减。说白了，就是每次采样后，系统的响应越来越小，最终趋于零。
如果某个极点的模等于1，系统就处于临界稳定——振荡不衰减也不发散。如果模大于1，那就完了，系统发散，输出越来越大，直到饱和或者炸掉。
注意： 临界稳定在实际工程中是不允许的。我曾经调试一个伺服系统，极点刚好落在单位圆上，仿真看着没问题，一上实物就抖得不行。后来才明白，模型误差和量化误差会把临界极点推到圆外。
4.3 映射关系可视化：一张图说清楚
为了让大家更直观地理解S平面到Z平面的映射，我画了一张图：
这张图左边是S平面，右边是Z平面。箭头表示映射方向。你看，S平面左半平面（绿色区域）映射到Z平面单位圆内部（也是绿色），右半平面（红色区域）映射到单位圆外部（红色）。
这个映射关系，说白了就是稳定性的「翻译器」。你只要记住：S平面的左半平面 = Z平面的单位圆内部，就够了。
4.4 实际分析中的坑：我踩过的那些雷
理论说完了，聊聊实战。我在做数字电源控制的时候，遇到过这么几个问题：
坑一：采样周期选不好，系统稳不住
有一次我设计一个电流环，连续域下相位裕度有60°，觉得稳得很。结果一离散化，系统振荡了。查了半天才发现，采样频率选得太低，导致高频极点映射到了单位圆附近。
经验公式：采样频率至少是系统带宽的10倍以上。低于这个值，你就要小心了。
坑二：零阶保持器带来的相位滞后
很多人做离散化的时候，直接拿连续系统的传递函数套Z变换，忽略了零阶保持器的影响。零阶保持器会引入额外的相位滞后，相当于在系统中加了一个延迟环节。
我建议：做离散化的时候，一定要把零阶保持器考虑进去。用下面的公式：
G(z) = (1 - z^(-1)) · Z[G(s)/s]这个公式看着麻烦，但它是准确的。偷懒用近似方法，早晚要还的。
坑三：双线性变换的频率畸变
双线性变换（Tustin变换）是常用的离散化方法，但它有个问题——频率会畸变。高频段的频率会被压缩。如果你用双线性变换设计数字滤波器，一定要做频率预畸变。
小技巧： 做预畸变时，把目标频率ω代入公式：ω_pre = (2/T) · tan(ωT/2)。这样双线性变换后的实际频率才能对上。
4.5 稳定性判据：不止一种方法
除了直接看极点位置，还有几种常用的稳定性判据：
| 判据名称 | 适用场景 | 核心思想 |
|---|---|---|
| 朱利判据 | 特征方程已知 | 构造朱利表，检查第一列符号 |
| 双线性变换+劳斯判据 | 系统阶数较高 | 将Z平面映射回W平面，用劳斯判据 |
| 奈奎斯特判据（离散版） | 频域分析 | 看开环频率特性包围(-1, j0)的情况 |
| 根轨迹法（离散版） | 参数变化分析 | 看根轨迹是否穿越单位圆 |
我个人最常用的是双线性变换+劳斯判据。为什么？因为劳斯判据我熟，不用学新东西。把Z平面的特征方程做个变换：
z = (1 + w) / (1 - w)然后对w的多项式用劳斯判据，判断是否所有根都在左半平面。这个方法简单粗暴，适合工程快速判断。
4.6 一个完整的分析例子
咱们看个具体的例子。假设一个采样系统的闭环传递函数为：
G(z) = K / (z^2 - 1.5z + 0.7)问：K取何值时系统稳定？
第一步，写出特征方程：
z^2 - 1.5z + 0.7 + K = 0第二步，用朱利判据。对于二阶系统，稳定条件很简单：
- |a_0| < a_2，即 |0.7 + K| < 1
- 特征多项式在z=1处为正：1 - 1.5 + 0.7 + K > 0 → K > -0.2
- 特征多项式在z=-1处为正：1 + 1.5 + 0.7 + K > 0 → K > -3.2
综合起来：-0.2 < K < 0.3
你看，K的范围其实很窄。如果K取0.4，系统就不稳定了。我在项目中就遇到过类似的情况——增益稍微调大一点，系统就开始抖。后来加了限幅和补偿才解决。
核心总结：
- S平面左半平面 ↔ Z平面单位圆内部
- 所有极点模小于1，系统稳定
- 采样频率要足够高，至少10倍带宽
- 离散化时别忘了零阶保持器
- 双线性变换要做频率预畸变
好了，采样系统稳定性分析就讲到这里。这些内容看着多，其实核心就一个映射关系。把这个搞明白了，后面不管是设计控制器还是分析性能，都会顺手很多。