制御工学I 第10回
安定性
ラウス,フルビッツの安定判別
2010/6/21 2
授業の予定
• 制御工学概論(1回) – 制御技術は現在様々な工学分野において重要な基本技術となっている。工学における制御工学 の位置づけと歴史について説明する。さらに、制御システムの基本構成と種類を紹介する。 • ラプラス変換(1回) – 制御工学、特に古典制御ではラプラス変換が重要な役割を果たしている。ラプラス変換と逆ラプラ ス変換の定義を紹介し、微分方程式のラプラス変換について解説する。 • 制御システムのモデリングと伝達関数(3回) – システムの相似性について概説し、システムの入出力特性を表す手法である伝達関数について詳 述する。システムの図的表現であるブロック線図とその等価変換について解説する。 • 過渡特性(3回) – システムの過渡状態を評価する方法であるインパルス応答とインディシャル応答について解説する。 システムの速応性や安定性の指標である整定時間、立ち上がり量、行き過ぎ量について述べる。 • 安定性(2回) – システムの安定性の概念を述べ、安定性を判定する代数的方法であるラウス-フルビッツの方法に ついて説明する。 • 周波数特性(4回) – 周波数領域におけるシステムの特性を周波数特性という。周波数特性と伝達関数との関係を説明 し、ベクトル軌跡とボード線図の作成方法を説明する。リアプノフ安定性
•
(微分量が0)とする状態
→平衡状態
– 平衡点• 平衡点に摂動を与えた系の解
の時間経過
を考える。ただし,
は初期値。
• 定義
ならば,全てのt>=t0に対して
となるような実数aが存在すれば,
平衡点
は安定。存在しなければ不安定。
– 平衡点 が安定かつ, のとき は漸近安定。 0 = x& xe e x(
t, x0,t0)
x 0 0, t xδ
≤ − x x0(
t x t)
x a x , 0, 0 − e ≤ e x e xlim
(
, 0, 0)
− = 0 ∞ → e t x t x t x xe2010/6/21 4
漸近安定
e x 0 xδ
a 不安定 安定 漸近安定内部安定性
• 零入力応答
– 入力が零で,初期条件のみに対する応答
– 初期状態
• 内部安定(漸近安定:リアプノフの安定性)
– どのような初期状態x0に対しても,
t→∞でx(t)=0となる
• システムが内部安定となる必要十分条件
– 特性方程式の根の実部が全て負⇔漸近安定
( )
0 x0 x =( )
kt e x t x = 02010/6/21 6
外部安定性
• 零状態応答
– 初期状態が零で,入力のみに対する応答• 有界な入力u(t)に対する零状態応答
– 有界な入力とは,下式に適当な正数k1が存在する• システムが外部安定・入出力安定・有界入力有界
出力(BIBO:Bounded Input Bounded Output)安定
– 有界な出力とは,下式において適当な正数k2が存在する – 有界な入力u(t)に対してk2が存在する場合BIBO安定
( )
t ≤ k < ∞ ≤ t < ∞ u 1 ,0( )
t ≤ k < ∞ ≤ t < ∞ y 2 ,0BIBO安定
• システムのBIBO安定と等価な条件
– 伝達関数の全ての極(特性方程式の根)が負の
実部を持つ。 → 漸近安定
– インパルス応答g(t)に対して,次式を満たす正数
kが存在
• 有界入力u(t)に対する出力の応答 任意のtに対して有界( )
≤ < ∞∫
∞ g t dt k 0( )
(
) ( )
(
) ( )
(
−)
≤ < ∞ ≤ − ≤ − =∫
∫
∫
k k d t g k d u t g d u t g t y t t t 1 0 1 0 0τ
τ
τ
τ
τ
τ
τ
τ
2010/6/21 8
BIBO安定
• インパルス応答g(t)が有界でない場合
– 有界な入力u(τ)
– 出力y(t)は有界ではない
• 変数変換が( )
(
)
( )
(
)
⎩ ⎨ ⎧ < − − = ≥ − = 0 1 0 1τ
τ
τ
τ
t g for u t g for u( )
t =∫
t g(
t −) ( )
u d =∫
t g(
t −)
d y 0 0τ
τ
τ
τ
τ
τ
σ
= t − dσ
= −dτ
( )
=∫
( )
− =∫
( )
− t t g d g d t y 0 0σ
σ
σ
σ
出力が有界でないのでBIBO安定が成立しない制御システムの安定性判別方法
(特性方程式の根を調べる方法)
• 特性方程式の特性根を求めて調べる
– 特性方程式を解く
– 根軌跡による方法
• 特性方程式の係数を用いて調べる
– ラウス-フルビッツの方法
• ベクトル軌跡による方法
– ナイキストの方法
• ボード線図を用いた方法
– ゲイン余裕,位相余裕
2010/6/21 10
特性方程式の係数を用いた
安定判別法
• 伝達関数の特性方程式の根の実部
– 高次の代数方程式の求解は困難
– 根の実部の正負判別で十分
– 代数方程式を解かないで安定判別
• 特性方程式の解(因数分解)
– pi:実数根,σi:複素根実部,ωi:複素根虚部
(
)
∏
{
(
)
}
∏
= − = − − + + + = + + + + = μ μω
σ
1 2 2 2 1 1 1 1 0 0 i i i n j j n n n n s p s s a s a s a a L特性方程式の係数を用いた
安定判別法
• 全ての根の実部が負
• s
kの係数は全て正
– 全ての根の実部が負となる必要条件
(
)
{
(
)
}
(
)
∏
{
}
∏
∏
∏
= − = = − = + + + + = + + + = μ μ μ μω
σ
σ
ω
σ
1 2 2 2 2 1 1 2 2 2 1 2 0 i i i i n j j i i i n j j s s p s s p s(
1L 2μ
)
,σ
0(
1Lμ
)
0 = − > = > j n i pj i 0 , 0 2 , 0 > 2 + 2 > > i i i j pσ
σ
ω
(
0 1)
0 = − > k n ak L 十分条件は?2010/6/21 12
ラウスの安定判別法
• 特性多項式の係数からラウス表を作成
– 係数を変換
– 変換した係数の変換
n n n s s a s a a0 + 1 +L+ −1 −1 + ⎪ ⎪ ⎩ ⎪ ⎪ ⎨ ⎧ = = = ⎪ ⎪ ⎩ ⎪ ⎪ ⎨ ⎧ = = = − − − − − M M 5 1 2 3 1 1 1 1 1 4 0 2 2 0 1 0 0 1 n n n n n a a a a aα
α
α
α
α
α
1 1 1 1 1 0 1 1 0 1 1 0 2 + + + + + + + + + + = − = + k i k k i k k i k k i k k iα
γ
α
α
α
α
α
α
α
ただし 1 , 00 1 0 0 1 = + = +α
α
α
γ
k kkラウスの安定判別法
• ラウス表
– 安定判別法(根の実部が全て負となる必要十分条件) n n n n m m m m m m m m n n n n s s s s s s s 0 1 0 2 1 2 0 3 2 3 2 3 1 3 0 2 1 2 2 2 2 2 1 2 0 1 1 1 2 1 2 1 1 1 0 0 0 1 0 2 0 2 0 1 0 0 0 1 2 3 2 1α
α
α
α
α
α
α
α
α
α
α
α
α
α
α
α
α
α
α
α
α
α
α
α
− − − − − − − − − − − − − M M L L L L M(
k n)
k > 0 =1,2,L,γ
(
i n)
i , , 2 , 1 0 0 > = Lα
または2010/6/21 14
フルビッツの安定判別法
• 特性多項式の係数からフルビッツ行列を作成
– nxn行列
– 要素an-1よりはじめる
• 一つ右に移る毎に 係数を2減らす • 一つ下に移る毎に 係数を1増やすす • an=1 • ak=0 k>n or k<0 ⎥⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎦ ⎤ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎣ ⎡ = − − − − − − − − − 0 1 2 3 1 4 2 5 3 1 0 0 0 0 0 1 0 0 0 1 0 0 0 1 0 a a a a a a a a a a H n n n n n n n n n M L L L Lフルビッツの安定判別法
• フルビッツ行列
⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎦ ⎤ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎣ ⎡ = − − − − − − − − − − − − − − − − 0 2 4 1 6 2 8 3 1 2 0 2 3 1 1 4 2 0 5 3 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 a a a a a a a a a a a a a a a a a a a a a a a a a H n n n n n n n n n n n n n n n n n n n n L M M L L M M L L L L L L L2010/6/21 16
フルビッツの安定判別法
• フルビッツ行列式の作成
– 左上kxk小行列の行列式⊿k(k=1,2,・・・,n)
(主座小行列式)
• 安定判別法
– 特性方程式の根の実部が全て負となる
H det = Δn ⎥ ⎦ ⎤ ⎢ ⎣ ⎡ = Δ − − − 2 3 1 2 1 n n n a a a ⎥ ⎥ ⎥ ⎦ ⎤ ⎢ ⎢ ⎢ ⎣ ⎡ = Δ − − − − − − − 3 1 4 2 5 3 1 3 0 1 n n n n n n n a a a a a a a 1 1 = − Δ an 0 > Δkラウス法とフルビッツ法の等価性
• フルビッツ行列のラウス配列要素表示
⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎦ ⎤ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎣ ⎡ = ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎦ ⎤ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎣ ⎡ = − − − − − − − − − 0 0 0 0 0 0 0 0 0 0 1 0 1 0 1 0 0 1 0 0 1 0 0 1 1 1 0 0 2 0 1 0 0 1 2 1 1 1 0 1 2 3 1 4 2 5 3 1 a a a a a a a a a a a a a a a a a a a a a H n n n n n n n n n L L L L L L L L2010/6/21 18
ラウス法とフルビッツ法の等価性
• 係数行列をつくる
⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎦ ⎤ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎣ ⎡ − − − = Γ 1 0 0 0 0 0 1 0 0 1 1 0 0 0 1 0 1 1 1 1 1 Oγ
γ
γ
ラウス法とフルビッツ法の等価性
• フルビッツ行列に係数行列をかける
⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎦ ⎤ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎣ ⎡ = ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎦ ⎤ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎣ ⎡ − − − − − − = Γ = 0 0 0 0 0 0 0 0 0 0 0 1 0 2 0 1 1 1 0 2 1 2 0 1 2 1 1 1 0 1 0 1 0 0 1 0 1 1 1 0 1 1 0 1 1 0 1 2 1 0 1 2 1 0 2 1 1 1 1 1 1 0 1 0 0 1 2 1 1 1 0 1 1 a a a a a a a a a a a a a a a a a a a a a a a a a a a H H L L L L L L L L γ γ γ γ γ γ2010/6/21 20
ラウス法とフルビッツ法の等価性
• 係数行列2をつくる
⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎦ ⎤ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎣ ⎡ − − = Γ 1 0 0 0 0 0 1 1 0 1 0 0 1 0 0 1 2 2 2 Oγ
γ
ラウス法とフルビッツ法の等価性
• フルビッツ行列に係数行列をかける
⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎦ ⎤ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎣ ⎡ = ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎦ ⎤ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎣ ⎡ − = Γ = 0 0 0 0 0 0 0 0 0 0 0 0 2 0 3 0 2 1 2 0 1 2 1 1 1 0 2 0 2 1 2 1 1 2 1 2 0 1 2 1 1 1 0 1 2 2 M L L L L M L L L L a a a a a a a a a a a a a a a H H γ2010/6/21 22
ラウス法とフルビッツ法の等価性
• 同様にして
⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎦ ⎤ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎣ ⎡ − − = Γ 1 0 1 0 1 0 1 0 1 0 0 1 0 0 0 1 O O k k k γ γ k-1個 1 − Γ = k k k H Hラウス法とフルビッツ法の等価性
• 変形フルビッツ行列
–
の左上kxk小行列の行列式(主座小行列式)
⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎦ ⎤ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎣ ⎡ = − − − − n n n n n a a a a a a a a a a a a a a a H 0 4 3 4 0 3 3 3 1 3 0 2 2 2 2 2 1 2 0 1 1 1 3 1 2 1 1 1 0 0 ~ M L L L L L H H H~ = Γn n−1 = ΓnΓn−1LΓ1 n k a a a k k , 1,2, , ~ 0 2 0 1 0 L = L = Δ H~2010/6/21 24
ラウス法とフルビッツ法の等価性
• 係数行列について
–
は下三角行列,対角要素は全て1
–
Γ~n = ΓnΓn−1LΓ1も対角要素が全て1の下三角行列
k Γ(
)
(
)
(
)
(
)
{
}
(
)
(
)
(
)
k k k k H k k H k k k k H k k k k H k k H Δ = × = × × Γ = × × × Γ = × Γ = × = Δ 小行列 の左上 小行列 の左上 小行列 の左上 小行列 の左上 小行列 の左上 小行列 の左上 小行列 の左上 1 det ~ det ~ det ~ det ~ det ~ n k a a a k k 0 , 1,2, , 2 0 1 0 L = L = Δ n k a0k > 0, =1,2,L, Δk > 0,k =1,2,L,n ラウスの安定判別法 フルビッツの安定判別法根軌跡法
• システムの過渡応答は極の位置(値)に関係
する
• 制御ゲインにより極の場所(値)が変わる
– どのように変化するのか注目
• 設計
– 適切なゲインの設定により,極を所期の場所に
配置
– システムの極や零の追加の影響の定性的な評
価
2010/6/21 26
根軌跡
• ゲインKを持つフィードバック制御システム
– 閉ループ伝達関数
• sの多項式で表す • ac+kbd=0→根は伝達関数の極 • bc=0 →根は伝達関数の零点 • K:0→∞で伝達関数の極(特性方程式の根)は変化 K G(s) + -R(s) C(s) H(s)( )
( )
KG( ) ( )
s( )
H s s KG s R s C + = 1( ) ( )
( ) ( )
( )
( )
s c s d s H s a s b s G = , =( )
( )
( )
( )
( )
( )
( )
( )
( )
( )
( ) ( )
( ) ( )
( ) ( )
a( ) ( )
s c s( ) ( )
Kb( ) ( )
s d s s c s b K s d s Kb s c s a s c s a s a s b K s c s d s a s b K s a s b K s R s C + = + = + = 1 根軌跡 根のs平面上の変化根軌跡の例
• 二次遅れシステム
– 特性方程式
• K=0 • 0≦K≦a2/4 • a2/4<K • 固有角周波数 • 減衰定数( )
( )
s as K K s R s C + + = 2 0 2 + + = K as s a s s1 = 0, 2 = − 2 4 2 , 2 2 1 K a a s s = − ± − 2 4 2 , 2 2 1 a K j a s s = − ± − K n = ω K a 2 = ζ -a 0 -a/2 x x2010/6/21 28
根軌跡の性質1
• 伝達関数
– d(s):n次,c(s):m次, m<n• 始点 K=0:
→伝達関数の極に対応
• 終点 K=∞:
→m個は零点に対応
– 極は零点より多い(m<n)ので残りは無限遠点に終わる• 実軸対称 特性方程式の根は実数または共役複素
数なので,根の配置は実軸対称
• 枝数 特性方程式の根はn個→根軌跡の枝数n本
( )
s = 0 d( ) ( )
1( )
( )
0 1+ = + = s d s c K s H s KG( )
s = 0 Kc2010/6/21 29
根軌跡
• s(s+2)(s+5)+K=0の根軌跡
– 漸近線
• 極の数n=3,零の数m=0,(r=n-m=3)– 無限遠点に向かう軌跡の漸近線
• 漸近線の実軸上の交点座標=(0-2-5)/3=-7/3 • 漸近線の方向角=π/3,π,5/3π K G(s) + -R(s) C(s) H(s)( )
( ) (
s s)(
s)
K K s R s C + + + = 5 2( ) ( )( ) ( )
, 1 5 2 1 = + + = H s s s s s G(
) (
)
r sa = 極の座標の総和 − 零点の座標の総和(
)
1 , , 2 , 1 , 0 , 1 2 − = + = k π k L r θ 0x
x
x
-2 -5●
-7/32010/6/21 30