§E13.16相平面

最終更新

線形系

x′=(1221)xx'= \begin{pmatrix}1&2\\2&1\end{pmatrix} x

の解のうち、直線x1+x2=0x_1+x_2=0上の解は原点へ近づくが、その直線からわずかに外れた解は原点から離れる。したがって、一つの特別な解が平衡点へ収束することだけでは、その平衡点の安定性を判断することができない。そのため、平衡点の安定性を判断するには、平衡点の近くのすべての初期状態から出る解の振る舞いを調べる必要がある。

相平面では、状態を平面上の点として表し、各点での変化の向きと解の軌道を同時に見る。鞍点では、原点へ近づく解の軌道はすべて一本の直線の上にある。平衡点の近くに始まる解が近くに留まるか、さらに平衡点へ収束するかという安定性の概念を用いると、軌道の長時間挙動を明確に区別することができる。双曲型平衡点の安定性は、Jacobi 行列の固有値の実部の符号から判定することができる。

本記事では、相線と相平面を用いて平衡点の安定性を調べ、二次元線形系の分類と、線形化による非線形系の平衡点の安定性の判定を扱う。

1 相線と相平面

定義 1.1.U⊆RU\subseteq\Rを開集合、f:U→Rf:U\to\Rを局所 Lipschitz 連続な関数とし、一階自律方程式x′=f(x)x'=f(x)を考える。UUの各点xxにf(x)f(x)の符号を記したものを、ffの相線 (phase line) という。

  1. ffの流れ、軌道、平衡点は§E13.15 定義 4.1のn=1n=1の場合である。

補題 1.2.U⊆RU\subseteq\Rを開集合、f:U→Rf:U\to\Rを局所 Lipschitz 連続な関数、φ\varphiをffの流れとする。x∈Ux\in Uはf(x)≠0f(x)\ne0を満たすとし、xxを含みffの零点を含まないUUの部分区間すべての和集合をCCとする。このとき、次が成り立つ。

  1. CCは開区間であり、ffはCC上でf(x)f(x)と同じ符号をとる。さらにγ(x)⊆C\gamma(x)\subseteq Cである。
  2. f(x)>0f(x)>0ならばt↦φt(x)t\mapsto\varphi_t(x)はI(x)I(x)上で狭義単調増加であり、f(x)<0f(x)<0ならばI(x)I(x)上で狭義単調減少である。
  3. コンパクト集合K⊆UK\subseteq Uが存在して、すべてのt∈I(x)t\in I(x)、t≥0t\ge0についてφt(x)∈K\varphi_t(x)\in Kであるならば、[0,∞)⊆I(x)[0,\infty)\subseteq I(x)であり、極限ℓ=lim⁡t→∞φt(x)\ell=\lim_{t\to\infty}\varphi_t(x)が存在し、ℓ∈K\ell\in Kはffの平衡点である。

証明.主張 (1)を示す。xxを含みffの零点を含まないUUの部分区間はすべて点xxを共有するから、その和集合CCは区間である。y∈Cy\in Cとすると、f(y)≠0f(y)\ne0であり、UUは開集合でffは連続であるから、ε>0\varepsilon>0が存在して、(y−ε,y+ε)⊆U(y-\varepsilon,y+\varepsilon)\subseteq Uかつこの区間上でf≠0f\ne0である。この区間と、xxとyyを含みffの零点を含まないUUの部分区間との和は、xxを含みffの零点を含まないUUの部分区間であるから、CCに含まれる。したがってCCは開区間である。a,b∈Ca,b\in Cについてf(a)f(a)とf(b)f(b)の符号が異なるならば、aaとbbを端点とする閉区間はCCに含まれ、§D1.12 系 1.2によりその閉区間にffの零点が存在する。これはCCがffの零点を含まないことに反するから、ffはCC上でf(x)f(x)と同じ符号をとる。

t∈I(x)t\in I(x)とし、y=φt(x)y=\varphi_t(x)とおく。f(y)=0f(y)=0ならば、§E13.15 命題 4.3 (2)によりγ(y)={y}\gamma(y)=\{y\}であり、y∈γ(x)y\in\gamma(x)であるから§E13.15 命題 4.3 (1)によりγ(x)=γ(y)={y}\gamma(x)=\gamma(y)=\{y\}である。このときx=yx=yであり、これはf(x)≠0f(x)\ne0に反する。したがってγ(x)\gamma(x)はffの零点を含まない。t1<t2t_1<t_2をI(x)I(x)の点とすると、§D1.12 系 1.2を[t1,t2][t_1,t_2]上の連続関数s↦φs(x)s\mapsto\varphi_s(x)に適用することにより、φt1(x)\varphi_{t_1}(x)とφt2(x)\varphi_{t_2}(x)の間の値はすべてγ(x)\gamma(x)に属する。したがってγ(x)\gamma(x)はxxを含みffの零点を含まないUUの部分区間であり、γ(x)⊆C\gamma(x)\subseteq Cである。

主張 (2)を示す。t1<t2t_1<t_2をI(x)I(x)の点とする。s↦φs(x)s\mapsto\varphi_s(x)はI(x)I(x)上で微分可能であり、その導関数はs↦f(φs(x))s\mapsto f(\varphi_s(x))である。§D1.14 定理 3.1により、c∈(t1,t2)c\in(t_1,t_2)が存在してφt2(x)−φt1(x)=f(φc(x))(t2−t1)\varphi_{t_2}(x)-\varphi_{t_1}(x)=f(\varphi_c(x))(t_2-t_1)である。主張 (1)によりφc(x)∈C\varphi_c(x)\in Cであり、f(φc(x))f(\varphi_c(x))はf(x)f(x)と同じ符号をもつ。

主張 (3)を示す。§E13.15 定理 5.3 (2)により、[0,∞)⊆I(x)[0,\infty)\subseteq I(x)であり、ω(x)\omega(x)は空でなく、ω(x)⊆K\omega(x)\subseteq Kである。f(x)>0f(x)>0ならばσ=1\sigma=1、f(x)<0f(x)<0ならばσ=−1\sigma=-1とおくと、主張 (2)によりs↦σφs(x)s\mapsto\sigma\varphi_s(x)はI(x)I(x)上で狭義単調増加である。ℓ∈ω(x)\ell\in\omega(x)と、φtm(x)→ℓ\varphi_{t_m}(x)\to\ellを満たす数列tm→∞t_m\to\inftyを取る。s≥0s\ge0について、tm>st_m>sならばσφtm(x)>σφs(x)\sigma\varphi_{t_m}(x)>\sigma\varphi_s(x)であるから、m→∞m\to\inftyとしてσℓ≥σφs(x)\sigma\ell\ge\sigma\varphi_s(x)である。ε>0\varepsilon>0に対してmmをσℓ−σφtm(x)<ε\sigma\ell-\sigma\varphi_{t_m}(x)<\varepsilonを満たすように取ると、t≥tmt\ge t_mについてσφtm(x)≤σφt(x)≤σℓ\sigma\varphi_{t_m}(x)\le\sigma\varphi_t(x)\le\sigma\ellであるから、∣φt(x)−ℓ∣<ε\lvert\varphi_t(x)-\ell\rvert<\varepsilonである。したがってφt(x)→ℓ\varphi_t(x)\to\ell(t→∞)(t\to\infty)である。ω(x)\omega(x)の任意の元はある数列tm′→∞t'_m\to\inftyに沿うφtm′(x)\varphi_{t'_m}(x)の極限であるからℓ\ellに等しく、ω(x)={ℓ}\omega(x)=\{\ell\}である。§E13.15 定理 5.3 (1)によりω(x)\omega(x)は不変集合であるからγ(ℓ)⊆{ℓ}\gamma(\ell)\subseteq\{\ell\}であり、§E13.15 命題 4.3 (2)によりf(ℓ)=0f(\ell)=0である。ℓ∈ω(x)⊆K\ell\in\omega(x)\subseteq Kである。▨

定義 1.3.U⊆R2U\subseteq\R^2を開集合、F:U→R2F:U\to\R^2を局所 Lipschitz 連続な写像とし、自律系x′=F(x)x'=F(x)を考える。状態x=(x1,x2)x=(x_1,x_2)を点として表す平面を相平面 (phase plane) という。各点x∈Ux\in Uに速度ベクトルF(x)F(x)を対応させたものを、FFの定めるベクトル場 (vector field) という。

  1. FFの流れ、軌道、平衡点は§E13.15 定義 4.1のn=2n=2の場合であり、相平面にはx∈Ux\in Uを通る軌道γ(x)\gamma(x)をttが増加する向きを付けて描く。

2 平衡点の安定性

定義 2.1.n∈N≥1n\in\NNとし、U⊆RnU\subseteq\R^nを開集合、F:U→RnF:U\to\R^nを局所 Lipschitz 連続な写像とする。

  1. FFの平衡点は§E13.15 定義 4.1の平衡点である。
  2. x∗∈Ux_*\in UをFFの平衡点とする。x∗x_*の開近傍W⊆UW\subseteq Uで、WWに含まれるFFの平衡点がx∗x_*だけであるものが存在するとき、x∗x_*をFFの孤立平衡点 (isolated equilibrium point) という。

定義 2.2.n∈N≥1n\in\NNとし、U⊆RnU\subseteq\R^nを開集合、F:U→RnF:U\to\R^nを局所 Lipschitz 連続な写像、φ\varphiをFFの流れ、x∗∈Ux_*\in UをFFの平衡点とする。Rn\R^nのノルム∥⋅∥\lVert\cdot\rVertを一つ固定する。

x∗x_*がLyapunov 安定 (Lyapunov stable) であるとは、任意のε>0\varepsilon>0に対してδ>0\delta>0が存在し、∥x−x∗∥<δ\lVert x-x_*\rVert<\deltaを満たす任意のx∈Ux\in Uについて、[0,∞)⊆I(x)[0,\infty)\subseteq I(x)であり、すべてのt≥0t\ge0で∥φt(x)−x∗∥<ε\lVert\varphi_t(x)-x_*\rVert<\varepsilonが成り立つことをいう。

x∗x_*が漸近安定 (asymptotically stable) であるとは、x∗x_*が Lyapunov 安定であり、さらにρ>0\rho>0が存在して、∥x−x∗∥<ρ\lVert x-x_*\rVert<\rhoを満たす任意のx∈Ux\in Uについて、[0,∞)⊆I(x)[0,\infty)\subseteq I(x)かつφt(x)→x∗\varphi_t(x)\to x_*(t→∞)(t\to\infty)が成り立つことをいう。Lyapunov 安定でない平衡点を不安定 (unstable) という。

補題 2.3.n∈N≥1n\in\NNとし、U⊆RnU\subseteq\R^nを開集合、F:U→RnF:U\to\R^nを局所 Lipschitz 連続な写像、x∗∈Ux_*\in UをFFの平衡点とする。∥⋅∥1\lVert\cdot\rVert_1、∥⋅∥2\lVert\cdot\rVert_2をRn\R^nの二つのノルムとする。このとき、x∗x_*が∥⋅∥1\lVert\cdot\rVert_1について Lyapunov 安定であることは、∥⋅∥2\lVert\cdot\rVert_2について Lyapunov 安定であるための必要十分条件であり、x∗x_*が∥⋅∥1\lVert\cdot\rVert_1について漸近安定であることは、∥⋅∥2\lVert\cdot\rVert_2について漸近安定であるための必要十分条件である。特に、x∗x_*が不安定であるかどうかはノルムの選び方によらない。

証明.φ\varphiをFFの流れとする。§E10.1 命題 2.2により、定数0<c≤C0<c\le Cが存在して、すべてのv∈Rnv\in\R^nでc∥v∥1≤∥v∥2≤C∥v∥1c\lVert v\rVert_1\le\lVert v\rVert_2\le C\lVert v\rVert_1が成り立つ。x∗x_*が∥⋅∥1\lVert\cdot\rVert_1について Lyapunov 安定であるとし、ε>0\varepsilon>0を取る。ε/C\varepsilon/Cに対して∥⋅∥1\lVert\cdot\rVert_1についてのδ>0\delta>0を取る。x∈Ux\in Uが∥x−x∗∥2<cδ\lVert x-x_*\rVert_2<c\deltaを満たすならば∥x−x∗∥1<δ\lVert x-x_*\rVert_1<\deltaであるから、[0,∞)⊆I(x)[0,\infty)\subseteq I(x)であり、すべてのt≥0t\ge0で∥φt(x)−x∗∥2≤C∥φt(x)−x∗∥1<ε\lVert\varphi_t(x)-x_*\rVert_2\le C\lVert\varphi_t(x)-x_*\rVert_1<\varepsilonである。したがってx∗x_*は∥⋅∥2\lVert\cdot\rVert_2について Lyapunov 安定である。さらにx∗x_*が∥⋅∥1\lVert\cdot\rVert_1について漸近安定であり、ρ>0\rho>0が漸近安定性の定義の条件を∥⋅∥1\lVert\cdot\rVert_1について満たすとする。x∈Ux\in Uが∥x−x∗∥2<cρ\lVert x-x_*\rVert_2<c\rhoを満たすならば∥x−x∗∥1<ρ\lVert x-x_*\rVert_1<\rhoであるから、[0,∞)⊆I(x)[0,\infty)\subseteq I(x)であり、∥φt(x)−x∗∥2≤C∥φt(x)−x∗∥1→0\lVert\varphi_t(x)-x_*\rVert_2\le C\lVert\varphi_t(x)-x_*\rVert_1\to0(t→∞)(t\to\infty)である。したがってx∗x_*は∥⋅∥2\lVert\cdot\rVert_2について漸近安定である。∥⋅∥1\lVert\cdot\rVert_1と∥⋅∥2\lVert\cdot\rVert_2の役割を入れ替えると、逆の含意を得る。不安定性は Lyapunov 安定性の否定である。▨

補題 2.4.n∈N≥1n\in\NNとし、U,V⊆RnU,V\subseteq\R^nを開集合、F:U→RnF:U\to\R^n、G:V→RnG:V\to\R^nを局所 Lipschitz 連続な写像とし、それらの流れをφ\varphi、ψ\psiとする。x∗∈Ux_*\in UをFFの平衡点、y∗∈Vy_*\in VをGGの平衡点とし、FFはx∗x_*でGGのy∗y_*と局所位相共役であるとする。このとき、次が成り立つ。

  1. x∗x_*が Lyapunov 安定ならば、y∗y_*は Lyapunov 安定である。
  2. x∗x_*が漸近安定ならば、y∗y_*は漸近安定である。

証明. 局所位相共役の定義の開近傍W⊆UW\subseteq U、W′⊆VW'\subseteq Vと同相写像h:W→W′h:W\to W'を取り、Rn\R^nのノルムを一つ固定する。p∈Rnp\in\R^n、r>0r>0についてB(p,r)={x∈Rn:∥x−p∥<r}B(p,r)=\{x\in\R^n:\lVert x-p\rVert<r\}と書く。

主張 2.4.1.ε>0\varepsilon>0とする。このとき、ε1>0\varepsilon_1>0が存在して、次が成り立つ。x∈Wx\in Wが[0,∞)⊆IF(x)[0,\infty)\subseteq I_F(x)を満たし、すべてのs≥0s\ge0でφs(x)∈B(x∗,ε1)\varphi_s(x)\in B(x_*,\varepsilon_1)であるならば、[0,∞)⊆IG(h(x))[0,\infty)\subseteq I_G(h(x))であり、すべてのt≥0t\ge0でψt(h(x))=h(φt(x))∈B(y∗,ε)\psi_t(h(x))=h(\varphi_t(x))\in B(y_*,\varepsilon)である。

証明.WWは開集合で、hhはx∗x_*で連続でありh(x∗)=y∗h(x_*)=y_*であるから、ε1>0\varepsilon_1>0が存在してB(x∗,ε1)⊆WB(x_*,\varepsilon_1)\subseteq Wかつh(B(x∗,ε1))⊆B(y∗,ε)h(B(x_*,\varepsilon_1))\subseteq B(y_*,\varepsilon)である。xxを仮定のとおりに取り、t≥0t\ge0とすると、00とttの間のすべてのssでφs(x)∈B(x∗,ε1)⊆W\varphi_s(x)\in B(x_*,\varepsilon_1)\subseteq Wである。局所位相共役の定義により、t∈IG(h(x))t\in I_G(h(x))であり、ψt(h(x))=h(φt(x))∈B(y∗,ε)\psi_t(h(x))=h(\varphi_t(x))\in B(y_*,\varepsilon)である。▨

主張 (1)を示す。x∗x_*が Lyapunov 安定であるとし、ε>0\varepsilon>0を取る。ε\varepsilonに対する主張 2.4.1のε1\varepsilon_1を取り、Lyapunov 安定性の定義によりε1\varepsilon_1に対するδ1>0\delta_1>0を取る。W′W'は開集合で、h−1h^{-1}はy∗y_*で連続でありh−1(y∗)=x∗h^{-1}(y_*)=x_*であるから、δ>0\delta>0が存在してB(y∗,δ)⊆W′B(y_*,\delta)\subseteq W'かつh−1(B(y∗,δ))⊆B(x∗,δ1)h^{-1}(B(y_*,\delta))\subseteq B(x_*,\delta_1)である。y∈Vy\in V、∥y−y∗∥<δ\lVert y-y_*\rVert<\deltaとし、x=h−1(y)x=h^{-1}(y)とおく。∥x−x∗∥<δ1\lVert x-x_*\rVert<\delta_1であるから、[0,∞)⊆IF(x)[0,\infty)\subseteq I_F(x)であり、すべてのs≥0s\ge0でφs(x)∈B(x∗,ε1)\varphi_s(x)\in B(x_*,\varepsilon_1)である。主張 2.4.1により、[0,∞)⊆IG(y)[0,\infty)\subseteq I_G(y)であり、すべてのt≥0t\ge0で∥ψt(y)−y∗∥<ε\lVert\psi_t(y)-y_*\rVert<\varepsilonである。

主張 (2)を示す。x∗x_*が漸近安定であるとし、ρ>0\rho>0を漸近安定性の定義の条件を満たすように取る。主張 (1)によりy∗y_*は Lyapunov 安定である。ε=1\varepsilon=1に対する主張 2.4.1のε1\varepsilon_1を取り、ε1\varepsilon_1に対する Lyapunov 安定性の定義のδ1>0\delta_1>0を取って、ρ1=min⁡{δ1,ρ}\rho_1=\min\{\delta_1,\rho\}とおく。h−1h^{-1}はy∗y_*で連続であるから、ρ′>0\rho'>0が存在してB(y∗,ρ′)⊆W′B(y_*,\rho')\subseteq W'かつh−1(B(y∗,ρ′))⊆B(x∗,ρ1)h^{-1}(B(y_*,\rho'))\subseteq B(x_*,\rho_1)である。y∈Vy\in V、∥y−y∗∥<ρ′\lVert y-y_*\rVert<\rho'とし、x=h−1(y)x=h^{-1}(y)とおく。∥x−x∗∥<δ1\lVert x-x_*\rVert<\delta_1であるから、主張 2.4.1により[0,∞)⊆IG(y)[0,\infty)\subseteq I_G(y)であり、すべてのt≥0t\ge0でψt(y)=h(φt(x))\psi_t(y)=h(\varphi_t(x))である。∥x−x∗∥<ρ\lVert x-x_*\rVert<\rhoであるからφt(x)→x∗\varphi_t(x)\to x_*であり、hhはx∗x_*で連続であるからψt(y)→h(x∗)=y∗\psi_t(y)\to h(x_*)=y_*(t→∞)(t\to\infty)である。▨

系 2.5.n∈N≥1n\in\NNとし、U,V⊆RnU,V\subseteq\R^nを開集合、F:U→RnF:U\to\R^n、G:V→RnG:V\to\R^nを局所 Lipschitz 連続な写像、h:U→Vh:U\to VをFFとGGの流れの位相共役とする。x∗∈Ux_*\in UをFFの平衡点とする。このとき、h(x∗)h(x_*)はGGの平衡点である。さらに、x∗x_*が Lyapunov 安定であることは、h(x∗)h(x_*)が Lyapunov 安定であるための必要十分条件であり、x∗x_*が漸近安定であることは、h(x∗)h(x_*)が漸近安定であるための必要十分条件である。

証明.FFとGGの流れをφ\varphi、ψ\psiとする。§E13.15 命題 7.2 (1)によりh(x∗)h(x_*)はGGの平衡点である。x∈Ux\in U、t∈IF(x)t\in I_F(x)ならば、位相共役の定義によりt∈IG(h(x))t\in I_G(h(x))かつh(φt(x))=ψt(h(x))h(\varphi_t(x))=\psi_t(h(x))であるから、W=UW=U、W′=VW'=VとしてFFはx∗x_*でGGのh(x∗)h(x_*)と局所位相共役である。y∈Vy\in V、t∈IG(y)t\in I_G(y)とし、x=h−1(y)x=h^{-1}(y)とおくと、IF(x)=IG(y)I_F(x)=I_G(y)であるからt∈IF(x)t\in I_F(x)であり、h(φt(x))=ψt(y)h(\varphi_t(x))=\psi_t(y)によりh−1(ψt(y))=φt(h−1(y))h^{-1}(\psi_t(y))=\varphi_t(h^{-1}(y))である。したがって、W=VW=V、W′=UW'=Uとして、GGはh(x∗)h(x_*)でFFのx∗x_*と局所位相共役である。この二つの局所位相共役に補題 2.4を適用すると、結論を得る。▨

例 2.6.U=RU=\R、f(x)=x(1−x)f(x)=x(1-x)とする。ffはR\R上でC1C^1級であるから、§E4.3 定理 2.1により局所 Lipschitz 連続である。ffの平衡点は00と11であり、ffは(0,1)(0,1)上で正、(−∞,0)(-\infty,0)と(1,∞)(1,\infty)上で負である。ffの流れをφ\varphiとし、R\Rには絶対値をノルムとして用いる。

0<x<10<x<1とする。補題 1.2のCCは(0,1)(0,1)であるから、補題 1.2 (1)と補題 1.2 (2)により、t∈I(x)t\in I(x)、t≥0t\ge0についてφt(x)∈[x,1)\varphi_t(x)\in[x,1)である。補題 1.2 (3)をK=[x,1]K=[x,1]に適用すると、[0,∞)⊆I(x)[0,\infty)\subseteq I(x)であり、φt(x)\varphi_t(x)は[x,1][x,1]に属する平衡点に収束する。その平衡点は11である。x>1x>1とすると、C=(1,∞)C=(1,\infty)であり、補題 1.2 (1)と補題 1.2 (2)によりt∈I(x)t\in I(x)、t≥0t\ge0についてφt(x)∈(1,x]\varphi_t(x)\in(1,x]である。補題 1.2 (3)をK=[1,x]K=[1,x]に適用すると、[0,∞)⊆I(x)[0,\infty)\subseteq I(x)であり、φt(x)→1\varphi_t(x)\to1(t→∞)(t\to\infty)である。x=1x=1ならば、§E13.15 命題 4.3 (2)によりすべてのt∈Rt\in\Rでφt(x)=1\varphi_t(x)=1である。以上により、x>0x>0ならば、すべてのt≥0t\ge0で∣φt(x)−1∣≤∣x−1∣\lvert\varphi_t(x)-1\rvert\le\lvert x-1\rvertであり、φt(x)→1\varphi_t(x)\to1(t→∞)(t\to\infty)である。

ε>0\varepsilon>0に対してδ=min⁡{ε,1}\delta=\min\{\varepsilon,1\}とおくと、∣x−1∣<δ\lvert x-1\rvert<\deltaならばx>0x>0であるから、[0,∞)⊆I(x)[0,\infty)\subseteq I(x)であり、すべてのt≥0t\ge0で∣φt(x)−1∣<ε\lvert\varphi_t(x)-1\rvert<\varepsilonである。また∣x−1∣<1\lvert x-1\rvert<1ならばφt(x)→1\varphi_t(x)\to1である。したがって11は漸近安定である。一方、任意のδ>0\delta>0について、0<x<min⁡{δ,1/2}0<x<\min\{\delta,1/2\}を満たすxxはφt(x)→1\varphi_t(x)\to1を満たすから、t≥0t\ge0が存在して∣φt(x)−0∣>1/2\lvert\varphi_t(x)-0\rvert>1/2である。したがって00はε=1/2\varepsilon=1/2について Lyapunov 安定性の条件を満たさず、不安定である。

3 線形系の安定性と二次元線形系の分類

命題 3.1.n∈N≥1n\in\NN、A∈Rn×nA\in\R^{n\times n}とし、F(x)=AxF(x)=Axで定まるF:Rn→RnF:\R^n\to\R^nの流れをφ\varphiとする。このとき、次が成り立つ。

  1. φ\varphiは完備な流れであり、すべてのt∈Rt\in\R、x∈Rnx\in\R^nについてφt(x)=etAx\varphi_t(x)=e^{tA}xである。
  2. S∈Rn×nS\in\R^{n\times n}を可逆行列とし、G(y)=SAS−1yG(y)=SAS^{-1}yでG:Rn→RnG:\R^n\to\R^nを定める。このとき、h(x)=Sxh(x)=Sxで定まるh:Rn→Rnh:\R^n\to\R^nはFFとGGの流れの位相共役である。

証明.主張 (1)を示す。Rn\R^nのノルムと対応する作用素ノルムについて∥Ax−Ay∥≤∥A∥∥x−y∥\lVert Ax-Ay\rVert\le\lVert A\rVert\lVert x-y\rVertであるから、FFは Lipschitz 連続である。x∈Rnx\in\R^nとし、z(t)=etAxz(t)=e^{tA}x(t∈Rt\in\R)とおく。§E13.12 定義 1.2によりetAe^{tA}は実行列であり、§E13.12 定理 2.1によりzzはR\R上でC1C^1級で、z′(t)=AetAx=F(z(t))z'(t)=Ae^{tA}x=F(z(t))、z(0)=xz(0)=xである。§E13.15 補題 1.2をJ=RJ=\R、τ=0\tau=0に適用すると、R⊆I(x)\R\subseteq I(x)であり、R\R上でz=φ(⋅,x)z=\varphi(\cdot,x)である。

主張 (2)を示す。hhは連続な逆写像y↦S−1yy\mapsto S^{-1}yをもつ同相写像である。GGの流れをψ\psiとする。主張 (1)をAAとSAS−1SAS^{-1}に適用すると、FFとGGの流れはともに完備であり、ψt(y)=etSAS−1y\psi_t(y)=e^{tSAS^{-1}}yである。§E13.12 命題 4.1を行列tAtAとSSに適用するとetSAS−1=SetAS−1e^{tSAS^{-1}}=Se^{tA}S^{-1}であるから、すべてのt∈Rt\in\R、x∈Rnx\in\R^nについてIG(h(x))=R=IF(x)I_G(h(x))=\R=I_F(x)であり、ψt(h(x))=SetAx=h(φt(x))\psi_t(h(x))=Se^{tA}x=h(\varphi_t(x))である。▨

補題 3.2.d∈N≥1d\in\NNとし、A∈Rd×dA\in\R^{d\times d}の複素固有値の実部がすべて負であるとして、その最大値をα\alphaとする。Rd\R^dのノルムを一つ固定し、行列には対応する作用素ノルムを用いる。このとき、0<γ<−α0<\gamma<-\alphaを満たす任意のγ\gammaについて、M≥1M\ge1が存在して

∥etA∥≤Me−γt(t≥0)\lVert e^{tA}\rVert\le Me^{-\gamma t}\qquad(t\ge0)(3.2.1)\text{(3.2.1)}

が成り立つ。特に、式 (3.2.1)を満たすM≥1M\ge1とγ>0\gamma>0が存在する。

証明. 作用素ノルムはRd×d\R^{d\times d}上のノルムであるから、§E13.12 系 6.6をK=R\mathbb K=\R、B=AB=Aに適用して結論を得る。▨

命題 3.3.n∈N≥1n\in\NN、A∈Rn×nA\in\R^{n\times n}とし、AAの相異なる複素固有値全体の集合をσ(A)\sigma(A)と書き、F(x)=AxF(x)=AxでF:Rn→RnF:\R^n\to\R^nを定める。Rn\R^nのノルムを一つ固定し、行列には対応する作用素ノルムを用いる。このとき、原点はFFの平衡点であり、次が成り立つ。

  1. 原点が Lyapunov 安定であることは、sup⁡t≥0∥etA∥<∞\sup_{t\ge0}\lVert e^{tA}\rVert<\inftyであるための必要十分条件である。さらに、sup⁡t≥0∥etA∥<∞\sup_{t\ge0}\lVert e^{tA}\rVert<\inftyであることは、すべてのβ∈σ(A)\beta\in\sigma(A)がℜβ≤0\Re\beta\le0を満たし、ℜβ=0\Re\beta=0を満たすすべてのβ∈σ(A)\beta\in\sigma(A)についてCn\C^nの部分空間としてGβ(A)=ker⁡(A−βI)G_\beta(A)=\ker(A-\beta I)が成り立つための必要十分条件である。
  2. 原点が漸近安定であることは、すべてのβ∈σ(A)\beta\in\sigma(A)がℜβ<0\Re\beta<0を満たすための必要十分条件である。このとき、すべてのx∈Rnx\in\R^nについてetAx→0e^{tA}x\to0(t→∞)(t\to\infty)である。

証明.命題 3.1 (1)により、FFの流れφ\varphiは完備であり、すべてのt∈Rt\in\R、x∈Rnx\in\R^nについてφt(x)=etAx\varphi_t(x)=e^{tA}xである。

主張 (1)を示す。C=sup⁡t≥0∥etA∥<∞C=\sup_{t\ge0}\lVert e^{tA}\rVert<\inftyとする。e0A=Ie^{0A}=IであるからC≥1C\ge1である。ε>0\varepsilon>0に対してδ=ε/C\delta=\varepsilon/Cとおくと、∥x∥<δ\lVert x\rVert<\deltaならば、すべてのt≥0t\ge0で∥φt(x)∥≤C∥x∥<ε\lVert\varphi_t(x)\rVert\le C\lVert x\rVert<\varepsilonである。したがって原点は Lyapunov 安定である。逆に、原点が Lyapunov 安定であるとし、ε=1\varepsilon=1に対するδ>0\delta>0を取る。x≠0x\ne0についてy=δx/(2∥x∥)y=\delta x/(2\lVert x\rVert)とおくと、∥y∥<δ\lVert y\rVert<\deltaであるから、すべてのt≥0t\ge0で∥etAy∥<1\lVert e^{tA}y\rVert<1であり、∥etAx∥<2∥x∥/δ\lVert e^{tA}x\rVert<2\lVert x\rVert/\deltaである。したがってsup⁡t≥0∥etA∥≤2/δ\sup_{t\ge0}\lVert e^{tA}\rVert\le2/\deltaである。作用素ノルムはRn×n\R^{n\times n}上のノルムであるから、§E13.12 補題 6.4 (2)をK=R\mathbb K=\R、B=AB=Aに適用すると、残りの同値を得る。

主張 (2)を示す。すべてのβ∈σ(A)\beta\in\sigma(A)がℜβ<0\Re\beta<0を満たすとする。補題 3.2により、M≥1M\ge1とγ>0\gamma>0が存在して、すべてのt≥0t\ge0で∥etA∥≤Me−γt\lVert e^{tA}\rVert\le Me^{-\gamma t}である。したがってsup⁡t≥0∥etA∥≤M\sup_{t\ge0}\lVert e^{tA}\rVert\le Mであり、主張 (1)により原点は Lyapunov 安定である。さらに、すべてのx∈Rnx\in\R^nで∥etAx∥≤Me−γt∥x∥→0\lVert e^{tA}x\rVert\le Me^{-\gamma t}\lVert x\rVert\to0(t→∞)(t\to\infty)であるから、原点は漸近安定である。

逆に、原点が漸近安定であるとし、ρ>0\rho>0を漸近安定性の定義の条件を満たすように取る。β∈σ(A)\beta\in\sigma(A)がℜβ≥0\Re\beta\ge0を満たすと仮定する。β\betaに対するAAの固有ベクトルv∈Cnv\in\C^nを取り、v=p+iqv=p+iq(p,q∈Rnp,q\in\R^n)と書く。§E13.12 補題 6.2 (3)によりetAv=etβve^{tA}v=e^{t\beta}vであり、etAe^{tA}は実行列であるから、etApe^{tA}pとetAqe^{tA}qはそれぞれetβve^{t\beta}vの実部と虚部である。vj≠0v_j\ne0となる成分vjv_jを取ると、すべてのt≥0t\ge0で∣etβvj∣=etℜβ∣vj∣≥∣vj∣>0\lvert e^{t\beta}v_j\rvert=e^{t\Re\beta}\lvert v_j\rvert\ge\lvert v_j\rvert>0である。etAp→0e^{tA}p\to0かつetAq→0e^{tA}q\to0(t→∞)(t\to\infty)ならば、§E10.1 命題 2.2により両者の各成分が00に収束し、etβvj→0e^{t\beta}v_j\to0となって前の不等式に反する。したがって、w=pw=pまたはw=qw=qのいずれかについて、etAwe^{tA}wは00に収束しない。このときw≠0w\ne0であり、w′=ρw/(2∥w∥)w'=\rho w/(2\lVert w\rVert)とおくと、∥w′∥<ρ\lVert w'\rVert<\rhoであるが、etAw′=ρetAw/(2∥w∥)e^{tA}w'=\rho e^{tA}w/(2\lVert w\rVert)は00に収束しない。これはρ\rhoの取り方に反する。したがって、すべてのβ∈σ(A)\beta\in\sigma(A)はℜβ<0\Re\beta<0を満たす。▨

定義 3.4.A∈R2×2A\in\R^{2\times2}とし、線形系x′=Axx'=Axの平衡点である原点を考える。AAの複素固有値を重複度を込めてλ1,λ2\lambda_1,\lambda_2とする。λ1,λ2\lambda_1,\lambda_2が相異なる実数で同じ符号をもつとき、原点を結節点 (node) という。実数λ≠0\lambda\ne0についてA=λIA=\lambda Iであるとき、原点を星形結節点 (star node) という。λ1=λ2\lambda_1=\lambda_2が00でない実数であり、A≠λ1IA\ne\lambda_1Iであるとき、原点を非正則結節点 (improper node) という。λ1,λ2\lambda_1,\lambda_2が実数でλ1λ2<0\lambda_1\lambda_2<0を満たすとき、原点を鞍点 (saddle point) という。λ1,λ2\lambda_1,\lambda_2が実数でなく、その実部が00でないとき、原点を渦状点 (spiral point) という。λ1,λ2\lambda_1,\lambda_2が実数でなく、その実部が00であるとき、原点を中心 (center of a planar linear system) という。

定理 3.5.A∈R2×2A\in\R^{2\times2}とし、F(x)=AxF(x)=Axで定まるF:R2→R2F:\R^2\to\R^2の流れをφ\varphiとし、∥⋅∥\lVert\cdot\rVertをR2\R^2のノルムとする。このとき、原点はFFの平衡点であり、次が成り立つ。

  1. 原点が結節点であるとき、AAの二つの固有値がともに負ならば原点は漸近安定であり、ともに正ならば原点は不安定である。
  2. 原点が鞍点であるとき、原点は不安定である。さらに、AAの負の固有値をλs\lambda_sとすると、x∈R2x\in\R^2について、Ax=λsxAx=\lambda_sxならばφt(x)→0\varphi_t(x)\to0(t→∞)(t\to\infty)であり、Ax≠λsxAx\ne\lambda_sxならば∥φt(x)∥→∞\lVert\varphi_t(x)\rVert\to\infty(t→∞)(t\to\infty)である。
  3. 原点が星形結節点または非正則結節点であるとき、AAの重複固有値が負ならば原点は漸近安定であり、正ならば原点は不安定である。
  4. 原点が渦状点であるとき、AAの固有値の実部が負ならば原点は漸近安定であり、正ならば原点は不安定である。原点が中心であるとき、原点は Lyapunov 安定であるが漸近安定でない。
  5. AAが固有値00をもつとき、原点はFFの孤立平衡点でない。AAの固有値が00と実数λ≠0\lambda\ne0であるとき、λ<0\lambda<0ならば原点は Lyapunov 安定であるが漸近安定でなく、λ>0\lambda>0ならば原点は不安定である。AAの固有値が00だけであるとき、原点が Lyapunov 安定であることは、A=0A=0であるための必要十分条件であり、このとき原点は漸近安定でない。

証明.命題 3.1 (1)により、すべてのt∈Rt\in\R、x∈R2x\in\R^2についてφt(x)=etAx\varphi_t(x)=e^{tA}xである。AAの相異なる複素固有値全体の集合をσ(A)\sigma(A)とする。

主張 3.5.1. 次の三つが成り立つ。σ(A)\sigma(A)のすべての元の実部が負ならば、原点は漸近安定である。σ(A)\sigma(A)が実部00の元をもつならば、原点は漸近安定でない。σ(A)\sigma(A)が実部が正の元をもつならば、原点は不安定である。

証明. 前の二つは命題 3.3 (2)から従い、三つ目は命題 3.3 (1)から従う。▨

主張 3.5.2.AAが相異なる二つの複素固有値β1,β2\beta_1,\beta_2をもつならば、j=1,2j=1,2についてC2\C^2の部分空間としてGβj(A)=ker⁡(A−βjI)G_{\beta_j}(A)=\ker(A-\beta_jI)である。

証明.§E13.12 定義 6.1によりC2=Gβ1(A)⊕Gβ2(A)\C^2=G_{\beta_1}(A)\oplus G_{\beta_2}(A)であり、各Gβj(A)G_{\beta_j}(A)はβj\beta_jに対する固有ベクトルを含むから零でない。したがって各Gβj(A)G_{\beta_j}(A)は一次元であり、零でない部分空間ker⁡(A−βjI)⊆Gβj(A)\ker(A-\beta_jI)\subseteq G_{\beta_j}(A)と一致する。▨

主張 3.5.1により、主張 (1)、主張 (3)、主張 (4)の渦状点についての主張、および主張 (2)の不安定性が従う。原点が中心であるとき、ω≠0\omega\ne0が存在してσ(A)={iω,−iω}\sigma(A)=\{i\omega,-i\omega\}であるから、主張 3.5.2と命題 3.3 (1)により原点は Lyapunov 安定であり、主張 3.5.1により漸近安定でない。これで主張 (4)は示された。

主張 (2)の後半を示す。AAの正の固有値をλu\lambda_uとする。A−λsIA-\lambda_sIとA−λuIA-\lambda_uIは実行列で可逆でないから、実の固有ベクトルvs,vu∈R2∖{0}v_s,v_u\in\R^2\setminus\{0\}が存在してAvs=λsvsAv_s=\lambda_sv_s、Avu=λuvuAv_u=\lambda_uv_uである。vu=cvsv_u=cv_sを満たすc∈Rc\in\Rが存在するならばλuvu=Avu=cλsvs=λsvu\lambda_uv_u=Av_u=c\lambda_sv_s=\lambda_sv_uとなってλu≠λs\lambda_u\ne\lambda_sに反するから、vs,vuv_s,v_uはR2\R^2の基底である。x=avs+bvux=av_s+bv_u(a,b∈Ra,b\in\R)と書くと、(A−λsI)x=b(λu−λs)vu(A-\lambda_sI)x=b(\lambda_u-\lambda_s)v_uであるから、Ax=λsxAx=\lambda_sxであることとb=0b=0であることは同値である。§E13.12 補題 6.2 (3)により

φt(x)=aeλstvs+beλutvu(t∈R)\varphi_t(x)=ae^{\lambda_st}v_s+be^{\lambda_ut}v_u \qquad(t\in\R)

である。b=0b=0ならば∥φt(x)∥=∣a∣eλst∥vs∥→0\lVert\varphi_t(x)\rVert=\lvert a\rvert e^{\lambda_st}\lVert v_s\rVert\to0(t→∞)(t\to\infty)である。N(avs+bvu)=max⁡{∣a∣,∣b∣}N(av_s+bv_u)=\max\{\lvert a\rvert,\lvert b\rvert\}はR2\R^2のノルムであるから、§E10.1 命題 2.2により定数K>0K>0が存在して、すべてのy∈R2y\in\R^2でN(y)≤K∥y∥N(y)\le K\lVert y\rVertである。b≠0b\ne0ならば∥φt(x)∥≥K−1N(φt(x))≥K−1∣b∣eλut→∞\lVert\varphi_t(x)\rVert\ge K^{-1}N(\varphi_t(x))\ge K^{-1}\lvert b\rvert e^{\lambda_ut}\to\infty(t→∞)(t\to\infty)である。

主張 (5)を示す。AAは固有値00をもつ実行列であるから、v∈R2∖{0}v\in\R^2\setminus\{0\}が存在してAv=0Av=0である。原点の任意の開近傍WWについて、十分小さいs>0s>0に対してsv∈Wsv\in Wであり、sv≠0sv\ne0はFFの平衡点である。したがって原点は孤立平衡点でない。AAの固有値が00とλ≠0\lambda\ne0であるとき、λ\lambdaはAAのトレースに等しいから実数である。主張 3.5.2によりG0(A)=ker⁡AG_0(A)=\ker Aであるから、λ<0\lambda<0ならば命題 3.3 (1)により原点は Lyapunov 安定であり、0∈σ(A)0\in\sigma(A)であるから主張 3.5.1により漸近安定でない。λ>0\lambda>0ならば主張 3.5.1により原点は不安定である。AAの固有値が00だけであるとき、§E13.12 定義 6.1によりG0(A)=C2G_0(A)=\C^2であるから、命題 3.3 (1)により、原点が Lyapunov 安定であることはker⁡A=C2\ker A=\C^2、すなわちA=0A=0と同値である。このとき0∈σ(A)0\in\sigma(A)であるから、主張 3.5.1により原点は漸近安定でない。▨

例 3.6 (鞍点).

A=(1221)A=\begin{pmatrix}1&2\\2&1\end{pmatrix}

の特性多項式は(z−1)2−4=(z−3)(z+1)(z-1)^2-4=(z-3)(z+1)であるから、固有値は3,−13,-1であり、原点は鞍点である。A(1,1)T=3(1,1)TA(1,1)^{\mathsf T}=3(1,1)^{\mathsf T}、A(1,−1)T=−(1,−1)TA(1,-1)^{\mathsf T}=-(1,-1)^{\mathsf T}であるから、x=u(1,1)T+v(1,−1)Tx=u(1,1)^{\mathsf T}+v(1,-1)^{\mathsf T}と書くと、x′=Axx'=Axの流れは

φt(x)=ue3t(1,1)T+ve−t(1,−1)T\varphi_t(x)=ue^{3t}(1,1)^{\mathsf T}+ve^{-t}(1,-1)^{\mathsf T}

である。定理 3.5 (2)により、φt(x)→0\varphi_t(x)\to0(t→∞)(t\to\infty)となるのはxxが直線x1+x2=0x_1+x_2=0上にあるときに限り、それ以外のxxについては∥φt(x)∥→∞\lVert\varphi_t(x)\rVert\to\inftyである。

例 3.7 (安定な渦状点).

A=(−1−11−1),J=(0−110)A=\begin{pmatrix}-1&-1\\1&-1\end{pmatrix},\qquad J=\begin{pmatrix}0&-1\\1&0\end{pmatrix}

とおくと、A=−I+JA=-I+Jであり、AAの固有値は−1±i-1\pm iであるから、原点は渦状点である。−I-IとJJは可換であるから、§E13.12 命題 3.2と§E13.12 例 4.4により

etA=e−tetJ=e−t(cos⁡t−sin⁡tsin⁡tcos⁡t)e^{tA}=e^{-t}e^{tJ} =e^{-t} \begin{pmatrix} \cos t&-\sin t\\ \sin t&\cos t \end{pmatrix}

である。したがって、Euclid ノルム∥⋅∥2\lVert\cdot\rVert_2について∥etAx∥2=e−t∥x∥2\lVert e^{tA}x\rVert_2=e^{-t}\lVert x\rVert_2であり、原点は漸近安定である。x≠0x\ne0の軌道上の点etAxe^{tA}xは、原点のまわりの偏角がxxの偏角からttだけ増え、原点からの距離がe−te^{-t}倍になった点である。

4 非線形系の線形化

注意 4.1.n∈N≥1n\in\NNとし、U⊆RnU\subseteq\R^nを開集合、F:U→RnF:U\to\R^nを写像、x∗∈Ux_*\in UをF(x∗)=0F(x_*)=0を満たす点とし、FFはx∗x_*で微分可能であるとする。A=DF(x∗)A=DF(x_*)とおき、x∗+u∈Ux_*+u\in Uを満たすu∈Rnu\in\R^nについてr(u)=F(x∗+u)−Aur(u)=F(x_*+u)-Auと定める。F(x∗)=0F(x_*)=0であるから、微分可能性の定義により∥r(u)∥/∥u∥→0\lVert r(u)\rVert/\lVert u\rVert\to0(u→0)(u\to0)である。xxがx′=F(x)x'=F(x)の古典解ならば、u=x−x∗u=x-x_*は

u′=F(x∗+u)=Au+r(u),A=DF(x∗),∥r(u)∥∥u∥⟶0(u→0)u' =F(x_*+u) =Au+r(u), \qquad A=DF(x_*), \qquad \frac{\lVert r(u)\rVert} {\lVert u\rVert}\longrightarrow0 \quad(u\to0)(4.1.1)\text{(4.1.1)}

を満たす。

定義 4.2.n∈N≥1n\in\NNとし、U⊆RnU\subseteq\R^nを開集合、F:U→RnF:U\to\R^nを写像、x∗∈Ux_*\in UをF(x∗)=0F(x_*)=0を満たす点とし、FFはx∗x_*で微分可能であるとする。Jacobi 行列DF(x∗)DF(x_*)が実部00の複素固有値をもたないとき、x∗x_*をFFの双曲型平衡点 (hyperbolic equilibrium point) という。

補題 4.3.n∈N≥1n\in\NNとし、A∈Rn×nA\in\R^{n\times n}は実部が00の複素固有値をもたないとする。∥⋅∥\lVert\cdot\rVertをRn\R^nのノルムとし、行列には対応する作用素ノルムを用いる。Ps,Pu,Es,EuP_s,P_u,E^s,E^uを、AAに対する§E13.12 補題 6.7(K=R\mathbb K=\R、d=nd=n)の射影と部分空間とする。K≥1K\ge1、α>0\alpha>0は、すべてのt≥0t\ge0について∥etAPs∥≤Ke−αt\lVert e^{tA}P_s\rVert\le Ke^{-\alpha t}を満たし、すべてのt≤0t\le0について∥etAPu∥≤Keαt\lVert e^{tA}P_u\rVert\le Ke^{\alpha t}を満たすとし、0<β<α0<\beta<\alphaとする。

∥x∥s=sup⁡t≥0eβt∥etAx∥(x∈Es),∥x∥u=sup⁡t≤0e−βt∥etAx∥(x∈Eu)\lVert x\rVert_s=\sup_{t\ge0}e^{\beta t}\lVert e^{tA}x\rVert\quad(x\in E^s),\qquad \lVert x\rVert_u=\sup_{t\le0}e^{-\beta t}\lVert e^{tA}x\rVert\quad(x\in E^u)

とおき、∥x∥∗=max⁡{∥Psx∥s,∥Pux∥u}\lVert x\rVert_*=\max\{\lVert P_sx\rVert_s,\lVert P_ux\rVert_u\}(x∈Rnx\in\R^n)、p=max⁡{∥Ps∥,∥Pu∥}p=\max\{\lVert P_s\rVert,\lVert P_u\rVert\}とおく。このとき、次が成り立つ。

  1. ∥⋅∥s\lVert\cdot\rVert_s、∥⋅∥u\lVert\cdot\rVert_uはそれぞれEsE^s、EuE^uのノルムであり、y∈Esy\in E^sについて∥y∥≤∥y∥s≤K∥y∥\lVert y\rVert\le\lVert y\rVert_s\le K\lVert y\rVert、y∈Euy\in E^uについて∥y∥≤∥y∥u≤K∥y∥\lVert y\rVert\le\lVert y\rVert_u\le K\lVert y\rVertである。さらに、∥⋅∥∗\lVert\cdot\rVert_*はRn\R^nのノルムであり、すべてのx∈Rnx\in\R^nについて12∥x∥≤∥x∥∗≤Kp∥x∥\tfrac12\lVert x\rVert\le\lVert x\rVert_*\le Kp\lVert x\rVertである。
  2. すべてのt≥0t\ge0についてetAEs=Ese^{tA}E^s=E^s、etAEu=Eue^{tA}E^u=E^uであり、y∈Esy\in E^sならば∥etAy∥s≤e−βt∥y∥s\lVert e^{tA}y\rVert_s\le e^{-\beta t}\lVert y\rVert_s、y∈Euy\in E^uならば∥e−tAy∥u≤e−βt∥y∥u\lVert e^{-tA}y\rVert_u\le e^{-\beta t}\lVert y\rVert_uである。

証明.主張 (1)を示す。§E13.12 補題 6.7 (4)により、y∈Esy\in E^sならばPsy=yP_sy=yであり、y∈Euy\in E^uならばPuy=yP_uy=yである。y∈Esy\in E^s、t≥0t\ge0について

eβt∥etAy∥=eβt∥etAPsy∥≤Ke−(α−β)t∥y∥≤K∥y∥e^{\beta t}\lVert e^{tA}y\rVert=e^{\beta t}\lVert e^{tA}P_sy\rVert\le Ke^{-(\alpha-\beta)t}\lVert y\rVert\le K\lVert y\rVert

であり、t=0t=0での値は∥y∥\lVert y\rVertであるから、∥y∥≤∥y∥s≤K∥y∥\lVert y\rVert\le\lVert y\rVert_s\le K\lVert y\rVertである。y∈Euy\in E^u、t≤0t\le0について

e−βt∥etAy∥=e−βt∥etAPuy∥≤Ke(α−β)t∥y∥≤K∥y∥e^{-\beta t}\lVert e^{tA}y\rVert=e^{-\beta t}\lVert e^{tA}P_uy\rVert\le Ke^{(\alpha-\beta)t}\lVert y\rVert\le K\lVert y\rVert

であるから、∥y∥≤∥y∥u≤K∥y∥\lVert y\rVert\le\lVert y\rVert_u\le K\lVert y\rVertである。各ttについてy↦eβt∥etAy∥y\mapsto e^{\beta t}\lVert e^{tA}y\rVertとy↦e−βt∥etAy∥y\mapsto e^{-\beta t}\lVert e^{tA}y\rVertは半ノルムであるから、それらの有限な上限∥⋅∥s\lVert\cdot\rVert_s、∥⋅∥u\lVert\cdot\rVert_uも半ノルムであり、∥y∥\lVert y\rVert以上であることからノルムである。PsP_s、PuP_uは線形であるから∥⋅∥∗\lVert\cdot\rVert_*はRn\R^nの半ノルムである。§E13.12 補題 6.7 (2)によりPs+Pu=IP_s+P_u=Iであるから、x∈Rnx\in\R^nについて

∥x∥≤∥Psx∥+∥Pux∥≤∥Psx∥s+∥Pux∥u≤2∥x∥∗\lVert x\rVert\le\lVert P_sx\rVert+\lVert P_ux\rVert\le\lVert P_sx\rVert_s+\lVert P_ux\rVert_u\le2\lVert x\rVert_*

であり、∥x∥∗≤Kmax⁡{∥Psx∥,∥Pux∥}≤Kp∥x∥\lVert x\rVert_*\le K\max\{\lVert P_sx\rVert,\lVert P_ux\rVert\}\le Kp\lVert x\rVertである。前の不等式により∥⋅∥∗\lVert\cdot\rVert_*はノルムである。

主張 (2)を示す。t≥0t\ge0とする。§E13.12 補題 6.7 (4)によりetAEs=Ese^{tA}E^s=E^s、etAEu=Eue^{tA}E^u=E^u、e−tAEu=Eue^{-tA}E^u=E^uであり、§E13.12 定理 2.1によりs∈Rs\in\RについてesAetA=e(s+t)Ae^{sA}e^{tA}=e^{(s+t)A}である。y∈Esy\in E^sについて

∥etAy∥s=sup⁡s≥0eβs∥e(s+t)Ay∥=e−βtsup⁡s≥teβs∥esAy∥≤e−βt∥y∥s\lVert e^{tA}y\rVert_s=\sup_{s\ge0}e^{\beta s}\lVert e^{(s+t)A}y\rVert=e^{-\beta t}\sup_{s\ge t}e^{\beta s}\lVert e^{sA}y\rVert\le e^{-\beta t}\lVert y\rVert_s

であり、y∈Euy\in E^uについて

∥e−tAy∥u=sup⁡s≤0e−βs∥e(s−t)Ay∥=e−βtsup⁡s≤−te−βs∥esAy∥≤e−βt∥y∥u\lVert e^{-tA}y\rVert_u=\sup_{s\le0}e^{-\beta s}\lVert e^{(s-t)A}y\rVert=e^{-\beta t}\sup_{s\le-t}e^{-\beta s}\lVert e^{sA}y\rVert\le e^{-\beta t}\lVert y\rVert_u

である。▨

補題 4.4.EEを有限次元実ベクトル空間とし、E=Eu⊕EsE=E^u\oplus E^sをEu≠{0}E^u\ne\{0\}を満たす直和分解とする。EuE^uとEsE^sにノルム∥⋅∥\lVert\cdot\rVertを与え、zu∈Euz^u\in E^u、zs∈Esz^s\in E^sについて∥zu+zs∥=max⁡{∥zu∥,∥zs∥}\lVert z^u+z^s\rVert=\max\{\lVert z^u\rVert,\lVert z^s\rVert\}によりEEのノルムを定める。線形写像B:E→EB:E\to EはB(Eu)⊆EuB(E^u)\subseteq E^u、B(Es)⊆EsB(E^s)\subseteq E^sを満たし、Bu=B∣Eu:Eu→EuB_u=B|_{E^u}:E^u\to E^uは可逆であるとし、Bs=B∣Es:Es→EsB_s=B|_{E^s}:E^s\to E^sとおく。実数a,ba,bは0≤a<b0\le a<b、b>1b>1を満たし、

∥Bsv∥≤a∥v∥(v∈Es),∥Bu−1v∥≤b−1∥v∥(v∈Eu)\lVert B_sv\rVert\le a\lVert v\rVert\quad(v\in E^s),\qquad \lVert B_u^{-1}v\rVert\le b^{-1}\lVert v\rVert\quad(v\in E^u)(4.4.1)\text{(4.4.1)}

が成り立つとする。N⊆EN\subseteq Eを00の開近傍、G:N→EG:N\to EをG(z)=Bz+g(z)G(z)=Bz+g(z)の形の写像とし、g(0)=0g(0)=0かつ∥g(z)∥/∥z∥→0\lVert g(z)\rVert/\lVert z\rVert\to0(z→0)(z\to0)であるとする。このとき、r>0r>0が存在して、{z∈E:∥z∥<r}⊆N\{z\in E:\lVert z\rVert<r\}\subseteq Nであり、任意のδ>0\delta>0に対して、∥z0∥<δ\lVert z_0\rVert<\deltaを満たすz0∈Ez_0\in Eと整数k≥1k\ge1で、zj=G(zj−1)z_j=G(z_{j-1})(1≤j≤k1\le j\le k)により定まるz1,…,zkz_1,\ldots,z_kが、0≤j≤k−10\le j\le k-1で∥zj∥<r\lVert z_j\rVert<rを満たし、∥zk∥≥r\lVert z_k\rVert\ge rを満たすものが存在する。

証明.η>0\eta>0をa+η<b−ηa+\eta<b-\etaかつμ=b−η>1\mu=b-\eta>1を満たすように取る。g(0)=0g(0)=0と∥g(z)∥/∥z∥→0\lVert g(z)\rVert/\lVert z\rVert\to0により、r>0r>0が存在して、∥z∥<r\lVert z\rVert<rならばz∈Nz\in Nかつ∥g(z)∥≤η∥z∥\lVert g(z)\rVert\le\eta\lVert z\rVertである。z∈Ez\in Eをz=zu+zsz=z^u+z^s(zu∈Euz^u\in E^u、zs∈Esz^s\in E^s)と分解し、錐C={z∈E:∥zs∥≤∥zu∥}\mathcal C=\{z\in E:\lVert z^s\rVert\le\lVert z^u\rVert\}を考える。ノルムの定義により∥g(z)u∥≤∥g(z)∥\lVert g(z)^u\rVert\le\lVert g(z)\rVert、∥g(z)s∥≤∥g(z)∥\lVert g(z)^s\rVert\le\lVert g(z)\rVertである。

主張 4.4.1.z∈Cz\in\mathcal C、∥z∥<r\lVert z\rVert<rならば、G(z)∈CG(z)\in\mathcal Cであり、∥G(z)∥≥μ∥z∥\lVert G(z)\rVert\ge\mu\lVert z\rVertである。

証明.∥z∥=∥zu∥\lVert z\rVert=\lVert z^u\rVertであり、BBはEuE^uとEsE^sを不変にするからG(z)u=Buzu+g(z)uG(z)^u=B_uz^u+g(z)^u、G(z)s=Bszs+g(z)sG(z)^s=B_sz^s+g(z)^sである。式 (4.4.1)から∥zu∥=∥Bu−1Buzu∥≤b−1∥Buzu∥\lVert z^u\rVert=\lVert B_u^{-1}B_uz^u\rVert\le b^{-1}\lVert B_uz^u\rVertであるから、

∥G(z)u∥≥∥Buzu∥−∥g(z)u∥≥(b−η)∥zu∥,∥G(z)s∥≤∥Bszs∥+∥g(z)s∥≤a∥zs∥+η∥zu∥≤(a+η)∥zu∥\begin{aligned} \lVert G(z)^u\rVert &\ge\lVert B_uz^u\rVert-\lVert g(z)^u\rVert \ge(b-\eta)\lVert z^u\rVert,\\ \lVert G(z)^s\rVert &\le\lVert B_sz^s\rVert+\lVert g(z)^s\rVert \le a\lVert z^s\rVert+\eta\lVert z^u\rVert \le(a+\eta)\lVert z^u\rVert \end{aligned}

である。a+η<b−ηa+\eta<b-\etaであるからG(z)∈CG(z)\in\mathcal Cであり、∥G(z)∥=∥G(z)u∥≥μ∥z∥\lVert G(z)\rVert=\lVert G(z)^u\rVert\ge\mu\lVert z\rVertである。▨

δ>0\delta>0とし、z0∈Euz_0\in E^uを0<∥z0∥<min⁡{δ,r}0<\lVert z_0\rVert<\min\{\delta,r\}を満たすように取る。z0∈Cz_0\in\mathcal Cである。j≥1j\ge1について、z0,…,zj−1z_0,\ldots,z_{j-1}が定まってC\mathcal Cに属し、ノルムがrr未満であり、∥zj−1∥≥μj−1∥z0∥\lVert z_{j-1}\rVert\ge\mu^{j-1}\lVert z_0\rVertであるならば、主張 4.4.1によりzj=G(zj−1)z_j=G(z_{j-1})が定まってC\mathcal Cに属し、∥zj∥≥μj∥z0∥\lVert z_j\rVert\ge\mu^j\lVert z_0\rVertである。μ>1\mu>1であるからμj∥z0∥≥r\mu^j\lVert z_0\rVert\ge rを満たすjjが存在し、したがって∥zk∥≥r\lVert z_k\rVert\ge rとなる最小のk≥1k\ge1が存在する。このkkについて、0≤j≤k−10\le j\le k-1で∥zj∥<r\lVert z_j\rVert<rである。▨

定理 4.5 (Lyapunov の間接法).n∈N≥1n\in\NNとし、U⊆RnU\subseteq\R^nを開集合、F:U→RnF:U\to\R^nを局所 Lipschitz 連続な写像、φ\varphiをFFの流れ、x∗∈Ux_*\in UをFFの平衡点とする。x∗x_*の開近傍W⊆UW\subseteq Uが存在してF∣WF|_WはC1C^1級であるとし、A=DF(x∗)A=DF(x_*)とおく。このとき、次が成り立つ。

  1. AAの複素固有値の実部がすべて負ならば、x∗x_*は漸近安定である。さらに、その実部の最大値をα\alphaとすると、Rn\R^nの任意のノルム∥⋅∥\lVert\cdot\rVertと0<γ<−α0<\gamma<-\alphaを満たす任意のγ\gammaについて、δ>0\delta>0とM≥1M\ge1が存在して、∥x−x∗∥<δ\lVert x-x_*\rVert<\deltaを満たす任意のx∈Ux\in Uについて[0,∞)⊆I(x)[0,\infty)\subseteq I(x)であり、すべてのt≥0t\ge0で∥φt(x)−x∗∥≤Me−γt∥x−x∗∥\lVert\varphi_t(x)-x_*\rVert\le Me^{-\gamma t}\lVert x-x_*\rVertである。
  2. AAが実部が正の複素固有値を少なくとも一つもつならば、x∗x_*は不安定である。
  3. x∗x_*が双曲型平衡点ならば、x∗x_*が漸近安定であることは、AAの複素固有値の実部がすべて負であるための必要十分条件であり、そうでないときx∗x_*は不安定である。

証明.主張 (1)を示す。∥⋅∥\lVert\cdot\rVertをRn\R^nのノルム、0<γ<−α0<\gamma<-\alphaとし、行列には対応する作用素ノルムを用いる。γ<γ′<−α\gamma<\gamma'<-\alphaを満たすγ′\gamma'を取る。補題 3.2により、M≥1M\ge1が存在して、すべてのt≥0t\ge0で∥etA∥≤Me−γ′t\lVert e^{tA}\rVert\le Me^{-\gamma't}である。η>0\eta>0をMη≤γ′−γM\eta\le\gamma'-\gammaを満たすように取る。注意 4.1のrrについて、ρ0>0\rho_0>0が存在して、閉球B‾={y∈Rn:∥y−x∗∥≤ρ0}\overline B=\{y\in\R^n:\lVert y-x_*\rVert\le\rho_0\}はWWに含まれ、∥u∥≤ρ0\lVert u\rVert\le\rho_0ならば∥r(u)∥≤η∥u∥\lVert r(u)\rVert\le\eta\lVert u\rVertである。δ=ρ0/M\delta=\rho_0/Mとおく。

x∈Ux\in U、∥x−x∗∥<δ\lVert x-x_*\rVert<\deltaとし、u(t)=φt(x)−x∗u(t)=\varphi_t(x)-x_*(t∈I(x)t\in I(x))、β+=sup⁡I(x)\beta_+=\sup I(x)とおく。[0,t][0,t]上で∥u∥<ρ0\lVert u\rVert<\rho_0が成り立つt∈[0,β+)t\in[0,\beta_+)全体の上限をτ\tauとする。∥u(0)∥<δ≤ρ0\lVert u(0)\rVert<\delta\le\rho_0とuuの連続性によりτ>0\tau>0である。[0,τ)[0,\tau)上でuuはC1C^1級であり、式 (4.1.1)によりu′=Au+bu'=Au+b、b(t)=r(u(t))b(t)=r(u(t))を満たし、bbは[0,τ)[0,\tau)上で連続である。§E13.12 系 2.5を区間[0,τ)[0,\tau)とt0=0t_0=0に適用すると、t∈[0,τ)t\in[0,\tau)について

u(t)=etAu(0)+∫0te(t−s)Ar(u(s)) dsu(t)=e^{tA}u(0)+\int_0^te^{(t-s)A}r(u(s))\,ds

である。t∈[0,τ)t\in[0,\tau)で∥u(t)∥<ρ0\lVert u(t)\rVert<\rho_0であるから、w(t)=eγ′t∥u(t)∥w(t)=e^{\gamma't}\lVert u(t)\rVertは

w(t)≤M∥u(0)∥+∫0tMη w(s) ds(0≤t<τ)w(t)\le M\lVert u(0)\rVert+\int_0^tM\eta\,w(s)\,ds \qquad(0\le t<\tau)

を満たす。0<t1<τ0<t_1<\tauを取り、§E13.7 系 1.2を区間[0,t1][0,t_1]、定数M∥u(0)∥M\lVert u(0)\rVertとMηM\eta、零関数に適用すると、0≤t≤t10\le t\le t_1でw(t)≤M∥u(0)∥eMηtw(t)\le M\lVert u(0)\rVert e^{M\eta t}である。t1t_1は任意であるから、t∈[0,τ)t\in[0,\tau)について

∥u(t)∥≤M∥u(0)∥e−(γ′−Mη)t≤Me−γt∥u(0)∥\lVert u(t)\rVert \le M\lVert u(0)\rVert e^{-(\gamma'-M\eta)t} \le Me^{-\gamma t}\lVert u(0)\rVert(4.5.1)\text{(4.5.1)}

である。τ<β+\tau<\beta_+と仮定すると、uuの連続性と式 (4.5.1)により∥u(τ)∥≤M∥u(0)∥<Mδ=ρ0\lVert u(\tau)\rVert\le M\lVert u(0)\rVert<M\delta=\rho_0であり、τ\tauを含む開区間上で∥u∥<ρ0\lVert u\rVert<\rho_0であるから、τ\tauより大きいt∈[0,β+)t\in[0,\beta_+)で[0,t][0,t]上∥u∥<ρ0\lVert u\rVert<\rho_0が成り立つ。これはτ\tauが上限であることに反するから、τ=β+\tau=\beta_+である。したがって、すべてのt∈I(x)t\in I(x)、t≥0t\ge0についてφt(x)∈B‾\varphi_t(x)\in\overline Bである。B‾\overline BはUUに含まれるコンパクト集合であるから、§E13.15 補題 3.1 (1)により[0,∞)⊆I(x)[0,\infty)\subseteq I(x)であり、式 (4.5.1)はすべてのt≥0t\ge0で成り立つ。ε>0\varepsilon>0について、∥x−x∗∥<min⁡{δ,ε/M}\lVert x-x_*\rVert<\min\{\delta,\varepsilon/M\}ならば、すべてのt≥0t\ge0で∥φt(x)−x∗∥<ε\lVert\varphi_t(x)-x_*\rVert<\varepsilonである。また∥x−x∗∥<δ\lVert x-x_*\rVert<\deltaならばφt(x)→x∗\varphi_t(x)\to x_*(t→∞)(t\to\infty)である。したがってx∗x_*は漸近安定である。

主張 (2)を示す。AAの複素固有値で実部が正のもの全体にわたる実部の最小値をλ+\lambda_+とし、c=λ+/2c=\lambda_+/2、A′=A−cIA'=A-cIとおく。μ∈C\mu\in\CがAAの固有値であることとμ−c\mu-cがA′A'の固有値であることは同値であるから、A′A'の固有値の実部は、実部が正のAAの固有値に対してはcc以上であり、それ以外のAAの固有値に対しては−c-c以下である。特にA′A'は実部00の固有値をもたず、実部が正の固有値をもつ。Rn\R^nのノルム∥⋅∥0\lVert\cdot\rVert_0を一つ固定し、行列には対応する作用素ノルムを用いる。§E13.12 補題 6.7をK=R\mathbb K=\RとA′A'に適用してPs,Pu,Es,EuP_s,P_u,E^s,E^uを取り、§E13.12 補題 6.7 (5)によりK≥1K\ge1とκ>0\kappa>0を取って、t≥0t\ge0で∥etA′Ps∥0≤Ke−κt\lVert e^{tA'}P_s\rVert_0\le Ke^{-\kappa t}、t≤0t\le0で∥etA′Pu∥0≤Keκt\lVert e^{tA'}P_u\rVert_0\le Ke^{\kappa t}とする。§E13.12 補題 6.7 (3)により、Rn=Es⊕Eu\R^n=E^s\oplus E^uであり、Eu≠{0}E^u\ne\{0\}である。補題 4.3を、A′A'、∥⋅∥0\lVert\cdot\rVert_0、KKと、同補題のα,β\alpha,\betaをそれぞれκ,κ/2\kappa,\kappa/2として適用し、EsE^sのノルム∥⋅∥s\lVert\cdot\rVert_s、EuE^uのノルム∥⋅∥u\lVert\cdot\rVert_uとRn\R^nのノルム∥⋅∥∗\lVert\cdot\rVert_*を取る。

cIcIとA′A'は可換であるから、§E13.12 命題 3.2により、すべてのt∈Rt\in\RでetA=ectetA′e^{tA}=e^{ct}e^{tA'}である。B=eAB=e^{A}とおく。補題 4.3 (2)をt=1t=1に適用すると、BBはEsE^sとEuE^uを不変にし、v∈Esv\in E^sならば

∥Bv∥s=ec∥eA′v∥s≤ec−κ/2∥v∥s\lVert Bv\rVert_s=e^{c}\lVert e^{A'}v\rVert_s\le e^{c-\kappa/2}\lVert v\rVert_s

であり、v∈Euv\in E^uならば

∥e−Av∥u=e−c∥e−A′v∥u≤e−(c+κ/2)∥v∥u\lVert e^{-A}v\rVert_u=e^{-c}\lVert e^{-A'}v\rVert_u\le e^{-(c+\kappa/2)}\lVert v\rVert_u

である。§E13.12 定理 2.1によりB∣EuB|_{E^u}は逆写像e−A∣Eue^{-A}|_{E^u}をもつ。a=ec−κ/2a=e^{c-\kappa/2}、b=ec+κ/2b=e^{c+\kappa/2}とおくと、c>0c>0かつκ>0\kappa>0であるから0≤a<b0\le a<bかつb>1b>1である。§E13.12 補題 6.7 (4)により、zs∈Esz^s\in E^s、zu∈Euz^u\in E^uについてPs(zu+zs)=zsP_s(z^u+z^s)=z^s、Pu(zu+zs)=zuP_u(z^u+z^s)=z^uであるから、∥zu+zs∥∗=max⁡{∥zu∥u,∥zs∥s}\lVert z^u+z^s\rVert_*=\max\{\lVert z^u\rVert_u,\lVert z^s\rVert_s\}である。したがって、E=RnE=\R^nに∥⋅∥∗\lVert\cdot\rVert_*を、EsE^sとEuE^uにそれぞれ∥⋅∥s\lVert\cdot\rVert_sと∥⋅∥u\lVert\cdot\rVert_uを与えると、BBとa,ba,bは式 (4.4.1)を満たす。

F∣WF|_Wの流れをψ\psi、その最大存在区間をIW(y)I_W(y)とし、W1={y∈W:1∈IW(y)}W_1=\{y\in W:1\in I_W(y)\}とおく。y∈Wy\in Wについて、すべてのs∈IW(y)s\in I_W(y)でψs(y)∈W\psi_s(y)\in Wであるから、§E13.15 補題 7.3 (2)をF∣WF|_W、N=WN=W、FFに適用すると、IW(y)⊆I(y)I_W(y)\subseteq I(y)であり、IW(y)I_W(y)上でψs(y)=φs(y)\psi_s(y)=\varphi_s(y)である。§E13.15 定理 2.1 (2)によりW1W_1は開集合であり、§E13.15 系 2.2 (4)によりIW(x∗)=RI_W(x_*)=\Rであるからx∗∈W1x_*\in W_1である。N0={z∈Rn:x∗+z∈W1}N_0=\{z\in\R^n:x_*+z\in W_1\}とおき、G:N0→RnG:N_0\to\R^nをG(z)=ψ1(x∗+z)−x∗G(z)=\psi_1(x_*+z)-x_*で定める。§E13.15 系 2.2 (3)によりGGはC1C^1級であり、§E13.15 系 2.2 (4)によりG(0)=0G(0)=0、DG(0)=eADG(0)=e^{A}である。g(z)=G(z)−Bzg(z)=G(z)-Bzとおくと、g(0)=0g(0)=0であり、微分可能性の定義により∥g(z)∥∗/∥z∥∗→0\lVert g(z)\rVert_*/\lVert z\rVert_*\to0(z→0)(z\to0)である。

主張 4.5.1.k≥1k\ge1とし、z0,…,zk−1∈N0z_0,\ldots,z_{k-1}\in N_0がzj=G(zj−1)z_j=G(z_{j-1})(1≤j≤k1\le j\le k)を満たすとする。このとき、y=x∗+z0y=x_*+z_0についてk∈I(y)k\in I(y)かつφk(y)=x∗+zk\varphi_k(y)=x_*+z_kである。

証明.y=x∗+z0y=x_*+z_0とおく。0≤j≤k−10\le j\le k-1についてj∈I(y)j\in I(y)かつφj(y)=x∗+zj\varphi_j(y)=x_*+z_jであるとすると、x∗+zj∈W1x_*+z_j\in W_1であるから1∈IW(x∗+zj)⊆I(x∗+zj)1\in I_W(x_*+z_j)\subseteq I(x_*+z_j)であり、φ1(x∗+zj)=ψ1(x∗+zj)=x∗+zj+1\varphi_1(x_*+z_j)=\psi_1(x_*+z_j)=x_*+z_{j+1}である。§E13.15 定理 1.3によりI(x∗+zj)=I(y)−jI(x_*+z_j)=I(y)-jであるからj+1∈I(y)j+1\in I(y)であり、φj+1(y)=φ1(φj(y))=x∗+zj+1\varphi_{j+1}(y)=\varphi_1(\varphi_j(y))=x_*+z_{j+1}である。j=0j=0では仮定は成り立つから、帰納法によりk∈I(y)k\in I(y)かつφk(y)=x∗+zk\varphi_k(y)=x_*+z_kである。▨

x∗x_*が Lyapunov 安定であると仮定する。補題 2.3により、x∗x_*はノルム∥⋅∥∗\lVert\cdot\rVert_*について Lyapunov 安定である。補題 4.4をE=RnE=\R^n、BB、N0N_0、GGに適用してr>0r>0を取り、Lyapunov 安定性の定義によりε=r\varepsilon=rに対するδ>0\delta>0を取る。補題 4.4により、∥z0∥∗<δ\lVert z_0\rVert_*<\deltaを満たすz0z_0とk≥1k\ge1が存在して、0≤j≤k−10\le j\le k-1で∥zj∥∗<r\lVert z_j\rVert_*<rであり、∥zk∥∗≥r\lVert z_k\rVert_*\ge rである。{z:∥z∥∗<r}⊆N0\{z:\lVert z\rVert_*<r\}\subseteq N_0であるからz0,…,zk−1∈N0z_0,\ldots,z_{k-1}\in N_0であり、主張 4.5.1によりy=x∗+z0∈Uy=x_*+z_0\in Uについてk∈I(y)k\in I(y)かつ∥φk(y)−x∗∥∗=∥zk∥∗≥r\lVert\varphi_k(y)-x_*\rVert_*=\lVert z_k\rVert_*\ge rである。一方、∥y−x∗∥∗<δ\lVert y-x_*\rVert_*<\deltaであるから、Lyapunov 安定性により∥φk(y)−x∗∥∗<r\lVert\varphi_k(y)-x_*\rVert_*<rである。この二つは両立しないから、x∗x_*は Lyapunov 安定でない。

主張 (3)を示す。x∗x_*が双曲型平衡点であるとする。AAの複素固有値の実部がすべて負ならば、主張 (1)によりx∗x_*は漸近安定である。そうでないならば、AAは実部00の固有値をもたないから実部が正の固有値をもち、主張 (2)によりx∗x_*は不安定であって、特に漸近安定でない。▨

注意 4.6.定理 4.5の設定で、A=DF(x∗)A=DF(x_*)の複素固有値の実部がすべて00以下であり、実部が00の固有値が存在する場合には、定理 4.5 (1)と定理 4.5 (2)のいずれの仮定も満たされない。この場合には、AAが同じでもx∗x_*の安定性が異なる系がある。たとえば、問題 5.5 (I)、問題 5.5 (II)、問題 5.5 (III)の三つの系は、いずれも原点を平衡点とし、原点で同じ Jacobi 行列(0−110)\begin{pmatrix}0&-1\\1&0\end{pmatrix}をもつが、それぞれの原点は、Lyapunov 安定であって漸近安定でない平衡点、漸近安定な平衡点、不安定な平衡点である。これらの判定は演習とする(問題 5.5)。

5 演習

問題 5.1 (零固有値と Jordan 構造). 二つの行列

A0=(0000),N=(0100)A_0=\begin{pmatrix}0&0\\0&0\end{pmatrix}, \qquad N=\begin{pmatrix}0&1\\0&0\end{pmatrix}

について、線形系x′=A0xx'=A_0xとx′=Nxx'=Nxの解を求めよ。両行列の固有値はともに0,00,0であるが、前者の原点は Lyapunov 安定であって漸近安定ではなく、後者の原点は不安定であることを証明せよ。

解答.

A0A_0とNNの特性多項式はともにz2z^2である。初期値を(u0,v0)T(u_0,v_0)^{\mathsf T}とすると、命題 3.1 (1)により、A0A_0に対する解はすべての実数ttに対して

u(t)=u0,v(t)=v0u(t)=u_0,\qquad v(t)=v_0

である。Euclid ノルムを用いる。任意のε>0\varepsilon>0に対してδ=ε\delta=\varepsilonと取れば、原点からの距離はε\varepsilon未満に保たれるため、原点は Lyapunov 安定である。非零の初期値から出る解は原点へ収束しないため、原点は漸近安定ではない。

N2=0N^2=0であるから、§E13.12 定理 5.1によりetN=I+tNe^{tN}=I+tNであり、解は

u(t)=u0+tv0,v(t)=v0(t∈R)u(t)=u_0+tv_0,\qquad v(t)=v_0 \qquad(t\in\R)

となる。ε=1\varepsilon=1と固定する。任意のδ>0\delta>0に対して0<a<min⁡{δ,1}0<a<\min\{\delta,1\}を取り、(u0,v0)=(0,a)(u_0,v_0)=(0,a)と置くと、初期値のノルムはδ\delta未満であるが、t=2/at=2/aではu(t)=2u(t)=2となり、解のノルムは11より大きい。したがって、NNに対する原点は不安定である。▨

問題 5.2 (重複固有値をもつ安定な線形系の軌道).λ<0\lambda<0とし、

J=(λ10λ)J=\begin{pmatrix}\lambda&1\\0&\lambda\end{pmatrix}

と置く。JJに対する線形系u′=λu+vu'=\lambda u+v、v′=λvv'=\lambda vを考える。初期値(u0,v0)(u_0,v_0)に対する解を求め、v0≠0v_0\ne0ならば十分大きいttでu(t)≠0u(t)\ne0となり、v(t)/u(t)→0v(t)/u(t)\to0が成り立つことを示せ。また、行列λI\lambda Iに対する系と比較し、両系の原点が漸近安定であることと、λI\lambda Iに対する非定数解の軌道が原点を通る直線上にあることを示せ。

解答.

§E13.12 定理 5.1から

u(t)=eλt(u0+tv0),v(t)=eλtv0(t∈R)u(t)=e^{\lambda t}(u_0+tv_0),\qquad v(t)=e^{\lambda t}v_0 \qquad(t\in\R)

を得る。v0≠0v_0\ne0のとき、t>∣u0/v0∣t>|u_0/v_0|ならばt+u0/v0>0t+u_0/v_0>0であるからu(t)≠0u(t)\ne0であり、

v(t)u(t)=1t+u0/v0⟶0(t→∞)\frac{v(t)}{u(t)}=\frac{1}{t+u_0/v_0}\longrightarrow0 \qquad(t\to\infty)

となる。v0=0v_0=0の場合はv(t)=0v(t)=0であり、解は固有直線v=0v=0上にある。

行列λI\lambda Iに対する解は

(u(t),v(t))=eλt(u0,v0)(t∈R)(u(t),v(t))=e^{\lambda t}(u_0,v_0) \qquad(t\in\R)

である。非零の初期値に対して、軌道は初期値と原点を結ぶ直線上にあり、方向は一定である。一方、JJに対するv0≠0v_0\ne0の解では、比v(t)/u(t)v(t)/u(t)が変化する。この解は、原点へ近づくにつれて固有直線v=0v=0の方向へ近づく。λI\lambda Iに対する原点は星形結節点であり、JJに対する原点は非正則結節点である。両行列の重複固有値λ\lambdaは負であるから、定理 3.5 (3)により両系の原点は漸近安定である。▨

問題 5.3 (Jacobi 行列による局所安定性の判定). 非線形系

x′=−x+xy,y′=−2y+x2\begin{aligned} x'&=-x+xy,\\ y'&=-2y+x^2 \end{aligned}

について、原点が平衡点であることを確認し、Jacobi 行列の固有値から原点の局所安定性を判定せよ。また、原点以外の平衡点をすべて求め、原点がすべての初期値から出る解を吸引するわけではないことを示せ。

解答.

右辺をF(x,y)F(x,y)と書くとF(0,0)=0F(0,0)=0であるから、原点は平衡点である。FFはR2\R^2上でC1C^1級であり、各点を中心とする閉球の上でDFDFは有界であるから、§E4.3 定理 2.1によりFFは局所 Lipschitz 連続である。Jacobi 行列は

DF(x,y)=(−1+yx2x−2)DF(x,y) = \begin{pmatrix} -1+y&x\\ 2x&-2 \end{pmatrix}

であり、

DF(0,0)=(−100−2)DF(0,0) = \begin{pmatrix} -1&0\\ 0&-2 \end{pmatrix}

の固有値は−1,−2-1,-2である。二つの固有値の実部がともに負であるから、定理 4.5 (1)により原点は漸近安定である。

平衡点はx(−1+y)=0x(-1+y)=0と−2y+x2=0-2y+x^{2}=0の共通解である。x=0x=0ならば第二式からy=0y=0である。x≠0x\ne0ならば第一式からy=1y=1であり、第二式からx2=2x^{2}=2、すなわちx=±2x=\pm\sqrt2である。したがって原点以外の平衡点は(2,1)(\sqrt2,1)と(−2,1)(-\sqrt2,1)だけである。これらの点から出る解は定数解であり、原点へ収束しない。ゆえに原点はすべての初期値から出る解を吸引するわけではない。▨

問題 5.4.問題 5.3の系の平衡点(2,1)(\sqrt2,1)と(−2,1)(-\sqrt2,1)が不安定であることを証明せよ。また、系

x′=x+y2,y′=−y3x'=x+y^2,\qquad y'=-y^3

について、原点がただ一つの平衡点であり、双曲型平衡点ではないが不安定であることを証明せよ。

解答.

問題 5.3の解答のとおり、その系の右辺FFは局所 Lipschitz 連続であり、R2\R^2上でC1C^1級であって、

DF(±2,1)=(0±2±22−2)DF(\pm\sqrt2,1) = \begin{pmatrix} 0&\pm\sqrt2\\ \pm2\sqrt2&-2 \end{pmatrix}

(複号同順)である。その特性多項式はz(z+2)−4=z2+2z−4z(z+2)-4=z^2+2z-4であるから、固有値は−1±5-1\pm\sqrt5であり、−1+5>0-1+\sqrt5>0である。定理 4.5 (2)により、(2,1)(\sqrt2,1)と(−2,1)(-\sqrt2,1)はいずれも不安定である。

後半の系の右辺をH(x,y)=(x+y2,−y3)H(x,y)=(x+y^2,-y^3)と書く。HHはR2\R^2上の多項式写像であるからC1C^1級であり、各点を中心とする閉球の上でDHDHは有界であるから、§E4.3 定理 2.1により局所 Lipschitz 連続である。H(x,y)=0H(x,y)=0ならばy3=0y^3=0からy=0y=0であり、x+y2=0x+y^2=0からx=0x=0である。したがって原点がただ一つの平衡点である。

DH(x,y)=(12y0−3y2),DH(0,0)=(1000)DH(x,y)= \begin{pmatrix} 1&2y\\ 0&-3y^2 \end{pmatrix}, \qquad DH(0,0)= \begin{pmatrix} 1&0\\ 0&0 \end{pmatrix}

の固有値は11と00である。固有値00の実部は00であるから、原点は双曲型平衡点ではなく、定理 4.5 (3)の仮定は満たされない。一方、DH(0,0)DH(0,0)は実部が正の固有値11をもつから、定理 4.5 (2)により原点は不安定である。▨

問題 5.5 (非双曲型平衡点で線形化だけでは安定性が定まらないことの証明).r2=x2+y2r^2=x^2+y^2と書き、次の三つの系を考える。

  1. x′=−y,y′=x,x'=-y,\qquad y'=x,
  2. x′=−y−xr2,y′=x−yr2,x'=-y-xr^2,\qquad y'=x-yr^2,
  3. x′=−y+xr2,y′=x+yr2.x'=-y+xr^2,\qquad y'=x+yr^2.

三つの系が原点で同じ Jacobi 行列をもち、その固有値の実部が零であることを示せ。そのうえで、(I)の系の原点は Lyapunov 安定であるが漸近安定ではなく、(II)の系の原点は漸近安定であり、(III)の系の原点は不安定であることを証明せよ。

解答.

三つの系の右辺はR2\R^2上の多項式写像であり、局所 Lipschitz 連続である。いずれの系でも原点は平衡点であり、xr2xr^2、yr2yr^2の偏導関数は原点で00であるから、原点における Jacobi 行列は

A=(0−110)A=\begin{pmatrix}0&-1\\1&0\end{pmatrix}

である。固有値は±i\pm iであり、実部は零である。Euclid ノルムを用いる。

各系の流れをφ\varphiとし、z∈R2∖{0}z\in\R^2\setminus\{0\}についてs(t)=∥φt(z)∥22s(t)=\lVert\varphi_t(z)\rVert_2^2(t∈I(z)t\in I(z))とおく。φt(z)=0\varphi_t(z)=0となるttが存在するならば、§E13.15 命題 4.3 (2)と§E13.15 命題 4.3 (1)によりγ(z)={0}\gamma(z)=\{0\}となってz≠0z\ne0に反するから、I(z)I(z)上でs>0s>0である。φt(z)=(x,y)\varphi_t(z)=(x,y)と書くと、s′=2xx′+2yy′s'=2xx'+2yy'である。

(I)の系ではs′=2x(−y)+2yx=0s'=2x(-y)+2yx=0であるから、§D1.14 定理 3.1によりssはI(z)I(z)上で定数であり、前方軌道は半径∥z∥2\lVert z\rVert_2の円周に含まれる。§E13.15 補題 3.1 (1)により[0,∞)⊆I(z)[0,\infty)\subseteq I(z)であり、すべてのt≥0t\ge0で∥φt(z)∥2=∥z∥2\lVert\varphi_t(z)\rVert_2=\lVert z\rVert_2である。任意のε>0\varepsilon>0に対してδ=ε\delta=\varepsilonと取ると Lyapunov 安定性の条件が成り立つ。z≠0z\ne0から出る解は原点へ収束しないから、原点は漸近安定でない。

(II)の系では

s′=2x(−y−xs)+2y(x−ys)=−2s2s'=2x(-y-xs)+2y(x-ys)=-2s^2

であるから、(1/s)′=−s′/s2=2(1/s)'=-s'/s^2=2であり、t∈I(z)t\in I(z)について1/s(t)=1/s(0)+2t1/s(t)=1/s(0)+2tである。したがってt∈I(z)t\in I(z)、t≥0t\ge0についてs(t)≤s(0)s(t)\le s(0)であり、前方軌道は半径∥z∥2\lVert z\rVert_2の閉円板に含まれる。§E13.15 補題 3.1 (1)により[0,∞)⊆I(z)[0,\infty)\subseteq I(z)であり、t≥0t\ge0について

s(t)=s(0)1+2s(0)t≤s(0),s(t)→0(t→∞)s(t)=\frac{s(0)}{1+2s(0)t}\le s(0),\qquad s(t)\to0\quad(t\to\infty)

である。任意のε>0\varepsilon>0に対してδ=ε\delta=\varepsilonと取ると Lyapunov 安定性の条件が成り立ち、すべてのzzについてφt(z)→0\varphi_t(z)\to0であるから、原点は漸近安定である。

(III)の系では

s′=2x(−y+xs)+2y(x+ys)=2s2s'=2x(-y+xs)+2y(x+ys)=2s^2

であるから、(1/s)′=−2(1/s)'=-2であり、t∈I(z)t\in I(z)について1/s(t)=1/s(0)−2t1/s(t)=1/s(0)-2tである。ε=1\varepsilon=1とし、δ>0\delta>0を任意に取る。0<∥z∥2<min⁡{δ,1}0<\lVert z\rVert_2<\min\{\delta,1\}を満たすzzを取り、t1=(1/s(0)−1)/2t_1=(1/s(0)-1)/2とおくと、s(0)<1s(0)<1であるからt1>0t_1>0である。t1∉I(z)t_1\notin I(z)ならば[0,∞)⊈I(z)[0,\infty)\not\subseteq I(z)である。t1∈I(z)t_1\in I(z)ならば1/s(t1)=11/s(t_1)=1であり、∥φt1(z)∥2=1\lVert\varphi_{t_1}(z)\rVert_2=1である。いずれの場合も、zzについて Lyapunov 安定性の条件はε=1\varepsilon=1に対して成り立たない。δ\deltaは任意であったから、原点は不安定である。▨

問題 5.6.

A=(−1−41−1),S=(1002)A=\begin{pmatrix}-1&-4\\1&-1\end{pmatrix},\qquad S=\begin{pmatrix}1&0\\0&2\end{pmatrix}

とし、F(x)=AxF(x)=Axで定まるF:R2→R2F:\R^2\to\R^2の流れをφ\varphiとする。SAS−1SAS^{-1}を求め、命題 3.1 (2)と系 2.5を用いて、原点がFFの漸近安定な平衡点であることを証明せよ。また、Euclid ノルムを∥⋅∥2\lVert\cdot\rVert_2とし、p=(1,−1)p=(1,-1)について、関数t↦∥φt(p)∥2t\mapsto\lVert\varphi_t(p)\rVert_2が00を含むある開区間で狭義単調増加であることを示せ。

解答.

S−1=(1001/2)S^{-1}=\begin{pmatrix}1&0\\0&1/2\end{pmatrix}であるから

SAS−1=(1002)(−1−21−1/2)=(−1−22−1)SAS^{-1} =\begin{pmatrix}1&0\\0&2\end{pmatrix} \begin{pmatrix}-1&-2\\1&-1/2\end{pmatrix} =\begin{pmatrix}-1&-2\\2&-1\end{pmatrix}

である。B=SAS−1B=SAS^{-1}、J=(0−110)J=\begin{pmatrix}0&-1\\1&0\end{pmatrix}とおくとB=−I+2JB=-I+2Jであり、−I-Iと2J2Jは可換であるから、§E13.12 命題 3.2と§E13.12 例 4.4(ω=2\omega=2)によりetB=e−te2tJe^{tB}=e^{-t}e^{2tJ}はe−te^{-t}と角2t2tの回転行列の積である。命題 3.1 (1)によりG(y)=ByG(y)=Byの流れはψt(y)=etBy\psi_t(y)=e^{tB}yであり、すべてのt∈Rt\in\Rで∥ψt(y)∥2=e−t∥y∥2\lVert\psi_t(y)\rVert_2=e^{-t}\lVert y\rVert_2である。任意のε>0\varepsilon>0に対してδ=ε\delta=\varepsilonと取ると Lyapunov 安定性の条件が成り立ち、すべてのyyについてψt(y)→0\psi_t(y)\to0であるから、原点はGGの漸近安定な平衡点である。命題 3.1 (2)によりh(x)=Sxh(x)=SxはFFとGGの流れの位相共役であり、h(0)=0h(0)=0であるから、系 2.5により原点はFFの漸近安定な平衡点である。

q(t)=∥φt(p)∥22q(t)=\lVert\varphi_t(p)\rVert_2^2とおく。φt(p)=(x1,x2)\varphi_t(p)=(x_1,x_2)と書くと、qqはR\R上でC1C^1級であり、

q′(t)=2⟨φt(p),Aφt(p)⟩=2(−x12−3x1x2−x22)q'(t)=2\langle\varphi_t(p),A\varphi_t(p)\rangle =2(-x_1^2-3x_1x_2-x_2^2)

である。t=0t=0では(x1,x2)=(1,−1)(x_1,x_2)=(1,-1)であるからq′(0)=2(−1+3−1)=2>0q'(0)=2(-1+3-1)=2>0である。q′q'は連続であるから、00を含む開区間J0J_0が存在してJ0J_0上でq′>0q'>0であり、§D1.14 定理 3.1によりqqはJ0J_0上で狭義単調増加である。∥φt(p)∥2=q(t)\lVert\varphi_t(p)\rVert_2=\sqrt{q(t)}であるから、t↦∥φt(p)∥2t\mapsto\lVert\varphi_t(p)\rVert_2もJ0J_0上で狭義単調増加である。一方、Sφt(p)=h(φt(p))=ψt(Sp)S\varphi_t(p)=h(\varphi_t(p))=\psi_t(Sp)であるから、∥Sφt(p)∥2=e−t∥Sp∥2\lVert S\varphi_t(p)\rVert_2=e^{-t}\lVert Sp\rVert_2は狭義単調減少である。▨

参考文献

  1. Lawrence Perko, Differential Equations and Dynamical Systems, 3rd ed., Texts in Applied Mathematics 7, Springer, New York, 2001.平面線形系の分類と双曲型平衡点の線形化による安定性判定を参考にした。
  2. Gerald Teschl, Ordinary Differential Equations and Dynamical Systems, American Mathematical Society, Providence, Rhode Island, 2012.相線、相平面、Lyapunov 安定性、線形化された安定性を参考にした。

前提記事

10 本の記事・単元を表示