113 年 國立臺灣聯合大學系統(清華、政治、陽明交通、中央四校聯招)研究所電機類《控制系統》

📄 試題原卷 免費註冊後即可對照原始考卷 PDF免費註冊
📄 以下 4 題共用同一段題幹
  1. (24%) Let G(s)=1P(s)G(s) = \frac{1}{P(s)} be the transfer function of a linear time-invariant system, where
P(s)=s4+(1−k1−k2)s3+2k1s2+(2k1−1)(1−k1−k2)s+2k1−1P(s) = s^4 + (1-k_1-k_2)s^3 + 2k_1s^2 + (2k_1-1)(1-k_1-k_2)s + 2k_1 - 1

and k1,k2>0k_1, k_2 > 0. Find all conditions on k1k_1 and k2k_2 for each of following cases. If no k1k_1 and k2k_2 can meet the case, explain the reason.

第 1-(a) 題6 分

(a) (6%) G(s)G(s) is stable.

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

這一題的完整詳解

核心觀念

  • 系統穩定性:對於線性時不變 (LTI) 系統,若所有特徵根(即傳遞函數分母的根)皆位於左半平面(實部 < 0),則系統 漸近穩定(在控制學教材裡常簡稱為 stable)。
  • Hurwitz 判準:對於實係數多項式
    P(s)=s4+a3s3+a2s2+a1s+a0P(s)=s^{4}+a_{3}s^{3}+a_{2}s^{2}+a_{1}s+a_{0}
    只要構成的 Routh 表 第一列的所有元素皆為正,就保證所有根位於左半平面。
  • 零列或零元素:若 Routh 表某一列的首項為 0,必須以 輔助多項式(該列上一列的多項式)取導數的係數取代,以免出現「除以 0」的情況。若最終仍有虛軸根,系統僅為 邊界穩定(不算漸近穩定)。

解題方法

  1. 寫出多項式的係數
a3=1−k1−k2,a2=2k1,a1=(2k1−1)(1−k1−k2),a0=2k1−1.\begin{aligned} a_{3}&=1-k_{1}-k_{2},\\ a_{2}&=2k_{1},\\ a_{1}&=(2k_{1}-1)(1-k_{1}-k_{2}),\\ a_{0}&=2k_{1}-1. \end{aligned}
  1. 建立 Routh 表(四階)
s第一列第二列第三列
s⁴1a₂a₀
s³a₃a₁0
s²b₁b₂0
s¹c₁00
s⁰d₁ (=a₀)

其中

b1=a3a2−a1a3,b2=a0,b_{1}=\frac{a_{3}a_{2}-a_{1}}{a_{3}},\qquad b_{2}=a_{0}, c1=b1a1−a3b2b1,d1=a0.c_{1}=\frac{b_{1}a_{1}-a_{3}b_{2}}{b_{1}},\qquad d_{1}=a_{0}.
  1. 逐項化簡

    • a0>0a_{0}>0(最後一列首項正)

2k1−1>0  ⟹  k1>12.2k_{1}-1>0\;\Longrightarrow\;k_{1}> \tfrac12 .

  • a3>0a_{3}>0(第二列首項正)

1−k1−k2>0  ⟹  k1+k2<1.1-k_{1}-k_{2}>0\;\Longrightarrow\;k_{1}+k_{2}<1 .

  • b1b_{1}
b1=(1−k1−k2)  2k1−(2k1−1)(1−k1−k2)1−k1−k2=(1−k1−k2)[2k1−(2k1−1)]1−k1−k2=1−k1−k21−k1−k2=1.\begin{aligned} b_{1}&=\frac{(1-k_{1}-k_{2})\;2k_{1}-(2k_{1}-1)(1-k_{1}-k_{2})}{1-k_{1}-k_{2}}\\ &=\frac{(1-k_{1}-k_{2})\bigl[2k_{1}-(2k_{1}-1)\bigr]}{1-k_{1}-k_{2}} =\frac{1-k_{1}-k_{2}}{1-k_{1}-k_{2}}=1 . \end{aligned}
 因而 $b_{1}=1>0$ 自動滿足。
  • c1c_{1}
🔒

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

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

免費註冊

第 1-(b) 題6 分

(b) (6%) G(s)G(s) has exactly one pair of pure imaginary poles.

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

這一題的完整詳解

核心觀念

  • 系統的極點即特徵多項式 P(s)=0P(s)=0 的根。
  • 若有純虛數極點,則必須滿足 s=jω  (ω≠0)s=j\omega\;(\omega\neq0) 使得 P(jω)=0P(j\omega)=0。
  • 透過把 s=jωs=j\omega 代入,分離實部與虛部,可得到 ω\omega 與參數之關係。
  • 多項式可因已知的純虛根而因式分解,剩餘二次因子若再產生純虛根則會得到第二對極點。

解題方法

  1. 寫出係數
a=1−k1−k2,b=2k1,c=(2k1−1)(1−k1−k2)=(2k1−1)a,d=2k1−1.\begin{aligned} a &=1-k_{1}-k_{2},\\ b &=2k_{1},\\ c &=(2k_{1}-1)(1-k_{1}-k_{2})=(2k_{1}-1)a,\\ d &=2k_{1}-1 . \end{aligned}
  1. 代入 s=jωs=j\omega

P(jω)=ω4−bω2+d  +  j(cω−aω3)=0.P(j\omega)=\omega^{4}-b\omega^{2}+d\;+\;j\bigl(c\omega-a\omega^{3}\bigr)=0 .

必須同時滿足

{ω4−bω2+d=0,c−aω2=0.\begin{cases} \omega^{4}-b\omega^{2}+d=0,\\[2mm] c-a\omega^{2}=0 . \end{cases}
  1. 由虛部得 ω2=ca=2k1−1\displaystyle \omega^{2}= \frac{c}{a}=2k_{1}-1
    (前提是 a≠0a\neq0。)
    因 ω2>0\omega^{2}>0,得到

k1>12.\boxed{k_{1}> \tfrac12 } .

  1. 檢查實部
    把 ω2=2k1−1\omega^{2}=2k_{1}-1 代入實部方程式

(2k1−1)2−b(2k1−1)+d=0.(2k_{1}-1)^{2}-b(2k_{1}-1)+d=0 .

代入 b=2k1,  d=2k1−1b=2k_{1},\;d=2k_{1}-1 可化簡為恆等式 0=00=0。
因此只要滿足第 3 步的條件,s=±j2k1−1s=\pm j\sqrt{2k_{1}-1} 就必是 P(s)P(s) 的根。

  1. 因式分解
    由上一步可寫出
🔒

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

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

免費註冊

第 1-(c) 題6 分

(c) (6%) G(s)G(s) has four poles on the imaginary axis.

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

這一題的完整詳解

核心觀念

  1. 轉移函數的極點(Poles):
    G(s)=1P(s)G(s) = \frac{1}{P(s)} 的極點即為特徵方程式 P(s)=0P(s) = 0 的根。
  2. 虛軸上的極點(Poles on the Imaginary Axis):
    若極點落在虛軸(jωj\omega 軸)上,表示極點形式為 s=±jωs = \pm j\omega(實部為零,Re(s)=0\text{Re}(s) = 0)。
  3. 因式分解與輔助多項式(Auxiliary Polynomial):
    若四階實係數多項式的所有四個根皆在虛軸上,其方程式必可因式分解為純偶數次方多項式(無奇數次方項),即 P(s)=(s2+ω12)(s2+ω22)=s4+(ω12+ω22)s2+ω12ω22P(s) = (s^2 + \omega_1^2)(s^2 + \omega_2^2) = s^4 + (\omega_1^2 + \omega_2^2)s^2 + \omega_1^2\omega_2^2(其中 ω1,ω2≥0\omega_1, \omega_2 \ge 0),奇次方項係數必須全部為零。

解題方法

方法一:多項式直接因式分解法(最簡捷)

觀察特徵多項式 P(s)P(s) 的各項係數:

P(s)=s4+(1−k1−k2)s3+2k1s2+(2k1−1)(1−k1−k2)s+2k1−1P(s) = s^4 + (1-k_1-k_2)s^3 + 2k_1s^2 + (2k_1-1)(1-k_1-k_2)s + 2k_1 - 1

將包含 (1−k1−k2)(1-k_1-k_2) 的奇次方項與其餘項分組:

P(s)=[s4+2k1s2+(2k1−1)]+(1−k1−k2)s(s2+2k1−1)P(s) = \left[ s^4 + 2k_1s^2 + (2k_1-1) \right] + (1-k_1-k_2)s \left( s^2 + 2k_1 - 1 \right)

其中第一組可十字交乘分解為:

s4+2k1s2+(2k1−1)=(s2+1)(s2+2k1−1)s^4 + 2k_1s^2 + (2k_1-1) = (s^2 + 1)(s^2 + 2k_1 - 1)

提出公因式 (s2+2k1−1)(s^2 + 2k_1 - 1):

P(s)=(s2+2k1−1)[s2+(1−k1−k2)s+1]P(s) = (s^2 + 2k_1 - 1)\left[ s^2 + (1-k_1-k_2)s + 1 \right]

令 P(s)=0P(s) = 0,極點由以下兩個二次方程式決定:

  1. s2+2k1−1=0  ⟹  s2=1−2k1s^2 + 2k_1 - 1 = 0 \implies s^2 = 1 - 2k_1
  2. s2+(1−k1−k2)s+1=0s^2 + (1-k_1-k_2)s + 1 = 0

要使 G(s)G(s) 的四個極點全部落在虛軸上,兩個二次因式的根都必須落在虛軸上:

  1. 針對因式一:s2=1−2k1s^2 = 1 - 2k_1

    • 欲使兩根為虛軸上的共軛根 s=±j2k1−1s = \pm j\sqrt{2k_1 - 1},必須滿足: 2k1−1>0  ⟹  k1>122k_1 - 1 > 0 \implies k_1 > \frac{1}{2}
    • (註:若包含原點 s=0s = 0 亦視為虛軸,則 2k1−1≥0  ⟹  k1≥122k_1 - 1 \ge 0 \implies k_1 \ge \frac{1}{2})
  2. 針對因式二:s2+(1−k1−k2)s+1=0s^2 + (1-k_1-k_2)s + 1 = 0

    • 其根為 s=−(1−k1−k2)±(1−k1−k2)2−42s = \frac{-(1-k_1-k_2) \pm \sqrt{(1-k_1-k_2)^2 - 4}}{2}。
    • 實部為零的充要條件為一次項係數為零: 1−k1−k2=0  ⟹  k1+k2=11 - k_1 - k_2 = 0 \implies k_1 + k_2 = 1
    • 此時因式化為 s2+1=0s^2 + 1 = 0,兩根為 s=±j1s = \pm j1,確實皆在虛軸上。
🔒

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

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

免費註冊

第 1-(d) 題6 分

(d) (6%) G(s)G(s) has one pole on the left-half plane and three poles on the right half plane. (Note: the imaginary axis is neither a part of the left-half plane, nor a part of the right-half plane.)

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

這一題的完整詳解

核心觀念
本題要求判斷 傳遞函數

G(s)=1P(s),P(s)=s4+(1−k1−k2)s3+2k1s2+(2k1−1)(1−k1−k2)s+2k1−1G(s)=\frac{1}{P(s)},\qquad P(s)=s^{4}+(1-k_{1}-k_{2})s^{3}+2k_{1}s^{2}+(2k_{1}-1)(1-k_{1}-k_{2})s+2k_{1}-1

的極點分布,使得 左半平面 (LHP) 只有 1 個極點,右半平面 (RHP) 有 3 個極點,且 不可有極點落在虛軸。
判斷極點分布最直接的工具是 Routh–Hurwitz 判別法:
第一列的符號變化次數 = 系統右半平面極點的個數。


解題方法

  1. 寫出 Routh 表(四次多項式)
    多項式係數
次方係數
s4s^{4}11
s3s^{3}a3=1−k1−k2a_{3}=1-k_{1}-k_{2}
s2s^{2}a2=2k1a_{2}=2k_{1}
s1s^{1}a1=(2k1−1)(1−k1−k2)a_{1}=(2k_{1}-1)(1-k_{1}-k_{2})
s0s^{0}a0=2k1−1a_{0}=2k_{1}-1
  1. 構造 Routh 表(只保留首列)
s412k12k1−1s31−k1−k2(2k1−1)(1−k1−k2)0s2(1−k1−k2)2k1−(2k1−1)(1−k1−k2)1−k1−k22k1−10s1000s02k1−1\begin{array}{c|ccc} s^{4} & 1 & 2k_{1} & 2k_{1}-1 \\[2mm] s^{3} & 1-k_{1}-k_{2} & (2k_{1}-1)(1-k_{1}-k_{2}) & 0 \\[2mm] s^{2} & \displaystyle\frac{(1-k_{1}-k_{2})2k_{1}- (2k_{1}-1)(1-k_{1}-k_{2})}{1-k_{1}-k_{2}} & 2k_{1}-1 & 0 \\[4mm] s^{1} & 0 & 0 & 0 \\[2mm] s^{0} & 2k_{1}-1 & & \end{array}
  1. 化簡
  • 第三列第一項

b1=(1−k1−k2)2k1−(2k1−1)(1−k1−k2)1−k1−k2=1b_{1}= \frac{(1-k_{1}-k_{2})2k_{1}-(2k_{1}-1)(1-k_{1}-k_{2})}{1-k_{1}-k_{2}} =1

  • 第三列第二項

b2=(1−k1−k2)(2k1−1)1−k1−k2=2k1−1b_{2}= \frac{(1-k_{1}-k_{2})(2k_{1}-1)}{1-k_{1}-k_{2}} = 2k_{1}-1

故

s412k12k1−1s31−k1−k2(2k1−1)(1−k1−k2)0s212k1−10s1000s02k1−1\begin{array}{c|ccc} s^{4} & 1 & 2k_{1} & 2k_{1}-1 \\[2mm] s^{3} & 1-k_{1}-k_{2} & (2k_{1}-1)(1-k_{1}-k_{2}) & 0 \\[2mm] s^{2} & 1 & 2k_{1}-1 & 0 \\[2mm] s^{1} & 0 & 0 & 0 \\[2mm] s^{0} & 2k_{1}-1 & \end{array}

第四列全為 0,需利用 輔助多項式 來處理。
輔助多項式取自上一步的非零列(s2s^{2} 行):

A(s)=1  s2+(2k1−1)A(s)=1\;s^{2}+(2k_{1}-1)

對其求導得到

A′(s)=2sA'(s)=2s

把導數的係數填入 s1s^{1} 行,得到

s1 行首項=2(正號)s^{1}\text{ 行首項}=2\quad(\text{正號})

最終第一列(即判斷極點數量所需的符號)為

  1,    1−k1−k2,    1,    2,    2k1−1  \boxed{ \; 1,\; \; 1-k_{1}-k_{2},\; \; 1,\; \; 2,\; \; 2k_{1}-1 \;}

  1. 計算符號變化次數
🔒

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

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

免費註冊
📄 以下 2 題共用同一段題幹
  1. (13%) Consider a second order system with a transfer function G(s)=ωn2s2+2ξωns+ωn2G(s) = \frac{\omega_n^2}{s^2+2\xi\omega_ns+\omega_n^2}, where ξ>0\xi > 0 and ωn>0\omega_n > 0 are the damping ratio and the undamped natural frequency of G(s)G(s), respectively. Let y(t)y(t) be the unit-step response of G(s)G(s). Suppose that y1y_1 is the maximum overshoot, which takes place at time t1t_1, while 1−y21-y_2 is the first local minimum of y(t)y(t) for t>t1t > t_1 (see Figure 1). Suppose y1y2=2.6858\frac{y_1}{y_2} = 2.6858 and t1=6.5866t_1 = 6.5866 sec.

🖼️【此處有附圖,請對照原卷 Figure 1: The unit-step response of Problem 2、Figure 2: The unity feedback system】

第 2-(a) 題8 分

(a) (8%) Find the values of ξ\xi and ωn\omega_n.

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

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

這一題的完整詳解

核心觀念

此題利用欠阻尼二階系統步階響應的兩個特性:相鄰峰值偏離穩態值的比例,以及第一次最大超越時間。

欠阻尼時 0<ξ<10<\xi<1,阻尼自然頻率為

ωd=ωn1−ξ2.\omega_d=\omega_n\sqrt{1-\xi^2}.

步階響應相鄰兩個極值與穩態值 11 的偏差,大小相差一個因子 e−ξπ/1−ξ2e^{-\xi\pi/\sqrt{1-\xi^2}}。因此,圖中的第一次超越量 y1y_1 與下一個谷值低於穩態值的量 y2y_2 滿足

y1y2=eξπ1−ξ2.\frac{y_1}{y_2} = e^{\frac{\xi\pi}{\sqrt{1-\xi^2}}}.

第一次最大超越發生時間則為

t1=πωd.t_1=\frac{\pi}{\omega_d}.

解題方法

先由超越量比值求阻尼比。令 R=2.6858R=2.6858,取自然對數:

ln⁡R=ξπ1−ξ2.\ln R=\frac{\xi\pi}{\sqrt{1-\xi^2}}.

整理可得

ξ=ln⁡Rπ2+(ln⁡R)2.\xi=\frac{\ln R}{\sqrt{\pi^2+(\ln R)^2}}.

代入 R=2.6858R=2.6858:

🔒

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

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

免費註冊

第 2-(b) 題5 分

(b) (5%) Consider the unity feedback system in Figure 2, where K>0K > 0. Let rr in Figure 2 be the unit-step input. Find the range of KK such that the percent maximum overshoot of yy is less than 50% and the steady-state error is less than 0.6.

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

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

這一題的完整詳解

核心觀念

這題先由題幹給的超越量與波谷資料,找出原系統的阻尼比 ξ\xi;再分析比例增益 KK 加入單位回授後的閉迴路阻尼比與穩態誤差。

標準二階系統的最大超越量比例為

Mp=exp⁡(−ξπ1−ξ2).M_p=\exp\left(-\frac{\xi\pi}{\sqrt{1-\xi^2}}\right).

第一個波谷相對於穩態值的偏差量為 Mp2M_p^2,因此題目所給的 y1/y2y_1/y_2 是 Mp/Mp2=1/MpM_p/M_p^2=1/M_p。峰值時間則為

t1=πωn1−ξ2.t_1=\frac{\pi}{\omega_n\sqrt{1-\xi^2}}.

解題方法

由 y1/y2=2.6858y_1/y_2=2.6858,可得

Mp=12.6858.M_p=\frac{1}{2.6858}.

令 L=ln⁡(2.6858)L=\ln(2.6858),由超越量公式可求阻尼比:

L=ξπ1−ξ2⟹ξ=Lπ2+L2≈0.3.L=\frac{\xi\pi}{\sqrt{1-\xi^2}} \quad\Longrightarrow\quad \xi=\frac{L}{\sqrt{\pi^2+L^2}} \approx 0.3.

再代入 t1=6.5866t_1=6.5866 秒:

ωn=πt11−ξ2≈0.5 rad/s.\omega_n =\frac{\pi}{t_1\sqrt{1-\xi^2}} \approx 0.5\ \text{rad/s}.

比例增益 KK 串接於原系統後,閉迴路傳遞函數為

T(s)=Kωn2s2+2ξωns+(1+K)ωn2.T(s)=\frac{K\omega_n^2} {s^2+2\xi\omega_n s+(1+K)\omega_n^2}.

因此閉迴路阻尼比為

🔒

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

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

免費註冊
📄 以下 2 題共用同一段題幹
  1. (13%) Consider the unity feedback control system in Figure 2, where K>0K > 0 and
G(s)=s2−s+4.25s(s2−0.8s+13.12)G(s) = \frac{s^2-s+4.25}{s(s^2-0.8s+13.12)}

第 3-(a) 題8 分

(a) (8%) Draw the root locus of the system for K>0K > 0. Find the arrival angles of the zeros, the departure angles of the complex poles, and the intersections of the root locus with the imaginary axis, if there is any.

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

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

這一題的完整詳解

核心觀念

根軌跡由閉迴路特徵方程 1+KG(s)=01+KG(s)=0 決定。分支從開迴路極點出發,終止於開迴路零點或無限遠;離開角、到達角由角度條件求得,虛軸交點則令 s=jωs=j\omega 代入特徵方程。

解題方法

圖中的方塊圖是單位負回授,前向路徑依序為增益 KK 與 G(s)G(s),且 K>0K>0。本題開迴路傳遞函數為

KG(s)=Ks2−s+4.25s(s2−0.8s+13.12).KG(s)=K\frac{s^2-s+4.25}{s(s^2-0.8s+13.12)}.

先求開迴路極點與零點:

p1=0,p2,3=0.4±j3.6,p_1=0,\qquad p_{2,3}=0.4\pm j3.6, z1,2=0.5±j2.z_{1,2}=0.5\pm j2.

因此共有三條根軌跡分支:兩條由共軛複數極點出發,終止於共軛零點;另一條由 s=0s=0 出發,走向無限遠。

實軸上只有一個極點 s=0s=0。依實軸根軌跡奇偶規則,根軌跡位於 (−∞,0)(-\infty,0),沒有正實軸根軌跡。無限遠漸近線數為 3−2=13-2=1,角度與形心為

θasym=180∘,σasym=∑pi−∑zi3−2=0.8−11=−0.2.\theta_{\mathrm{asym}}=180^\circ,\qquad \sigma_{\mathrm{asym}}= \frac{\sum p_i-\sum z_i}{3-2} =\frac{0.8-1}{1}=-0.2.

所以從原點出發的實軸分支沿負實軸走向 −∞-\infty。

複數極點的離開角

令角度由正實軸起算,使用

θd=180∘+∑i∠(pk−zi)−∑j≠k∠(pk−pj).\theta_d =180^\circ+\sum_i\angle(p_k-z_i) -\sum_{j\ne k}\angle(p_k-p_j).

對上方極點 p+=0.4+j3.6p_+=0.4+j3.6,各角度和為

∑i∠(p+−zi)=∠(−0.1+j1.6)+∠(−0.1+j5.6)=93.576∘+91.023∘=184.599∘,∑j≠+∠(p+−pj)=∠(0.4+j3.6)+∠(j7.2)=83.660∘+90∘=173.660∘.\begin{aligned} \sum_i\angle(p_+-z_i) &=\angle(-0.1+j1.6)+\angle(-0.1+j5.6)\\ &=93.576^\circ+91.023^\circ =184.599^\circ,\\ \sum_{j\ne +}\angle(p_+-p_j) &=\angle(0.4+j3.6)+\angle(j7.2)\\ &=83.660^\circ+90^\circ =173.660^\circ. \end{aligned}

因此

θd(p+)=180∘+184.599∘−173.660∘=190.939∘.\theta_d(p_+)=180^\circ+184.599^\circ-173.660^\circ =190.939^\circ.

下方極點由共軛對稱得到

θd(p−)=169.061∘.\theta_d(p_-)=169.061^\circ.

複數零點的到達角

令零點附近根軌跡方向以 s−zks-z_k 表示,使用

θa=180∘−∑i≠k∠(zk−zi)+∑j∠(zk−pj).\theta_a =180^\circ-\sum_{i\ne k}\angle(z_k-z_i) +\sum_j\angle(z_k-p_j).

對上方零點 z+=0.5+j2z_+=0.5+j2:

∑i≠+∠(z+−zi)=90∘,∑j∠(z+−pj)=75.964∘−86.424∘+88.977∘=78.517∘.\begin{aligned} \sum_{i\ne +}\angle(z_+-z_i)&=90^\circ,\\ \sum_j\angle(z_+-p_j) &=75.964^\circ-86.424^\circ+88.977^\circ\\ &=78.517^\circ. \end{aligned}

故

θa(z+)=180∘−90∘+78.517∘=168.517∘.\theta_a(z_+)=180^\circ-90^\circ+78.517^\circ =168.517^\circ.

由共軛對稱,下方零點的到達角為

θa(z−)=191.483∘.\theta_a(z_-)=191.483^\circ.

虛軸交點

閉迴路特徵方程為

s3+(K−0.8)s2+(13.12−K)s+4.25K=0.s^3+(K-0.8)s^2+(13.12-K)s+4.25K=0.

令 s=jωs=j\omega,分別比較實部與虛部:

(0.8−K)ω2+4.25K=0,(0.8-K)\omega^2+4.25K=0,
🔒

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

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

免費註冊

第 3-(b) 題5 分

(b) (5%) Find the range of KK for the closed-loop system to be stable.

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

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

這一題的完整詳解

核心觀念

單位負回授系統的閉迴路特徵方程為 1+KG(s)=01+KG(s)=0。三次多項式

s3+a2s2+a1s+a0s^3+a_2s^2+a_1s+a_0

依 Routh–Hurwitz 判據穩定的條件為 a2>0a_2>0、a1>0a_1>0、a0>0a_0>0,且 a2a1>a0a_2a_1>a_0。

解題方法

代入題目給定的 G(s)G(s),閉迴路特徵方程為

s(s2−0.8s+13.12)+K(s2−s+4.25)=0s(s^2-0.8s+13.12)+K(s^2-s+4.25)=0

整理得

s3+(K−0.8)s2+(13.12−K)s+4.25K=0s^3+(K-0.8)s^2+(13.12-K)s+4.25K=0

由於題目已給 K>0K>0,常數項 4.25K4.25K 為正。其餘穩定條件為

🔒

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

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

免費註冊
📄 以下 6 題共用同一段題幹
  1. (38%) Consider a unity-feedback system with the open-loop transfer function
G(s)=K(s+1)s(s−1),where G(jω)=−2K1+ω2+jK(1−ω2)ω(1+ω2)G(s) = \frac{K(s+1)}{s(s-1)}, \quad \text{where } G(j\omega) = \frac{-2K}{1+\omega^2} + j\frac{K(1-\omega^2)}{\omega(1+\omega^2)}

第 4-(a) 題10 分

(a) (10%) Show that part of the root locus is a circle. Write down the circle equation by letting s=x+jys = x + jy.

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

這一題的完整詳解

核心觀念

  • 根軌跡 (Root Locus):描述閉迴路系統特徵根隨增益 KK 變化的軌跡。
  • 根軌跡特徵:
    1. 必須滿足特徵方程 1+KG(s)=01+KG(s)=0。
    2. 任意點 ss 必須同時滿足 相角條件 ∠KG(s)=(2k+1)π\angle KG(s)= (2k+1)\pi。
    3. 若 KG(s)KG(s) 為實數且為負,則該點必在根軌跡上。
  • 虛部與實部分離:將 s=x+jys=x+jy 代入 KG(s)KG(s),分離實部與虛部,可直接求得滿足實部條件的幾何形狀。

解題方法

  1. 寫出開環傳遞函數

G(s)=K(s+1)s(s−1)G(s)=\frac{K(s+1)}{s(s-1)}

根軌跡由方程 1+G(s)=01+G(s)=0 或等價的 KG(s)=−1KG(s) = -1 決定。

  1. 代入 s=x+jys=x+jy

G(s)=K((x+jy)+1)(x+jy)((x+jy)−1)G(s)=\frac{K((x+jy)+1)}{(x+jy)\big((x+jy)-1\big)}

分子:(x+1)+jy(x+1)+jy
分母:(x+jy)(x−1+jy)=(x(x−1)−y2)+j y(2x−1)(x+jy)(x-1+jy)= (x(x-1)-y^{2})+j\,y(2x-1)

  1. 求 G(s)K\displaystyle \frac{G(s)}{K} 的實部與虛部

G(s)K=(x+1)+jy(x(x−1)−y2)+j y(2x−1)\frac{G(s)}{K}= \frac{(x+1)+jy}{(x(x-1)-y^{2})+j\,y(2x-1)}

乘以共軛分母:

G(s)K=[(x+1)+jy][(x(x−1)−y2)−j y(2x−1)](x(x−1)−y2)2+y2(2x−1)2\frac{G(s)}{K}= \frac{\big[(x+1)+jy\big]\big[(x(x-1)-y^{2})-j\,y(2x-1)\big]} {(x(x-1)-y^{2})^{2}+y^{2}(2x-1)^{2}}

只取分子實部(因根軌跡要求 KG(s)=−1KG(s)=-1 為實數):

Re ⁣[G(s)K]=(x+1)[x(x−1)−y2]+y2(2x−1)(x(x−1)−y2)2+y2(2x−1)2.\text{Re}\!\left[\frac{G(s)}{K}\right]= \frac{(x+1)\big[x(x-1)-y^{2}\big] + y^{2}(2x-1)} {(x(x-1)-y^{2})^{2}+y^{2}(2x-1)^{2}} .
  1. 相角條件
    為使 KG(s)=−1KG(s)=-1,必須同時滿足
Im ⁣[G(s)K]=0,Re ⁣[G(s)K]<0.\text{Im}\!\left[\frac{G(s)}{K}\right]=0 ,\qquad \text{Re}\!\left[\frac{G(s)}{K}\right] <0 .
🔒

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

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

免費註冊

第 4-(b) 題10 分

(b) (10%) Sketch the Nyquist plot for K>0K>0 as detail as possible.

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

這一題的完整詳解

核心觀念

  1. 奈氏軌跡(Nyquist Plot)與映射原理:
    • 封閉的奈氏路徑(Nyquist Contour)沿著 ss 平面的虛軸順時針包圍整個右半平面(RHP)。
    • 針對開迴路極點位於原點 s=0s = 0 的情況,奈氏路徑必須在原點處以半徑 ϵ→0+\epsilon \to 0^+ 的小半圓(通常繞過右半平面,即由 −j0+-j0^+ 經 0+0^+ 繞至 +j0++j0^+)進行迴避(indentation)。
    • 奈氏軌跡為 ss 沿著該路徑時,開迴路轉移函數 G(s)G(s) 在複數平面上的對應映射軌跡。
  2. 複數分解與關鍵幾何特徵:
    • 頻率響應的實部與虛部: G(jω)=Re[G(jω)]+j Im[G(jω)]G(j\omega) = \text{Re}[G(j\omega)] + j\,\text{Im}[G(j\omega)]
    • 漸近線與交點:
      • 與實軸交點:令 Im[G(jω)]=0\text{Im}[G(j\omega)] = 0,解出交點頻率與對應實部數值。
      • 與虛軸交點:令 Re[G(jω)]=0\text{Re}[G(j\omega)] = 0,確認是否與虛軸相交。
      • 高頻行為(ω→∞\omega \to \infty):幅值趨近於 00,確認進入原點的角度。
      • 低頻小半圓映射(ϵ→0+\epsilon \to 0^+):分析無窮大圓弧的半徑、起點與終點相位、轉向。

解題方法

第一步:定義標準 Nyquist Contour

在 ss 平面上選取包圍整個右半平面的順時針路徑:

  1. 路徑 I(正虛軸):s=jωs = j\omega,ω\omega 由 0+0^+ 增加到 +∞+\infty。
  2. 路徑 II(無窮大半圓):s=Rejθs = R e^{j\theta},其中 R→∞R \to \infty,θ\theta 由 +90∘+90^\circ 順時針轉至 −90∘-90^\circ。
  3. 路徑 III(負虛軸):s=jωs = j\omega,ω\omega 由 −∞-\infty 增加到 0−0^-。
  4. 路徑 IV(原點小半圓):s=ϵejϕs = \epsilon e^{j\phi},其中 ϵ→0+\epsilon \to 0^+,ϕ\phi 由 −90∘-90^\circ 順時針轉至 +90∘+90^\circ(繞過右側以排除原點極點)。

第二步:各段路徑在 G(s)G(s) 平面上的映射分析

1. 路徑 I:正虛軸映射(ω:0+→+∞\omega: 0^+ \to +\infty)

題目已給出:

G(jω)=−2K1+ω2+j K(1−ω2)ω(1+ω2)G(j\omega) = \frac{-2K}{1+\omega^2} + j\,\frac{K(1-\omega^2)}{\omega(1+\omega^2)}
  • 實部與虛部表達式:

    Re[G(jω)]=−2K1+ω2<0(∀ω≥0,K>0)\text{Re}[G(j\omega)] = \frac{-2K}{1+\omega^2} < 0 \quad (\forall \omega \ge 0, K > 0) Im[G(jω)]=K(1−ω2)ω(1+ω2)\text{Im}[G(j\omega)] = \frac{K(1-\omega^2)}{\omega(1+\omega^2)}

    實部恆為負值,表示正虛軸映射之軌跡完全落在第二與第三象限,永不跨越虛軸。

  • 起點(ω→0+\omega \to 0^+):

    lim⁡ω→0+Re[G(jω)]=−2K\lim_{\omega \to 0^+} \text{Re}[G(j\omega)] = -2K lim⁡ω→0+Im[G(jω)]=+∞\lim_{\omega \to 0^+} \text{Im}[G(j\omega)] = +\infty

    軌跡發端於第二象限遠處,垂直漸近線為 Re=−2K\text{Re} = -2K。

  • 與實軸交點:
    令 Im[G(jω)]=0\text{Im}[G(j\omega)] = 0:

    1−ω2=0  ⟹  ω=1 rad/s1 - \omega^2 = 0 \implies \omega = 1\text{ rad/s}

    將 ω=1\omega = 1 代入實部:

    G(j1)=−2K1+12=−KG(j1) = \frac{-2K}{1+1^2} = -K

    軌跡在 ω=1\omega = 1 處由第二象限跨越負實軸進入第三象限,交點座標為 (−K,0)(-K, 0)。

  • 終點(ω→+∞\omega \to +\infty):

    lim⁡ω→+∞∣G(jω)∣=lim⁡ω→+∞Kω=0\lim_{\omega \to +\infty} |G(j\omega)| = \lim_{\omega \to +\infty} \frac{K}{\omega} = 0

    角度分析:

    G(s)≈Ks  ⟹  ∠G(jω)→−90∘G(s) \approx \frac{K}{s} \implies \angle G(j\omega) \to -90^\circ

    軌跡從第三象限沿著切線方向 −90∘-90^\circ(負虛軸方向)收斂至原點 (0,0)(0, 0)。


2. 路徑 II:無窮大半圓映射(R→∞R \to \infty)

令 s=Rejθs = R e^{j\theta},其中 R→∞R \to \infty,θ\theta 由 +90∘→−90∘+90^\circ \to -90^\circ:

G(s)≈Kss2=Ks=KRe−jθ→0⋅e−jθG(s) \approx \frac{K s}{s^2} = \frac{K}{s} = \frac{K}{R} e^{-j\theta} \to 0 \cdot e^{-j\theta}

在 G(s)G(s) 平面上,該段映射緊縮在原點 (0,0)(0, 0) 處。


3. 路徑 III:負虛軸映射(ω:−∞→0−\omega: -\infty \to 0^-)

根據複數對稱性,負虛軸映射與正虛軸映射關於實軸互相對稱:

  • ω→−∞\omega \to -\infty:自原點以 +90∘+90^\circ 方向出發進入第二象限。
  • ω=−1\omega = -1:穿過實軸交點 (−K,0)(-K, 0)。
  • ω→0−\omega \to 0^-:實部趨近於 −2K-2K,虛部趨近於 −∞-\infty(位於第三象限)。

🔒

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

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

免費註冊

第 4-(c) 題4 分

(c) (4%) Analyze the stability from the Nyquist plot for all KK.

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

這一題的完整詳解

核心觀念

  • Nyquist 穩定判據:對於單回授系統,若開迴路傳遞函數 G(s)G(s) 在右半平面(RHP)有 PP 個極點,則其 Nyquist 曲線須順時針圍繞點 −1+0j-1+0j NN 圈,使得

N=P−ZN = P - Z

其中 ZZ 為閉迴路特徵方程 1+G(s)=01+G(s)=0 在 RHP 的根數。系統穩定⇔Z=0Z=0。

  • 極點與零點分布:G(s)=K(s+1)s(s−1)G(s)=\dfrac{K(s+1)}{s(s-1)} 有兩個實數極點 s=0s=0(左半平面)與 s=1s=1(右半平面),因此 P=1P=1。
  • Nyquist 曲線的端點與對稱性:由題目給出的頻率響應

G(jω)=−2K1+ω2+jK(1−ω2)ω(1+ω2)G(j\omega)=\frac{-2K}{1+\omega^{2}}+j\frac{K(1-\omega^{2})}{\omega(1+\omega^{2})}

可直接觀察 ω→0+\omega\to0^{+} 與 ω→∞\omega\to\infty 時的極限,且因係數實數,負頻率會產生鏡像曲線。


解題方法

  1. 計算 Nyquist 曲線的端點
ω→0+ ⁣:  ℜ{G(jω)}→−2K,ℑ{G(jω)}=K(1−ω2)ω(1+ω2)→{+∞(K>0)−∞(K<0)\begin{aligned} \omega\to0^{+}\! &: \; \Re\{G(j\omega)\}\to -2K,\\ &\Im\{G(j\omega)\}= \frac{K(1-\omega^{2})}{\omega(1+\omega^{2})}\to \begin{cases} +\infty & (K>0)\\ -\infty & (K<0) \end{cases} \end{aligned}

因此曲線從正(或負)無窮大的虛軸點 (0,±∞)(0,\pm\infty) 出發,橫向靠近實軸的 (−2K,0)(-2K,0)。

ω→∞ ⁣:  ℜ{G(jω)}→0,ℑ{G(jω)}→0−.\omega\to\infty\! :\; \Re\{G(j\omega)\}\to 0,\qquad \Im\{G(j\omega)\}\to 0^{-} .

所以曲線最終收斂到原點,且是從 下方 逼近。

  1. 找出實軸交點

    ℑ{G(jω)}=0\Im\{G(j\omega)\}=0 時,必須滿足 1−ω2=01-\omega^{2}=0 → ω=1\omega=1。
    代入得到實部

ℜ{G(j1)}=−2K1+1=−K.\Re\{G(j1)\}= \frac{-2K}{1+1}= -K .

因此 Nyquist 曲線在 ω=1\omega=1 處穿過實軸於點 (−K,0)(-K,0)。

  1. 判斷曲線是否包圍 −1-1

    • K>0K>0 時:曲線由 (0,+∞j)(0,+\infty j) 向下移動,先經過左上象限,於 ω=1\omega=1 交於 (−K,0)(-K,0),再向原點下降。
🔒

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

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

免費註冊

第 4-(d) 題3 分

(d) (3%) When K=1/3K=1/3 and r(t)=1r(t)=1, what is the steady state output?

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

這一題的完整詳解

核心觀念

  • 單位負回授系統:閉迴路轉移函數為
    T(s)=G(s)1+G(s)T(s)=\frac{G(s)}{1+G(s)}
  • 穩態誤差 (Steady‑state error):對於單位階躍輸入 r(t)=1r(t)=1,若系統在 s=0s=0 具有至少一個極點(type‑1 系統),則 lim⁡t→∞e(t)=0\displaystyle \lim_{t\to\infty}e(t)=0,即輸出最終會等於參考   y∞=1\;y_{\infty}=1。
  • 最終值定理 (Final‑value theorem)
y∞=lim⁡t→∞y(t)=lim⁡s→0sY(s)Y(s)=T(s)R(s),  R(s)=1sy_{\infty}= \lim_{t\to\infty}y(t)=\lim_{s\to0}sY(s) \qquad Y(s)=T(s)R(s),\; R(s)=\frac1s

前提:閉迴路系統的所有極點必須位於左半平面(或僅有一個簡單的 s=0s=0 極點)。在本題中,儘管閉迴路極點會因 KK 而移動,但在 K=13K=\frac13 時仍屬於 type‑1,可直接使用極限方式得到 y∞=1\displaystyle y_{\infty}=1(考試常以此簡化)。


解題方法

  1. 寫出閉迴路轉移函數
T(s)=G(s)1+G(s)=K(s+1)s(s−1)1+K(s+1)s(s−1)=K(s+1)s(s−1)+K(s+1)T(s)=\frac{G(s)}{1+G(s)}= \frac{\displaystyle \frac{K(s+1)}{s(s-1)}} {1+\displaystyle\frac{K(s+1)}{s(s-1)}} =\frac{K(s+1)}{s(s-1)+K(s+1)}
  1. 代入 K=13K=\frac13
🔒

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

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

免費註冊

第 4-(e) 題3 分

(e) (3%) When K=3K=3 and r(t)=cos⁡(t)r(t)=\cos(t), what is the steady state output?

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

這一題的完整詳解

核心觀念

  • 單位回授系統:閉迴路傳遞函數 T(s)=G(s)1+G(s)T(s)=\dfrac{G(s)}{1+G(s)}。
  • 正弦輸入的穩態響應:若 r(t)=cos⁡(ωt)r(t)=\cos(\omega t),則 yss(t)=∣T(jω)∣cos⁡(ωt+∠T(jω))y_{\text{ss}}(t)=|T(j\omega)|\cos(\omega t+\angle T(j\omega))。
  • 頻率響應:只需要求 G(jω)G(j\omega) 在 ω=1 \omega=1\,rad/s(因為 r(t)=cos⁡tr(t)=\cos t)的值,進而求出 T(jω)T(j\omega)。
  • 系統穩定性:閉迴路特徵方程 1+G(s)=01+G(s)=0 必須有左半平面根,才能保證存在穩態解。

解題方法

  1. 求 G(jω)G(j\omega) 的實虛部
    題目已給出

G(jω)=−2K1+ω2+jK(1−ω2)ω(1+ω2).G(j\omega)=\frac{-2K}{1+\omega^{2}}+j\frac{K(1-\omega^{2})}{\omega(1+\omega^{2})}.

  1. 代入 K=3K=3、ω=1\omega=1
Re⁡{G(j1)}=−2⋅31+12=−62=−3,Im⁡{G(j1)}=3(1−12)1 (1+12)=0.\begin{aligned} \operatorname{Re}\{G(j1)\}&=\frac{-2\cdot 3}{1+1^{2}}=\frac{-6}{2}=-3,\\[2pt] \operatorname{Im}\{G(j1)\}&=\frac{3(1-1^{2})}{1\,(1+1^{2})}=0. \end{aligned}

因此

G(j1)=−3(純實數).G(j1)=-3\quad(\text{純實數}).

  1. 求閉迴路頻率響應 T(jω)T(j\omega)

T(jω)=G(jω)1+G(jω).T(j\omega)=\frac{G(j\omega)}{1+G(j\omega)}.

🔒

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

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

免費註冊

第 4-(f) 題8 分

(f) (8%) If the gain crossover frequency =3=\sqrt{3} rad/sec, what are the corresponding KK, the Gain Margin, Phase Margin and the phase crossover frequency?

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

這一題的完整詳解

核心觀念

本題評量頻域穩定度分析之四大核心指標的定義與計算:

  1. 增益交越頻率(Gain Crossover Frequency, ωgc\omega_{gc}):
    開環路弦波轉移函數之振幅大小等於 11(即 0 dB0\text{ dB})所對應的角頻率,滿足: ∣G(jωgc)∣=1|G(j\omega_{gc})| = 1
  2. 相位交越頻率(Phase Crossover Frequency, ωpc\omega_{pc}):
    開環路弦波轉移函數之相位角為 −180∘-180^\circ(負實軸交點)所對應的角頻率,滿足: ∠G(jωpc)=−180∘⟺Im[G(jωpc)]=0且Re[G(jωpc)]<0\angle G(j\omega_{pc}) = -180^\circ \quad \Longleftrightarrow \quad \text{Im}[G(j\omega_{pc})] = 0 \quad \text{且} \quad \text{Re}[G(j\omega_{pc})] < 0
  3. 增益邊限(Gain Margin, GMGM):
    在相位交越頻率 ωpc\omega_{pc} 處,開環振幅倒數所代表的裕度量: GM=1∣G(jωpc)∣(或 GMdB=20log⁡101∣G(jωpc)∣ dB)GM = \frac{1}{|G(j\omega_{pc})|} \quad \left(\text{或 } GM_{\text{dB}} = 20\log_{10}\frac{1}{|G(j\omega_{pc})|} \text{ dB}\right)
  4. 相位邊限(Phase Margin, PMPM):
    在增益交越頻率 ωgc\omega_{gc} 處,系統相位距離臨界相位 −180∘-180^\circ 之差值: PM=180∘+∠G(jωgc)PM = 180^\circ + \angle G(j\omega_{gc})

解題方法

題目已將 G(jω)G(j\omega) 拆解為實部與虛部之代數形式:

G(jω)=−2K1+ω2+jK(1−ω2)ω(1+ω2)G(j\omega) = \frac{-2K}{1+\omega^2} + j\frac{K(1-\omega^2)}{\omega(1+\omega^2)}

步驟一:求開環頻率響應之振幅式 ∣G(jω)∣|G(j\omega)|

利用實部與虛部計算振幅平方:

∣G(jω)∣2=[Re(G(jω))]2+[Im(G(jω))]2=4K2(1+ω2)2+K2(1−ω2)2ω2(1+ω2)2=K2ω2(1+ω2)2[4ω2+(1−2ω2+ω4)]=K2(1+2ω2+ω4)ω2(1+ω2)2=K2(1+ω2)2ω2(1+ω2)2=K2ω2\begin{aligned} |G(j\omega)|^2 &= \left[\text{Re}(G(j\omega))\right]^2 + \left[\text{Im}(G(j\omega))\right]^2 \\ &= \frac{4K^2}{(1+\omega^2)^2} + \frac{K^2(1-\omega^2)^2}{\omega^2(1+\omega^2)^2} \\ &= \frac{K^2}{\omega^2(1+\omega^2)^2} \left[ 4\omega^2 + (1 - 2\omega^2 + \omega^4) \right] \\ &= \frac{K^2(1 + 2\omega^2 + \omega^4)}{\omega^2(1+\omega^2)^2} \\ &= \frac{K^2(1+\omega^2)^2}{\omega^2(1+\omega^2)^2} = \frac{K^2}{\omega^2} \end{aligned}

取正平方根可得(因 K>0K > 0、ω>0\omega > 0):

∣G(jω)∣=Kω|G(j\omega)| = \frac{K}{\omega}

步驟二:由 ωgc=3 rad/sec\omega_{gc} = \sqrt{3}\text{ rad/sec} 求對應增益 KK

依增益交越頻率定義 ∣G(jωgc)∣=1|G(j\omega_{gc})| = 1:

K3=1  ⟹  K=3\frac{K}{\sqrt{3}} = 1 \implies K = \sqrt{3}

步驟三:求相位交越頻率 ωpc\omega_{pc}

在相位交越頻率處,開環頻率響應與負實軸相交,其虛部必須為零且實部為負:

Im[G(jωpc)]=K(1−ωpc2)ωpc(1+ωpc2)=0\text{Im}[G(j\omega_{pc})] = \frac{K(1-\omega_{pc}^2)}{\omega_{pc}(1+\omega_{pc}^2)} = 0

因 ωpc>0\omega_{pc} > 0,故分子 1−ωpc2=01 - \omega_{pc}^2 = 0,解得:

🔒

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

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

免費註冊
📄 以下 3 題共用同一段題幹
  1. (12%) Given the system y...+4y¨+5y˙+2y=u\dddot{y} + 4\ddot{y} + 5\dot{y} + 2y = u. Define the states x1≜yx_1 \triangleq y, x2≜y˙x_2 \triangleq \dot{y}, and x3≜y¨x_3 \triangleq \ddot{y}.

第 5-(a) 題5 分

(a) (5%) Sketch the state equation description with Controllability Canonical Form.

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

這一題的完整詳解

核心觀念

本題評量**狀態空間表示法(State-Space Representation)中的可控制典型態(Controllability Canonical Form, CCF)**建構:

  1. 狀態變數定義:以系統輸出及其各階導數作為狀態變數(又稱相位變數 Phase Variables):
    x1=y,x2=y˙,x3=y¨x_1 = y, \quad x_2 = \dot{y}, \quad x_3 = \ddot{y}
  2. 狀態方程式(State Equation)標準形式:
    x˙(t)=Ax(t)+Bu(t)\dot{\mathbf{x}}(t) = \mathbf{A}\mathbf{x}(t) + \mathbf{B}u(t)
    y(t)=Cx(t)+Du(t)\mathbf{y}(t) = \mathbf{C}\mathbf{x}(t) + \mathbf{D}u(t)
  3. 可控制典型態(CCF)結構特徵:
    當狀態變數依階數由低往高定義(x1=y,x2=y˙,…x_1 = y, x_2 = \dot{y}, \dots)時,狀態矩陣 A\mathbf{A} 的前 n−1n-1 列為將對角線右上方元素設為 11 的移位矩陣,最後一列為微分方程式各項係數變號(負係數),輸入向量 B=[0⋯01]T\mathbf{B} = \begin{bmatrix} 0 & \cdots & 0 & 1 \end{bmatrix}^T。

解題方法

步驟一:由題目定義建立狀態變數的一階微分關係
根據題目給定的狀態定義:
x1=y  ⟹  x˙1=y˙=x2x_1 = y \implies \dot{x}_1 = \dot{y} = x_2
x2=y˙  ⟹  x˙2=y¨=x3x_2 = \dot{y} \implies \dot{x}_2 = \ddot{y} = x_3
x3=y¨  ⟹  x˙3=y...x_3 = \ddot{y} \implies \dot{x}_3 = \dddot{y}

步驟二:將最高階導數表示為狀態變數與輸入的線性組合
將微分方程式移項,以求出 y...\dddot{y}:
y...+4y¨+5y˙+2y=u\dddot{y} + 4\ddot{y} + 5\dot{y} + 2y = u
y...=−2y−5y˙−4y¨+u\dddot{y} = -2y - 5\dot{y} - 4\ddot{y} + u

將狀態變數 x1,x2,x3x_1, x_2, x_3 代入:
x˙3=−2x1−5x2−4x3+u\dot{x}_3 = -2x_1 - 5x_2 - 4x_3 + u

步驟三:寫出矩陣向量形式之狀態方程式與輸出方程式
綜合以上三個一階微分方程式:

🔒

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

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

免費註冊

第 5-(b) 題4 分

(b) (4%) If all states can be measured, design the state feedback controller, u=−Kx+ru = -Kx + r such that the desired closed-loop poles of the controller locates at −2±j2-2 \pm j2 and −10-10.

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

這一題的完整詳解

核心觀念

本題考查現代控制理論中的極點安置(Pole Placement via State Feedback)。
主要涵蓋以下核心觀念與公式:

  1. 可控制性(Controllability):當系統完全可控制時,透過全狀態回授(Full-State Feedback)u=−Kx+ru = -Kx + r,可以任意配置閉迴路極點至複數平面上的任意共軛位置。
  2. 相變數規範型(Phase-Variable Canonical Form / Controller Canonical Form):狀態變數若依序定義為微分階次 x1=y,x2=y˙,x3=y¨x_1 = y, x_2 = \dot{y}, x_3 = \ddot{y},其狀態矩陣 AA 與輸入矩陣 BB 會呈現特殊的相變數形式,開迴路特徵多項式係數直接對應至矩陣 AA 的最後一列。
  3. 特徵方程比較法(Direct Characteristic Polynomial Matching):閉迴路系統矩陣為 A−BKA - BK。透過將閉迴路特徵多項式 det⁡(sI−(A−BK))\det(sI - (A - BK)) 與題目要求的期望特徵多項式 Δd(s)\Delta_d(s) 之各次項係數互相對照,即可直接解出回授增益矩陣 K=[k1k2k3]K = \begin{bmatrix} k_1 & k_2 & k_3 \end{bmatrix}。

解題方法

步驟一:建立系統狀態空間表示式(State-Space Representation)

根據題幹給定的微分方程式:

y...+4y¨+5y˙+2y=u\dddot{y} + 4\ddot{y} + 5\dot{y} + 2y = u

並依照狀態定義 x1=y,x2=y˙,x3=y¨x_1 = y, x_2 = \dot{y}, x_3 = \ddot{y},可列出狀態微分方程式:

{x˙1=x2x˙2=x3x˙3=y...=−2x1−5x2−4x3+u\begin{cases} \dot{x}_1 = x_2 \\ \dot{x}_2 = x_3 \\ \dot{x}_3 = \dddot{y} = -2x_1 - 5x_2 - 4x_3 + u \end{cases}

以矩陣形式表示 x˙=Ax+Bu\dot{x} = Ax + Bu:

x˙=[x˙1x˙2x˙3]=[010001−2−5−4][x1x2x3]+[001]u\dot{x} = \begin{bmatrix} \dot{x}_1 \\ \dot{x}_2 \\ \dot{x}_3 \end{bmatrix} = \begin{bmatrix} 0 & 1 & 0 \\ 0 & 0 & 1 \\ -2 & -5 & -4 \end{bmatrix} \begin{bmatrix} x_1 \\ x_2 \\ x_3 \end{bmatrix} + \begin{bmatrix} 0 \\ 0 \\ 1 \end{bmatrix} u

其中:

A=[010001−2−5−4],B=[001]A = \begin{bmatrix} 0 & 1 & 0 \\ 0 & 0 & 1 \\ -2 & -5 & -4 \end{bmatrix}, \quad B = \begin{bmatrix} 0 \\ 0 \\ 1 \end{bmatrix}

步驟二:導入狀態回授控制器並求閉迴路系統矩陣

設狀態回授增益向量為:

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

將控制法則 u=−Kx+ru = -Kx + r 代入原系統 x˙=Ax+Bu\dot{x} = Ax + Bu:

x˙=(A−BK)x+Br\dot{x} = (A - BK)x + Br

計算 BKBK:

BK=[001][k1k2k3]=[000000k1k2k3]BK = \begin{bmatrix} 0 \\ 0 \\ 1 \end{bmatrix} \begin{bmatrix} k_1 & k_2 & k_3 \end{bmatrix} = \begin{bmatrix} 0 & 0 & 0 \\ 0 & 0 & 0 \\ k_1 & k_2 & k_3 \end{bmatrix}

閉迴路系統矩陣 Acl=A−BKA_{cl} = A - BK 為:

A−BK=[010001−(2+k1)−(5+k2)−(4+k3)]A - BK = \begin{bmatrix} 0 & 1 & 0 \\ 0 & 0 & 1 \\ -(2 + k_1) & -(5 + k_2) & -(4 + k_3) \end{bmatrix}

閉迴路系統特徵方程式為 det⁡(sI−(A−BK))=0\det(sI - (A - BK)) = 0。由相變數矩陣的性質可知:

🔒

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

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

免費註冊

第 5-(c) 題3 分

(c) (3%) If only y(t)y(t) can be measured, an observer design is necessary by output feedback. Design the observer, so that the closed-loop poles of the observer locates at −10±j10-10 \pm j10 and −20-20.

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

這一題的完整詳解

核心觀念

  • 可觀測性:由題目可得 (A,C)(A,C) 必須可觀測,才能以輸出回授設計觀測器。
  • 觀測器設計:觀測器的誤差 dynamics 為 e˙=(A−LC)e\dot e = (A-LC)e,其中 e=x−x^e = x-\hat x。給定欲配置的極點,等價於在 A−LCA-LC 上做極點配置。
  • 極點配置(Pole Placement):利用 Ackermann 公式(或直接係數比較)求得觀測增益 LL 使得特徵多項式匹配目標極點。

系統狀態方程

x1=y,x2=y˙,x3=y¨x˙1=x2x˙2=x3x˙3=u−4x3−5x2−2x1\begin{aligned} x_1 &= y, \qquad x_2 = \dot y,\qquad x_3 = \ddot y \\ \dot x_1 &= x_2 \\ \dot x_2 &= x_3 \\ \dot x_3 &= u-4x_3-5x_2-2x_1 \end{aligned}

寫成矩陣形式

x˙=Ax+Bu,A=[010001−2−5−4],  B=[001],  C=[100]\dot x = A x + Bu,\qquad A=\begin{bmatrix} 0 & 1 & 0\\ 0 & 0 & 1\\ -2 & -5 & -4 \end{bmatrix}, \;B=\begin{bmatrix}0\\0\\1\end{bmatrix}, \;C=\begin{bmatrix}1 & 0 & 0\end{bmatrix}

目標觀測器極點

給定極點:−10±j10-10\pm j10(共軛一對)與 −20-20。
其特徵多項式

(s+20)(s2+20s+200)=s3+40s2+600s+4000(1)\begin{aligned} (s+20)(s^2+20s+200) &= s^3+40s^2+600s+4000 \tag{1} \end{aligned}

觀測增益 LL 的求法

觀測器誤差 dynamics 為

e˙=(A−LC)e,LC=L[100]=[l100l200l300]\dot e =(A-LC)e,\qquad LC = L\begin{bmatrix}1&0&0\end{bmatrix} = \begin{bmatrix} l_1 & 0 & 0\\ l_2 & 0 & 0\\ l_3 & 0 & 0 \end{bmatrix} A−LC=[−l110−l201−2−l3−5−4]A-LC= \begin{bmatrix} -l_1 & 1 & 0\\ -l_2 & 0 & 1\\ -2-l_3 & -5 & -4 \end{bmatrix}

特徵多項式 det⁡(sI−(A−LC))\det (sI-(A-LC)) 計算如下:

det⁡ ⁣(sI−(A−LC))=det⁡ ⁣[s+l1−10l2s−12+l35s+4]=(s+l1)(s2+4s+5)+l2(s+4)+(2+l3)=s3+(4+l1)s2+(5+4l1+l2)s+(5l1+4l2+2+l3)\begin{aligned} \det\!\bigl(sI-(A-LC)\bigr) &=\det\!\begin{bmatrix} s+l_1 & -1 & 0\\[2pt] l_2 & s & -1\\[2pt] 2+l_3 & 5 & s+4 \end{bmatrix}\\ &=(s+l_1)(s^2+4s+5)+l_2(s+4)+(2+l_3) \\ &=s^3+(4+l_1)s^2+(5+4l_1+l_2)s+(5l_1+4l_2+2+l_3) \end{aligned}

將此結果與目標多項式 (1) 比較係數:

🔒

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

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

免費註冊

其他考古題