線形力学系

提供: testwiki
2021年11月12日 (金) 08:33時点におけるimported>Nnhによる版 (曖昧さ回避ページスカラーへのリンクを解消、リンク先をスカラー (物理学)に変更(DisamAssist使用))
(差分) ← 古い版 | 最新版 (差分) | 新しい版 → (差分)
ナビゲーションに移動 検索に移動

線形力学系(せんけいりきがくけい、テンプレート:Lang-en-short)とは、行列で定義され、線形性を持つ力学系である。

定義

一般に テンプレート:Math における線形力学系は、ベクトル値関数 テンプレート:Math と、テンプレート:Mvar 次の正方行列 テンプレート:Mvar により、次のような微分方程式で表される。

ddt𝐱(t)=A𝐱(t)

ただしこれは、テンプレート:Math連続的に変化する場合であり、離散系の場合には、

𝐱m+1=A𝐱m

で表される。

これが線形であるとは、テンプレート:Mathテンプレート:Math が解ならば、任意のスカラー テンプレート:Mvar について、線形結合 テンプレート:Math も解である、ということを意味している。

線形力学系は、多くの非線形の場合と異なり、完全に解くことができる。このとき、解は行列の指数 テンプレート:Math(連続系)、もしくは累乗 テンプレート:Mvar(離散系)によって表現され、その振る舞いは一般的に行列 テンプレート:Mvar固有値固有ベクトルによって理解できる。 非線形のときでも、変数変換により線型化して解くことができることもある。また、不動点の周りでの線形近似は、非線形系を理解するのに役立つ(ハートマン=グロブマンの定理)。

線形力学系の解

初期値 テンプレート:Math が、行列 テンプレート:Mvar固有ベクトル テンプレート:Math ならば、初期条件は

ddt𝐱(t)|t=0=A𝐯k=λk𝐯k

となる。ただし、テンプレート:Mvar は、固有ベクトル テンプレート:Math に対応する固有値である。このとき、解は、

𝐱(t)=𝐯keλkt

となる。

もし テンプレート:Mvar対角化可能ならば、任意の初期値 テンプレート:Math は、固有ベクトルの線形結合で一意に表される。つまり、次のような係数 テンプレート:Mvar が一意に存在する。

𝐱0=k=1nak𝐯k

このとき解は、

𝐱(t)=k=1nak𝐯keλkt

となる。

対角化不可能な場合でも一般に行列の指数関数を用いて

𝐱(t)=etA𝐱0(etA=n=0tnn!An)

と、解を導くことができる。

二次元の場合

二次元の線形力学系は、

ddt(xy)=A(xy)

で表される。この系では、テンプレート:Mvarテンプレート:Math正方行列である。テンプレート:Mvar の固有値は、行列式 テンプレート:Math と、トレース テンプレート:Mvar を用いて、

λ1=τ+τ24Δ2
λ2=ττ24Δ2

のように書くことができる。

また、Δ=λ1λ2 であり、τ=λ1+λ2 である。

もし、Δ<0 ならば、固有値の符号が異なり原点は、鞍点 (テンプレート:En) となる。

Δ=0 ならば、原点は孤立した平衡点ではない。

Δ>0 ならば、固有値の符号が同じになり、τ<0 ならば(漸近)安定、τ=0 ならば中立安定τ>0 ならば不安定になる。また固有値が実数ならば節点 (テンプレート:En) となる。ただし、二つの固有値が同じときには対角化可能なときスター、不可能なとき退化節点 (テンプレート:En) となる。最後に複素数のときは、渦状点 (テンプレート:En) となる。

参考文献

関連項目