射影力学系のソースを表示

[[数学]]における'''射影力学系'''（しゃえいりきがくけい、{{Lang-en-short|projected dynamical system}}）とは、解がある制約集合に制限された[[力学系]]の挙動を調べる数学理論である。この学問では、[[数理最適化|最適化]]や[[平衡点]]の問題などの静的な分野と、[[常微分方程式]]の動的な分野との関連や応用が示されている。

射影力学系は、次の'''射影微分方程式'''（projected differential equation）の[[フロー (数学)|フロー]]として与えられる：

:<math>
\frac{dx(t)}{dt} = \Pi_K(x(t),-F(x(t))).
</math>

ここで ''K'' は制約集合である。この形状の微分方程式は、不連続なベクトル場を持つという点において注目すべきものである。

== 射影力学系の歴史 ==

射影力学系は、特に時間などのあるパラメータに関する平衡点問題における非静的な解の挙動を力学的にモデルするという狙いの下で発展してきた。この系の挙動は常微分方程式のそれと異なり、例えば[[金融]]モデルにおける投資の非負性や、[[オペレーションズ・リサーチ]]における[[凸集合|凸]][[多面体]]集合など、根底にある平衡点問題が有効となるようなある制約集合に制限された解を扱うものである。射影力学系の発展による恩恵を受けた特に重要な平衡点問題のクラスの一つに、{{仮リンク|変分不等式|en|variational inequality}}が挙げられる。

射影力学系の理論の構築は1990年代に行われた。しかし、それ以前の数学の文献、特に変分不等式と微分包含式との関連に関する文献にも、同様の概念は見つけられている。

== 射影と錐 ==

射影微分方程式の任意の解は、すべての時間に対して制約集合 ''K'' の内側に留まる必要がある。このようなことは、射影作用素と二つの特に重要な[[凸錐]]のクラスを用いることで、達成される。ここで ''K'' はある[[ヒルベルト空間]] ''X'' の[[閉集合|閉]][[凸集合|凸]]部分集合である。

集合 ''K'' の点 ''x'' での'''法錐'''（normal cone）は、次で与えられる。

:<math>
N_K(x) = \{ p \in V | \langle p, x - x^* \rangle \geq 0, \forall x^* \in K \}.
</math>

集合 ''K'' の点 ''x'' での{{仮リンク|接錐|en|tangent cone}}は、次で与えられる。

:<math>
T_K(x) = \overline{\bigcup_{h>0} \frac{1}{h} (K-x)}.
</math>

''X'' の点 ''x'' での、''K'' への'''射影作用素'''（projection operator）は、''K'' 内の点 <math>P_K(x)</math> で

:<math>
\| x-P_K(x) \| \leq \| x-y \|
</math>

をすべての ''K'' 内の ''y'' に対して満たすものとして与えられる。

''X'' 内のあるベクトル ''v'' の、''K'' の点 ''x'' での'''ベクトル射影作用素'''は、次で与えられる。

:<math>
\Pi_K(x,v)=\lim_{\delta \to 0^+} \frac{P_K(x+\delta v)-x}{\delta}.
</math>

== 射影微分方程式 ==

ヒルベルト空間 ''X'' の閉凸部分集合 ''K'' と、''K'' から ''X'' への元を取るベクトル場 ''-F'' が与えられたとき、''K'' と ''-F'' に関する射影微分方程式は次で定義される。

:<math>
\frac{dx(t)}{dt} = \Pi_K(x(t),-F(x(t))).
</math>

''K'' の[[内部 (位相空間論)|内部]]で、解は制限のない常微分方程式におけるものと同じように振る舞う。しかし、ベクトル場はその集合の境界に沿って不連続なので、射影部分方程式は不連続な微分方程式のクラスに属する。この事実によって、常微分方程式の多くの理論が適用できなくなるが、''-F'' が[[リプシッツ連続]]なベクトル場であるときは、''K'' 内の各初期点 ''x(0)=x<sub>0</sub>'' を通る[[絶対連続]]な解が区間 <math>[0,\infty)</math> 上で唯一つ存在する。

この微分方程式は、代替的に

:<math>
\frac{dx(t)}{dt} = P_{T_K(x(t))}(-F(x(t)))
</math>

や

:<math>
\frac{dx(t)}{dt} = -F(x(t))-P_{N_K(x(t))}(-F(x(t)))
</math>

で特徴付けることも出来る。負号を使ってベクトル場を ''-F'' と記述する慣習は、射影力学系が変分不等式と共有する特定の関係性によるものである。文献においてそのような慣習は、変分不等式においてはベクトル場は正で、対応する射影力学系においては負であると述べられている。

== 関連項目 ==
* {{仮リンク|微分変分不等式|en|Differential variational inequality}}
* [[力学系|力学系理論]]
* [[常微分方程式]]
* {{仮リンク|変分不等式|en|Variational inequality}}
* [[微分包含式]]
* {{仮リンク|相補性理論|en|complementarity theory}}

== 参考文献 ==

* Aubin, J.P. and Cellina, A., ''Differential Inclusions'', Springer-Verlag, Berlin (1984).
* Nagurney, A. and Zhang, D., ''Projected Dynamical Systems and Variational Inequalities with Applications'', Kluwer Academic Publishers (1996). 
* Cojocaru, M., and  Jonker L., ''Existence of solutions to projected differential equations on Hilbert spaces'', Proc. Amer. Math. Soc., 132(1), 183-193 (2004).
* Brogliato, B., and Daniilidis, A., and [[:en:Claude Lemaréchal|Lemaréchal, C.]], and Acary, V., "On the equivalence between complementarity systems, projected systems and differential inclusions", ''Systems and Control Letters'', vol.55, pp.45-51 (2006)

{{DEFAULTSORT:しやえいりきかくけい}}
[[Category:微分方程式]]
[[Category:力学系]]
[[Category:数学に関する記事]]