题目
3.15已知单位反馈系统的开环传递函数分别为:
(1)G(s)=100(0.1s+1)(s+5)
(2)G(s)=50s(0.1s+1)(s+5)
(3)G(s)=10(2s+1)s2(s2+6s+100)
试求输入分别为r(1)=2t和r(t)=2+2t+t2时,系统的稳态误差。
知识点
知识点窍:稳态误差、静态误差系统法、终值定理
逻辑推理:根据线性叠加原理,用静态误差系数法或终值定理来求系统的稳态误差。
解题过程
开环传递函数时间常数形式
G(s)H(s)=K(1+τ1s)(1+τ2s)⋯(1+τms)sv(1+T1s)(1+T2s)⋯(1+Tn−vs)=Km∏j=1(τjs+1)svn−v∏i=1(Tis+1)
当v=0时,称该系统为0型系统
当v=1时,称该系统为1型系统
当v=2时,称该系统为2型系统
(1)
G(s)=100(0.1s+1)(s+5)=20(0.1s+1)(0.2s+1)
该形式已经是时间常数形式了
由上式可知,该系统是0型系统,且K=20。
查表可得,0型系统在1(t),t,12t2信号作用下的稳态误差分别为:
11+K,∞,∞
根据线性叠加原理,该系统带输入为r(t)=2t时的稳态误差为
ess1=2×∞=∞
该系统在输入
r(t)=2+2t+t2=2+2t+2×12t2
时的稳态误差为:
ess2=2×11+K+2×∞+2×∞=21+K+∞+∞=∞
(2)
首先要化为时间常数形式:
G(s)H(s)=K(1+τ1s)(1+τ2s)⋯(1+τms)sv(1+T1s)(1+T2s)⋯(1+Tn−vs)=Km∏j=1(τjs+1)svn−v∏i=1(Tis+1)
注意:分子分母的常数项都是1
G(s)=50s(0.1s+1)(s+5)=10s(0.1s+1)(0.2s+1)
sv=s1,v=1,所以该系统为1型系统
开环增益K=10,
查表可得,1型系统在输入信号1(t),t,12t2的作用下的稳态误差分别为:
0,1K,∞
根据线性叠加原理,该系统带输入为r(t)=2t时的稳态误差为:
ess1=2×1K=2K=210=0.2
该系统在输入r(t)=2+2t+2×12t2是的稳态误差为:
ess2=2×0+2K+2×∞=∞
(3)
G(s)=10(2s+1)s2(s2+6s+100)
只有稳定的系统才有稳态误差,所以首先需要判定此系统的稳定性对于单位负反馈系统有H(s)=1,所以系统的
所以闭环传递函数:
Φ(s)=G(s)1+G(s)H(s)=G(s)1+G(s)×1=10(2s+1)s2(s2+6s+100)1+10(2s+1)s2(s2+6s+100)=10(2s+1)s2(s2+6s+100)+10(2s+1)
令闭环传递函数分母为0即可得到闭环特征方程(其实开环传递函数分母加上分子):
闭环特征方程为:
D(s)=10(2s+1)s2(s2+6s+100)+10(2s+1)=s4+6s3+100s2+20s+10=0
用劳思稳定判据来确定此系统的稳定性,列劳思表如下
s4 | 1 | 100 | 10 |
s3 | 6 | 20 | |
s2 | 6×6×100−1×206=6×5806=580 | 6×6×10−1×06=60 | |
s1 | 580×580×20−6×60580=11240 | ||
s0 | 11240×60−580×011240=60 |
劳斯表第一列全部为正数,所以系统稳定
用终值定理来求系统的稳态误差,公式如下:
ess3=lim
代入闭环传递函数
G(s)=\dfrac{10(2s+1)}{s^2(s^2+6s+100)},反馈通路传递函数H(s)=1
得到:
\begin{aligned} e_{ss3}=&\lim\limits_{s\to 0}s\dfrac{1}{1+\dfrac{10(2s+1)}{s^2(s^2+6s+100)}}R(s) \\ =&\lim\limits_{s\to 0}s\dfrac{s^2(s^2+6s+100)}{s^2(s^2+6s+100)+10(2s+1)}R(s) \\ \end{aligned}
当r(t)=2t时,R(s)=\dfrac{2}{s^2},代入稳态误差表达式,求极限就可以得到稳态误差:
\begin{aligned} e_{ss3}=&\lim\limits_{s\to 0}s\times\dfrac{1}{1+\dfrac{10(2s+1)}{s^2(s^2+6s+100)}}\times R(s) \\ =&\lim\limits_{s\to 0}s\times\dfrac{s^2(s^2+6s+100)}{s^2(s^2+6s+100)+10(2s+1)}\times\dfrac{2}{s^2} \\ =&\lim\limits_{s\to 0}s\times\dfrac{2(s^2+6s+100)}{s^2(s^2+6s+100)+10(2s+1)} \\ =&0 \\ \end{aligned}
当r(t)=2+2t+2\times\dfrac{1}{2}t^2,
\begin{aligned} R(s)=&\dfrac{2}{s}+\dfrac{2}{s^2}+\dfrac{2}{s^3} \\ =&\dfrac{2s^2+2s+2}{s^3} \\ \end{aligned}
代入稳态误差表达式得到:
\begin{aligned} &e_{ss3} \\ &=\lim\limits_{s\to 0}s\dfrac{s^2(s^2+6s+100)}{s^2(s^2+6s+100)+10(2s+1)}R(s) \\ &=\lim\limits_{s\to 0}s\dfrac{s^2(s^2+6s+100)}{s^2(s^2+6s+100)+10(2s+1)}\dfrac{2s^2+2s+2}{s^3} \\ &=\lim\limits_{s\to 0}\dfrac{s^2(s^2+6s+100)}{s^2(s^2+6s+100)+20s+10)}\dfrac{2s^2+2s+2}{s^2} \\ &=\lim\limits_{s\to 0}\dfrac{(s^2+6s+100)(2s^2+2s+2)}{s^2(s^2+6s+100)+20s+10)} \\ &=\dfrac{100\times2}{10} \\ &=20 \\ \end{aligned}
原文链接: 自动控制原理 课后习题 3-15