如何求解带延迟系统的最大增益极限?
这个问题问到点子上了!带纯延迟项 ( e^{-Ls} ) 的反馈系统,没法直接用有理多项式求根的方法找最大增益极限,因为延迟项是非有理函数。咱们得切换到频率域的临界稳定分析思路,具体操作步骤如下:
1. 先推导临界稳定的特征方程
对于单位反馈系统,闭环特征方程是 ( 1 + K G(s) = 0 )。把 ( G(s) = \frac{num(s)}{den(s)}e^{-Ls} ) 代进去,整理一下得到:
$$den(s) + K \cdot num(s) \cdot e^{-Ls} = 0$$
系统临界稳定时,特征根刚好落在虚轴上(也就是 ( s = jω ),( ω > 0 )),把这个代入上面的方程:
$$den(jω) + K \cdot num(jω) \cdot e^{-jLω} = 0$$
再整理成更方便分析的形式:
$$K \cdot \frac{num(jω)}{den(jω)} = -e^{jLω}$$
2. 拆分幅值和相位两个核心条件
这个复方程可以拆成幅值相等和相位匹配两个独立条件,咱们分别处理:
幅值条件
两边取幅值(注意 ( |-e^{jLω}| = 1 ),因为指数项的幅值永远是1):
$$K \cdot \left| \frac{num(jω)}{den(jω)} \right| = 1$$
直接就能解出增益和频率的关系:
$$K = \left| \frac{den(jω)}{num(jω)} \right|$$
相位条件
两边取相位,这里要注意 ( -1 ) 的相位是 ( π + 2πn )(( n ) 是任意整数,因为相位是周期性的):
$$\angle num(jω) - \angle den(jω) = π + 2πn + Lω$$
3. 找到临界频率,计算最大增益极限
接下来的步骤就是找满足相位条件的所有临界频率 ( ω_c ),然后代入幅值条件算出对应的增益 ( K_c ):
- 先遍历所有正的 ( ω ),找到所有满足相位条件的 ( ω_c )(可能有多个,尤其是延迟 ( L ) 较大时)
- 对每个 ( ω_c ),用幅值条件算出对应的 ( K_c )
- 这些 ( K_c ) 里的最小值就是系统稳定的最大增益极限——超过这个值,系统就会不稳定
4. 实用工具和技巧
如果手动计算太麻烦,这些工具能帮到你:
- 奈奎斯特图法:先画有理部分 ( \frac{num(s)}{den(s)} ) 的奈奎斯特曲线,延迟项会让这条曲线随着频率增加不断顺时针旋转(每次旋转角度是 ( -Lω ))。当曲线刚好穿过负实轴上的 ( (-1, 0) ) 点时,对应的 ( K ) 就是临界增益。
- 数值求解工具:比如用MATLAB的
margin函数(它支持带延迟的系统),或者自己写个小脚本遍历 ( ω ),先找到满足相位条件的 ( ω_c ),再代入算 ( K )。
举个简单例子帮你理解:假设 ( G(s) = \frac{1}{s+1}e^{-Ls} ),特征方程是 ( s+1 + K e^{-Ls} = 0 )。代入 ( s=jω ) 后,幅值条件是 ( \sqrt{ω²+1} = K ),相位条件是 ( arctan(ω) = π + Lω )。先解相位条件找到 ( ω_c ),再代入幅值条件就得到临界 ( K ) 了。
内容的提问来源于stack exchange,提问作者euraad

