2018南京航天航空大学820自动控制原理第九题

2018南京航天航空大学820自动控制原理第九题_第1张图片
1.系统 M ( s ) M(s) M(s)的状态空间方程:
{ x ˙ 1 = − 4 x ˙ 1 + u y = x ˙ 1 \left\{\begin{array}{l} \boldsymbol{\dot x_1} = -4 \boldsymbol{\dot x_1} + \boldsymbol u \\ \boldsymbol{y} = \boldsymbol{\dot x_1} \end{array}\right. {x˙1=4x˙1+uy=x˙1
系统 N ( s ) N(s) N(s)的状态空间方程:
{ [ x ˙ 2 x ˙ 3 ] = [ − 3 0 0 − 2 ] [ x 2 x 3 ] + [ 1 1 ] u y = [ − 2 4 ] x \left\{\begin{array}{l} \begin{bmatrix}\boldsymbol{\dot x_2} \\ \boldsymbol{\dot x_3} \end{bmatrix} = \begin{bmatrix} -3 & 0 \\ 0 &-2\end{bmatrix} \begin{bmatrix}\boldsymbol{x_2} \\ \boldsymbol{x_3}\end{bmatrix} + \begin{bmatrix} 1 \\ 1 \end{bmatrix} \boldsymbol u \\ \boldsymbol y = \begin{bmatrix} -2 & 4 \end{bmatrix} \boldsymbol x \end{array}\right. [x˙2x˙3]=[3002][x2x3]+[11]uy=[24]x
两个子系统串联,存在零极点对消,系统不能控。系统串联后的状态空间方程
{ x ˙ = [ − 4 − 2 4 0 − 3 0 0 0 − 2 ] x + [ 0 1 1 ] u y = [ 1 0 0 ] x \left\{\begin{array}{l} \boldsymbol{\dot x} = \begin{bmatrix} -4 & -2 & 4 \\ 0 & -3 & 0 \\ 0 & 0 & -2 \end{bmatrix} \boldsymbol x + \begin{bmatrix} 0 \\ 1 \\ 1 \end{bmatrix} \boldsymbol u \\ \boldsymbol y= \begin{bmatrix} 1 & 0 & 0 \end{bmatrix} \boldsymbol x \end{array}\right. x˙=400230402x+011uy=[100]x
2.系统的单位阶跃响应近似为一个调节时间 t s t_s ts 4.75 s 4.75s 4.75s的临界阻尼系统,即系统存在两个相同的极点,可得
t s = 4 ξ ω n = 4.75 s t_s=\tfrac{4}{\xi\omega_n}=4.75s ts=ξωn4=4.75s
解得
s 1 , 2 = 0.8421 s_{1,2}=0.8421 s1,2=0.8421
由1可知,系统状态 x 1 x_1 x1不能控,设状态反馈后系统的特征方程为
Δ ‘ = ( s + 0.8421 ) 2 = s 2 + 1.6842 s + 0.70913241 = 0 \begin{aligned} \Delta^` & =(s+0.8421)^2 \\ & = s^2 + 1.6842s + 0.70913241 \\ & = 0 \end{aligned} Δ=(s+0.8421)2=s2+1.6842s+0.70913241=0
设状态反馈矩阵为 K = [ K 1 K 2 ] \boldsymbol K = \begin{bmatrix}K_1 & K_2 \end{bmatrix} K=[K1K2],对应系统的特征多项式为
Δ = d e t [ s I − ( A − B K ) ] = ( s + 3 + K 1 ) ( s + 2 + K 2 ) − K 1 K 2 = s 2 + ( 5 + K 1 + K 2 ) s + 6 + 2 K 1 + 3 K 2 = 0 \begin{aligned} \Delta &= det[s\boldsymbol{I} - (\boldsymbol{A} - \boldsymbol{B}\boldsymbol{K})] \\ & = (s+3+K_1)(s+2+K_2) - K_1 K_2 \\ & = s^2 + (5 + K_1 + K_2)s + 6 + 2K_1 + 3K_2 \\ & =0 \end{aligned} Δ=det[sI(ABK)]=(s+3+K1)(s+2+K2)K1K2=s2+(5+K1+K2)s+6+2K1+3K2=0
对比两式可解得
{ K 1 = − 4.6565321 K 2 = 1.34073241 \left\{\begin{array}{l} K_1 = -4.6565321 \\ K_2 = 1.34073241 \end{array}\right. {K1=4.6565321K2=1.34073241
3.系统的能观性矩阵为
Q o = [ 1 0 0 − 4 − 2 4 16 14 − 24 ] \boldsymbol Q_o = \begin{bmatrix} 1 & 0 & 0 \\ -4 & -2 & 4 \\ 16 & 14 & -24 \end{bmatrix} Qo=141602140424
r a n k Q o = 3 rank\boldsymbol Q_o =3 rankQo=3,系统完全能观,可以任意配置极点调节状态观测器的动态性能。 r a n k C = 1 rank\boldsymbol C =1 rankC=1,可以设计 n − m = 2 n - m = 2 nm=2维观测器,将极点配置到 s 1 , 2 = 0.8421 s_{1,2}=0.8421 s1,2=0.8421来实现状态反馈控制。

你可能感兴趣的:(现代控制理论)