此系列属于胡寿松《自动控制原理题海与考研指导》(第三版)习题精选,仅包含部分经典习题,需要完整版习题答案请自行查找,本系列属于知识点巩固部分,搭配如下几个系列进行学习,可用于期末考试和考研复习。
自动控制原理(第七版)知识提炼
自动控制原理(第七版)课后习题精选
自动控制原理(第七版)附录MATLAB基础
第三章:时域分析法
Example 3.31
设单位反馈系统的开环传递函数为: G ( s ) = K ( 0.2 s + 1 ) ( 0.4 s + 1 ) G(s)=\displaystyle\frac{K}{(0.2s+1)(0.4s+1)} G(s)=(0.2s+1)(0.4s+1)K,当 K K K分别取 1 , 3 , 7 1,3,7 1,3,7时,列表比较系统的自然频率、阻尼比、上升时间、峰值时间、调节时间、超调量和在单位阶跃输入下的稳态误差,并指出 K K K值变化对系统产生的影响.
解:
依题意可得,系统的闭环传递函数为:
Φ ( s ) = G ( s ) 1 + G ( s ) = K 0.08 s 2 + 0.6 s + 1 + K = 12.5 K s 2 + 7.5 s + 12.5 ( 1 + K ) \Phi(s)=\frac{G(s)}{1+G(s)}=\frac{K}{0.08s^2+0.6s+1+K}=\frac{12.5K}{s^2+7.5s+12.5(1+K)} Φ(s)=1+G(s)G(s)=0.08s2+0.6s+1+KK=s2+7.5s+12.5(1+K)12.5K
可得系统的自然频率和阻尼比:
ω n = 12.5 ( 1 + K ) , ζ = 3.75 / ω n \omega_n=\sqrt{12.5(1+K)},\zeta=3.75/\omega_n ωn=12.5(1+K),ζ=3.75/ωn
系统动态性能指标计算公式:
上升时间: t r = π − arccos ζ ω n 1 − ζ 2 ;峰值时间: t p = π ω n 1 − ζ 2 调节时间: t s = 4.4 ζ ω n ( Δ = 2 % ) ;超调量: σ % = e − π ζ / 1 − ζ 2 × 100 % \begin{aligned} &上升时间:t_r=\frac{\pi-\arccos\zeta}{\omega_n\sqrt{1-\zeta^2}};峰值时间:t_p=\frac{\pi}{\omega_n\sqrt{1-\zeta^2}}\\\\ &调节时间:t_s=\frac{4.4}{\zeta\omega_n}(\Delta=2\%);超调量:\sigma\%={\rm e}^{-\pi\zeta/\sqrt{1-\zeta^2}}\times100\% \end{aligned} 上升时间:tr=ωn1−ζ2π−arccosζ;峰值时间:tp=ωn1−ζ2π调节时间:ts=ζωn4.4(Δ=2%);超调量:σ%=e−πζ/1−ζ2×100%
单位阶跃输入下的稳态误差可由终值定理可得:
e s s ( ∞ ) = 1 1 + K e_{ss}(\infty)=\frac{1}{1+K} ess(∞)=1+K1
【性能比较】
K K K | 1 1 1 | 3 3 3 | 7 7 7 |
---|---|---|---|
自然频率 ω n \omega_n ωn | 5 5 5 | 7.07 7.07 7.07 | 10 10 10 |
阻尼比 ζ \zeta ζ | 0.75 0.75 0.75 | 0.53 0.53 0.53 | 0.375 0.375 0.375 |
上升时间 t r / s t_r/s tr/s | 0.73 0.73 0.73 | 0.36 0.36 0.36 | 0.21 0.21 0.21 |
峰值时间 t p / s t_p/s tp/s | 0.95 0.95 0.95 | 0.52 0.52 0.52 | 0.34 0.34 0.34 |
调节时间 t s ( Δ = 2 % ) / s t_s(\Delta=2\%)/s ts(Δ=2%)/s | 1.17 1.17 1.17 | 1.17 1.17 1.17 | 1.17 1.17 1.17 |
超调量 σ % \sigma\% σ% | 2.8 2.8 2.8 | 14.0 14.0 14.0 | 28.1 28.1 28.1 |
稳态误差 e s s ( ∞ ) e_{ss}(\infty) ess(∞) | 0.5 0.5 0.5 | 0.25 0.25 0.25 | 0.125 0.125 0.125 |
当 K K K增加时,系统的自然频率变大,阻尼比变小,系统的快速性增强,但超调量增大,稳态误差越来越小.
【单位阶跃响应】
Example 3.32
系统结构图如下图所示,已知 K 1 > 0 , K 2 > 0 K_1>0,K_2>0 K1>0,K2>0,确定闭环系统稳定时 K 1 K_1 K1和 K 2 K_2 K2的可调范围,并在 K 1 − K 2 K_1-K_2 K1−K2平面上画出参数稳定域的概略图形.
解:
由图可得,系统开环传递函数为:
G ( s ) = K 1 s ( 0.05 s + 1 ) ( 0.1 s + 1 ) + K 1 K 2 s G(s)=\frac{K_1}{s(0.05s+1)(0.1s+1)+K_1K_2s} G(s)=s(0.05s+1)(0.1s+1)+K1K2sK1
系统的闭环特征方程为:
D ( s ) = s ( 0.05 s + 1 ) ( 0.1 s + 1 ) + K 1 K 2 s + K 1 = 0.005 s 3 + 0.15 s 2 + ( 1 + K 1 K 2 ) s + K 1 = 0 D(s)=s(0.05s+1)(0.1s+1)+K_1K_2s+K_1=0.005s^3+0.15s^2+(1+K_1K_2)s+K_1=0 D(s)=s(0.05s+1)(0.1s+1)+K1K2s+K1=0.005s3+0.15s2+(1+K1K2)s+K1=0
劳斯表如下:
s 3 s^3 s3 | 0.005 0.005 0.005 | 1 + K 1 K 2 1+K_1K_2 1+K1K2 |
---|---|---|
s 2 s^2 s2 | 0.15 0.15 0.15 | K 1 K_1 K1 |
s 1 s^1 s1 | 0.15 ( 1 + K 1 K 2 ) − 0.005 K 1 0.15 \displaystyle\frac{0.15(1+K_1K_2)-0.005K_1}{0.15} 0.150.15(1+K1K2)−0.005K1 | |
s 0 s^0 s0 | K 1 K_1 K1 |
要使闭环系统稳定,则有:
{ 0.15 ( 1 + K 1 K 2 ) − 0.005 K 1 > 0 K 1 > 0 K 2 > 0 ⇒ { K 1 > 30 K 2 > 1 / 30 − 1 / K 1 \begin{cases} &0.15(1+K_1K_2)-0.005K_1>0\\ &K_1>0\\ &K_2>0 \end{cases}\Rightarrow \begin{cases} &K_1>30\\ &K_2>1/30-1/K_1 \end{cases} ⎩
⎨
⎧0.15(1+K1K2)−0.005K1>0K1>0K2>0⇒{
K1>30K2>1/30−1/K1
【 K 1 − K 2 K_1-K_2 K1−K2平面上参数稳定域概略图】
【不同 K 1 、 K 2 K_1、K_2 K1、K2的单位阶跃响应】
Example 3.33
已知系统微分方程为: m y ¨ ( t ) + c y ˙ ( t ) + k y ( t ) = k x ( t ) m\ddot{y}(t)+c\dot{y}(t)+ky(t)=kx(t) my¨(t)+cy˙(t)+ky(t)=kx(t),其初始条件全部为零.
求:
- x ( t ) = 1 ( t ) x(t)=1(t) x(t)=1(t)时的输出响应 y ( t ) y(t) y(t);
- 输出 y ( t ) y(t) y(t)无振荡的条件;
解:
依题意可得:
Y ( s ) X ( s ) = k m s 2 + c s + k = k / m s 2 + ( c / m ) s + k / m \frac{Y(s)}{X(s)}=\frac{k}{ms^2+cs+k}=\frac{k/m}{s^2+(c/m)s+k/m} X(s)Y(s)=ms2+cs+kk=s2+(c/m)s+k/mk/m
可得系统的自然频率和阻尼比:
ω n = k m , ζ = c 2 m k \omega_n=\sqrt{\displaystyle\frac{k}{m}},\zeta=\frac{c}{2\sqrt{mk}} ωn=mk,ζ=2mkc
-
x ( t ) = 1 ( t ) x(t)=1(t) x(t)=1(t)时的输出响应 y ( t ) y(t) y(t);
-
当 ζ = c 2 m k ∈ [ 0 , 1 ) \zeta=\displaystyle\frac{c}{2\sqrt{mk}}\in[0,1) ζ=2mkc∈[0,1)时,即欠阻尼条件下,有:
y ( t ) = 1 − 1 1 − ζ 2 e − ζ ω n t sin ( ω n 1 − ζ 2 t + arccos ζ ) y(t)=1-\frac{1}{\sqrt{1-\zeta^2}}{\rm e}^{-\zeta\omega_nt}\sin(\omega_n\sqrt{1-\zeta^2}t+\arccos\zeta) y(t)=1−1−ζ21e−ζωntsin(ωn1−ζ2t+arccosζ) -
当 ζ = c 2 m k = 1 \zeta=\displaystyle\frac{c}{2\sqrt{mk}}=1 ζ=2mkc=1时,即临界阻尼条件下,有:
y ( t ) = 1 − e − ω n t sin ( 1 + ω n t ) y(t)=1-{\rm e}^{-\omega_nt}\sin(1+\omega_nt) y(t)=1−e−ωntsin(1+ωnt) -
当 ζ = c 2 m k > 1 \zeta=\displaystyle\frac{c}{2\sqrt{mk}}>1 ζ=2mkc>1时,即过阻尼条件下,有:
y ( t ) = 1 − T 1 T 1 − T 2 e − t / T 1 + T 2 T 1 − T 2 e − t / T 2 y(t)=1-\frac{T_1}{T_1-T_2}{\rm e}^{-t/T_1}+\frac{T_2}{T_1-T_2}{\rm e}^{-t/T_2} y(t)=1−T1−T2T1e−t/T1+T1−T2T2e−t/T2
其中:
T 1 = 1 ω ( ζ − ζ 2 − 1 ) , T 2 = 1 ω n ( ζ + ζ 2 − 1 ) T_1=\frac{1}{\omega(\zeta-\sqrt{\zeta^2-1})},T_2=\frac{1}{\omega_n(\zeta+\sqrt{\zeta^2-1})} T1=ω(ζ−ζ2−1)1,T2=ωn(ζ+ζ2−1)1
-
-
输出 y ( t ) y(t) y(t)无振荡的条件;
要使输出 y ( t ) y(t) y(t)无振荡,则要求 ζ ≥ 1 \zeta≥1 ζ≥1,即当 c 2 m k ≥ 1 \displaystyle\frac{c}{2\sqrt{mk}}≥1 2mkc≥1时,输出 y ( t ) y(t) y(t)无振荡.
Example 3.34
设符合控制系统如下图所示:
要求:
- 计算扰动 n ( t ) = t n(t)=t n(t)=t引起的稳态误差;
- 设计 K c K_c Kc,使系统在 r ( t ) = t r(t)=t r(t)=t作用下无稳态误差.
解:
-
计算扰动 n ( t ) = t n(t)=t n(t)=t引起的稳态误差;
仅考虑扰动 n ( t ) = t n(t)=t n(t)=t,即 N ( s ) = 1 s 2 N(s)=\displaystyle\frac{1}{s^2} N(s)=s21时,由图可得:
C ( s ) N ( s ) = s ( s + K 2 K 3 ) ( T s + 1 ) s ( s + K 2 K 3 ) ( T s + 1 ) + K 1 K 2 K 4 \frac{C(s)}{N(s)}=\frac{s(s+K_2K_3)(Ts+1)}{s(s+K_2K_3)(Ts+1)+K_1K_2K_4} N(s)C(s)=s(s+K2K3)(Ts+1)+K1K2K4s(s+K2K3)(Ts+1)
此时系统由扰动引起的误差函数为:
E n ( s ) = − C ( s ) = − s ( s + K 2 K 3 ) ( T s + 1 ) s ( s + K 2 K 3 ) ( T s + 1 ) + K 1 K 2 K 4 N ( s ) E_n(s)=-C(s)=-\frac{s(s+K_2K_3)(Ts+1)}{s(s+K_2K_3)(Ts+1)+K_1K_2K_4}N(s) En(s)=−C(s)=−s(s+K2K3)(Ts+1)+K1K2K4s(s+K2K3)(Ts+1)N(s)
利用终值定理求系统稳态误差:
e s s n ( ∞ ) = lim s → 0 s E n ( s ) = − lim s → 0 s ⋅ s ( s + K 2 K 3 ) ( T s + 1 ) s ( s + K 2 K 3 ) ( T s + 1 ) + K 1 K 2 K 4 ⋅ 1 s 2 e_{ssn}(\infty)=\lim_{s\to0}sE_n(s)=-\lim_{s\to0}s·\frac{s(s+K_2K_3)(Ts+1)}{s(s+K_2K_3)(Ts+1)+K_1K_2K_4}·\frac{1}{s^2} essn(∞)=s→0limsEn(s)=−s→0lims⋅s(s+K2K3)(Ts+1)+K1K2K4s(s+K2K3)(Ts+1)⋅s21
则由扰动 n ( t ) = t n(t)=t n(t)=t引起的稳态误差为:
e s s n ( ∞ ) = − K 3 K 1 K 4 e_{ssn}(\infty)=-\frac{K_3}{K_1K_4} essn(∞)=−K1K4K3 -
设计 K c K_c Kc,使系统在 r ( t ) = t r(t)=t r(t)=t作用下无稳态误差.
仅考虑输入作用,即 R ( s ) = 1 s 2 R(s)=\displaystyle\frac{1}{s^2} R(s)=s21时,由图可得:
Φ ( s ) = C ( s ) R ( s ) = K c K 4 s + K 1 K 2 K 4 s 2 ( T s + 1 ) + K 2 K 3 s ( T s + 1 ) + K 1 K 2 K 4 \Phi(s)=\frac{C(s)}{R(s)}=\frac{K_cK_4s+K_1K_2K_4}{s^2(Ts+1)+K_2K_3s(Ts+1)+K_1K_2K_4} Φ(s)=R(s)C(s)=s2(Ts+1)+K2K3s(Ts+1)+K1K2K4KcK4s+K1K2K4
系统的误差函数为:
E r ( s ) = R ( s ) − C ( s ) = [ 1 − Φ ( s ) ] ⋅ R ( s ) = s 2 ( T s + 1 ) + K 2 K 3 s ( T s + 1 ) − K c K 4 s s 2 ( T s + 1 ) + K 2 K 3 s ( T s + 1 ) + K 1 K 2 K 4 ⋅ R ( s ) E_r(s)=R(s)-C(s)=[1-\Phi(s)]·R(s)=\frac{s^2(Ts+1)+K_2K_3s(Ts+1)-K_cK_4s}{s^2(Ts+1)+K_2K_3s(Ts+1)+K_1K_2K_4}·R(s) Er(s)=R(s)−C(s)=[1−Φ(s)]⋅R(s)=s2(Ts+1)+K2K3s(Ts+1)+K1K2K4s2(Ts+1)+K2K3s(Ts+1)−KcK4s⋅R(s)
由终值定理求系统的稳态误差,有:
e s s r ( ∞ ) = lim s → 0 s E r ( s ) = lim s → 0 s ⋅ s 2 ( T s + 1 ) + K 2 K 3 s ( T s + 1 ) − K c K 4 s s 2 ( T s + 1 ) + K 2 K 3 s ( T s + 1 ) + K 1 K 2 K 4 ⋅ 1 s 2 = lim s → 0 s ( T s + 1 ) + K 2 K 3 T s + ( K 2 K 3 − K c K 4 ) s 2 ( T s + 1 ) + K 2 K 3 s ( T s + 1 ) + K 1 K 2 K 4 \begin{aligned} e_{ssr}(\infty)&=\lim_{s\to0}sE_r(s)=\lim_{s\to0}s·\frac{s^2(Ts+1)+K_2K_3s(Ts+1)-K_cK_4s}{s^2(Ts+1)+K_2K_3s(Ts+1)+K_1K_2K_4}·\frac{1}{s^2}\\\\ &=\lim_{s\to0}\frac{s(Ts+1)+K_2K_3Ts+(K_2K_3-K_cK_4)}{s^2(Ts+1)+K_2K_3s(Ts+1)+K_1K_2K_4} \end{aligned} essr(∞)=s→0limsEr(s)=s→0lims⋅s2(Ts+1)+K2K3s(Ts+1)+K1K2K4s2(Ts+1)+K2K3s(Ts+1)−KcK4s⋅s21=s→0lims2(Ts+1)+K2K3s(Ts+1)+K1K2K4s(Ts+1)+K2K3Ts+(K2K3−KcK4)
要使系统在 r ( t ) = t r(t)=t r(t)=t作用下无稳态误差,有:
K 2 K 3 − K c K 4 = 0 ⇒ K c = K 2 K 3 K 4 K_2K_3-K_cK_4=0\Rightarrow{K_c}=\frac{K_2K_3}{K_4} K2K3−KcK4=0⇒Kc=K4K2K3
即当 K c = K 2 K 3 K 4 K_c=\displaystyle\frac{K_2K_3}{K_4} Kc=K4K2K3时,系统在 r ( t ) = t r(t)=t r(t)=t作用下无稳态误差.
Example 3.35
已知单位反馈的角度随动系统的开环传递函数: G ( s ) = K 1 s ( J s + F ) G(s)=\displaystyle\frac{K_1}{s(Js+F)} G(s)=s(Js+F)K1,其中: J = 200 , F = 1600 J=200,F=1600 J=200,F=1600.假设系统处于临界阻尼状态,输入 r ( t ) = a m sin ω t , a m = 34 ° , ω = 1 r(t)=a_m\sin\omega{t},a_m=34°,\omega=1 r(t)=amsinωt,am=34°,ω=1.求:
- 系统稳态输出 c s ( t ) c_s(t) cs(t)的摆动幅度 c m c_m cm;
- 系统稳态输出 c s ( t ) c_s(t) cs(t)对输入 r ( t ) r(t) r(t)的滞后时间 τ \tau τ;
解:
根据题意可得,系统开环传递函数为:
G ( s ) = K 1 s ( J s + F ) = K 1 / 200 s ( s + 8 ) G(s)=\frac{K_1}{s(Js+F)}=\frac{K_1/200}{s(s+8)} G(s)=s(Js+F)K1=s(s+8)K1/200
即:
ω n 2 = K 1 / 200 , 2 ζ ω n = 8 \omega_n^2=K_1/200,2\zeta\omega_n=8 ωn2=K1/200,2ζωn=8
由于系统处于临界阻尼状态,则有:
ζ = 1 , ω n = 4 , K 1 = 3200 \zeta=1,\omega_n=4,K_1=3200 ζ=1,ωn=4,K1=3200
则系统闭环传递函数为:
Φ ( s ) = G ( s ) 1 + G ( s ) = 16 s 2 + 8 s + 16 \Phi(s)=\frac{G(s)}{1+G(s)}=\frac{16}{s^2+8s+16} Φ(s)=1+G(s)G(s)=s2+8s+1616
-
系统稳态输出 c s ( t ) c_s(t) cs(t)的摆动幅度 c m c_m cm;
当输入为: r ( t ) = a m sin ω t r(t)=a_m\sin\omega{t} r(t)=amsinωt时,
Φ ( j ω ) = 16 ( 16 − ω 2 ) + j 8 ω \Phi({\rm j}\omega)=\frac{16}{(16-\omega^2)+{\rm j}8\omega} Φ(jω)=(16−ω2)+j8ω16
由于: a m = 34 ° , ω = 1 a_m=34°,\omega=1 am=34°,ω=1.有:
∣ Φ ( j 1 ) ∣ = 16 1 5 2 + 8 2 = 16 17 |\Phi({\rm j}1)|=\frac{16}{\sqrt{15^2+8^2}}=\frac{16}{17} ∣Φ(j1)∣=152+8216=1716c m = 34 ° × 16 17 = 32 ° c_m=34°\times\frac{16}{17}=32° cm=34°×1716=32°
故系统稳态输出 c s ( t ) c_s(t) cs(t)的摆动幅度 c m = 32 ° c_m=32° cm=32°.
-
系统稳态输出 c s ( t ) c_s(t) cs(t)对输入 r ( t ) r(t) r(t)的滞后时间 τ \tau τ;
由于
Φ ( j ω ) ∣ ω = 1 = 16 15 + j 8 \left.\Phi({\rm j}\omega)\right|_{\omega=1}=\frac{16}{15+{\rm j}8} Φ(jω)∣ω=1=15+j816
则
φ = arctan ( 8 15 ) = 28.1 ° = 0.49 r a d \varphi=\arctan\left(\frac{8}{15}\right)=28.1°=0.49{\rm rad} φ=arctan(158)=28.1°=0.49rad
此时系统稳态输出为:
c s ( t ) = c m sin ( t − φ ) = 32 ° sin ( t − 28.1 ° ) c_s(t)=c_m\sin(t-\varphi)=32°\sin(t-28.1°) cs(t)=cmsin(t−φ)=32°sin(t−28.1°)
故系统稳态输出 c s ( t ) c_s(t) cs(t)对输入 r ( t ) r(t) r(t)的滞后时间 τ = φ / ω = 0.49 s \tau=\varphi/\omega=0.49{\rm s} τ=φ/ω=0.49s. -
【时间响应曲线】
Example 3.36
设控制系统如下图所示,系统输入端除有用信号 r ( t ) r(t) r(t)以外,还夹杂有扰动 n ( t ) n(t) n(t).已知 r ( t ) = 10 t , n ( t ) = 0.1 sin 10 t r(t)=10t,n(t)=0.1\sin10t r(t)=10t,n(t)=0.1sin10t.
计算系统稳态误差的最大值,并概略画出初始状态为零时的输出响应 c ( t ) c(t) c(t)曲线.
解:
由图可得,系统开环传递函数为:
G ( s ) = 10 s ( 0.1 s + 1 ) ( 0.4 s + 1 ) G(s)=\frac{10}{s(0.1s+1)(0.4s+1)} G(s)=s(0.1s+1)(0.4s+1)10
仅考虑有用信号 r ( t ) = 10 t r(t)=10t r(t)=10t作用,即 R ( s ) = 10 s 2 R(s)=\displaystyle\frac{10}{s^2} R(s)=s210时,由系统的开环传递函数可知,该系统为Ⅰ型系统,且 K v = 10 K_v=10 Kv=10,则在 r ( t ) = 10 t r(t)=10t r(t)=10t作用下,系统的稳态误差为:
e s s r ( ∞ ) = R K v = 1 e_{ssr}(\infty)=\frac{R}{K_v}=1 essr(∞)=KvR=1
仅考虑扰动信号 n ( t ) = 0.1 sin 10 t n(t)=0.1\sin10t n(t)=0.1sin10t作用,系统的误差为:
E n ( s ) = − C ( s ) = − G ( s ) 1 + G ( s ) N ( s ) = − 10 s ( 0.1 s + 1 ) ( 0.4 s + 1 ) + 10 N ( s ) E_n(s)=-C(s)=-\frac{G(s)}{1+G(s)}N(s)=-\frac{10}{s(0.1s+1)(0.4s+1)+10}N(s) En(s)=−C(s)=−1+G(s)G(s)N(s)=−s(0.1s+1)(0.4s+1)+1010N(s)
则误差传递函数为:
Φ e n ( s ) = − 10 0.04 s 3 + 0.5 s 2 + s + 10 \Phi_{en}(s)=-\frac{10}{0.04s^3+0.5s^2+s+10} Φen(s)=−0.04s3+0.5s2+s+1010
由于 n ( t ) = 0.1 sin 10 t n(t)=0.1\sin10t n(t)=0.1sin10t,即 ω = 10 \omega=10 ω=10,则:
A [ Φ e n ( j ω ) ] = ∣ Φ e n ( j 10 ) ∣ = ∣ 10 − j 0.04 ω 3 + j ω − 0.5 ω 2 + 10 ∣ ω = 10 = 10 ( 10 − 50 ) 2 + ( 10 − 40 ) 2 = 10 50 = 0.2 ∠ [ Φ e n ( j ω ) ] = ∠ ( 10 − j 0.04 ω 3 + j ω − 0.5 ω 2 + 10 ∣ ω = 10 ) = ∠ ( 10 − 40 − j 30 ) = − 36.87 ° \begin{aligned} A[\Phi_{en}({\rm j}\omega)]&=|\Phi_{en}({\rm j}10)|=\left|\frac{10}{-{\rm j}0.04\omega^3+{\rm j}\omega-0.5\omega^2+10}\right|_{\omega=10}\\\\ &=\frac{10}{\sqrt{(10-50)^2+(10-40)^2}}=\frac{10}{50}=0.2\\\\ \angle[\Phi_{en}({\rm j}\omega)]&=\angle\left(\left.\displaystyle\frac{10}{-{\rm j}0.04\omega^3+{\rm j}\omega-0.5\omega^2+10}\right|_{\omega=10}\right)\\\\ &=\angle\left(\displaystyle\frac{10}{-40-{\rm j}30}\right)=-36.87° \end{aligned} A[Φen(jω)]∠[Φen(jω)]=∣Φen(j10)∣=∣
∣−j0.04ω3+jω−0.5ω2+1010∣
∣ω=10=(10−50)2+(10−40)210=5010=0.2=∠(−j0.04ω3+jω−0.5ω2+1010∣
∣ω=10)=∠(−40−j3010)=−36.87°
则扰动信号作用下,系统输出为:
c n ( t ) = 0.02 sin ( 10 t − 36.87 ° ) c_n(t)=0.02\sin(10t-36.87°) cn(t)=0.02sin(10t−36.87°)
故稳态误差的最大值为:
max ( e s s ( ∞ ) ) = 1 + 0.2 × 0.1 = 1.02 \max(e_{ss}(\infty))=1+0.2\times0.1=1.02 max(ess(∞))=1+0.2×0.1=1.02
【初始状态为零的输出响应曲线】
Example 3.37
已知单位反馈控制系统的开环传递函数为: G ( s ) = 22 ( s + 1 ) ( s + 3 ) G(s)=\displaystyle\frac{22}{(s+1)(s+3)} G(s)=(s+1)(s+3)22,系统输入量为: r ( t ) r(t) r(t),输出量为: c ( t ) c(t) c(t).
要求:
- 当 r ( t ) = 1 ( t ) r(t)=1(t) r(t)=1(t)时, c ( t ) c(t) c(t)的最大值和稳态值;
- 当 r ( t ) = A sin ω t r(t)=A\sin\omega{t} r(t)=Asinωt时,选取 ω \omega ω值使系统稳态输出的振幅最大,并求此最大振幅.
解:
-
当 r ( t ) = 1 ( t ) r(t)=1(t) r(t)=1(t)时, c ( t ) c(t) c(t)的最大值和稳态值;
根据题意可得,系统闭环传递函数为:
Φ ( s ) = G ( s ) 1 + G ( s ) = 22 s 2 + 4 s + 25 \Phi(s)=\frac{G(s)}{1+G(s)}=\frac{22}{s^2+4s+25} Φ(s)=1+G(s)G(s)=s2+4s+2522
当输入信号为: r ( t ) = 1 ( t ) r(t)=1(t) r(t)=1(t),即 R ( s ) = 1 s R(s)=\displaystyle\frac{1}{s} R(s)=s1时,
C ( s ) = Φ ( s ) R ( s ) = 22 s ( s 2 + 4 s + 25 ) C(s)=\Phi(s)R(s)=\frac{22}{s(s^2+4s+25)} C(s)=Φ(s)R(s)=s(s2+4s+25)22
根据终值定理可得:
c ( ∞ ) = lim s → 0 s C ( s ) = lim s → 0 s ⋅ 22 s ( s 2 + 4 s + 25 ) = 0.88 c(\infty)=\lim_{s\to0}sC(s)=\lim_{s\to0}s·\frac{22}{s(s^2+4s+25)}=0.88 c(∞)=s→0limsC(s)=s→0lims⋅s(s2+4s+25)22=0.88
由系统闭环传递函数可得,系统自然频率和阻尼比为:
ω n = 5 , ζ = 0.4 \omega_n=5,\zeta=0.4 ωn=5,ζ=0.4
则系统的超调量为:
σ % = e − π ζ / 1 − ζ 2 × 100 % = 25.4 % \sigma\%={\rm e}^{-\pi\zeta/\sqrt{1-\zeta^2}}\times100\%=25.4\% σ%=e−πζ/1−ζ2×100%=25.4%
系统的输出最大值为:
c max = c ( ∞ ) × ( 1 + σ % ) = 1.1 c_{\max}=c(\infty)\times(1+\sigma\%)=1.1 cmax=c(∞)×(1+σ%)=1.1
故 r ( t ) = 1 ( t ) r(t)=1(t) r(t)=1(t)时, c ( t ) c(t) c(t)的最大值为 1.1 1.1 1.1,稳态值为 0.88 0.88 0.88。 -
当 r ( t ) = A sin ω t r(t)=A\sin\omega{t} r(t)=Asinωt时,选取 ω \omega ω值使系统稳态输出的振幅最大,并求此最大振幅.
当 r ( t ) = A sin ω t r(t)=A\sin\omega{t} r(t)=Asinωt时,
∣ Φ ( j ω ) ∣ = ∣ 22 ( 25 − ω 2 ) + j 4 ω ∣ = 22 ( 25 − ω 2 ) 2 + 16 ω 2 = 22 ( ω 2 − 17 ) 2 + 336 |\Phi({\rm j}\omega)|=\left|\frac{22}{(25-\omega^2)+{\rm j}4\omega}\right|=\frac{22}{\sqrt{(25-\omega^2)^2+16\omega^2}}=\frac{22}{\sqrt{(\omega^2-17)^2+336}} ∣Φ(jω)∣=∣ ∣(25−ω2)+j4ω22∣ ∣=(25−ω2)2+16ω222=(ω2−17)2+33622
即当 ω 2 = 17 \omega^2=17 ω2=17时,系统稳态输出的振幅最大,该最大振幅为:
A max = A × 22 336 = 1.2 A A_{\max}=A\times\frac{22}{\sqrt{336}}=1.2A Amax=A×33622=1.2A -
【系统的时间响应】
Example 3.38
一复合控制系统如下图所示,图中 G 0 ( s ) = 10 s ( 0.1 s + 1 ) ( 0.5 s + 1 ) G_0(s)=\displaystyle\frac{10}{s(0.1s+1)(0.5s+1)} G0(s)=s(0.1s+1)(0.5s+1)10,如果要求系统在加速度输入信号 r ( t ) = 1 2 t 2 r(t)=\displaystyle\frac{1}{2}t^2 r(t)=21t2作用下无稳态误差,如何设计物理上可实现的 G r ( s ) G_r(s) Gr(s)?
解:
根据梅森增益公式,可求系统的闭环传递函数为:
Φ ( s ) = G r ( s ) G 0 ( s ) + G 0 ( s ) 1 + G 0 ( s ) \Phi(s)=\frac{G_r(s)G_0(s)+G_0(s)}{1+G_0(s)} Φ(s)=1+G0(s)Gr(s)G0(s)+G0(s)
则系统的误差传递函数为:
Φ e ( s ) = 1 − Φ ( s ) = 1 − G r ( s ) G 0 ( s ) 1 + G 0 ( s ) = s ( 0.1 s + 1 ) ( 0.5 s + 1 ) − 10 G r ( s ) s ( 0.1 s + 1 ) ( 0.5 s + 1 ) + 10 \Phi_e(s)=1-\Phi(s)=\frac{1-G_r(s)G_0(s)}{1+G_0(s)}=\frac{s(0.1s+1)(0.5s+1)-10G_r(s)}{s(0.1s+1)(0.5s+1)+10} Φe(s)=1−Φ(s)=1+G0(s)1−Gr(s)G0(s)=s(0.1s+1)(0.5s+1)+10s(0.1s+1)(0.5s+1)−10Gr(s)
当输入信号为加速度信号 r ( t ) = 1 2 t 2 r(t)=\displaystyle\frac{1}{2}t^2 r(t)=21t2,即 R ( s ) = 1 s 3 R(s)=\displaystyle\frac{1}{s^3} R(s)=s31时,利用终值定理求解系统稳态误差,有:
e s s ( ∞ ) = lim s → 0 s Φ e ( s ) R ( s ) = lim s → 0 s ( 0.1 s + 1 ) ( 0.5 s + 1 ) − 10 G r ( s ) s ( 0.1 s + 1 ) ( 0.5 s + 1 ) + 10 ⋅ 1 s 2 e_{ss}(\infty)=\lim_{s\to0}s\Phi_e(s)R(s)=\lim_{s\to0}\frac{s(0.1s+1)(0.5s+1)-10G_r(s)}{s(0.1s+1)(0.5s+1)+10}·\frac{1}{s^2} ess(∞)=s→0limsΦe(s)R(s)=s→0lims(0.1s+1)(0.5s+1)+10s(0.1s+1)(0.5s+1)−10Gr(s)⋅s21
设 F ( s ) = s ( 0.1 s + 1 ) ( 0.5 s + 1 ) − 10 G r ( s ) F(s)=s(0.1s+1)(0.5s+1)-10G_r(s) F(s)=s(0.1s+1)(0.5s+1)−10Gr(s),若要求系统在加速度输入信号 r ( t ) = 1 2 t 2 r(t)=\displaystyle\frac{1}{2}t^2 r(t)=21t2作用下无稳态误差,则有 F ( s ) F(s) F(s)不含低于 s s s三次幂的项.由于:
F ( s ) = s ( 0.1 s + 1 ) ( 0.5 s + 1 ) − 10 G r ( s ) = 0.05 s 3 + 0.6 s 2 + s − 10 G r ( s ) F(s)=s(0.1s+1)(0.5s+1)-10G_r(s)=0.05s^3+0.6s^2+s-10G_r(s) F(s)=s(0.1s+1)(0.5s+1)−10Gr(s)=0.05s3+0.6s2+s−10Gr(s)
故选择
G r ( s ) = 0.6 s 2 + s 10 = 0.06 s 2 + 0.1 s G_r(s)=\frac{0.6s^2+s}{10}=0.06s^2+0.1s Gr(s)=100.6s2+s=0.06s2+0.1s
可使系统在加速度输入信号 r ( t ) = 1 2 t 2 r(t)=\displaystyle\frac{1}{2}t^2 r(t)=21t2作用下无稳态误差.
由于无源网络无法提供纯微分信号,于是需增加一项因式 T s + 1 Ts+1 Ts+1.取 G r ( s ) = a s 2 + b s T s + 1 G_r(s)=\displaystyle\frac{as^2+bs}{Ts+1} Gr(s)=Ts+1as2+bs,则在物理上可考虑采用测速发电机或无源网络组合线路近似实现上述补偿方案.
【系统时间响应曲线】
Example 3.39
某非单位反馈控制系统如下图所示,图中 G 1 ( s ) G_1(s) G1(s)的单位阶跃响应为: 8 5 ( 1 − e − 5 t ) \displaystyle\frac{8}{5}(1-{\rm e}^{-5t}) 58(1−e−5t).
- 若 r ( t ) = 20 ⋅ 1 ( t ) r(t)=20·1(t) r(t)=20⋅1(t),求系统稳态输出 c ( ∞ ) c(\infty) c(∞);
- 若 r ( t ) = 20 ⋅ 1 ( t ) r(t)=20·1(t) r(t)=20⋅1(t),求系统的超调量 σ % \sigma\% σ%、调节时间 t s t_s ts和稳态误差 e s s ( ∞ ) e_{ss}(\infty) ess(∞),并绘制输出响应曲线;
- 若 n ( t ) n(t) n(t)为可测量的阶跃扰动信号,为消除扰动对稳态输出的影响,设计顺馈补偿装置 G n ( s ) G_n(s) Gn(s),并画出相应的系统结构图;
解:
由于 G 1 ( s ) G_1(s) G1(s)的单位阶跃响应为: 8 5 ( 1 − e − 5 t ) \displaystyle\frac{8}{5}(1-{\rm e}^{-5t}) 58(1−e−5t),则:
1 s G 1 ( s ) = 8 5 ( 1 s − 1 s + 5 ) \frac{1}{s}G_1(s)=\frac{8}{5}\left(\frac{1}{s}-\frac{1}{s+5}\right) s1G1(s)=58(s1−s+51)
即:
G 1 ( s ) = 8 s + 5 G_1(s)=\frac{8}{s+5} G1(s)=s+58
-
若 r ( t ) = 20 ⋅ 1 ( t ) r(t)=20·1(t) r(t)=20⋅1(t),求系统稳态输出 c ( ∞ ) c(\infty) c(∞);
由图可得,系统的开环传递函数为:
G ( s ) H ( s ) = 8 s + 5 ⋅ 2 s + 4 ⋅ 5 = 50 ( s + 5 ) ( s + 4 ) G(s)H(s)=\frac{8}{s+5}·\frac{2}{s+4}·5=\frac{50}{(s+5)(s+4)} G(s)H(s)=s+58⋅s+42⋅5=(s+5)(s+4)50
则系统的闭环传递函数为:
Φ e ( s ) = 16 s 2 + 9 s + 100 \Phi_e(s)=\frac{16}{s^2+9s+100} Φe(s)=s2+9s+10016
若 r ( t ) = 20 ⋅ 1 ( t ) r(t)=20·1(t) r(t)=20⋅1(t),即 R ( s ) = 20 s R(s)=\displaystyle\frac{20}{s} R(s)=s20,由终值定理可得:
c ( ∞ ) = lim s → 0 s Φ r ( s ) R ( s ) = lim s → 0 s ⋅ 16 s 2 + 9 s + 100 ⋅ 20 s = 3.2 c(\infty)=\lim_{s\to0}s\Phi_r(s)R(s)=\lim_{s\to0}s·\frac{16}{s^2+9s+100}·\frac{20}{s}=3.2 c(∞)=s→0limsΦr(s)R(s)=s→0lims⋅s2+9s+10016⋅s20=3.2 -
若 r ( t ) = 20 ⋅ 1 ( t ) r(t)=20·1(t) r(t)=20⋅1(t),求系统的超调量 σ % \sigma\% σ%、调节时间 t s t_s ts和稳态误差 e s s ( ∞ ) e_{ss}(\infty) ess(∞),并绘制输出响应曲线;
由 Φ r ( s ) = 16 s 2 + 9 s + 100 \Phi_r(s)=\displaystyle\frac{16}{s^2+9s+100} Φr(s)=s2+9s+10016可知,系统的自然频率和阻尼比为:
ω n = 10 , ζ = 0.45 \omega_n=10,\zeta=0.45 ωn=10,ζ=0.45
则系统的动态性能指标如下:
超调量: σ % = e − π ζ / 1 − ζ 2 × 100 % = 20.5 % 峰值时间: t p = π ω n 1 − ζ 2 = 0.352 s 调节时间: t s = 4.4 ζ ω n = 0.98 s ( Δ = 2 % ) \begin{aligned} &超调量:\sigma\%={\rm e}^{-\pi\zeta/\sqrt{1-\zeta^2}}\times100\%=20.5\%\\\\ &峰值时间:t_p=\frac{\pi}{\omega_n\sqrt{1-\zeta^2}}=0.352{\rm s}\\\\ &调节时间:t_s=\frac{4.4}{\zeta\omega_n}=0.98{\rm s}(\Delta=2\%) \end{aligned} 超调量:σ%=e−πζ/1−ζ2×100%=20.5%峰值时间:tp=ωn1−ζ2π=0.352s调节时间:ts=ζωn4.4=0.98s(Δ=2%)
由 G ( s ) H ( s ) = 80 ( s + 5 ) ( s + 4 ) G(s)H(s)=\displaystyle\frac{80}{(s+5)(s+4)} G(s)H(s)=(s+5)(s+4)80可知,该系统为 0 0 0型系统,且 K p = 4 K_p=4 Kp=4,等效单位反馈系统如下图所示:
则系统在输入 r ( t ) = 20 ⋅ 1 ( t ) r(t)=20·1(t) r(t)=20⋅1(t)下的稳态误差为:
e s s ( ∞ ) = R ′ 1 + K p = 0.8 ( R ′ = 4 ) e_{ss}(\infty)=\frac{R'}{1+K_p}=0.8(R'=4) ess(∞)=1+KpR′=0.8(R′=4)
【系统输出响应曲线】
-
若 n ( t ) n(t) n(t)为可测量的阶跃扰动信号,为消除扰动对稳态输出的影响,设计顺馈补偿装置 G n ( s ) G_n(s) Gn(s),并画出相应的系统结构图;
在扰动端加入顺馈补偿装置 G n ( s ) G_n(s) Gn(s)的系统结构图如下图所示:
根据梅森增益公式可得系统的传递函数为:
C ( s ) N ( s ) = G 2 ( s ) + G n ( s ) G 1 ( s ) G 2 ( s ) 1 + G 1 G 2 ( s ) H ( s ) \frac{C(s)}{N(s)}=\frac{G_2(s)+G_n(s)G_1(s)G_2(s)}{1+G_1G_2(s)H(s)} N(s)C(s)=1+G1G2(s)H(s)G2(s)+Gn(s)G1(s)G2(s)
由于 n ( t ) n(t) n(t)为可量测的阶跃扰动信号,为消除扰动对稳态输出的影响,有:
G n ( s ) = − 1 G 1 ( s ) = − s + 5 8 G_n(s)=-\frac{1}{G_1(s)}=-\frac{s+5}{8} Gn(s)=−G1(s)1=−8s+5
故当引入顺馈补偿装置: G n ( s ) = − s + 5 8 G_n(s)=-\displaystyle\frac{s+5}{8} Gn(s)=−8s+5时可消除扰动对稳态输出的影响.
Example 3.40
一单位反馈系统,其开环传递函数为: G ( s ) = 100 s ( 0.1 s + 1 ) G(s)=\displaystyle\frac{100}{s(0.1s+1)} G(s)=s(0.1s+1)100,输入信号为: r ( t ) = 3.5 sin 2.8 t r(t)=3.5\sin2.8t r(t)=3.5sin2.8t,求系统的稳态误差.
解:
由系统的开环传递函数可得系统的闭环传递函数为:
Φ ( s ) = G ( s ) 1 + G ( s ) = 100 0.1 s 2 + s + 100 \Phi(s)=\frac{G(s)}{1+G(s)}=\frac{100}{0.1s^2+s+100} Φ(s)=1+G(s)G(s)=0.1s2+s+100100
则系统的误差传递函数为:
Φ e ( s ) = 1 − Φ ( s ) = 0.1 s 2 + s 0.1 s 2 + s + 100 \Phi_e(s)=1-\Phi(s)=\frac{0.1s^2+s}{0.1s^2+s+100} Φe(s)=1−Φ(s)=0.1s2+s+1000.1s2+s
由于输入信号为: r ( t ) = 3.5 sin 2.8 t r(t)=3.5\sin2.8t r(t)=3.5sin2.8t,则:
∣ Φ e ( j 2.8 ) ∣ = ∣ − 0.1 ω 2 + j ω 100 − 0.1 ω 2 + j ω ∣ ω = 2.8 = ( 0.1 × 2. 8 2 ) + 2. 8 2 ( 100 − 0.1 × 2. 8 2 ) 2 + 2. 8 2 = 0.03 ∠ [ Φ e ( j 2.8 ) ] = ∠ ( − 0.1 ω 2 + j ω 2 100 − 0.1 ω 2 + j ω ∣ ω = 2.8 ) = ∠ ( − 0.784 + j 2.8 99.216 + j 2.8 ) = − 2.3 ° \begin{aligned} &|\Phi_e({\rm j}2.8)|=\left|\displaystyle\frac{-0.1\omega^2+{\rm j}\omega}{100-0.1\omega^2+{\rm j}\omega}\right|_{\omega=2.8}=\frac{\sqrt{(0.1\times2.8^2)+2.8^2}}{\sqrt{(100-0.1\times2.8^2)^2+2.8^2}}=0.03\\\\ &\angle[\Phi_e({\rm j}2.8)]=\angle\left(\left.\displaystyle\frac{-0.1\omega^2+{\rm j}\omega^2}{100-0.1\omega^2+{\rm j}\omega}\right|_{\omega=2.8}\right)=\angle\left(\displaystyle\frac{-0.784+{\rm j}2.8}{99.216+{\rm j}2.8}\right)=-2.3° \end{aligned} ∣Φe(j2.8)∣=∣
∣100−0.1ω2+jω−0.1ω2+jω∣
∣ω=2.8=(100−0.1×2.82)2+2.82(0.1×2.82)+2.82=0.03∠[Φe(j2.8)]=∠(100−0.1ω2+jω−0.1ω2+jω2∣
∣ω=2.8)=∠(99.216+j2.8−0.784+j2.8)=−2.3°
故系统的稳态误差为: e s s ( ∞ ) = 0.105 sin ( 2.8 t − 2.3 ° ) e_{ss}(\infty)=0.105\sin(2.8t-2.3°) ess(∞)=0.105sin(2.8t−2.3°).
【输出响应曲线】
Example 3.41
设位置随动系统如下图所示,图中: θ i \theta_i θi为输入轴转角, θ o \theta_o θo为输出轴转角, K 1 = 500 , K 2 = 0.1 , K_1=500,K_2=0.1, K1=500,K2=0.1, K 3 = 1 , T = 0.003 K_3=1,T=0.003 K3=1,T=0.003.
求:
- 闭环传递函数 Θ o ( s ) / Θ i ( s ) \Theta_o(s)/\Theta_i(s) Θo(s)/Θi(s);
- 当输入轴以角速度 ω r = 10 \omega_r=10 ωr=10转动时,输出转角的稳态误差 Δ θ s s ( ∞ ) \Delta\theta_{ss}(\infty) Δθss(∞);
解:
-
闭环传递函数 Θ o ( s ) / Θ i ( s ) \Theta_o(s)/\Theta_i(s) Θo(s)/Θi(s);
由图可得,系统的开环传递函数为:
G ( s ) = K 1 K 2 T s 2 + s + K 3 s G(s)=\frac{K_1K_2}{Ts^2+s+K_3s} G(s)=Ts2+s+K3sK1K2
则系统的闭环传递函数为:
Θ o ( s ) Θ i ( s ) = K 1 K 2 T s 2 + ( 1 + K 3 ) s + K 1 K 2 = 50 0.003 s 2 + 2 s + 50 \frac{\Theta_o(s)}{\Theta_i(s)}=\frac{K_1K_2}{Ts^2+(1+K_3)s+K_1K_2}=\frac{50}{0.003s^2+2s+50} Θi(s)Θo(s)=Ts2+(1+K3)s+K1K2K1K2=0.003s2+2s+5050 -
当输入轴以角速度 ω r = 10 \omega_r=10 ωr=10转动时,输出转角的稳态误差 Δ θ s s ( ∞ ) \Delta\theta_{ss}(\infty) Δθss(∞);
当输入轴以角速度 ω r = 10 \omega_r=10 ωr=10转动时,输入为: Θ i ( s ) = 10 s 2 \Theta_i(s)=\displaystyle\frac{10}{s^2} Θi(s)=s210;系统的误差函数为:
Δ Θ ( s ) = Θ o ( s ) − Θ i ( s ) = 0.003 s 2 + 2 s 0.003 s 2 + 2 s + 50 Θ i ( s ) \Delta\Theta(s)=\Theta_o(s)-\Theta_i(s)=\frac{0.003s^2+2s}{0.003s^2+2s+50}\Theta_i(s) ΔΘ(s)=Θo(s)−Θi(s)=0.003s2+2s+500.003s2+2sΘi(s)
由终值定理可得,系统的稳态误差为:
Δ θ s s ( ∞ ) = lim s → 0 s Δ Θ ( s ) = lim s → 0 s ⋅ 0.003 s 2 + 2 s 0.003 s 2 + 2 s + 50 ⋅ 10 s 2 = 0.4 \Delta\theta_{ss}(\infty)=\lim_{s\to0}s\Delta\Theta(s)=\lim_{s\to0}s·\frac{0.003s^2+2s}{0.003s^2+2s+50}·\frac{10}{s^2}=0.4 Δθss(∞)=s→0limsΔΘ(s)=s→0lims⋅0.003s2+2s+500.003s2+2s⋅s210=0.4
故当输入轴以角速度 ω r = 10 \omega_r=10 ωr=10转动时,输出转角的稳态误差 Δ θ s s ( ∞ ) = 0.4 \Delta\theta_{ss}(\infty)=0.4 Δθss(∞)=0.4。 -
【输出响应曲线】