リャプノフ関数
📂動力学リャプノフ関数
定義
空間 X と関数 f:X→X に対して次のような ベクトルフィールド が 微分方程式 によって与えられているとする。
x˙=f(x)
上記のような 自律システム の一点 x0∈X が与えられているとき、x0 の近傍 N(x0) で定義された スカラー関数 V∈C1(N(x0),R) が次の条件を満たすと リアプノフ関数Liapunov function と呼ばれる。
- (i): V(x0)=0 が成立し、x=x0 の場合は V(x)>0
- (ii): x∈N(x0)∖{x0} において V′(x)≤0
- C1(A,B) は定義域が A で、値域が B であり、微分可能でその導関数が連続な関数の集合を意味する。
- V が C1(N(x0),R) に属するということは、x0 近傍で定義されたスカラー関数であり微分可能でv′ が連続であることを意味する。
説明
リアプノフ関数は、与えられたシステム x˙=f(x) に依存して存在する可能性のある関数であり、特に X=Rn の場合に固定点 x0=x の 安定性 を確認するために考慮される。リアプノフ関数の存在性は即ち安定性を意味し、V はシステム f に応じて適切に定義されるべき関数である。v′ は文字通り時間 t に関する微分であり、V を微分すると x˙ に対する項が生じるため、システム f との関係が明らかになる。
この説明によると、一見リアプノフ関数が 非線形システム を理解する万能ツールのように思えるが、そもそも非線形システムが難しいため、このリアプノフ関数を見つけることは簡単なことではない。リアプノフ関数を一般的に見つける方法はなく、一つの重要なシステムに対してでもこのリアプノフ関数を見つけること自体が研究テーマになるほど難しいのである。
例
非常に簡単な例としてリアプノフ関数を見つける過程を見てみよう:
x˙=y˙=−x+4y−x−y3
このシステムには固定点 (0,0) が与えられている。言及したように、一般的にリアプノフ関数を見つける方法はないため、直感が必要となる。リアプノフ関数を見つけることに慣れれば、その関数を見つけるのも早くなるだろう。ここでは V(x,y)=x2+ay2 がリアプノフ関数であると仮定し、a≥0 の値を具体化する方法で見つけてみよう。
Part 1.
V(0,0)=0 であり、(x,y)=(0,0) の場合は V(x,y)>0 である。(0,0) で V=0 であり、a が負でないと仮定したので V>0 である。
Part 2.
(x,y)∈N((0,0))∖{(0,0)} において V′(x,y)≤0 である。V を t によって微分すると
dtdV===2xx˙+2ayy˙2x(−x+4y)+2ay(−x−y3)−2x2+(8−2a)xy−2ay2
ここで a=4 ならば V′<0 である。
したがって、固定点 (0,0) が与えられると、リアプノフ関数は V(x,y)=x2+4y2 として存在することが保証され、その存在性によって (0,0) が リアプノフ安定性 を持つこともわかる。