MASTERY · CONTROL ENGINEERING

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

2.3k WORDS11 MIN READ#考研

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

来源:1996-2023年上海交通大学答案.pdf 第 41-47 页 说明:逐字录入原件文字,原扫描件斜向红色水印已剔除。 注:原答案第一题对于系统(2)求 C2(s)C_2(s) 时,分子印为 k+5k+5(上一式分子为 k(s+5)k(s+5)),正文照录。 注:原答案第八题【标注】第 2 条中”化标准型的维多一般不超过 3 维”印为”维多”,正文照录。


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

入学考试试题参考答案

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

一、【分析】本题考查单位阶跃输入条件下,两系统的输出稳态值的计算方法

【解】对于系统(1)来说,

C1(s)R1(s)=k(s3)s(s+5)1+ks(s+5)=k(s3)s2+5s+k\frac{C_1(s)}{R_1(s)}=\frac{\dfrac{k(s-3)}{s(s+5)}}{1+\dfrac{k}{s(s+5)}}=\frac{k(s-3)}{s^2+5s+k}

C1(s)=1sk(s3)s2+5s+kC_1(s)=\frac{1}{s}\cdot\frac{k(s-3)}{s^2+5s+k}

系统要有稳态输出,则系统必须 稳定· 经判断,要使系统稳定,则必须 k>0k>0

c1()=lims0sC1(s)=lims0k(s3)s2+5s+k=3 (其中k>0)c_1(\infty)=\lim_{s\to 0}sC_1(s)=\lim_{s\to 0}\frac{k(s-3)}{s^2+5s+k}=-3\ (\text{其中}k>0)

对于系统(2)来说,

C2(s)R2(s)=ks(s3)1+ks(s+5)=k(s+5)s(s+5)(s3)+k(s3)\frac{C_2(s)}{R_2(s)}=\frac{\dfrac{k}{s(s-3)}}{1+\dfrac{k}{s(s+5)}}=\frac{k(s+5)}{s(s+5)(s-3)+k(s-3)}

C2(s)=1sk+5s(s+5)(s3)+k(s3)C_2(s)=\frac{1}{s}\cdot\frac{k+5}{s(s+5)(s-3)+k(s-3)}

系统的闭环特征方程为:D(s)=s(s+5)(s3)+k(s3)=s3+2s2+(k15)s3kD(s)=s(s+5)(s-3)+k(s-3)=s^3+2s^2+(k-15)s-3k

可列劳斯表

s31k15s223ks15k302s03k\begin{array}{ccc} s^3 & 1 & k-15 \\ s^2 & 2 & -3k \\ s^1 & \dfrac{5k-30}{2} & \\ s^0 & -3k & \end{array}

系统要有稳态输出,则系统必须稳定; 经过劳斯表可以得出,不存在 kk 使系统稳定,即系统不稳定。

c2()=\therefore c_2(\infty)=\infty

【评注】本题是一道较为基础的题型,难度系数极低,只需按常规的方法算出传递函数,对输入进行拉氏变换,使用终值定理 lims0sc(s)=c()\lim\limits_{s\to 0}sc(s)=c(\infty) 即可,但需要在说明终值定理前,对系

统的稳定性进行判定。这类题属于送分题,同学们只需认真算数即可。

二、【分析】本题考查传递函数绘制根轨迹以及不同条件下,参数不同取值范围。

【解】(1)开环传递函数 G(s)=Kv(s2)s(s+2)D(s)=s2+(2Kv)s+2KvG(s)=\dfrac{-K_v(s-2)}{s(s+2)}\qquad D(s)=s^2+(2-K_v)s+2K_v

①两个极点 p1=0, p2=2p_1=0,\ p_2=-2 一个零点 z1=2z_1=2

②实轴上的根轨迹 [2,0], [2,+][-2,0],\ [2,+\infty]

③与虚轴的交点 s=jωs=j\omega 带代入 D(s)=0D(s)=0

④分离点,M(s)=s+2, N(s)=s(s+2), M(s)N(s)M(s)N(s)=0M(s)=s+2,\ N(s)=s(s+2),\ M'(s)N(s)-M(s)N'(s)=0

得出 s1=22+2, s2=222s_1=2\sqrt{2}+2,\ s_2=2-2\sqrt{2}

可做 00^\circ 根轨迹如右图所示

图

无超调时,ξ1\xi\geqslant 1,那么

s=4.8s=4.8 时,代入 D(s)Kv=11.7D(s)\quad K_v=11.7

s=0.8s=-0.8 时,代入 D(s)Kv=0.34D(s)\quad K_v=0.34\qquad \therefore 无超调 0<Kv0.30<K_v\leqslant 0.3

(2)G(s)=2(2s)s(s+a)D(s)=s2+(a2)s+4G(s)=\dfrac{2(2-s)}{s(s+a)}\qquad D(s)=s^2+(a-2)s+4

等效开环传递函数 G(s)=ass22s+4G'(s)=\dfrac{as}{s^2-2s+4}

①系统存在两个极点 p1,2=1±j3p_{1,2}=1\pm j\sqrt{3} 一个零点 z1=0z_1=0

②实轴上的根轨迹 (,0](-\infty,0]

③与虚轴交点 s=jωs=j\omega 代入 D(s)=0D(s)=0

④系统的分离点:M(s)=s, N(s)=s22s+4, M(s)N(s)M(s)N(s)=0M(s)=s,\ N(s)=s^2-2s+4,\ M'(s)N(s)-M(s)N'(s)=0

可得 s=2, s=2s=-2,\ s=2(舍) 可做如下图根轨迹

图

由图可知,当 ξ\xi 为 0.707 时,对应根轨迹上的点 (2,j2)(-\sqrt{2},j\sqrt{2})

s=2+j2s=-\sqrt{2}+j\sqrt{2} 代入 D(s)D(s) 令其等于零,{4j+2(a2)j=042(a2)=0a=2+22\begin{cases}-4j+\sqrt{2}(a-2)j=0\\ 4-\sqrt{2}(a-2)=0\end{cases}\Rightarrow a=2+2\sqrt{2}

【评注】1、本题考查的根轨迹:(1)问用常规方法;(2)问用闭环特征多项式得到等效开环传递函数来做根轨迹。同学们应注意 00^\circ 根轨迹和 180180^\circ 根轨迹的区别、具体请看胡寿松第六版 P150 和 P157 注意区别记忆。 2、对于根轨迹为圆的情况,交大教材也作了具体说明。详情请看交大教材 P107 图 4-13 以及图 4-14。

三、【分析】本题考查系统开环传递函数的求解以及奈氏判据应用于伯德图问题。 【解】(1)由系统的伯德图可知,系统属于非最小相位系统

260lgω1lg5=40260lg0.5lgω2=20\frac{26-0}{\lg\omega_1-\lg\sqrt{5}}=-40\qquad \frac{26-0}{\lg 0.5-\lg\omega_2}=-20

(系统初始频段的延长线与实轴的交点 ωc=kv\omega_c=k_v

ω1=0.5ω2=10=kv\omega_1=0.5\qquad\qquad \omega_2=10=k_v

G(s)=10(10.1s)s(2s+1)G(s)=\dfrac{10(1-0.1s)}{s(2s+1)}\qquad 验证 φ(0.5)=138φ(10)=222\varphi(0.5)=-138^\circ\qquad \varphi(10)=-222^\circ

(2) ωg=5 rads1φ(ωg)=180\omega_g=\sqrt{5}\ rad\cdot s^{-1}\qquad \varphi(\omega_g)=-180^\circ

L(ωc)>0L(\omega_c)>0 时,N+=N=0N=N+NN_+=N_-=0\qquad N=N_+-N_-

注意相频曲线应补线,初始向上补到 00^\circz=P2NP=0z=P-2N\qquad P=0

ωc=ωg=5\omega_c=\omega_g=\sqrt{5},系统临界稳定。

【分析】1、在之前的题目中,我们已经提到过,当伯德图同时给出幅频曲线和相频曲线,系统一般为非最小相位系统。当仅给出幅频曲线,那么系统很大可能为最小相位系统。一旦给出相频曲线,我们应当考虑到非最小相位系统或具有延迟环节。 2、将奈氏判据应用于伯德图中,我们 1999 年上海交通大学研究生入学考试真题第三题中已有详细解释。请同学们结合两道题,同时注意补线。

四、【分析】本题考查给定条件下,确定校正 环节和各个参数,用到了长除法。

【解】由 σp=eπξ1ξ2=16%ξ=0.5\sigma_p=e^{-\frac{\pi\xi}{\sqrt{1-\xi^2}}}=16\%\to \xi=0.5

ts=3ξωn=3ωn=2t_s=\dfrac{3}{\xi\omega_n}=3\to \omega_n=2

所以系统的主导极点位置为 s1,2=ξωn±jωn1ξ2=1±j3s_{1,2}=-\xi\omega_n\pm j\omega_n\sqrt{1-\xi^2}=-1\pm j\sqrt{3}

通过系统方块图,闭环系统特征方程为 D(s)=0.1s3+s2+KDs+Kp=0D(s)=0.1s^3+s^2+K_Ds+K_p=0

其中 Φ(s)=s2+2s+4\Phi(s)=s^2+2s+4

Φ(s)\Phi(s) 可以整除 D(s)D(s),即使用长除法,得 D(s)=Φ(s)(0.1s+0.8)D(s)=\Phi(s)(0.1s+0.8)

且有 {KD0.4=1.6Kp3.2=0{KD=2Kp=3.2\begin{cases}K_D-0.4=1.6\\ K_p-3.2=0\end{cases}\Rightarrow \begin{cases}K_D=2\\ K_p=3.2\end{cases}

五、【分析】本题考查非线性系统的负倒数描述函数和线性部分自振问题。

【解】1N(M)=πM4K111(DM)2-\dfrac{1}{N(M)}=-\dfrac{\pi M}{4K_1}\dfrac{1}{\sqrt{1-\left(\dfrac{D}{M}\right)^2}}

绘制图可得:

图

M0M\to 0 时,1N(M)-\dfrac{1}{N(M)}\to -\infty

MM\to \infty 时,1N(M)-\dfrac{1}{N(M)}\to -\infty

M=2DM=\sqrt{2}D 时,1N(M)M=2D=πD2K1-\left.\dfrac{1}{N(M)}\right|_{M=\sqrt{2}D}=-\dfrac{\pi D}{2K_1}

线性部分 G(s)=2s(s+1)(0.5s+1)G(s)=\dfrac{2}{s(s+1)(0.5s+1)}

ωx=1T1T2=2\omega_x=\dfrac{1}{\sqrt{T_1T_2}}=\sqrt{2}

与实轴交点 kT1T2T1+T2=23-k\dfrac{T_1T_2}{T_1+T_2}=-\dfrac{2}{3}\qquad 渐近线:k(T1+T2)=3-k(T_1+T_2)=-3

由图可知 πD2K1<23-\dfrac{\pi D}{2K_1}<-\dfrac{2}{3} 时系统稳定,也即当 D>4K13πD>\dfrac{4K_1}{3\pi} 时,系统稳定。

(2)若 πD2K1>23-\dfrac{\pi D}{2K_1}>-\dfrac{2}{3} 时,系统存在自振,此时自振角频率为 ωx=1T1T2=2\omega_x=\dfrac{1}{\sqrt{T_1T_2}}=\sqrt{2}

1N(M)=23-\dfrac{1}{N(M)}=-\dfrac{2}{3} 得振幅 M=32π[1+14π2D236K12]12M=\dfrac{3\sqrt{2}}{\pi}\left[1+\sqrt{1-\dfrac{4\pi^2 D^2}{36K_1^2}}\right]^\frac{1}{2}

六、【分析】本题考查单位反馈下,稳态系统的输出函数式。

【解】ess()lims0sE(s)=lims0sR(s)[1GB(s)]e_{ss}(\infty)\lim\limits_{s\to 0}sE(s)=\lim\limits_{s\to 0}s\cdot R(s)[1-G_B(s)],其中 GB(s)G_B(s) 为闭环传递函数

=lims0s2s3s3+ω2s2s3+ω2s2+k1s+k1ω1=\lim_{s\to 0}s\cdot\frac{2}{s^3}\cdot\frac{s^3+\omega_2 s^2}{s^3+\omega_2 s^2+k_1s+k_1\omega_1}

=2ω2k1ω1=\frac{2\omega_2}{k_1\omega_1}

因此 c(t)=r(t)ess=t22ω2k1ω1c(t)=r(t)-e_{ss}=t^2-\dfrac{2\omega_2}{k_1\omega_1}

【评注】1、系统稳定后,系统的输出 c(t)=r(t)essc(t)=r(t)-e_{ss}

2、同学们应当学会使用终值定理 ess=lims0sE(s)e_{ss}=\lim\limits_{s\to 0}s\cdot E(s),同时牢记系统单位负反馈下,

E(s)=R(s)[1GB(s)]E(s)=R(s)[1-G_B(s)],这个公式的使用是具有前提的。

七、【分析】本题考查现代控制理论中可观测性、可控性的判别。

【解】(1)①可控性 rankPc=rank(b,Ab)=rank(1a1+11a2)=2=n\text{rank}P_c=\text{rank}(b,Ab)=\text{rank}\begin{pmatrix}1 & a_1+1\\ 1 & a_2\end{pmatrix}=2=n\quada2a110a_2-a_1-1\ne 0

②可观性 rankPo=rank(CCA)=rank(10a11)=2=n\text{rank}P_o=\text{rank}\begin{pmatrix}C\\ CA\end{pmatrix}=\text{rank}\begin{pmatrix}1 & 0\\ a_1 & 1\end{pmatrix}=2=n 无论如何系统必定可观。

(2)根据约当规范型判据,此线性定常系统完全可控的充分必要条件

B^σ1=(b^r11b^r12)\hat{B}_{\sigma 1}=\begin{pmatrix}\hat{b}_{r11}\\ \hat{b}_{r12}\end{pmatrix}B^σ2=b^r21\hat{B}_{\sigma 2}=\hat{b}_{r21}B^σ3=b^r31\hat{B}_{\sigma 3}=\hat{b}_{r31} 均为线性无关的,维数至少为 3 维

八、【分析】本题考查可控矩阵化为可控标准型以及状态反馈下的极点配置。

【解】(1)可控矩阵 rankPc=rank(b,Ab)=rank(1001)=2=n\text{rank}P_c=\text{rank}(b,Ab)=\text{rank}\begin{pmatrix}1 & 0\\ 0 & 1\end{pmatrix}=2=n 系统可控

特征多项式 sIA=s11s+1=s2+s+1α2=1α1=1|sI-A|=\begin{vmatrix}s & 1\\ -1 & s+1\end{vmatrix}=s^2+s+1\qquad \alpha_2=1\qquad \alpha_1=1

Aˉ=[0111]Bˉ=[01]Cˉ=C[Ab,b][10α11]C存在时应该这样求变换矩\bar{A}=\begin{bmatrix}0 & 1\\ -1 & -1\end{bmatrix}\qquad \bar{B}=\begin{bmatrix}0\\ 1\end{bmatrix}\qquad \bar{C}=C[Ab,b]\begin{bmatrix}1 & 0\\ \alpha_1 & 1\end{bmatrix}\Rightarrow \text{当} C \text{存在时应该这样求变换矩}

p1=[Ab,b][10α11]=[0110][1011]=[1110]p^{-1}=[Ab,b]\begin{bmatrix}1 & 0\\ \alpha_1 & 1\end{bmatrix}=\begin{bmatrix}0 & 1\\ 1 & 0\end{bmatrix}\begin{bmatrix}1 & 0\\ 1 & 1\end{bmatrix}=\begin{bmatrix}1 & 1\\ 1 & 0\end{bmatrix}

(2)可控性判定 rankPc=rank[b,Ab]=rank(10200101)=2=n\text{rank}P_c=\text{rank}[b,Ab]=\text{rank}\begin{pmatrix}1 & 0 & 2 & 0\\ 0 & 1 & 0 & 1\end{pmatrix}=2=n 系统可控。

做到该步后不必向下做。由于上海交通大学电气工程学院官网的考纲公示中已写明在现代控制部分中只考查单变量系统的极点配置和观测器设计,本部分为超纲内容,感兴趣的同学可以查阅施颂椒编写的现代控制理论中 P211 多输入系统的极点配算法。

【标注】1、可观性、可控性的内容已经反复提到,这里不再赘述。

2、化为可控标准型、可观标准型是我们同学必须掌握的重点内容,其中最关键的是构造状态转移矩阵。请同学们参考施颂椒编写的现代控制理论中的可控标准型和可观标准型的化法部分,结合具体例题 P151 例 4-16 和 P153 例 4-17。化标准型的维多一般不超过 3 维,同学们在记忆时只需记下 2 维,3 维的变换公式。

Discussion

Comments

Post

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