歴史 リャプノフ安定性は、1892年にハリコフ大学(現在の VNカラジン・ハリコフ国立大学 ) で「運動の安定性の一般問題」という 論文を提出したロシアの数学者アレクサンドル・ミハイロヴィチ・リャプノフにちなんで名付けられました。 [ 1 ] A.M.リャプノフは、広く普及していた平衡点の周りで線形化する局所的方法と比較して、非線形力学系の安定性の解析に対するグローバルなアプローチを開発することに成功した先駆者でした。彼の研究は、最初はロシア語で発表され、その後フランス語に翻訳されましたが、長年にわたりほとんど注目されませんでした。A.M.リャプノフによって確立された運動の安定性の数学理論は、科学技術へのその応用の時代をかなり先取りしていました。さらに、リャプノフ自身はこの分野に応用せず、彼自身の関心は天文学的応用を伴う回転流体質量の安定性にありました。安定性の分野で研究を続ける博士課程の学生は彼にはおらず、1918年に自殺したため、彼自身の運命はひどく悲劇的だった。 数十年間、安定性の理論は完全に忘れ去られていた。1930年代にカザン航空研究所で働いていたロシア・ソビエトの数学者で機械学者のニコライ・グルイェヴィチ・チェタエフは、 AM リャプノフが行った発見の驚くべき大きさを最初に認識した人物だった。NG チェタエフ[ 3 ] による理論への貢献は非常に重要で、多くの数学者、物理学者、エンジニアは彼をリャプノフの直接の後継者であり、安定性の数学理論の創造と発展における次の科学的後継者とみなしている。
冷戦 期には、いわゆる「リアプノフの第2の方法」(下記参照)が、他の方法では扱えない強い非線形性を含む航空宇宙誘導システム の安定性に適用できることがわかったため、この分野への関心が急激に高まりました。当時から現在に至るまで、制御およびシステム分野の文献に多数の論文が発表されています。[ 4 ] [ 5 ] [ 6 ] [ 7 ] [ 8 ] 最近では、リアプノフ指数 (安定性を議論するリアプノフの第1の方法に関連) の概念が、カオス理論 に関連して広く関心を集めています。リアプノフ安定性法は、交通割り当て問題における平衡解を見つけるためにも適用されています。[ 9 ]
離散時間システムの定義 離散時間 システムの定義は、連続時間システムの定義とほぼ同じです。以下では、より数学的な文献でよく用いられる別の表現を用いて、この定義を示します。
( X , d )を距離空間 とし、f : X → X を 連続関数 とする。X内の点x は、 以下 の条件を満たす場合にリアプノフ安定で あると言われる。
∀ ϵ > 0 ∃ δ > 0 ∀ y ∈ X [ d ( x 、 y ) < δ ⇒ ∀ n ∈ N d ( f n ( x ) 、 f n ( y ) ) < ϵ ] 。 {\displaystyle \forall \epsilon >0\ \exists \delta >0\ \forall y\in X\ \left[d(x,y)<\delta \Rightarrow \forall n\in \mathbf {N} \ d\left(f^{n}(x),f^{n}(y)\right)<\epsilon \right].} xが 漸近的に安定で あるとは、 x がその安定集合 の内部に属する場合、すなわち 、
∃ δ > 0 [ d ( x 、 y ) < δ ⇒ リム n → ∞ d ( f n ( x ) 、 f n ( y ) ) = 0 ] 。 {\displaystyle \exists \delta >0\left[d(x,y)<\delta \Rightarrow \lim _{n\to \infty }d\left(f^{n}(x),f^{n}(y)\right)=0\right].}
線形状態空間モデルの安定性 線形状態空間 モデル
x ˙ = A x {\displaystyle {\dot {\textbf {x}}}=A{\textbf {x}}} 、どこA {\displaystyle A} は有限行列であり、の固有値 のすべての実部が であれば漸近安定(実際には指数安定)である。 A {\displaystyle A} 負の値をとる。この条件は、次の条件と同等である。[ 13 ]
A T M + M A {\displaystyle A^{\textsf {T}}M+MA} ある正定値 行列に対して負定値である。M = M T {\displaystyle M=M^{\textsf {T}}} (関連するリアプノフ関数はV ( x ) = x T M x {\displaystyle V(x)=x^{\textsf {T}}Mx} )
それに対応して、時間離散線形状態空間 モデル
x t + 1 = A x t {\displaystyle {\textbf {x}}_{t+1}=A{\textbf {x}}_{t}} は、すべての固有値がA {\displaystyle A} 絶対値 が1より小さい。
この後者の条件は、スイッチドシステムに一般化されています。線形スイッチド離散時間システム(行列の集合によって支配される){ A 1 、 … 、 A m } {\displaystyle \{A_{1},\dots ,A_{m}\}} )
x t + 1 = A 私 t x t 、 A 私 t ∈ { A 1 、 … 、 A m } {\displaystyle {{\textbf {x}}_{t+1}}=A_{i_{t}}{\textbf {x}}_{t},\quad A_{i_{t}}\in \{A_{1},\dots ,A_{m}\}} 集合の結合スペクトル半径が 漸近的に安定(実際には指数的に安定)である場合{ A 1 、 … 、 A m } {\displaystyle \{A_{1},\dots ,A_{m}\}} 1より小さい。
入力(または制御)を持つシステムは次の形式をとる。
x ˙ = f ( x 、 u ) {\displaystyle {\dot {\textbf {x}}}={\textbf {f}}({\textbf {x}},{\textbf {u}})} ここで、(一般的に時間依存の)入力 u(t) は、制御 、外部入力 、 刺激 、外乱 、または強制関数と見なすことができます。 [ 14 ] は、リアプノフ安定な平衡点の近傍では、小さな外乱の下でもシステムが安定していることが示されています。より大きな入力外乱の場合、このようなシステムの研究は制御理論の対象であり、 制御工学 に応用されています。入力のあるシステムの場合、システムの安定性に対する入力の影響を定量化する必要があります。この分析の主な 2 つのアプローチは、BIBO 安定性 (線形システム の場合)と入力状態安定性 (ISS)(非線形システム の場合)です。
無制限の摂動を伴う時間変動システムのリアプノフ安定性 システムについて考えてみましょう 。ζ ˙ ( t ) = g ( t 、 ζ ( t ) ) ; t ≥ t 0 、 {\displaystyle {\dot {\zeta }}(t)=g{\big (}t,\zeta (t){\big )};t\geq t_{0},} ζ ( t 0 ) = ζ 0 、 {\displaystyle \zeta (t_{0})=\zeta _{0},} どこ( t 0 、 ζ 0 ) ∈ R × R m {\displaystyle (t_{0},\zeta _{0})\in \mathbb {R} \times \mathbb {R} ^{m}} 、 解決ζ ( t ) {\displaystyle \zeta (t)} でR m {\displaystyle \mathbb {R} ^{m}} (m {\displaystyle m} は厳密に正の整数であり、かつ明確に定義された関数である。g : [ t 0 、 ∞ ) × R m → R m {\displaystyle g:[t_{0},\infty )\times \mathbb {R} ^{m}\rightarrow \mathbb {R} ^{m}} とg ( t 、 0 ) = 0 {\displaystyle g(t,0)=0} 、∀ t ≥ t 0 {\displaystyle \forall t\geq t_{0}} 。
システムがカラテオドリ条件を満たすと仮定します。つまり、マッピングはt ↦ g ( t 、 ζ ) {\displaystyle t\mapsto g(t,\zeta )} 局所的に本質的に有界である[ t 0 、 ∞ ) × R m {\displaystyle [t_{0},\infty )\times \mathbb {R} ^{m}} は、すべてのζ ∈ R m {\displaystyle \zeta \in \mathbb {R} ^{m}} そしてほぼすべてのt ≥ t 0 {\displaystyle t\geq t_{0}} このシステムは、最大区間上で定義される局所的に絶対連続な局所カラテオドリ解を許容する。[ t 0 、 ω ) {\displaystyle [t_{0},\omega )} 。
定数が存在すると仮定するα > 0 {\displaystyle \alpha >0} 、β > 0 {\displaystyle \beta >0} とα < β {\displaystyle \alpha <\beta } 局所的に絶対連続な関数 r 1 ∈ C 0 ( R 、 R ) {\displaystyle r_{1}\in C^{0}(\mathbb {R} ,\mathbb {R} )} 、r 2 ∈ C 0 ( R 、 R ) {\displaystyle r_{2}\in C^{0}(\mathbb {R} ,\mathbb {R} )} そしてルベーグ可測関数h : R → R {\displaystyle h:\mathbb {R} \rightarrow \mathbb {R} } 以下の条件を満たすこと:
(私) ( − 1 ) α = − 1 {\displaystyle (-1)^{\alpha }=-1} そして( − 1 ) β {\displaystyle (-1)^{\beta }} 明確に定義されています(( − 1 ) β {\displaystyle (-1)^{\beta }} 解の正値が保証されている場合は、明確に定義される必要がある。
(ii)r 1 ( t ) > 0 {\displaystyle r_{1}(t)>0} 、r 2 ( t ) > 0 {\displaystyle r_{2}(t)>0} 、∀ t > t 0 {\displaystyle \forall t>t_{0}} そしてh ( t ) > 0 {\displaystyle h(t)>0} ほぼすべてのt > t 0 {\displaystyle t>t_{0}} 。
(iii)リム t → ∞ r 2 ( t ) r 1 ( t ) = ∞ {\displaystyle \lim _{t\rightarrow \infty }{\frac {r_{2}(t)}{r_{1}(t)}}=\infty } (これから、摂動は無限大であることが導かれる。これは後述するリアプノフの不等式によって示される。)
(iv) リム t → ∞ Λ ( t ) = 0 {\displaystyle \lim _{t\rightarrow \infty }\Lambda (t)=0} 、 どこ Λ ( t ) := r 1 ( t ) r ˙ 2 ( t ) − r ˙ 1 ( t ) r 2 ( t ) h ( t ) ( r 1 ( t ) ) 2 β − α − 1 β − α ( r 2 ( t ) ) β − 2 α + 1 β − α 、 ほぼすべての t > t 0 。 {\displaystyle \Lambda (t):={\frac {r_{1}(t){\dot {r}}_{2}(t)-{\dot {r}}_{1}(t)r_{2}(t)}{h(t)(r_{1}(t))^{\frac {2\beta -\alpha -1}{\beta -\alpha }}(r_{2}(t))^{\frac {\beta -2\alpha +1}{\beta -\alpha }}}},{\mbox{ for almost all }}t>t_{0}.}
(v)各解についてζ ( t ) {\displaystyle \zeta (t)} 最大存在間隔を持つシステムの[ t 0 、 ω ) {\displaystyle [t_{0},\omega )} 正の定数が存在するδ {\displaystyle \delta } 、σ {\displaystyle \sigma } 、c 1 {\displaystyle c_{1}} 、c 2 {\displaystyle c_{2}} 、そしてリアプノフ関数 V ∈ C 1 ( R × R m 、 R + ) {\displaystyle V\in C^{1}(\mathbb {R} \times \mathbb {R} ^{m},\mathbb {R} _{+})} 満足(ほぼすべてのt ∈ ( t 0 、 ω ) {\displaystyle t\in (t_{0},\omega )} 満たすV ( t 、 ζ ( t ) ) < δ {\displaystyle V(t,\zeta (t))<\delta } ):
c 1 | κ | σ ≤ V ( t 、 κ ) ≤ c 2 | κ | σ 、 ∀ t ∈ R 、 ∀ κ ∈ R m 、 {\displaystyle c_{1}|\kappa |^{\sigma }\leq V(t,\kappa )\leq c_{2}|\kappa |^{\sigma },\forall t\in \mathbb {R} ,\forall \kappa \in \mathbb {R} ^{m},}
∇ V ( t 、 ζ ( t ) ) ⋅ g ( t 、 ζ ( t ) ) ≤ ( − r 1 ( t ) V α ( t 、 ζ ( t ) ) + r 2 ( t ) V β ( t 、 ζ ( t ) ) ) h ( t ) 。 {\displaystyle \nabla V(t,\zeta (t))\cdot g(t,\zeta (t))\leq {\big (}-r_{1}(t)V^{\alpha }(t,\zeta (t))+r_{2}(t)V^{\beta }(t,\zeta (t)){\big )}h(t).}
そして、存在するc 3 > 0 {\displaystyle c_{3}>0} 各| ζ 0 | < c 3 {\displaystyle |\zeta _{0}|<c_{3}} 1つ取得ω = ∞ {\displaystyle \omega =\infty } そして| ζ ( t ) | ≤ c 2 c 1 σ | ζ 0 | 、 ∀ t ≥ t 0 、 {\displaystyle |\zeta (t)|\leq {\sqrt[{\sigma }]{\frac {c_{2}}{c_{1}}}}|\zeta _{0}|,\forall t\geq t_{0},} となることによってζ = 0 {\displaystyle \zeta =0} は一様に安定である。さらに、原点は漸近的に安定である。 [ 25 ]
例 アリー効果を伴う時間変動人口増加モデルは、 [ 25 ]で表すことができる。
N ˙ ( t ) = R ( t ) N ( t ) ( N ( t ) A ( t ) − 1 ) ( 1 − N ( t ) K ( t ) ) 、 {\displaystyle {\dot {N}}(t)=R(t)N(t)\left({\frac {N(t)}{A(t)}}-1\right)\left(1-{\frac {N(t)}{K(t)}}\right),} どこt ≥ t 0 {\displaystyle t\geq t_{0}} 、 州N ( t ) ∈ R {\displaystyle N(t)\in \mathbb {R} } ルベーグ可測関数 R : R → R {\displaystyle R:\mathbb {R} \rightarrow \mathbb {R} } とR ( t ) > 0 {\displaystyle R(t)>0} ほぼすべてのt > t 0 {\displaystyle t>t_{0}} 、局所的に絶対連続な関数 A ∈ C 0 ( R 、 R ) {\displaystyle A\in C^{0}(\mathbb {R} ,\mathbb {R} )} そして K ∈ C 0 ( R 、 R ) {\displaystyle K\in C^{0}(\mathbb {R} ,\mathbb {R} )} となることによって A ( t ) > 0 {\displaystyle A(t)>0} そしてK ( t ) > 0 {\displaystyle K(t)>0} すべての人々のためにt > t 0 {\displaystyle t>t_{0}} 。
等式の右辺は局所的にリプシッツ連続である。N {\displaystyle N} したがって、最大存在区間を持つ一意の解が存在する。[ t 0 、 ω ) {\displaystyle [t_{0},\omega )} 起源N = 0 {\displaystyle N=0} は平衡点です。私たちは、N = 0 {\displaystyle N=0} 一様安定かつ漸近安定な絶滅平衡。この目的のために、以下を仮定する。
∫ t 0 ∞ R ( t ) d t = ∞ 、 {\displaystyle \int _{t_{0}}^{\infty }R(t)dt=\infty ,} リム t → ∞ ( 1 A ( t ) + 1 K ( t ) ) = ∞ 、 {\displaystyle \lim _{t\rightarrow \infty }\left({\frac {1}{A(t)}}+{\frac {1}{K(t)}}\right)=\infty ,} そして リム t → ∞ ( A ˙ ( t ) ( K ( t ) ) 2 + K ˙ ( t ) ( A ( t ) ) 2 R ( t ) A ( t ) K ( t ) ( A ( t ) + K ( t ) ) ) = 0. {\displaystyle \lim _{t\rightarrow \infty }\left({\frac {{\dot {A}}(t)(K(t))^{2}+{\dot {K}}(t)(A(t))^{2}}{R(t)A(t)K(t)(A(t)+K(t))}}\right)=0.} させてV = N 2 {\displaystyle V=N^{2}} 不等式は満たされるc 1 = 1 {\displaystyle c_{1}=1} 、c 2 = 1 {\displaystyle c_{2}=1} そしてσ = 2 {\displaystyle \sigma =2} 簡単な計算で
V ˙ ( t ) ≤ 2 R ( t ) ( − V ( t ) + ( 1 A ( t ) + 1 K ( t ) ) V 3 2 ( t ) ) ほぼすべての t > t 0 、 {\displaystyle {\dot {V}}(t)\leq 2R(t)\left(-V(t)+\left({\frac {1}{A(t)}}+{\frac {1}{K(t)}}\right)V^{\frac {3}{2}}(t)\right){\mbox{ for almost all }}t>t_{0},} これは、不等式の形をとる。α = 1 {\displaystyle \alpha =1} 、β = 3 2 {\displaystyle \beta ={\frac {3}{2}}} 、h ( t ) = 2 R ( t ) {\displaystyle h(t)=2R(t)} 、r 1 ( t ) = 1 {\displaystyle r_{1}(t)=1} 、r 2 ( t ) = 1 A ( t ) + 1 K ( t ) {\displaystyle r_{2}(t)={\frac {1}{A(t)}}+{\frac {1}{K(t)}}} そしてδ {\displaystyle \delta } は任意です( t 0 、 ∞ ) {\displaystyle (t_{0},\infty )} さらに、関数を簡単に示すことができる。Λ ( t ) {\displaystyle \Lambda (t)} ゼロになるt {\displaystyle t} 無限大になる。したがって、すべての条件が満たされ、したがって存在する。c 3 > 0 {\displaystyle c_{3}>0} 各| N 0 | < c 3 {\displaystyle |N_{0}|<c_{3}} 1つ取得ω = ∞ {\displaystyle \omega =\infty } そして
| N ( t ) | ≤ c 2 c 1 σ | N 0 | 、 ∀ t ≥ t 0 、 {\displaystyle |N(t)|\leq {\sqrt[{\sigma }]{\frac {c_{2}}{c_{1}}}}|N_{0}|,\forall t\geq t_{0},} そしてそれゆえにN = 0 {\displaystyle N=0} は均一に安定している。実際、N = 0 {\displaystyle N=0} 一様安定であり、漸近安定である。
参考文献 1 2 リャプノフ、AM 運動の安定性の一般問題 (ロシア語)、博士論文、ハリコフ大学、1892年 英語訳:(1)運動の安定性 、アカデミックプレス、ニューヨーク&ロンドン、1966年 (2)運動の安定性の一般問題 、(AT フラー訳)、テイラー&フランシス、ロンドン、1992年。スミルノフによる伝記とリャプノフの著作の広範な参考文献が含まれています。↑ Chetaev、NG 力学の安定した軌道について、Kazan Univ Sci Notes、vol.4 no.1 1936; 『運動の安定性』、元々は ОГИЗ によって 1946 年にロシア語で出版されました。 Гос。 изд-во технико-теорет。 лит.、Москва-Ленинград.モートン・ナドラー訳、オックスフォード、1961年、200ページ。 ↑ レトフ、午前 (1955 年)。Устойчивость нелинейных регулируемых систем [ 非線形制御システムの安定性 ] (ロシア語)。モスクワ: ゴステヒズダット。 英語訳、プリンストン大学出版局、1961年↑ Kalman, RE ; Bertram, J. F (1960). "Lyapunovの「第2の方法」による制御システムの解析と設計:I—連続時間システム". Journal of Basic Engineering . 82 (2): 371– 393. doi : 10.1115/1.3662604 . ↑ LaSalle, JP ; Lefschetz, S. (1961). Stability by Lyapunov's Second Method with Applications . New York: Academic Press. ↑ Parks, PC (1962). 「自動制御理論におけるリアプノフの方法」. Control . I Nov 1962 II Dec 1962. ↑ Kalman, RE (1963). "自動制御におけるLur'e問題に対するLyapunov関数" . Proc Natl Acad Sci USA . 49 (2): 201–205 . Bibcode : 1963PNAS...49..201K . doi : 10.1073 / pnas.49.2.201 . PMC 299777. PMID 16591048 . ↑ Smith, MJ; Wisten, MB (1995). "連続的な日々の交通量割り当てモデルと連続的な動的ユーザー均衡の存在". Annals of Operations Research . 60 (1): 59– 79. doi : 10.1007/BF02031940 . S2CID 14034490 . ↑ ハーン、ヴォルフガング ( 1967)。 運動の安定性 。 シュプリンガー。pp. 191–194 、 セクション 40。doi : 10.1007/978-3-642-50085-5。ISBN 978-3-642-50087-9 。↑ Braun, Philipp; Grune, Lars; Kellett, Christopher M. (2021). (In-)Stability of Differential Inclusions: Notions, Equivalences, and Lyapunov-like Characterizations . Springer. pp. 19–20 , Example 2.18. doi : 10.1007/978-3-030-76317-6 . ISBN 978-3-030-76316-9 . S2CID 237964551 . ↑ Vinograd, RE (1957). 「非線形微分方程式の研究における特性指数法の不適切性」 . Doklady Akademii Nauk (ロシア語). 114 (2): 239– 240. ↑ Goh, BS (1977). "多種システムにおける地球規模の安定性". The American Naturalist . 111 (977): 135–143 . Bibcode : 1977ANat..111..135G . doi : 10.1086/283144 . S2CID 84826590 . ↑ マルキン IG 運動の安定性理論、モスクワ 1952 (ゴステヒズダット) 第 II 章 第 4 節 (ロシア語) 英語訳、言語サービス局、ワシントン AEC -tr-3352; 原題 On stability under constant acting disturbances Prikl Mat 1944, vol. 8 no.3 241-245 (ロシア語); Amer. Math. Soc. transl. no. 8 1 2 スロットイン、ジャン=ジャック E.リー・ウェイピン (1991)。 非線形制御を適用 。ニュージャージー州: プレンティス・ホール。 ↑ I. Barbălat、Systèmes d'équations différentielles d'oscillations non Linéaires、Rev. Math。 Pures Appl. 4 (1959) 267–270、p. 269. ↑ B. Farkas 他、「Barbălat の補題の変奏」、Amer. Math. Monthly (2016) 128、no. 8、825-830、DOI: 10.4169/amer.math.monthly.123.8.825、p. 827。 ↑ B. Farkas 他、「Barbălat の補題の変奏」、Amer. Math. Monthly (2016) 128、no. 8、825-830、DOI: 10.4169/amer.math.monthly.123.8.825、p. 826。 ↑ Naser, MFM (2020).ある種の時間変動微分方程式の状態収束 . IMA Journal of Mathematical Control and Information, 37 , 27–38. ↑ Naser, MFM & Ikhouane, F. (2019). Stability of time–-varying systems in the absence of strict Lyapunov functions . IMA Journal of Mathematical Control and Information . DOI: 10.1093/imamci/dnx056. ↑ Jiang, Z., Lin, Y. & Wang, Y. (2009). 非線形時変システムの安定化:制御リアプノフ関数アプローチ . Journal of Systems Science and Complexity , 22 , 683–696. ↑ Malisoff, M. & Mazenc, F. (2005). 時変システムに対する厳密な入力状態安定リアプノフ関数に関する追加考察 . Automatica , 41 , 1973–1978. ↑ Mazenc, F. (2003). 時間変動システムのための厳密なリアプノフ関数 . Automatica , 39 , 349–353. ↑ Mu, X. & Cheng, D. (2005). 時変非線形制御システムの安定性と安定化について . Asian Journal of Control , 1 2 Naser, MFM、「非自律システムの無限摂動を伴う解の時間無限近傍での挙動 」、IMA Journal of Mathematical Control & Information、39、pp. 60–79、定理6.2、2022年。
さらに読む Bhatia, Nam Parshad; Szegő, Giorgio P. (2002).力学系の安定性理論 . Springer. ISBN 978-3-540-42748-3 。 チェルビン、ロバート(1971)。2流体プラズマシステムのリアプノフ安定性とフィードバック制御 (博士論文)。コロンビア大学。 ガンドルフォ、ジャンカルロ(1996)。経済動学 (第3 版)。ベルリン:シュプリンガー。407-428頁 。ISBN 978-3-540-60988-9 。 Shcherbakov, Pavel S. (1992), 「アレクサンドル・ミハイロヴィチ・リャプノフ:運動の安定性に関する博士論文発表100周年について」、Automatica 、28 (5): 865–871 、doi : 10.1016/0005-1098(92)90140-B Parks, PC (1992). "AM Lyapunovの安定性理論 ― 100年後". IMA Journal of Mathematical Control & Information . 9 (4): 275–303 . doi : 10.1093/imamci/9.4.275 . スロットイン、ジャン=ジャック・E.リー・ウェイピン (1991)。非線形制御を適用 。ニュージャージー州: プレンティス・ホール。 Teschl, G. (2012).常微分方程式と力学系 .プロビデンス :アメリカ数学会 . ISBN 978-0-8218-8328-0 。ウィギンズ、S. (2003).応用非線形力学系とカオス入門 (第2 版). ニューヨーク:シュプリンガー・フェルラーク . ISBN 978-0-387-00177-7 。
この記事は、 PlanetMath の asymptotically stable の資料を組み込んでおり、Creative Commons Attribution/Share-Alike Licenseの下でライセンスされています。