典型二階系統的階躍響應,包括有過衝 (overshoot)、振鈴 (ringing),在安定時間 (settling time)後回到穩定值。
階躍響應 是指系統在其輸入為單位階躍函數 時,其輸出的變化[ 1] 。在電子工程 或自動控制 領域中,階躍響應是指系統 的輸入在很短時間由0變成1時,其輸出的時域 特性。此概念可以延伸到使用抽象數學概念的動力系統 ,以演化參數表示其特性。
分析系統的階躍響應有助於了解系統的特性,因為當輸入在長時間穩態後,有快速而大幅度的變化,可以看出系統各個部份的特性。而且也可以知道系統的穩定性。
階躍響應的時域相關特性
系統的階躍響應可以用與以下時域 特性的量來描述:
過衝 (overshoot)。
上昇時間 (rise time)。
安定時間 (settling time)。
振鈴 (ringing)。
對於線性 動態系統來說,從這些特徵中可以推知許多該系統的資訊。
一階線性電路的階躍響應
RC電路
一階RC電路的階躍響應,沒有過衝及振鈴,在三倍時間常數時輸出到達輸入的95%
考慮如右圖的RC電路 ,頻域 下輸出電壓Vc 和輸入電壓Vin 的關係可表示為下式:
V
c
(
s
)
=
1
/
C
s
R
+
1
/
C
s
V
i
n
(
s
)
=
1
1
+
R
C
s
V
i
n
(
s
)
=
1
1
+
τ
s
V
i
n
(
s
)
{\displaystyle V_{c}(s)={\frac {1/Cs}{R+1/Cs}}V_{in}(s)={\frac {1}{1+RCs}}V_{in}(s)={\frac {1}{1+\tau s}}V_{in}(s)}
其中
τ
=
R
C
{\displaystyle \tau =RC}
為此系統的時間常數
考慮以下形式的輸入電壓Vin (t):
V
i
n
(
t
)
=
0
,
t
≤
0
{\displaystyle V_{in}(t)=0,t\leq 0}
V
i
n
(
t
)
=
V
i
n
,
t
>
0
{\displaystyle V_{in}(t)=V_{in},t>0}
則輸出電壓Vc (t)可以表示為以下的形式:
V
c
(
t
)
=
V
i
n
(
1
−
e
−
t
τ
)
{\displaystyle V_{c}(t)=V_{in}(1-e^{-{\frac {t}{\tau }}})}
反饋放大器的階躍響應
圖1:理想的負反饋模型;開迴路增益為 A OL ,反饋係數為 β。
本節介紹了一個簡單的負反饋放大器 的階躍響應如圖1所示。反饋放大器由一個增益為 A OL 的主開迴路放大器和反饋因子為 β 的反饋迴路。以下會分析此回授放大器,確認其階躍響應和控制響應的時間常數之間的關係,也看階躍響應和回授量之間的關係。
負反饋放大器增益為(見負反饋放大器 ):
A
F
B
=
A
O
L
1
+
β
A
O
L
{\displaystyle A_{FB}={\frac {A_{OL}}{1+\beta A_{OL}}}}
,
其中
A OL = 開迴路增益
A FB = 閉迴路增益(存在負反饋時的增益)
β = 反饋因子。
有一個主導極點
在許多情況下,可以用時間常數 τ 的單一主導極點很好地模擬順向放大器,它的開迴路增益為:
A
O
L
=
A
0
1
+
j
ω
τ
{\displaystyle A_{OL}={\frac {A_{0}}{1+j\omega \tau }}}
,
零頻率增益為 A 0 ,角頻率 ω = 2πf 。這種順向放大器有單位階躍響應
S
O
L
(
t
)
=
A
0
(
1
−
e
−
t
/
τ
)
{\displaystyle S_{OL}(t)=A_{0}(1-e^{-t/\tau })}
,
是從零到新平衡值 A 0 的指數趨近。
單極點放大器的傳遞函數導出閉迴路增益:
A
F
B
=
A
0
1
+
β
A
0
×
1
1
+
j
ω
τ
1
+
β
A
0
{\displaystyle A_{FB}={\frac {A_{0}}{1+\beta A_{0}}}\times {\frac {1}{1+j\omega {\frac {\tau }{1+\beta A_{0}}}}}}
此閉迴路增益與開迴路增益是形式相同:均為單極濾波器。其階躍響應的形式相同:一個趨於新平衡值的指數衰減。但閉迴路階躍函數的時間常數為 τ / (1 + β A 0 ),因此因此它比前向放大器的響應快,是其 1 + β A 0 倍:
S
F
B
(
t
)
=
A
0
1
+
β
A
0
(
1
−
e
−
t
(
1
+
β
A
0
)
/
τ
)
{\displaystyle S_{FB}(t)={\frac {A_{0}}{1+\beta A_{0}}}(1-e^{-t(1+\beta A_{0})/\tau })}
,
由於反饋因子 β 的增加,階躍響應會更快,直到最初假設的一個主導極點不再準確。如果有第二個極點,則隨着閉迴路時間常數區域第二個極點的時間常數,需要進行雙極點分析。
雙極點放大器
在開迴路增益有兩個極點的情況下(兩個時間常數 ,τ1 和τ2 ),階躍響應更為複雜。開迴路增益為:
A
O
L
=
A
0
(
1
+
j
ω
τ
1
)
(
1
+
j
ω
τ
2
)
,
{\displaystyle A_{OL}={\frac {A_{0}}{(1+j\omega \tau _{1})(1+j\omega \tau _{2})}},}
零頻率增益為 A 0 ,角頻率 ω = 2πf 。
分析
雙極點放大器的傳遞函數可以導出閉迴路增益:
A
F
B
=
A
0
1
+
β
A
0
×
1
1
+
j
ω
τ
1
+
τ
2
1
+
β
A
0
+
(
j
ω
)
2
τ
1
τ
2
1
+
β
A
0
{\displaystyle A_{FB}={\frac {A_{0}}{1+\beta A_{0}}}\times {\frac {1}{1+j\omega {\frac {\tau _{1}+\tau _{2}}{1+\beta A_{0}}}+(j\omega )^{2}{\frac {\tau _{1}\tau _{2}}{1+\beta A_{0}}}}}}
圖2:雙極點反饋放大器的共軛極點位置;Re (s) = 實軸,Im (s) = 虛軸。
放大器的時間相關性通過切換變量 s = j ω 很容易發現,於是增益變為:
A
F
B
=
A
0
τ
1
τ
2
×
1
s
2
+
s
(
1
τ
1
+
1
τ
2
)
+
1
+
β
A
0
τ
1
τ
2
{\displaystyle A_{FB}={\frac {A_{0}}{\tau _{1}\tau _{2}}}\times {\frac {1}{s^{2}+s\left({\frac {1}{\tau _{1}}}+{\frac {1}{\tau _{2}}}\right)+{\frac {1+\beta A_{0}}{\tau _{1}\tau _{2}}}}}}
這個表達式的極點(即分母的零點)位於:
2
s
=
−
(
1
τ
1
+
1
τ
2
)
{\displaystyle 2s=-\left({\frac {1}{\tau _{1}}}+{\frac {1}{\tau _{2}}}\right)}
±
(
1
τ
1
−
1
τ
2
)
2
−
4
β
A
0
τ
1
τ
2
,
{\displaystyle \pm {\sqrt {\left({\frac {1}{\tau _{1}}}-{\frac {1}{\tau _{2}}}\right)^{2}-{\frac {4\beta A_{0}}{\tau _{1}\tau _{2}}}}},}
這說明對於足夠大 βA 0 平方根就會變成虛數,極點的位置是共軛複數,s + 與 s − ;參見圖2:
s
±
=
−
ρ
±
j
μ
,
{\displaystyle s_{\pm }=-\rho \pm j\mu ,\,}
其中
ρ
=
1
2
(
1
τ
1
+
1
τ
2
)
,
{\displaystyle \rho ={\frac {1}{2}}\left({\frac {1}{\tau _{1}}}+{\frac {1}{\tau _{2}}}\right),}
而
μ
=
1
2
4
β
A
0
τ
1
τ
2
−
(
1
τ
1
−
1
τ
2
)
2
{\displaystyle \mu ={\frac {1}{2}}{\sqrt {{\frac {4\beta A_{0}}{\tau _{1}\tau _{2}}}-\left({\frac {1}{\tau _{1}}}-{\frac {1}{\tau _{2}}}\right)^{2}}}}
使用極座標系,根的半徑的模為 |s |(圖2):
|
s
|
=
|
s
±
|
=
ρ
2
+
μ
2
,
{\displaystyle |s|=|s_{\pm }|={\sqrt {\rho ^{2}+\mu ^{2}}},}
而角座標 φ 為:
cos
ϕ
=
ρ
|
s
|
,
sin
ϕ
=
μ
|
s
|
{\displaystyle \cos \phi ={\frac {\rho }{|s|}},\sin \phi ={\frac {\mu }{|s|}}}
拉普拉斯轉換 表告訴我們這樣一個系統的時間響應是由兩個函數的組合而成的:
e
−
ρ
t
sin
(
μ
t
)
and
{\displaystyle e^{-\rho t}\sin(\mu t)\quad {\text{and}}\quad }
e
−
ρ
t
cos
(
μ
t
)
,
{\displaystyle e^{-\rho t}\cos(\mu t),}
也就是說,解在時間上阻尼振盪。特別是,該系統的單位階躍響應為:[ 2]
S
(
t
)
=
(
A
0
1
+
β
A
0
)
(
1
−
e
−
ρ
t
sin
(
μ
t
+
ϕ
)
sin
(
ϕ
)
)
,
{\displaystyle S(t)=\left({\frac {A_{0}}{1+\beta A_{0}}}\right)\left(1-e^{-\rho t}\ {\frac {\sin \left(\mu t+\phi \right)}{\sin(\phi )}}\right)\ ,}
化簡為
S
(
t
)
=
1
−
e
−
ρ
t
sin
(
μ
t
+
ϕ
)
sin
(
ϕ
)
{\displaystyle S(t)=1-e^{-\rho t}\ {\frac {\sin \left(\mu t+\phi \right)}{\sin(\phi )}}\ }
當 A 0 趨於無窮大時,反饋係數 β 為1。
注意到響應的阻尼是由 ρ 決定的,也就是由開迴路放大器的時間常數確定的。相反,振盪的頻率是由 μ 確定的,也就是由, 通過 βA 0 由反饋參數確定的。因為 ρ 涉及到時間常數的倒數之和,所以可以發現 ρ 主要受到兩個時間常數中較短 的那個影響。
結果
圖3:線性雙極點反饋放大器的階躍響應;時間的單位是 1/ρ,也就是以 A OL 的時間常數為變量;曲線是用 mu = μ 的三個不同取值畫成的, μ 由 β 控制。
圖3顯示了參數 μ 取3個不同值時,單元階躍輸入的時間響應。可以看出隨着 μ 增加,振盪頻率也會增加,但振盪被包含在由指數型函數 [ 1 − exp (−ρt) ] 和 [ 1 + exp(−ρt) ] 確定的兩條漸近線以內。這兩條漸近線是由 ρ 決定的,所以也就由開迴路放大器的時間常數決定,與反饋無關。
終值的振盪現象被稱為振鈴 。過衝 是指擺動最大值高於終值,顯然會隨着 μ 增加。同樣,下沖是指擺動最小值低於終值,同樣也會隨着 μ 增加。安定時間 是指從終值出發,降到低於某個特定水平(終值的10%)用的時間。
穩定時間對 μ 的依賴性不明顯,而雙極點系統的近似可能不能達到用穩定時間對反饋的依賴性作出現實中的結論的準確性。但漸近線 [ 1 − exp (−ρt) ] 與 [ 1 + exp (−ρt) ] 顯然影響穩定時間,它們被開迴路放大器的時間常數控制,特別是在2個時間常數中的時間較短的。這表明開迴路放大器的設計必須滿足穩定時間的規定。
此分析的有兩個主要結論:
反饋控制了給定開迴路放大器並給定開迴路時間常數 τ1 與 τ2 時終值上下振盪的幅度。
開迴路放大器決定了穩定時間。它確定了圖3中的時標,開迴路放大器越快,時標越快。
順便說一句,可以注意到實際中與線性雙極點模型的偏離主要來自兩個方面:其一,實際放大器的極點多於兩個,零點也是;其二,實際放大器是非線性的,所以它們的階躍響應會隨着訊號幅度轉換。
圖4:α 的三個不同取值的階躍響應。頂部:α = 4;中間:α = 2;底部:α = 0.5。隨着 α 減小,極點分離也會減小,而過衝增加。
控制過衝
以下會說明如何用適當的參數選擇來減少過沖。
利用以上的公式,可以將階躍響應微分找最大值來計算過沖量。其過沖量最大值S max 為 :[ 3]
S
max
=
1
+
exp
(
−
π
ρ
μ
)
{\displaystyle S_{\max }=1+\exp \left(-\pi {\frac {\rho }{\mu }}\right)}
階躍響應的終值為1,因此其指數即為過沖量。可以看出若μ = 0,其過沖量為0,也就是:
4
β
A
0
τ
1
τ
2
=
(
1
τ
1
−
1
τ
2
)
2
{\displaystyle {\frac {4\beta A_{0}}{\tau _{1}\tau _{2}}}=\left({\frac {1}{\tau _{1}}}-{\frac {1}{\tau _{2}}}\right)^{2}}
令x = ( τ1 / τ2 )1 / 2 ,可以求解二個時間常數之間的比例,結果為
x
=
β
A
0
+
β
A
0
+
1
{\displaystyle x={\sqrt {\beta A_{0}}}+{\sqrt {\beta A_{0}+1}}\,}
因為β A 0 >> 1,因此平方根中的1可以省略,得到
τ
1
τ
2
=
4
β
A
0
{\displaystyle {\frac {\tau _{1}}{\tau _{2}}}=4\beta A_{0}}
換句話說,第一個時間常數需遠大於第二個時間常數。有時系統為了一些特性,需要允許一些過沖量,以下的關係中引入一個因子α:
τ
1
τ
2
=
α
β
A
0
,
{\displaystyle {\frac {\tau _{1}}{\tau _{2}}}=\alpha \beta A_{0},}
α可以依允許的過沖量來設計。
圖4就是描述其程式。比較上圖(α = 4)及下圖(α = 0.5)可以看出α較小,可以加快響應的速度,但也讓過沖量變大。中間的圖α = 2為幅度最平坦的濾波器 ,在波德圖 上沒有尖點。此設計有經驗法則 內建的安全預度,可以處理像重零點、重極點、非線性(例如和訊號振幅相依的特性)及製造的變異,這些都可能造成過大的過沖量。極點擺放位置(也就是α)的調整是頻率補償 的主題,其中一個方式是極點分離 。
穩定時間控制
圖3中階躍響應中振鈴的幅度是由阻尼因數 exp ( −ρ t ) 決定的。也就是說,如果我們指定出可接受的階躍響應離終值的偏移量 Δ,即:
S
(
t
)
≤
1
+
Δ
,
{\displaystyle S(t)\leq 1+\Delta ,\,}
在時間長於穩定時間 t S 這個前提下,無論 β A OL 的值為多少時這個條件都能滿足。穩定時間 t S 為:[ 4]
Δ
=
e
−
ρ
t
S
or
t
S
=
ln
(
1
Δ
)
ρ
=
τ
2
2
ln
(
1
Δ
)
1
+
τ
2
τ
1
≈
2
τ
2
ln
(
1
Δ
)
,
{\displaystyle \Delta =e^{-\rho t_{S}}{\text{ or }}t_{S}={\frac {\ln \left({\frac {1}{\Delta }}\right)}{\rho }}=\tau _{2}{\frac {2\ln \left({\frac {1}{\Delta }}\right)}{1+{\frac {\tau _{2}}{\tau _{1}}}}}\approx 2\tau _{2}\ln \left({\frac {1}{\Delta }}\right),}
因為過沖的條件,τ1 = αβA OL τ2 ,因此τ1 >> τ2 成立。一般穩定時間條件是指穩定時間和其單位增益的頻寬成反比的情形,原因是1/(2π τ2 )接近放大器在典型主極點補償 下的頻寬。不過此結果比經驗法則 的結果更準確。例如,若Δ = 1/e4 = 1.8 %,其穩定時間為t S = 8 τ2 。
一般而言,對過沖量的控制會決定二個時間常數的比例,穩定時間t S 會決定 τ2 [ 5] [ 6] [ 7] 。
相位裕度
圖5:波德圖可以找出相位裕度,其尺度是對數的,因此二個刻度之間的間隔是其比例,例如f 0 dB = βA 0 × f 1 。
其次,極點比例τ1 /τ2 也和回授放大器的相位裕度有關[ 8] 。圖5是二個極點放大器的波德圖 ,頻率到第二個極點的位置。圖5的假設是頻率f 0 dB 在位在f 1 = 1/(2πτ1 )的最小極點及位在f 2 = 1/(2πτ2 )的第二極點之間。如圖5所示,若 α ≥ 1,此假設即成立。
利用圖5,頻率(用f 0 dB 表示)為迴路增益 βA 0 滿足單位增益或是0 dB條件的位置,可以定義為:
|
β
A
OL
(
f
0 db
)
|
=
1
{\displaystyle |\beta A_{\text{OL}}(f_{\text{0 db}})|=1\ }
波德增益圖中增益下降的斜率是 20 dB/decade,頻率每增加十倍,增益下降的比例相同:
f
0 dB
=
β
A
0
f
1
{\displaystyle f_{\text{0 dB}}=\beta A_{0}f_{1}\,}
相位裕度是頻率在f 0 dB 處,相位和−180°之間的距離,因此裕度為:
ϕ
m
=
180
∘
−
arctan
(
f
0 dB
/
f
1
)
−
arctan
(
f
0 dB
/
f
2
)
{\displaystyle \phi _{m}=180^{\circ }-\arctan(f_{\text{0 dB}}/f_{1})-\arctan(f_{\text{0 dB}}/f_{2})\ }
因為f 0 dB / f 1 = βA 0 >> 1,有關f 1 的項為90°,因此相位裕度為:
ϕ
m
=
90
∘
−
arctan
(
f
0 dB
/
f
2
)
{\displaystyle \phi _{m}=90^{\circ }-\arctan(f_{\text{0 dB}}/f_{2})\,}
=
90
∘
−
arctan
(
β
A
0
f
1
α
β
A
0
f
1
)
{\displaystyle =90^{\circ }-\arctan \left({\frac {\beta A_{0}f_{1}}{\alpha \beta A_{0}f_{1}}}\right)}
=
90
∘
−
arctan
(
1
α
)
=
arctan
(
α
)
{\displaystyle =90^{\circ }-\arctan \left({\frac {1}{\alpha }}\right)=\arctan \left(\alpha \right)}
若α = 1,則 φm = 45°,若α = 2,則φm = 63.4°. Sansen[ 9] 建議α = 3,對應的φm = 71.6°「是一種很好的啟始條件。」
若τ2 縮短,α會增加。穩定時間t S 也會減小。若τ1 變大,α也會增加,穩定時間會略有變動。若用到極點分離 技巧,τ1 和 τ2 都會變化。
若放大器有二個以上的極點,圖5的波德圖仍然可以計算相位裕度,只要將f 2 視為「等效的第二極點」位置即可[ 10] 。
數學定義
本節以抽象概念下的動態系統
S
{\displaystyle \textstyle {\mathfrak {S}}}
,來提供形式性的數學定義。所有符號及假設列在下方。
t
∈
T
{\displaystyle \textstyle t\in T}
是系統的演化參數 ,為方便說明,簡稱為時間 。
x
|
t
∈
M
{\displaystyle \textstyle {\boldsymbol {x}}|_{t}\in M}
是系統在時間
t
{\displaystyle t\,}
時的狀態 ,為方便說明,稱為輸出。
Φ
:
T
×
M
⟶
M
{\displaystyle \textstyle \Phi :T\times M\longrightarrow M}
是動態系統的演化函數 。
Φ
(
0
,
x
)
=
x
0
∈
M
{\displaystyle \textstyle \Phi (0,{\boldsymbol {x}})={\boldsymbol {x}}_{0}\in M}
是動態系統的初始狀態 。
H
(
t
)
{\displaystyle \textstyle H(t)\,}
是單位階躍函數 。
非線性動態系統
若針對一個一般的動態系統,其階躍響應可定義如下:
x
|
t
=
Φ
{
H
(
t
)
}
(
t
,
x
0
)
{\displaystyle {\boldsymbol {x}}|_{t}=\Phi _{\{H(t)\}\left(t,{{\boldsymbol {x}}_{0}}\right)}\,}
其階躍響應是系統輸入為單位階躍函數時的演化函數。表示式中H (t )為下標。
線性動態系統
對於一個線性非時變系統 ,令
S
≡
S
{\displaystyle \textstyle {\mathfrak {S}}\ \equiv \ S}
,其階躍響應可以用單位階躍函數
H
(
t
)
{\displaystyle \textstyle H(t)}
和系統衝激響應
h
(
t
)
{\displaystyle \textstyle h(t)}
的卷積 來表示:
a
(
t
)
=
h
∗
H
(
t
)
=
H
∗
h
(
t
)
=
∫
−
∞
+
∞
h
(
τ
)
H
(
t
−
τ
)
d
τ
=
∫
−
∞
t
h
(
τ
)
d
τ
{\displaystyle a(t)={h*H}(t)={H*h}(t)=\int _{-\infty }^{+\infty }h(\tau )H(t-\tau )\,d\tau =\int _{-\infty }^{t}h(\tau )\,d\tau }
對線性非時變系統而言就是將後面的式子積分。相對的,對於線性非時變系統,階躍響應的微分即為衝激響應:
h
(
t
)
=
d
d
t
a
(
t
)
{\displaystyle h(t)={\frac {d}{dt}}\,a(t)}
不過此關係在非線性系統或是時變系統 並不成立[ 1] 。
對於一個線性非時變系統 ,其階躍響應可以用單位階躍函數
H
(
t
)
{\displaystyle \textstyle H(t)}
和系統衝激響應
h
(
t
)
{\displaystyle \textstyle h(t)}
的卷積 來表示:
a
(
t
)
=
h
∗
H
(
t
)
=
H
∗
h
(
t
)
=
∫
−
∞
+
∞
h
(
τ
)
H
(
t
−
τ
)
d
τ
=
∫
−
∞
t
h
(
τ
)
d
τ
{\displaystyle a(t)={h*H}(t)={H*h}(t)=\int _{-\infty }^{+\infty }h(\tau )H(t-\tau )\,d\tau =\int _{-\infty }^{t}h(\tau )\,d\tau }
參照
參考文獻與註釋
延伸閱讀
Robert I. Demrow Settling time of operational amplifiers [1] (頁面存檔備份 ,存於互聯網檔案館 )
Cezmi Kayabasi Settling time measurement techniques achieving high precision at high speeds [2] (頁面存檔備份 ,存於互聯網檔案館 )
Vladimir Igorevic Arnol'd "Ordinary differential equations", various editions from MIT Press and from Springer Verlag, chapter 1 "Fundamental concepts"
外部連結