MASTERY · CONTROL ENGINEERING

1996 年上海交通大学 816 自动控制原理 参考答案

1.8k WORDS7 MIN READ#考研

1996 年上海交通大学 816 自动控制原理 参考答案

来源:1996-2023年上海交通大学答案.pdf 第 1-7 页 说明:逐字录入原件文字,原扫描件斜向红色水印已剔除。 注:第二题解中”由 k<0k<0 所作的奈氏图”原文如此(结合上下文应为 k>0k>0);第七题原真题试卷中输入矩阵 BB 遗漏第三行元素,参考答案按 B=(123)B=\begin{pmatrix} 1 \\ 2 \\ 3 \end{pmatrix} 求解。


上海交通大学一九九六年攻读硕士学位研究生

入学考试试题参考答案

注意:答案必须写在答题纸上,写在试卷或草稿纸上均无效。

一、【分析】此题主要考察状态空间表达式转化为传递函数的方法,理解阶跃响应的定义,学会使用终值定理。

【解】

(1)

Y(s)U(s)=C[sIA]1B=(10)(s12ks+1+k)1(0k)=ks2+(1+k)s+2k\dfrac{Y(s)}{U(s)}=C[sI-A]^{-1}B=\begin{pmatrix} 1 & 0 \end{pmatrix}\begin{pmatrix} s & -1 \\ 2k & s+1+k \end{pmatrix}^{-1}\begin{pmatrix} 0 \\ k \end{pmatrix}=\dfrac{k}{s^2+(1+k)s+2k}

(2)设 r(t)=M1(t)r(t)=M\cdot 1(t),则 R(s)=MsR(s)=\dfrac{M}{s}

Y(s)U(s)=C(s)R(s)=输出输入\dfrac{Y(s)}{U(s)}=\dfrac{C(s)}{R(s)}=\dfrac{\text{输出}}{\text{输入}}

c()=lims0sR(s)C(s)R(s)=lims0sR(s)ks2+(1+k)s+2kc(\infty)=\lim_{s\to 0}s\cdot R(s)\cdot\dfrac{C(s)}{R(s)}=\lim_{s\to 0}s\cdot R(s)\cdot\dfrac{k}{s^2+(1+k)s+2k}

=lims0sMsks2+(1+k)s+2k=lims0M12=2M=4r(t)=4=\lim_{s\to 0}s\cdot\dfrac{M}{s}\cdot\dfrac{k}{s^2+(1+k)s+2k}=\lim_{s\to 0}M\cdot\dfrac{1}{2}=2 \qquad M=4 \qquad r(t)=4

【评注】

  1. 本题考查现代控制理论中,由状态空间求传递函数

当系统的状态空间表达式为 x˙=Ax+Buy=Cx+Du\dot{x}=Ax+Bu \quad y=Cx+Du 时,在零初始条件下,则

G(s)=Y(s)U(s)=C[sIA]1B+DG(s)=\dfrac{Y(s)}{U(s)}=C[sI-A]^{-1}B+D

  1. 系统对阶跃信号的响应,lims0sC(s)=c()\lim\limits_{s\to 0}sC(s)=c(\infty) 应该使用该终值定理公式,由图可知,系统最终的响应数值为 2,因此可得 c()=2c(\infty)=2

二、【分析】本题主要考查求解等效开环传递函数的方法,以及用奈氏判据判断系统的稳定性。

【解】用到公式 Z=PRZ=P-R 其中 P=0P=0,为保证 Z=0Z=0,则必 R=0R=0

由系统框图,来得到系统的闭环特征多项式为:

D(s)=s3+3s2+3s+2+kD(s)=s^3+3s^2+3s+2+k

由闭环特征多项式得到系统等效开环传递函数为:

G(s)H(s)=ks3+3s2+3s+2G(s)H(s)=\dfrac{k}{s^3+3s^2+3s+2}

那么可得当 k>0k>0 时,有:

图

k7>1k<7-\dfrac{k}{7}>-1 \quad k<7

G(j0+)H(j0+)=k20G(j0^+)H(j0^+)=\dfrac{k}{2}\angle 0^\circ

G(j)H(j)=0270G(j\infty)H(j\infty)=0\angle -270^\circ

G(jω)H(jω)=k(23ω2)jk(3ωω3)(23ω2)2+(3ωω3)2G(j\omega)H(j\omega)=\dfrac{k(2-3\omega^2)-jk(3\omega-\omega^3)}{(2-3\omega^2)^2+(3\omega-\omega^3)^2}

与实轴交点:3ωω3=0ω=3Re=k73\omega-\omega^3=0 \qquad \omega=\sqrt{3} \quad Re=-\dfrac{k}{7}

与虚轴交点:3ω2+2=0ω=23Im=3372kj-3\omega^2+2=0 \qquad \omega=\sqrt{\dfrac{2}{3}} \quad \text{Im}=\dfrac{3\sqrt{3}}{7\sqrt{2}}kj

故而由 k<0k<0〔原文如此〕所作的奈氏图可知 kk 的范围为:0<k<70<k<7

【评注】1. 学会根据闭环特征多项式来得到等效开环传递函数,因为我们期待的可变参数 kk 一般希望它位于分子,充当增益。若 kk 本身位于分子,则无需求等效开环传递函数的步骤。

  1. 奈氏图的稳定性判断标准 Z=PRZ=P-R

ZZT0T_0 半开平面闭环极点的个数,PPT0T_0 半开平面开环极点的个数,RR 为包围 (1,j0)(-1, j0) 点的圈数,由等效开环传递函数可知 P=0P=0,那么 R=0R=0 时,才能保证 Z=0Z=0,即右半开平面无闭环极点,所以系统可保证稳定

  1. 学会画奈氏图,熟记其画图步骤,把握关键点,例如与实轴或虚轴的交点 渐近线

三、【分析】此题考查根轨迹的绘制,以及相关数据的计算。

【解】根据开环传递函数,应画出 00^\circ 根轨迹

①系统存在三个开环的极点 s1=0s_1=0s23=2±2js_{23}=-2\pm 2j

二个开环零点 s1=2s_1=-2s2=1s_2=1

②实轴根轨迹:[20][-2 \quad 0][1+)[1 \quad +\infty)

③系统极点 2+2j-2+2j 的出射角为:

θ=2kπ+146.3+9013590=11.3\theta=2k\pi+146.3^\circ+90^\circ-135^\circ-90^\circ=11.3^\circ

④系统与虚轴的交点 s=jωs=j\omega 代入

D(s)=s3+(4k1)s2+(8k1)s+2k1D(s)=s^3+(4-k_1)s^2+(8-k_1)s+2k_1

{2k1(4k1)ω2=0(8k1)ωω3=0{ω=0k1=0{ω=2.26k1=2.88\begin{cases} 2k_1-(4-k_1)\omega^2=0 \\ (8-k_1)\omega-\omega^3=0 \end{cases} \Rightarrow \begin{cases} \omega=0 \\ k_1=0 \end{cases} \text{或} \begin{cases} \omega=2.26 \\ k_1=2.88 \end{cases}

⑤与实轴的分离点

d((2+s)(1s)s(s2+4s+8))ds=0,s1=3.19,s2=3.8(舍),s3,4=0.69±j0.92(舍)\dfrac{d\left(\dfrac{(2+s)(1-s)}{s(s^2+4s+8)}\right)}{ds}=0, \quad s_1=3.19, \quad s_2=-3.8\text{(舍)}, \quad s_{3,4}=-0.69\pm j0.92\text{(舍)}

所以可画出系统的根轨迹为:

图

由图可得,当系统稳定时,0<k1<2.880<k_1<2.88,此时开环增益为 0<k<0.720<k<0.72

【评注】1. 学会区分画 00^\circ 根轨迹和 180180^\circ 根轨迹的判别标准,注意 00^\circ180180^\circ 根轨迹的不同大致有三点不同:①实轴上的根轨迹,②出射角以及入射角计算,③渐近线的角度问题

  1. 学会计算关键点①与虚轴交点,②分离点,出射角,入射角

四、【分析】由非线性系统所给出的相关条件,来画出系统相轨迹

【解】n=c¨n=\ddot{c}e=rc=Rce=r-c=-R-ce˙=c˙\dot{e}=-\dot{c}e¨=c¨\ddot{e}=-\ddot{c}

m=eβc˙=e+βe˙m=e-\beta\dot{c}=e+\beta\dot{e}

e(t)=r(t)c(t)\because e(t)=r(t)-c(t) 且原系统处于静止状态,即 c(0)=0c(0)=0

e(0)=r(0)=Re˙(0)=r˙(0)=0(R,0)为起始点\therefore e(0)=r(0)=-R \qquad \dot{e}(0)=\dot{r}(0)=0 \qquad (-R, 0)\text{为起始点}

故而可得系统的相轨迹

图

【评注】1. 由题目的条件,得到相轨迹的起始点。

  1. 注意题目的问法是画出 ee˙e-\dot{e} 还是 cc˙c-\dot{c} 的根轨迹。若不加说明,默认绘制 ee˙e-\dot{e} 相轨迹。相轨迹的画法极为简单,列出的式子只含有 e¨\ddot{e}e˙\dot{e}ee 即可,找出 ee˙e-\dot{e} 的函数关系。同样,应注明范围,进而寻找到相轨迹的开关线,关键点应相应的求出。画相轨迹应注明其走向,学会典型几种相轨迹的画法,另外对于复杂的方块图,应先学会化简,再进行画图,是最基本的思路方法。

五、【分析】此题主要考查伯德图及其相应的转折频率,相位裕量等基础知识

【解】(1)由图可知系统存在一个积分环节及两个惯性环节

可得 G(s)=ks(0.1s+1)(10s+1)G(s)=\dfrac{k}{s(0.1s+1)(10s+1)},由 I 型系统,初始频段和横坐标的交点可知 k=10k=10

G(s)=10s(0.1s+1)(10s+1)G(s)=\dfrac{10}{s(0.1s+1)(10s+1)}

ωc=1γ=180+φ(ωc)=18090arctan10ωcarctan0.1ωc=0\therefore \omega_c=1 \qquad \therefore \gamma=180^\circ+\varphi(\omega_c)=180^\circ-90^\circ-\arctan 10\omega_c-\arctan 0.1\omega_c=0^\circ

(2)现要求 γ=30\gamma'=30^\circ 也即取 ωc=ωm\omega_c'=\omega_m

使得 γ=18090arctan0.1ωmarctan10ωm=30ωm=0.167\gamma'=180^\circ-90^\circ-\arctan 0.1\omega_m-\arctan 10\omega_m=30^\circ \quad \therefore \omega_m=0.167

由(1)问的 L(ω)L(\omega) 的范围可知

k10ω2=1k=0.28\dfrac{k}{10\omega^2}=1 \qquad k=0.28

【评注】1. 伯德图是考研的重点题型,也是必考题型,通常较为基础,一般会和系统的校正相结合来考查。明确求出截止频率,相位裕度和幅值裕度的方法同事〔原文如此〕学会根据幅值裕度和相位裕度判断稳定性的方法,详情见交大教材 P148

  1. 必记伯德图的初始频段或初始频段的反向延长线必经过 (1,20lgk)(1, 20\lg k) 这一点,同时应理解,I 型或 II 型系统的初始段的延长线与实轴交点的意义,可见交大教材 P147 ωv=kv\omega_v=k_vωa=ka\omega_a=\sqrt{k_a} 的公式应熟记

  2. 学会列不同 ω\omega 范围取值下的 L(ω)L(\omega) 的表示方法,进而求出相对应的截止频率,求一个系统的传递函数,应判断该系统是否为最小相位系统。通常给出相频曲线时,应多考虑一步:该系统是否为最小相位系统,进而来求解题目

六、【分析】本题目考查能控性和能观性的判断方法及用状态反馈进行极点配置问题。

【解】(1)可控性矩阵 Pc=(b,Ab,A2b)=(111000124)rankPc=2<3P_c=(b, Ab, A^2b)=\begin{pmatrix} 1 & 1 & 1 \\ 0 & 0 & 0 \\ 1 & 2 & 4 \end{pmatrix} \quad \text{rank}P_c=2<3 所以不完全可控

可观性矩阵 Po=(CCACA2)=(111102114)rankPo=3=nP_o=\begin{pmatrix} C \\ CA \\ CA^2 \end{pmatrix}=\begin{pmatrix} 1 & 1 & 1 \\ 1 & 0 & 2 \\ 1 & 1 & 4 \end{pmatrix} \quad \text{rank}P_o=3=n 所以系统

可观测

(2)对系统进行能控性分解

P1=(110001120)P=(201101010)PAP1=(022131001)P^{-1}=\begin{pmatrix} 1 & 1 & 0 \\ 0 & 0 & 1 \\ 1 & 2 & 0 \end{pmatrix} \qquad P=\begin{pmatrix} 2 & 0 & -1 \\ -1 & 0 & 1 \\ 0 & 1 & 0 \end{pmatrix} \qquad PAP^{-1}=\begin{pmatrix} 0 & -2 & 2 \\ 1 & 3 & -1 \\ 0 & 0 & -1 \end{pmatrix}

可知不可控极点为 1-1,和所要配置的极点 1-1 相同,故存在配置的可能。

设状态反馈阵为 (k1,k2,k3)(k_1, k_2, k_3) 那么 sI(Abk)=(s+1)(s+2)2|sI-(A-bk)|=(s+1)(s+2)^2 也即

(s+1)[(s+k11)(s+k32)k1k3]=(s+1)(s+2)2(s+1)\left[(s+k_1-1)(s+k_3-2)-k_1 k_3\right]=(s+1)(s+2)^2

\begin{cases} k_1+k_3-3=4 \\ 2-k_3-2k_1=4 \end{cases} \Rightarrow \begin{cases} k_1=-9 \\ k_3=16 \end{cases} \qquad k=(-9, 0, 16)\text{ 其中 } k_2 \text{ 可任意值,\quad 这里 } k_2=0$$ 【评析】1. 对于系统的可观性和可控性的判断是最基本的题型,主要用当两个公式①可控性 $P_c=(b, Ab, \cdots, A^{n-1}b)$,②可观性 $P_o=\begin{pmatrix} C \\ CA \\ \vdots \\ CA^{n-1} \end{pmatrix}$ 系统一般最多为 3 维,也即 $n=3$ 2. 对于配置极点的问题一般分为系统完全可控和系统不完全可控两种情况,对于完全能控,只需按教材的操作方法做即可。这里重点讲一下不完全可控的系统,当不可控极点与期待配置的极点相同,那么系统存在进行极点配置的可能,96 年真题与 07 年真题中的现代控制理论部分具有异曲同工之妙。不难发现,同一个知识点,十一年后又出现在了考卷上,这也要求同学会把握该题型的做题思路,同时学会对系统进行能控性分解找到不可控极点。 七、【分析】本题考查约当标准型的基本代法。 【解】$$|\lambda I-A|=\begin{vmatrix} s-2 & -4 & -5 \\ 0 & s-1 & 0 \\ 0 & 0 & s-1 \end{vmatrix}=0 \Rightarrow s_1=2, \; s_2=s_3=1 \qquad \widetilde{A}=\begin{pmatrix} 1 & * & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 2 \end{pmatrix}

Rank(IA)=rank(145000000)=1\text{Rank}(I-A)=\text{rank}\begin{pmatrix} -1 & -4 & -5 \\ 0 & 0 & 0 \\ 0 & 0 & 0 \end{pmatrix}=1

J1=(101), 其中 u[nrank(IA)]=2(31)=0, 位置的数为 0J_1=\begin{pmatrix} 1 & * \\ 0 & 1 \end{pmatrix}, \text{ 其中 } u-[n-\text{rank}(I-A)]=2-(3-1)=0, \quad \therefore * \text{ 位置的数为 } 0

Jordan 标准型为 x~˙(t)=(100010002)x~(t)+(2324)u(t)\therefore \text{Jordan 标准型为 } \dot{\tilde{x}}(t)=\begin{pmatrix} 1 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 2 \end{pmatrix}\tilde{x}(t)+\begin{pmatrix} 2 \\ 3 \\ 24 \end{pmatrix}u(t)

【评析】本题为常规的约旦标准型的化简方法,是现代控制理论的基础部分,该部分和线性代数结合较多,因此也要求同学会拥有较好的数学基础和良好的计算能力,该题难度系数较低,只要做好相关习题便可以熟练掌握。

Discussion

Comments

Post

Thoughts, corrections, and follow-up notes are welcome here.