109 年 國立臺灣大學工程科學及海洋工程學系碩士班丙組《控制系統(B)》

📄 試題原卷 免費註冊後即可對照原始考卷 PDF免費註冊

第 1 題

A car active suspension control system with a pneumatic actuator in parallel with the passive suspension system is shown in Fig.1, where acceleration is measured by an accelerometer. The feedback signals of acceleration y¨m(t)\ddot{y}_m(t) and velocity y˙m(t)\dot{y}_m(t) are described in Fig.1. The transfer function of acceleration y¨m(s)\ddot{y}_m(s) to the ground displacement yg(s)y_g(s) can be obtained as
Ym(s)Yg(s)=s2(Ds+K)(Cg+M)s2+(Cr+D)s+K\frac{Y_m(s)}{Y_g(s)} = \frac{s^2(Ds+K)}{(C_g+M)s^2 + (C_r+D)s+K}
where M=1M=1 and D=K=Cr=2D=K=C_r=2 are given.

(a) Please sketch the root locus according to the variation of CgC_g and Cg≥0C_g \ge 0. (10%)
(b) Solve CgC_g to yield a damping ratio of 0.69 for closed-loop poles. (10%)

🖼️【此處有附圖,請對照原卷】(Fig.1)

🖼️ 本題含圖表,以下為原卷對應頁面:
原卷第 1 頁

登入後即可作答並保存紀錄。

這一題的完整詳解

核心觀念

本題考查二階系統的極點位置、根軌跡,以及阻尼比定義。

依原卷圖示,可變參數為加速度回授增益 CaC_a,且

M=1,D=K=Cv=2M=1,\qquad D=K=C_v=2

因此特徵方程式為分母:

(Ca+M)s2+(Cv+D)s+K=0(C_a+M)s^2+(C_v+D)s+K=0

代入已知數值得

(Ca+1)s2+4s+2=0(C_a+1)s^2+4s+2=0

令

A=Ca+1,A≥1A=C_a+1,\qquad A\geq 1

則特徵方程式為

As2+4s+2=0As^2+4s+2=0

(a)CaC_a 變化時的根軌跡

1. 起始極點

當 Ca=0C_a=0 時,A=1A=1,特徵方程式為

s2+4s+2=0s^2+4s+2=0

故兩個閉迴路極點為

s=−4±16−82=−2±2s=\frac{-4\pm\sqrt{16-8}}{2} =-2\pm\sqrt{2}

因此起始位置為

s1=−2+2≈−0.586s_1=-2+\sqrt{2}\approx -0.586 s2=−2−2≈−3.414s_2=-2-\sqrt{2}\approx -3.414

2. 極點進入重根的位置

特徵方程式的判別式為

Δ=42−4(A)(2)=16−8A\Delta=4^2-4(A)(2)=16-8A

當 Δ=0\Delta=0 時形成重根:

16−8A=0⇒A=216-8A=0 \quad\Rightarrow\quad A=2

因此

Ca=A−1=1C_a=A-1=1

此時重根為

s=−42A=−1s=-\frac{4}{2A}=-1

所以兩個實極點分別由 −0.586-0.586 與 −3.414-3.414 向 −1-1 移動,並在 s=−1s=-1 匯合。


3. 複數極點軌跡

令

s=σ+jωs=\sigma+j\omega

由特徵方程式

As2+4s+2=0As^2+4s+2=0

取虛部:

2Aσω+4ω=02A\sigma\omega+4\omega=0

對 ω≠0\omega\neq 0,可得

Aσ=−2⇒A=−2σA\sigma=-2 \quad\Rightarrow\quad A=-\frac{2}{\sigma}

再代入實部:

A(σ2−ω2)+4σ+2=0A(\sigma^2-\omega^2)+4\sigma+2=0

整理得

ω2=−σ(1+σ)\omega^2=-\sigma(1+\sigma)

因此複數根的軌跡為

σ2+σ+ω2=0\sigma^2+\sigma+\omega^2=0

配方後:

(σ+12)2+ω2=(12)2\left(\sigma+\frac12\right)^2+\omega^2 =\left(\frac12\right)^2

這是一個圓心在 (−1/2,0)(-1/2,0)、半徑為 1/21/2 的圓。由於系統係數為正,實部始終為負,因此只取左半平面內的上下兩段圓弧:

  • 從 s=−1s=-1 出發;
  • 上、下對稱地沿圓弧移動;
  • 當 Ca→∞C_a\to\infty 時,極點趨近原點 s=0s=0;
  • 不會穿越虛軸。

根軌跡示意如下:

🔒

後續完整解題步驟與【答案】

免費註冊,享三天全站完整詳解閱覽。

免費註冊

第 2 題

A flight control system is shown in Fig.2 with controller D(s) and controlled plant G(s).
(a) Plot the root locus if D(s)=KD(s) = K, including asymptotes, intersection points with imaginary axis and the corresponding value of K. (8%)
(b) If the controller is given as D(s)=Ks+2s+aD(s) = K \frac{s+2}{s+a}, please find the value of aa to make the root locus pass through the desired point (−1+j4)(-1+j4). (8%)
(c) Please find the value of K at the desired point (−1+j4)(-1+j4) in (b). (4%)

🖼️【此處有附圖,請對照原卷】(Fig.2)

🖼️ 本題含圖表,以下為原卷對應頁面:
原卷第 1 頁

登入後即可作答並保存紀錄。

這一題的完整詳解

核心觀念

此題考查單位負回授系統的根軌跡。閉迴路特徵方程式由

1+D(s)G(s)=01+D(s)G(s)=0

決定。

但原卷 Fig.2 僅標示方塊 G(s)G(s),未提供 G(s)G(s) 的實際傳遞函數。因此無法得到唯一的數值根軌跡、虛軸交點、KK 與 aa。以下以

G(s)=B(s)A(s)G(s)=\frac{B(s)}{A(s)}

表示一般形式,給出完整定案公式。


(a) D(s)=KD(s)=K 的根軌跡

此時開迴路傳遞函數為

L(s)=KB(s)A(s)L(s)=K\frac{B(s)}{A(s)}

閉迴路特徵方程式為

A(s)+KB(s)=0A(s)+KB(s)=0

根軌跡的基本規則如下:

  1. 根軌跡由開迴路極點出發,終止於開迴路零點或無限遠。
  2. 實軸上的根軌跡區段,其右側開迴路極點與零點總數必須為奇數。
  3. 若極點數為 nn、零點數為 mm,則無限遠漸近線數為
n−mn-m
  1. 漸近線交點為
σa=∑pi−∑zin−m\sigma_a= \frac{\displaystyle\sum p_i-\displaystyle\sum z_i}{n-m}
  1. 漸近線角度為
θk=(2k+1)180∘n−m,k=0,1,…,n−m−1\theta_k=\frac{(2k+1)180^\circ}{n-m}, \qquad k=0,1,\ldots,n-m-1

虛軸交點與對應增益

令

s=jωs=j\omega

代入特徵方程式:

A(jω)+KB(jω)=0A(j\omega)+KB(j\omega)=0

因此

K=−A(jω)B(jω)K=-\frac{A(j\omega)}{B(j\omega)}

必須同時滿足:

Im⁡[−A(jω)B(jω)]=0\operatorname{Im}\left[-\frac{A(j\omega)}{B(j\omega)}\right]=0

且

K>0K>0

所得的 ω\omega 即為虛軸交點 s=±jωs=\pm j\omega,所得實數即為相對應的 KK。

由於題目未給出 A(s)A(s) 與 B(s)B(s),無法繪製唯一的根軌跡,也無法列出數值漸近線與虛軸交點。


(b) 求使根軌跡通過 s0=−1+j4s_0=-1+j4 的 aa

控制器為

D(s)=Ks+2s+aD(s)=K\frac{s+2}{s+a}

閉迴路特徵方程式為

1+Ks+2s+aG(s)=01+K\frac{s+2}{s+a}G(s)=0

在指定點 s0=−1+j4s_0=-1+j4,令

Q=G(s0)(s0+2)Q=G(s_0)(s_0+2)

其中

s0+2=1+j4s_0+2=1+j4

設

Q=qr+jqiQ=q_r+jq_i

以及

🔒

後續完整解題步驟與【答案】

免費註冊,享三天全站完整詳解閱覽。

免費註冊

第 3 題

A hydraulic control system controlled by the state feedback with integral control is shown in Fig.3, where
A=[−313−5]A = \begin{bmatrix} -3 & 1 \\ 3 & -5 \end{bmatrix}, B=[10]B = \begin{bmatrix} 1 \\ 0 \end{bmatrix}, K=[k1k2]K = \begin{bmatrix} k_1 & k_2 \end{bmatrix}, and k1k_1 is the parameter of integral control.

(a) Please find the open loop transfer function. (7%)
(b) Solve KK and k1k_1 for the closed-loop poles at -2, -3 and -4. (13%)

🖼️【此處有附圖,請對照原卷】(Fig.3)

🖼️ 本題含圖表,以下為原卷對應頁面:
原卷第 2 頁

登入後即可作答並保存紀錄。

這一題的完整詳解

核心觀念

本題考查:

  • 狀態空間模型轉換為開迴路轉移函數:
    G(s)=C(sI−A)−1BG(s)=C(sI-A)^{-1}B
  • 積分控制器:
    kis\frac{k_i}{s}
  • 狀態回授控制律:
    u=kiz−Kx,K=[k1k2]u=k_i z-Kx,\qquad K=\begin{bmatrix}k_1&k_2\end{bmatrix}
  • 以閉迴路特徵多項式係數比對法求回授增益。

其中積分器狀態滿足:

z˙=r−y=r−Cx\dot z=r-y=r-Cx


(a)開迴路轉移函數

由

A=[−313−5],B=[10],C=[01]A=\begin{bmatrix}-3&1\\3&-5\end{bmatrix},\quad B=\begin{bmatrix}1\\0\end{bmatrix},\quad C=\begin{bmatrix}0&1\end{bmatrix}

可得

sI−A=[s+3−1−3s+5]sI-A= \begin{bmatrix} s+3&-1\\ -3&s+5 \end{bmatrix}

其行列式為

det⁡(sI−A)=(s+3)(s+5)−3=s2+8s+12\det(sI-A)=(s+3)(s+5)-3 =s^2+8s+12

因此

(sI−A)−1B=1s2+8s+12[s+53](sI-A)^{-1}B =\frac{1}{s^2+8s+12} \begin{bmatrix} s+5\\ 3 \end{bmatrix}

所以液壓系統的開迴路轉移函數為

G(s)=C(sI−A)−1B=3s2+8s+12G(s)=C(sI-A)^{-1}B =\frac{3}{s^2+8s+12}

圖中另有積分控制器 kis\frac{k_i}{s},因此包含積分器的前向迴路傳遞函數為

L(s)=kisG(s)=3kis(s2+8s+12)L(s)=\frac{k_i}{s}G(s) =\frac{3k_i}{s(s^2+8s+12)}

(b)求 KK 與積分增益 kik_i

控制律為

u=kiz−Kxu=k_i z-Kx

其中

K=[k1k2]K=\begin{bmatrix}k_1&k_2\end{bmatrix}

增廣狀態向量取為

xa=[xz]x_a=\begin{bmatrix}x\\z\end{bmatrix}

在計算閉迴路極點時令 r=0r=0,則

x˙=Ax+B(kiz−Kx)=(A−BK)x+Bkiz\dot x=Ax+B(k_i z-Kx) =(A-BK)x+B k_i z

且

z˙=−Cx\dot z=-Cx

因此增廣閉迴路系統矩陣為

🔒

後續完整解題步驟與【答案】

免費註冊,享三天全站完整詳解閱覽。

免費註冊

第 4 題

A frequency domain controller design problem is shown in Fig.4.
(a) Please find the controller KK and the gain crossover frequency to make the system have a phase margin of 50° (13%)
(b) Please find the gain margin in (a). (7%)

🖼️【此處有附圖,請對照原卷】(Fig.4)

🖼️ 本題含圖表,以下為原卷對應頁面:
原卷第 2 頁

登入後即可作答並保存紀錄。

這一題的完整詳解

核心觀念

開迴路傳遞函數由 Fig. 4 可讀為

L(s)=Ks(s2+s+4)L(s)=\frac{K}{s(s^2+s+4)}

單位負回授系統的定義如下:

  • 增益交越頻率 ωgc\omega_{gc}:∣L(jωgc)∣=1|L(j\omega_{gc})|=1
  • 相位裕度:
PM=180∘+∠L(jωgc)PM=180^\circ+\angle L(j\omega_{gc})
  • 增益裕度:
GM=1∣L(jωpc)∣GM=\frac{1}{|L(j\omega_{pc})|}

其中 ωpc\omega_{pc} 為相位交越頻率,即 ∠L(jωpc)=−180∘\angle L(j\omega_{pc})=-180^\circ。


(a) 求控制器增益 KK 與增益交越頻率

令 s=jωs=j\omega,則

L(jω)=Kjω[(4−ω2)+jω]L(j\omega)=\frac{K}{j\omega\left[(4-\omega^2)+j\omega\right]}

當 ω<2\omega<2 時,

∠L(jω)=−90∘−tan⁡−1(ω4−ω2)\angle L(j\omega) =-90^\circ-\tan^{-1}\left(\frac{\omega}{4-\omega^2}\right)

題目要求相位裕度為 50∘50^\circ,因此增益交越頻率處的相位必須為

∠L(jωgc)=−180∘+50∘=−130∘\angle L(j\omega_{gc}) =-180^\circ+50^\circ=-130^\circ

所以

−90∘−tan⁡−1(ωgc4−ωgc2)=−130∘-90^\circ-\tan^{-1}\left(\frac{\omega_{gc}}{4-\omega_{gc}^2}\right) =-130^\circ

得到

tan⁡−1(ωgc4−ωgc2)=40∘\tan^{-1}\left(\frac{\omega_{gc}}{4-\omega_{gc}^2}\right)=40^\circ

令 t=tan⁡40∘t=\tan40^\circ,則

ωgc4−ωgc2=t\frac{\omega_{gc}}{4-\omega_{gc}^2}=t

整理為

tωgc2+ωgc−4t=0t\omega_{gc}^2+\omega_{gc}-4t=0

取正根:

ωgc=−1+1+16t22t\omega_{gc} =\frac{-1+\sqrt{1+16t^2}}{2t}

代入 t=tan⁡40∘t=\tan40^\circ:

🔒

後續完整解題步驟與【答案】

免費註冊,享三天全站完整詳解閱覽。

免費註冊

第 5 題

Consider a unity negative feedback control system, as shown in Fig.5, which open loop transfer function is
G(s)=K(s+q)s(s+p)G(s) = \frac{K(s+q)}{s(s+p)}
In order to make the closed-loop control system achieve the specifications of the closed-loop poles at (−1±j1)(-1 \pm j1) and the steady state error of 0.1 for ramp input, please solve the control parameters KK, pp, and qq. (20%)

R(s) + E(s) C(s)
G(s)
Fig.5

🖼️ 本題含圖表,以下為原卷對應頁面:
原卷第 2 頁

登入後即可作答並保存紀錄。

這一題的完整詳解

本題為一標準的根軌跡與穩態誤差設計問題。首先,利用閉迴路極點的位置來確定系統的特徵方程式,進而求解參數 pp。接著,利用穩態誤差的要求,求解參數 qq。最後,綜合以上條件,求解控制器增益 KK。

系統為單位負回授系統,開迴路傳遞函數為 G(s)=K(s+q)s(s+p)G(s) = \frac{K(s+q)}{s(s+p)}。
閉迴路傳遞函數為 T(s)=G(s)1+G(s)=K(s+q)s(s+p)+K(s+q)T(s) = \frac{G(s)}{1+G(s)} = \frac{K(s+q)}{s(s+p) + K(s+q)}。
閉迴路系統的特徵方程式為 s(s+p)+K(s+q)=0s(s+p) + K(s+q) = 0。
s2+ps+Ks+Kq=0s^2 + ps + Ks + Kq = 0。
s2+(p+K)s+Kq=0s^2 + (p+K)s + Kq = 0。

閉迴路極點要求在 −1±j1-1 \pm j1。
這意味著特徵方程式的根為 s1=−1+j1s_1 = -1+j1 和 s2=−1−j1s_2 = -1-j1。
由極點可以寫出特徵方程式:(s−(−1+j1))(s−(−1−j1))=0(s - (-1+j1))(s - (-1-j1)) = 0。
(s+1−j1)(s+1+j1)=0(s+1-j1)(s+1+j1) = 0。
((s+1)−j1)((s+1)+j1)=0((s+1)-j1)((s+1)+j1) = 0。
(s+1)2−(j1)2=0(s+1)^2 - (j1)^2 = 0。
s2+2s+1−(−1)=0s^2 + 2s + 1 - (-1) = 0。
s2+2s+2=0s^2 + 2s + 2 = 0。

將此方程式與 s2+(p+K)s+Kq=0s^2 + (p+K)s + Kq = 0 比較係數:

  1. ss 項:p+K=2p+K = 2。
  2. 常數項:Kq=2Kq = 2。

穩態誤差要求為 0.1 對於斜坡輸入 (ramp input)。
斜坡輸入的拉氏轉換為 R(s)=1/s2R(s) = 1/s^2。
穩態誤差 esse_{ss} 對於斜坡輸入的公式為 ess=1Kve_{ss} = \frac{1}{K_v},其中 KvK_v 是速度常數。
Kv=lim⁡s→0sG(s)K_v = \lim_{s \to 0} s G(s)。
sG(s)=sK(s+q)s(s+p)=K(s+q)s+ps G(s) = s \frac{K(s+q)}{s(s+p)} = \frac{K(s+q)}{s+p}。

🔒

後續完整解題步驟與【答案】

免費註冊,享三天全站完整詳解閱覽。

免費註冊

其他考古題