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

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

第 1 題25 分

Consider the feedback control system in Figure 1. Suppose that
G(s)=1s(s+1),C(s)=K(s+b)s+aG(s) = \frac{1}{s(s+1)}, \quad C(s) = \frac{K(s+b)}{s+a}
where a,b,K>0a, b, K > 0. Answer the following questions.
🖼️【此處有附圖,請對照原卷】Figure 1: Feedback Control System

(a) Find the range of KK (in terms of aa and bb) in which the closed-loop system is stable.
(b) Let b=2(a+1)b = 2(a+1). Find KK (in terms of aa) such that the closed-loop system has two pure imaginary poles. What are these pure imaginary poles (in terms of aa)?
(c) Let b=0.5b = 0.5. Find the numerical range of KK such that the steady-state error of the closed-loop system with respect to the parabolic input, i.e. r(t)=t2r(t) = t^2 for t≥0t \ge 0, is less than 0.010.01.
(d) Find the values of a,ba, b and KK such that the closed-loop transfer function from rr to yy is a second-order system with damping ratio ζ=0.5\zeta = 0.5 and natural frequency ωn=4\omega_n = 4.

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

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

這一題的完整詳解

核心觀念

本題考查單位負回授的閉迴路特徵方程、三次系統的 Routh–Hurwitz 穩定判據、系統型別與穩態誤差,以及透過極零相消得到指定的二階閉迴路傳遞函數。

解題方法

依原圖,加總點是單位負回授:參考輸入 rr 接正端,輸出 yy 回授至負端;前向串接為 C(s)=K(s+b)s+aC(s)=\dfrac{K(s+b)}{s+a} 與 G(s)=1s(s+1)G(s)=\dfrac{1}{s(s+1)}。因此

L(s)=C(s)G(s)=K(s+b)(s+a)s(s+1),T(s)=L(s)1+L(s).L(s)=C(s)G(s)=\frac{K(s+b)}{(s+a)s(s+1)},\qquad T(s)=\frac{L(s)}{1+L(s)}.

閉迴路特徵多項式為

P(s)=(s+a)s(s+1)+K(s+b)=s3+(a+1)s2+(a+K)s+Kb.P(s)=(s+a)s(s+1)+K(s+b) =s^3+(a+1)s^2+(a+K)s+Kb.

(a) 閉迴路穩定的 KK 範圍

對 P(s)P(s) 建立 Routh 表:

s31a+Ks2a+1Kbs1(a+1)(a+K)−Kba+10s0Kb\begin{array}{c|cc} s^3 & 1 & a+K\\ s^2 & a+1 & Kb\\ s^1 & \dfrac{(a+1)(a+K)-Kb}{a+1} & 0\\ s^0 & Kb & \end{array}

由於 a,b,K>0a,b,K>0,第一欄除 s1s^1 列外皆為正。穩定條件是

(a+1)(a+K)−Kb>0⟺a(a+1)+K(a+1−b)>0.(a+1)(a+K)-Kb>0 \quad\Longleftrightarrow\quad a(a+1)+K(a+1-b)>0.

若 b≤a+1b\le a+1,此不等式對所有 K>0K>0 成立。若 b>a+1b>a+1,則須有 K<a(a+1)b−a−1K<\dfrac{a(a+1)}{b-a-1}。

(b) 具有兩個純虛數極點

令 b=2(a+1)b=2(a+1)。Routh 表中 s1s^1 列的分子變為

(a+1)(a+K)−2K(a+1)=(a+1)(a−K).(a+1)(a+K)-2K(a+1)=(a+1)(a-K).

純虛數極點出現在穩定邊界,因此令其為零,得 K=aK=a。代回特徵多項式:

P(s)=s3+(a+1)s2+2as+2a(a+1)=(s+a+1)(s2+2a).P(s)=s^3+(a+1)s^2+2as+2a(a+1) =(s+a+1)(s^2+2a).

故純虛數極點為 s=±j2as=\pm j\sqrt{2a}。

(c) b=0.5b=0.5 時拋物線輸入的穩態誤差

開迴路 L(s)L(s) 在原點有一個極點,因此是第一型系統。其加速度誤差常數為

Ka=lim⁡s→0s2L(s)=lim⁡s→0Ks(s+b)(s+a)(s+1)=0.K_a=\lim_{s\to 0}s^2L(s) =\lim_{s\to 0}\frac{Ks(s+b)}{(s+a)(s+1)}=0.
🔒

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

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

免費註冊

第 2 題25 分

Consider the feedback control system in Figure 1. Suppose that
G(s)=s2+3s+2.5s2(s2+3s+2)G(s) = \frac{s^2 + 3s + 2.5}{s^2(s^2 + 3s + 2)}
Answer the following questions.
🖼️【此處有附圖,請對照原卷】Figure 1: Feedback Control System

(a) Let C(s)=KC(s) = K. Show that the closed-loop system is unstable for any KK.
(b) Let C(s)=KC(s) = K. Sketch the root locus of the system for K>0K > 0. You need to specify the intersection and angles of the asymptotes, the departure angles of the loci starting from the poles at s=0s = 0, and the arrival angles of the loci entering the open-loop zeros.
(c) Let C(s)=Ks+bs+aC(s) = K\frac{s+b}{s+a}. Find the numerical values of aa and bb such that there exists positive gain KK that stabilizes the closed-loop system. Sketch the root locus to verify your design.

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

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

這一題的完整詳解

核心觀念

此題考查單位負回授系統的:

  • 閉迴路特徵方程式;
  • Routh–Hurwitz 穩定判據;
  • 根軌跡的實軸區段、漸近線、出發角與到達角;
  • 利用補償器 C(s)=Ks+bs+aC(s)=K\dfrac{s+b}{s+a} 改變根軌跡以達成穩定。

由圖為單位負回授,因此閉迴路特徵方程式為

1+C(s)G(s)=0.1+C(s)G(s)=0.

令

N(s)=s2+3s+2.5=(s+1.5−0.5j)(s+1.5+0.5j),N(s)=s^2+3s+2.5=(s+1.5-0.5j)(s+1.5+0.5j), D(s)=s2(s2+3s+2)=s2(s+1)(s+2).D(s)=s^2(s^2+3s+2)=s^2(s+1)(s+2).

(a) 證明任意 K>0K>0 時系統皆不穩定

令 C(s)=KC(s)=K,則

1+KN(s)D(s)=0,1+K\frac{N(s)}{D(s)}=0,

所以閉迴路特徵方程式為

D(s)+KN(s)=0.D(s)+KN(s)=0.

展開得

s2(s2+3s+2)+K(s2+3s+2.5)=0,s^2(s^2+3s+2)+K(s^2+3s+2.5)=0,

即

s4+3s3+(K+2)s2+3Ks+2.5K=0.s^4+3s^3+(K+2)s^2+3Ks+2.5K=0.

Routh 表為

s41K+22.5Ks333K0s23(K+2)−3K3=22.5K0s12(3K)−3(2.5K)2=−0.75K00s02.5K\begin{array}{c|ccc} s^4 & 1 & K+2 & 2.5K\\ s^3 & 3 & 3K & 0\\ s^2 & \dfrac{3(K+2)-3K}{3}=2 & 2.5K & 0\\ s^1 & \dfrac{2(3K)-3(2.5K)}{2} =-0.75K & 0 & 0\\ s^0 & 2.5K & & \end{array}

第一欄為

1,3,2,−0.75K,2.5K.1,\quad 3,\quad 2,\quad -0.75K,\quad 2.5K.

當 K>0K>0 時,符號為

+,+,+,−,+.+,\quad +,\quad +,\quad -,\quad +.

共有兩次符號變化,因此有兩個右半平面根。故閉迴路系統對任意 K>0K>0 皆不穩定。


(b) C(s)=KC(s)=K 時的根軌跡

開迴路極點與零點

開迴路傳遞函數為

L(s)=Ks2+3s+2.5s2(s+1)(s+2).L(s)=K\frac{s^2+3s+2.5}{s^2(s+1)(s+2)}.

極點為

s=0, 0, −1, −2,s=0,\ 0,\ -1,\ -2,

零點為

z1,2=−1.5±0.5j.z_{1,2}=-1.5\pm0.5j.

實軸根軌跡

實軸上的點若其右側的開迴路極點與零點總數為奇數,則位於根軌跡上。

因此實軸根軌跡只有

−2<s<−1.-2<s<-1.

漸近線

極點數 n=4n=4,零點數 m=2m=2,故有

n−m=2n-m=2

條漸近線。

漸近線交點為

σa=∑pi−∑zin−m.\sigma_a= \frac{\sum p_i-\sum z_i}{n-m}.

其中

∑pi=0+0−1−2=−3,\sum p_i=0+0-1-2=-3, ∑zi=(−1.5+0.5j)+(−1.5−0.5j)=−3.\sum z_i=(-1.5+0.5j)+(-1.5-0.5j)=-3.

所以

σa=−3−(−3)2=0.\sigma_a=\frac{-3-(-3)}{2}=0.

漸近線角度為

θq=(2q+1)180∘2,q=0,1.\theta_q=\frac{(2q+1)180^\circ}{2}, \qquad q=0,1.

因此

θ=90∘, 270∘.\boxed{\theta=90^\circ,\ 270^\circ}.

也就是通過原點的上下兩條垂直漸近線。


s=0s=0 雙重極點的出發角

在 s=0s=0 附近,其他極點為 −1,−2-1,-2,兩個零點為

−1.5±0.5j.-1.5\pm0.5j.

由角度條件可得

∑∠(s−zi)−∑∠(s−pi)=180∘.\sum \angle(s-z_i) -\sum \angle(s-p_i) =180^\circ.

兩個零點相對於原點的角度和為 360∘360^\circ;其他兩個極點的角度和亦為 360∘360^\circ。設出發角為 θ\theta,則雙重極點貢獻 2θ2\theta:

360∘−(360∘+2θ)=180∘(mod360∘).360^\circ-(360^\circ+2\theta) =180^\circ \pmod {360^\circ}.

因此

−2θ=180∘(mod360∘),-2\theta=180^\circ \pmod {360^\circ},

得到

θ=90∘, 270∘.\boxed{\theta=90^\circ,\ 270^\circ}.

兩條根軌跡由原點沿著上下方向離開。


複數零點的到達角

考慮上方零點

z1=−1.5+0.5j.z_1=-1.5+0.5j.

其他零點相對於 z1z_1 的角度為 90∘90^\circ。

各極點相對於 z1z_1 的角度為

∠(z1−0)=161.565∘,\angle(z_1-0)=161.565^\circ, ∠(z1+1)=135∘,∠(z1+2)=45∘.\angle(z_1+1)=135^\circ, \qquad \angle(z_1+2)=45^\circ.

由於 s=0s=0 為雙重極點,極點角度總和為

🔒

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

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

免費註冊

第 3 題25 分

Let the state equations of an LTI system be represented by
dx(t)dt=Ax(t)+Bu(t)\frac{dx(t)}{dt} = Ax(t) + Bu(t)
y(t)=Cx(t)y(t) = Cx(t)
where
A=[−10α0−2100γ],B=[β10],C=[111]A = \begin{bmatrix} -1 & 0 & \alpha \\ 0 & -2 & 1 \\ 0 & 0 & \gamma \end{bmatrix}, \quad B = \begin{bmatrix} \beta \\ 1 \\ 0 \end{bmatrix}, \quad C = \begin{bmatrix} 1 & 1 & 1 \end{bmatrix}
(a) Find the state transition matrix.
(b) With the initial condition x(0)=0x(0)=0 and β=1\beta=1, choose BB such that the transient response of y(t)y(t) has only the component of te−2tte^{-2t}.
(c) With B=[010]B = \begin{bmatrix} 0 \\ 1 \\ 0 \end{bmatrix} and u(t)u(t), a unit-step input, choose the initial condition x(0)=[x1(0) x2(0) x3(0)]Tx(0) = [x_1(0) \ x_2(0) \ x_3(0)]^T such that x(t)x(t) does NOT contain any transient parts, i.e., no e−te^{-t}, e−2te^{-2t} and te−2tte^{-2t}.

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

這一題的完整詳解

核心觀念

本題考查:

  1. 狀態轉移矩陣 Φ(t)=eAt\Phi(t)=e^{At}。
  2. 上三角系統矩陣指數的計算。
  3. 單位階輸入下的零狀態響應:
    x(t)=∫0teA(t−τ)B dτx(t)=\int_0^t e^{A(t-\tau)}B\,d\tau
  4. 穩態值與暫態項的分離。
  5. 重根 γ=−2\gamma=-2 所造成的 te−2tte^{-2t} 項。

(a) 狀態轉移矩陣

由於 AA 為上三角矩陣,先求各狀態間的傳遞關係。

令

Φ(t)=eAt=[ϕ11(t)ϕ12(t)ϕ13(t)ϕ21(t)ϕ22(t)ϕ23(t)ϕ31(t)ϕ32(t)ϕ33(t)]\Phi(t)=e^{At} = \begin{bmatrix} \phi_{11}(t)&\phi_{12}(t)&\phi_{13}(t)\\ \phi_{21}(t)&\phi_{22}(t)&\phi_{23}(t)\\ \phi_{31}(t)&\phi_{32}(t)&\phi_{33}(t) \end{bmatrix}

對角元素直接為

ϕ11(t)=e−t,ϕ22(t)=e−2t,ϕ33(t)=eγt\phi_{11}(t)=e^{-t},\qquad \phi_{22}(t)=e^{-2t},\qquad \phi_{33}(t)=e^{\gamma t}

因為 AA 的第 (1,2)(1,2)、(2,1)(2,1)、(3,1)(3,1)、(3,2)(3,2) 元素皆為零,所以

ϕ12(t)=ϕ21(t)=ϕ31(t)=ϕ32(t)=0\phi_{12}(t)=\phi_{21}(t)=\phi_{31}(t)=\phi_{32}(t)=0

對於 ϕ23(t)\phi_{23}(t):

x˙2+2x2=x3\dot{x}_2+2x_2=x_3

因此

ϕ23(t)=∫0te−2(t−τ)eγτ dτ=eγt−e−2tγ+2\phi_{23}(t) = \int_0^t e^{-2(t-\tau)}e^{\gamma\tau}\,d\tau = \frac{e^{\gamma t}-e^{-2t}}{\gamma+2}

當 γ=−2\gamma=-2 時取極限,得到

ϕ23(t)=te−2t\phi_{23}(t)=te^{-2t}

對於 ϕ13(t)\phi_{13}(t):

x˙1+x1=αx3\dot{x}_1+x_1=\alpha x_3

所以

ϕ13(t)=α∫0te−(t−τ)eγτ dτ=αeγt−e−tγ+1\phi_{13}(t) = \alpha\int_0^t e^{-(t-\tau)}e^{\gamma\tau}\,d\tau = \alpha\frac{e^{\gamma t}-e^{-t}}{\gamma+1}

當 γ=−1\gamma=-1 時取極限,得到

ϕ13(t)=αte−t\phi_{13}(t)=\alpha te^{-t}

因此,當 γ≠−1,−2\gamma\neq -1,-2 時,

Φ(t)=[e−t0αeγt−e−tγ+10e−2teγt−e−2tγ+200eγt]\boxed{ \Phi(t)= \begin{bmatrix} e^{-t} & 0 & \displaystyle \alpha\frac{e^{\gamma t}-e^{-t}}{\gamma+1}\\[8pt] 0 & e^{-2t} & \displaystyle \frac{e^{\gamma t}-e^{-2t}}{\gamma+2}\\[8pt] 0&0&e^{\gamma t} \end{bmatrix} }

特殊情形如下。

當 γ=−2\gamma=-2 時,

Φ(t)=[e−t0α(e−t−e−2t)0e−2tte−2t00e−2t]\boxed{ \Phi(t)= \begin{bmatrix} e^{-t} & 0 & \alpha(e^{-t}-e^{-2t})\\ 0 & e^{-2t} & te^{-2t}\\ 0&0&e^{-2t} \end{bmatrix} }

當 γ=−1\gamma=-1 時,

Φ(t)=[e−t0αte−t0e−2te−t−e−2t00e−t]\boxed{ \Phi(t)= \begin{bmatrix} e^{-t} & 0 & \alpha te^{-t}\\ 0&e^{-2t}&e^{-t}-e^{-2t}\\ 0&0&e^{-t} \end{bmatrix} }

(b) 零初始條件下的暫態輸出

題目給定

x(0)=0,β=1x(0)=0,\qquad \beta=1

因此

B=[110]B= \begin{bmatrix} 1\\1\\0 \end{bmatrix}

且輸入為單位階輸入 u(t)=1u(t)=1。

由於 BB 的第三個元素為零,且 x3(0)=0x_3(0)=0,有

x˙3=γx3\dot{x}_3=\gamma x_3

因此

x3(t)=0x_3(t)=0

此時前兩個狀態彼此獨立:

x˙1=−x1+1\dot{x}_1=-x_1+1

得到

x1(t)=1−e−tx_1(t)=1-e^{-t}

以及

x˙2=−2x2+1\dot{x}_2=-2x_2+1

得到

x2(t)=12(1−e−2t)x_2(t)=\frac{1}{2}\left(1-e^{-2t}\right)

輸出為

y(t)=x1(t)+x2(t)+x3(t)y(t)=x_1(t)+x_2(t)+x_3(t)

所以

y(t)=1−e−t+12(1−e−2t)y(t) = 1-e^{-t} +\frac{1}{2}(1-e^{-2t})

即

y(t)=32−e−t−12e−2ty(t)=\frac{3}{2}-e^{-t}-\frac{1}{2}e^{-2t}

暫態響應為

ytr(t)=−e−t−12e−2ty_{\mathrm{tr}}(t) = -e^{-t}-\frac{1}{2}e^{-2t}

其中沒有 te−2tte^{-2t} 項,而且 α\alpha、γ\gamma 均不影響此結果,原因是 x3(t)x_3(t) 永遠為零。

因此,依照題面所給的

B=[β10],β=1B= \begin{bmatrix} \beta\\1\\0 \end{bmatrix}, \qquad \beta=1
🔒

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

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

免費註冊

第 4 題25 分

A unity-feedback system with the open loop transfer function L(s)=K(s−z)s(s+p)L(s) = \frac{K(s-z)}{s(s+p)}, where pp and zz are all positive.
L(jω)=K(jω−z)jω(jω+p)=K(p+jω)jω(p+jω)=K(p+jω)pjω+(jω)2=K(p+jω)−ω2+pjωL(j\omega) = \frac{K(j\omega-z)}{j\omega(j\omega+p)} = \frac{K(p+j\omega)}{j\omega(p+j\omega)} = \frac{K(p+j\omega)}{p j\omega + (j\omega)^2} = \frac{K(p+j\omega)}{-\omega^2 + p j\omega}
, and pp, and zz are all positive.

(a) Sketch the Nyquist plot for K>0K>0.
(b) Find the stability region of KK and the oscillation frequency if exists from Nyquist plot.
(c) Let p=z=1p = z = 1 and K=1/3K = 1/\sqrt{3}. Find the Gain margin and Phase Margin and the corresponding gain crossover frequency and phase crossover frequency, if available.

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

這一題的完整詳解

核心觀念

本題考查:

  1. Nyquist 圖的繪製與關鍵頻率判斷。
  2. Nyquist 穩定判據。
  3. 閉迴路特徵方程式與穩定範圍。
  4. Gain Margin(增益裕度)與 Phase Margin(相位裕度)。

單位回授系統的閉迴路特徵方程式為

1+L(s)=0.1+L(s)=0.

其中

L(s)=K(s−z)s(s+p),p>0, z>0.L(s)=\frac{K(s-z)}{s(s+p)},\qquad p>0,\ z>0.

題目中所列的 L(jω)L(j\omega) 中間推導將 jω−zj\omega-z 寫成 p+jωp+j\omega,此等式並不成立。以下依原始傳遞函數正確計算。


(a)K>0K>0 時的 Nyquist 圖

令 s=jωs=j\omega:

L(jω)=K(jω−z)jω(p+jω).L(j\omega) =\frac{K(j\omega-z)}{j\omega(p+j\omega)}.

分母為

jω(p+jω)=−ω2+jpω.j\omega(p+j\omega)=-\omega^2+jp\omega.

乘以分母共軛,得

L(jω)=K(jω−z)(−ω2−jpω)ω4+p2ω2.L(j\omega) = K\frac{(j\omega-z)(-\omega^2-jp\omega)} {\omega^4+p^2\omega^2}.

整理後:

L(jω)=K[p+zω2+p2+jpz−ω2ω(ω2+p2)].L(j\omega) = K\left[ \frac{p+z}{\omega^2+p^2} +j\frac{pz-\omega^2}{\omega(\omega^2+p^2)} \right].

因此

Re⁡{L(jω)}=K(p+z)ω2+p2>0\operatorname{Re}\{L(j\omega)\} = \frac{K(p+z)}{\omega^2+p^2}>0

對所有 ω>0\omega>0 均成立。故正頻率部分的 Nyquist 軌跡完全位於右半平面。

關鍵位置

當 ω→0+\omega\to 0^+ 時:

Re⁡{L(jω)}→K(p+z)p2,\operatorname{Re}\{L(j\omega)\}\to \frac{K(p+z)}{p^2}, Im⁡{L(jω)}∼Kzpω→+∞.\operatorname{Im}\{L(j\omega)\} \sim \frac{Kz}{p\omega}\to+\infty.

所以軌跡由正虛軸無限遠附近開始。

當 ω→∞\omega\to\infty 時:

L(jω)→0,L(j\omega)\to 0,

且

Re⁡{L(jω)}→0+,Im⁡{L(jω)}→0−.\operatorname{Re}\{L(j\omega)\}\to 0^+,\qquad \operatorname{Im}\{L(j\omega)\}\to 0^-.

所以軌跡最後由第四象限趨近原點。

與實軸交點

令虛部為零:

K(pz−ω2)ω(ω2+p2)=0.\frac{K(pz-\omega^2)} {\omega(\omega^2+p^2)}=0.

因此

ω=pz.\omega=\sqrt{pz}.

此時實部為

Re⁡{L(jpz)}=K(p+z)pz+p2=Kp.\operatorname{Re}\{L(j\sqrt{pz})\} = \frac{K(p+z)}{pz+p^2} = \frac{K}{p}.

所以正頻率曲線通過實軸上的點

Kp.\boxed{\frac{K}{p}}.

負頻率部分為正頻率部分對實軸的鏡射。

概略形狀如下:

  • 從上方無限遠開始;
  • 位於右半平面;
  • 以 ω=pz\omega=\sqrt{pz} 穿越正實軸 K/pK/p;
  • 最後由下方趨近原點;
  • 再以實軸對稱形成完整 Nyquist 圖。

此外,由於開迴路含有原點極點,嚴格繪製 Nyquist 圖時須對原點作小半圓凹避;此凹避部分對穩定性判斷不可忽略。


(b)KK 的穩定範圍與振盪頻率

方法一:由閉迴路特徵方程式判斷

由

1+K(s−z)s(s+p)=01+\frac{K(s-z)}{s(s+p)}=0

可得

s(s+p)+K(s−z)=0.s(s+p)+K(s-z)=0.

因此閉迴路特徵方程式為

s2+(p+K)s−Kz=0.s^2+(p+K)s-Kz=0.

對二階系統

s2+a1s+a0=0s^2+a_1s+a_0=0

而言,漸近穩定的必要且充分條件為

a1>0,a0>0.a_1>0,\qquad a_0>0.

本題中

a1=p+K,a0=−Kz.a_1=p+K,\qquad a_0=-Kz.

由於 p,z>0p,z>0,故

p+K>0⇒K>−p,p+K>0\quad\Rightarrow\quad K>-p,

以及

−Kz>0⇒K<0.-Kz>0\quad\Rightarrow\quad K<0.

合併得

−p<K<0.\boxed{-p<K<0}.

在 K=0K=0 時,特徵方程式為

s(s+p)=0,s(s+p)=0,

含有原點極點,屬於臨界穩定而非漸近穩定。


方法二:由 Nyquist 圖判斷

當 K<0K<0 時,Nyquist 圖相當於將 K>0K>0 的圖形反轉,並位於左半平面。其與實軸交點為

−∣K∣p=Kp.-\frac{|K|}{p}=\frac{K}{p}.

當此交點尚未越過 −1-1 時,系統穩定:

Kp>−1.\frac{K}{p}>-1.

因此

K>−p.K>-p.

再配合 K<0K<0,得到

−p<K<0.\boxed{-p<K<0}.

當 K>0K>0 時,閉迴路特徵方程式的常數項 −Kz<0-Kz<0,兩根乘積為負,因此必有一根位於右半平面,系統必不穩定。


臨界振盪條件

令閉迴路極點位於虛軸上:

s=jω,ω>0.s=j\omega,\qquad \omega>0.

代入特徵方程式:

(jω)2+(p+K)jω−Kz=0.(j\omega)^2+(p+K)j\omega-Kz=0.

整理實部與虛部:

−ω2−Kz+jω(p+K)=0.-\omega^2-Kz+j\omega(p+K)=0.

因此必須同時滿足

−ω2−Kz=0,-\omega^2-Kz=0, ω(p+K)=0.\omega(p+K)=0.

因為 ω>0\omega>0,所以

K=−p.K=-p.

再代入實部方程式:

−ω2+pz=0.-\omega^2+pz=0.

故臨界振盪角頻率為

ωosc=pz rad/s.\boxed{\omega_{\mathrm{osc}}=\sqrt{pz}\ \text{rad/s}}.

此時特徵方程式為

s2+pz=0,s^2+pz=0,

極點為

🔒

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

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

免費註冊

第 3 題25 分

Let the state equations of an LTI system be represented by
dx(t)dt=Ax(t)+Bu(t)\frac{dx(t)}{dt} = Ax(t) + Bu(t)
y(t)=Cx(t)y(t) = Cx(t)
where
A=[−10α0−2100γ],B=[β10],C=[111]A = \begin{bmatrix} -1 & 0 & \alpha \\ 0 & -2 & 1 \\ 0 & 0 & \gamma \end{bmatrix}, \quad B = \begin{bmatrix} \beta \\ 1 \\ 0 \end{bmatrix}, \quad C = \begin{bmatrix} 1 & 1 & 1 \end{bmatrix}
(a) Find the state transition matrix.
(b) With the initial condition x(0)=0x(0)=0 and β=1\beta=1, choose BB such that the transient response of y(t)y(t) has only the component of te−2tte^{-2t}.
(c) With B=[010]B = \begin{bmatrix} 0 \\ 1 \\ 0 \end{bmatrix} and u(t)u(t), a unit-step input, choose the initial condition x(0)=[x1(0) x2(0) x3(0)]Tx(0) = [x_1(0) \ x_2(0) \ x_3(0)]^T such that x(t)x(t) does NOT contain any transient parts, i.e., no e−te^{-t}, e−2te^{-2t} and te−2tte^{-2t}.

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

這一題的完整詳解

核心觀念

本題考查:

  1. 狀態轉移矩陣 Φ(t)=eAt\Phi(t)=e^{At}。
  2. 上三角系統矩陣指數的計算。
  3. 單位階輸入下的零狀態響應:
    x(t)=∫0teA(t−τ)B dτx(t)=\int_0^t e^{A(t-\tau)}B\,d\tau
  4. 穩態值與暫態項的分離。
  5. 重根 γ=−2\gamma=-2 所造成的 te−2tte^{-2t} 項。

(a) 狀態轉移矩陣

由於 AA 為上三角矩陣,先求各狀態間的傳遞關係。

令

Φ(t)=eAt=[ϕ11(t)ϕ12(t)ϕ13(t)ϕ21(t)ϕ22(t)ϕ23(t)ϕ31(t)ϕ32(t)ϕ33(t)]\Phi(t)=e^{At} = \begin{bmatrix} \phi_{11}(t)&\phi_{12}(t)&\phi_{13}(t)\\ \phi_{21}(t)&\phi_{22}(t)&\phi_{23}(t)\\ \phi_{31}(t)&\phi_{32}(t)&\phi_{33}(t) \end{bmatrix}

對角元素直接為

ϕ11(t)=e−t,ϕ22(t)=e−2t,ϕ33(t)=eγt\phi_{11}(t)=e^{-t},\qquad \phi_{22}(t)=e^{-2t},\qquad \phi_{33}(t)=e^{\gamma t}

因為 AA 的第 (1,2)(1,2)、(2,1)(2,1)、(3,1)(3,1)、(3,2)(3,2) 元素皆為零,所以

ϕ12(t)=ϕ21(t)=ϕ31(t)=ϕ32(t)=0\phi_{12}(t)=\phi_{21}(t)=\phi_{31}(t)=\phi_{32}(t)=0

對於 ϕ23(t)\phi_{23}(t):

x˙2+2x2=x3\dot{x}_2+2x_2=x_3

因此

ϕ23(t)=∫0te−2(t−τ)eγτ dτ=eγt−e−2tγ+2\phi_{23}(t) = \int_0^t e^{-2(t-\tau)}e^{\gamma\tau}\,d\tau = \frac{e^{\gamma t}-e^{-2t}}{\gamma+2}

當 γ=−2\gamma=-2 時取極限,得到

ϕ23(t)=te−2t\phi_{23}(t)=te^{-2t}

對於 ϕ13(t)\phi_{13}(t):

x˙1+x1=αx3\dot{x}_1+x_1=\alpha x_3

所以

ϕ13(t)=α∫0te−(t−τ)eγτ dτ=αeγt−e−tγ+1\phi_{13}(t) = \alpha\int_0^t e^{-(t-\tau)}e^{\gamma\tau}\,d\tau = \alpha\frac{e^{\gamma t}-e^{-t}}{\gamma+1}

當 γ=−1\gamma=-1 時取極限,得到

ϕ13(t)=αte−t\phi_{13}(t)=\alpha te^{-t}

因此,當 γ≠−1,−2\gamma\neq -1,-2 時,

Φ(t)=[e−t0αeγt−e−tγ+10e−2teγt−e−2tγ+200eγt]\boxed{ \Phi(t)= \begin{bmatrix} e^{-t} & 0 & \displaystyle \alpha\frac{e^{\gamma t}-e^{-t}}{\gamma+1}\\[8pt] 0 & e^{-2t} & \displaystyle \frac{e^{\gamma t}-e^{-2t}}{\gamma+2}\\[8pt] 0&0&e^{\gamma t} \end{bmatrix} }

特殊情形如下。

當 γ=−2\gamma=-2 時,

Φ(t)=[e−t0α(e−t−e−2t)0e−2tte−2t00e−2t]\boxed{ \Phi(t)= \begin{bmatrix} e^{-t} & 0 & \alpha(e^{-t}-e^{-2t})\\ 0 & e^{-2t} & te^{-2t}\\ 0&0&e^{-2t} \end{bmatrix} }

當 γ=−1\gamma=-1 時,

Φ(t)=[e−t0αte−t0e−2te−t−e−2t00e−t]\boxed{ \Phi(t)= \begin{bmatrix} e^{-t} & 0 & \alpha te^{-t}\\ 0&e^{-2t}&e^{-t}-e^{-2t}\\ 0&0&e^{-t} \end{bmatrix} }

(b) 零初始條件下的暫態輸出

題目給定

x(0)=0,β=1x(0)=0,\qquad \beta=1

因此

B=[110]B= \begin{bmatrix} 1\\1\\0 \end{bmatrix}

且輸入為單位階輸入 u(t)=1u(t)=1。

由於 BB 的第三個元素為零,且 x3(0)=0x_3(0)=0,有

x˙3=γx3\dot{x}_3=\gamma x_3

因此

x3(t)=0x_3(t)=0

此時前兩個狀態彼此獨立:

x˙1=−x1+1\dot{x}_1=-x_1+1

得到

x1(t)=1−e−tx_1(t)=1-e^{-t}

以及

x˙2=−2x2+1\dot{x}_2=-2x_2+1

得到

x2(t)=12(1−e−2t)x_2(t)=\frac{1}{2}\left(1-e^{-2t}\right)

輸出為

y(t)=x1(t)+x2(t)+x3(t)y(t)=x_1(t)+x_2(t)+x_3(t)

所以

y(t)=1−e−t+12(1−e−2t)y(t) = 1-e^{-t} +\frac{1}{2}(1-e^{-2t})

即

y(t)=32−e−t−12e−2ty(t)=\frac{3}{2}-e^{-t}-\frac{1}{2}e^{-2t}

暫態響應為

ytr(t)=−e−t−12e−2ty_{\mathrm{tr}}(t) = -e^{-t}-\frac{1}{2}e^{-2t}

其中沒有 te−2tte^{-2t} 項,而且 α\alpha、γ\gamma 均不影響此結果,原因是 x3(t)x_3(t) 永遠為零。

因此,依照題面所給的

B=[β10],β=1B= \begin{bmatrix} \beta\\1\\0 \end{bmatrix}, \qquad \beta=1
🔒

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

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

免費註冊

第 4 題25 分

A unity-feedback system with the open loop transfer function L(s)=K(s−z)s(s+p)L(s) = \frac{K(s-z)}{s(s+p)}, where pp and zz are all positive.
(a) Sketch the Nyquist plot for K>0K>0.
(b) Find the stability region of KK and the oscillation frequency if exists from Nyquist plot.
(c) Let p=z=1p = z = 1 and K=1/3K = 1/\sqrt{3}. Find the Gain margin and Phase Margin and the corresponding gain crossover frequency and phase crossover frequency, if available.

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

這一題的完整詳解

此題考查奈奎斯特穩定判準、頻率響應分析、增益裕度與相位裕度。
開迴路傳遞函數為 L(s)=K(s−z)s(s+p)L(s) = \frac{K(s-z)}{s(s+p)},其中 p,z,K>0p, z, K > 0。

(a) 繪製奈奎斯特圖。
我們分析 L(jω)L(j\omega) 的幅值和相角。
L(jω)=K(jω−z)jω(jω+p)=K(−z+jω)−ω2+jpωL(j\omega) = \frac{K(j\omega-z)}{j\omega(j\omega+p)} = \frac{K(-z+j\omega)}{-\omega^2+jp\omega}.
當 ω→0+\omega \to 0^+:
L(jω)≈K(−z)jpω=−Kzp1jω=−Kzp(−j/ω)=jKzpωL(j\omega) \approx \frac{K(-z)}{jp\omega} = \frac{-Kz}{p} \frac{1}{j\omega} = \frac{-Kz}{p}(-j/\omega) = j \frac{Kz}{p\omega}.
幅值 ∣L(jω)∣→∞|L(j\omega)| \to \infty。相角 ∠L(jω)→π2\angle L(j\omega) \to \frac{\pi}{2}。
軌跡從虛軸正向上方無限遠處開始。

當 ω→∞\omega \to \infty:
L(jω)≈K(jω)jω(jω)=KjωL(j\omega) \approx \frac{K(j\omega)}{j\omega(j\omega)} = \frac{K}{j\omega}.
幅值 ∣L(jω)∣→0|L(j\omega)| \to 0。相角 ∠L(jω)→−π2\angle L(j\omega) \to -\frac{\pi}{2}。
軌跡趨近於原點,從虛軸下方進入。

相位角 ∠L(jω)=arg⁡(K)+arg⁡(jω−z)−arg⁡(jω)−arg⁡(jω+p)\angle L(j\omega) = \arg(K) + \arg(j\omega-z) - \arg(j\omega) - \arg(j\omega+p)。
arg⁡(jω−z)=π−arctan⁡(ωz)\arg(j\omega-z) = \pi - \arctan(\frac{\omega}{z}) (因為 −z-z 是負實部,ω\omega 是正虛部,在第二象限)。
arg⁡(jω)=π2\arg(j\omega) = \frac{\pi}{2}。
arg⁡(jω+p)=arctan⁡(ωp)\arg(j\omega+p) = \arctan(\frac{\omega}{p}) (因為 pp 是正實部,ω\omega 是正虛部,在第一象限)。
∠L(jω)=0+(π−arctan⁡(ωz))−π2−arctan⁡(ωp)=π2−arctan⁡(ωz)−arctan⁡(ωp)\angle L(j\omega) = 0 + (\pi - \arctan(\frac{\omega}{z})) - \frac{\pi}{2} - \arctan(\frac{\omega}{p}) = \frac{\pi}{2} - \arctan(\frac{\omega}{z}) - \arctan(\frac{\omega}{p})。
當 ω=zp\omega = \sqrt{zp} 時,arctan⁡(zpz)+arctan⁡(zpp)=arctan⁡(pz)+arctan⁡(zp)=π2\arctan(\frac{\sqrt{zp}}{z}) + \arctan(\frac{\sqrt{zp}}{p}) = \arctan(\sqrt{\frac{p}{z}}) + \arctan(\sqrt{\frac{z}{p}}) = \frac{\pi}{2}。
所以,當 ω=zp\omega = \sqrt{zp} 時,∠L(jω)=π2−π2=0\angle L(j\omega) = \frac{\pi}{2} - \frac{\pi}{2} = 0。
此時軌跡與實軸相交。
實軸交點的幅值 ∣L(jzp)∣=K(z+p)p2+zp=K(z+p)p(p+z)=Kp|L(j\sqrt{zp})| = \frac{K(z+p)}{p^2+zp} = \frac{K(z+p)}{p(p+z)} = \frac{K}{p}。
所以,軌跡與實軸的交點是 (Kp,0)(\frac{K}{p}, 0)。

奈奎斯特圖草圖:軌跡從 ω=0\omega=0 的 (0,∞)(0, \infty) 開始,沿正虛軸上升,然後繞過實軸上的點 (Kp,0)(\frac{K}{p}, 0), 最後沿負虛軸趨近於原點 (0,0)(0,0)。由於 L(s)L(s) 的極點在 s=0s=0 和 s=−ps=-p,以及零點在 s=z>0s=z>0,這條軌跡是 ω∈[0,∞)\omega \in [0, \infty) 的部分。奈奎斯特圖還包括 ω∈(−∞,0]\omega \in (-\infty, 0] 的軌跡,這是 ω∈[0,∞)\omega \in [0, \infty) 軌跡的共軛。由於實部 Re(L(jω))=K(z+p)p2+ω2>0Re(L(j\omega)) = \frac{K(z+p)}{p^2+\omega^2} > 0, 軌跡始終在右半平面。

(b) 穩定區域 K>0K>0。
振盪頻率:不存在。
根據奈奎斯特穩定判準,N=P−ZN = P-Z。
L(s)L(s) 的極點為 s=0s=0 和 s=−ps=-p。在右半平面 (RHP) 沒有極點,所以 P=0P=0。
為了使閉迴路系統穩定,我們需要 Z=0Z=0。
N=P−Z=0−0=0N=P-Z = 0-0 = 0。
這表示奈奎斯特圖(考慮 ω∈(−∞,∞)\omega \in (-\infty, \infty))逆時針繞過 −1+j0-1+j0 點的次數必須為 0。
由於 Re(L(jω))>0Re(L(j\omega)) > 0, 軌跡始終在右半平面,永遠不會繞過 −1+j0-1+j0 點。
因此,N=0N=0 對於所有 K>0K>0 都成立。
所以,系統對於所有 K>0K>0 都是穩定的。
穩定區域為 K>0K>0。

🔒

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

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

免費註冊

其他考古題