1 相線と相平面
定義 1.1.U⊆Rを開集合、f:U→Rを局所 Lipschitz 連続な関数とし、一階自律方程式x′=f(x)を考える。Uの各点xにf(x)の符号を記したものを、fの相線 (phase line) という。
補題 1.2.U⊆Rを開集合、f:U→Rを局所 Lipschitz 連続な関数、φをfの流れとする。x∈Uはf(x)=0を満たすとし、xを含みfの零点を含まないUの部分区間すべての和集合をCとする。このとき、次が成り立つ。
- Cは開区間であり、fはC上でf(x)と同じ符号をとる。さらにγ(x)⊆Cである。
- f(x)>0ならばt↦φt(x)はI(x)上で狭義単調増加であり、f(x)<0ならばI(x)上で狭義単調減少である。
- コンパクト集合K⊆Uが存在して、すべてのt∈I(x)、t≥0についてφt(x)∈Kであるならば、[0,∞)⊆I(x)であり、極限ℓ=limt→∞φt(x)が存在し、ℓ∈Kはfの平衡点である。
証明.主張 (1)を示す。xを含みfの零点を含まないUの部分区間はすべて点xを共有するから、その和集合Cは区間である。y∈Cとすると、f(y)=0であり、Uは開集合でfは連続であるから、ε>0が存在して、(y−ε,y+ε)⊆Uかつこの区間上でf=0である。この区間と、xとyを含みfの零点を含まないUの部分区間との和は、xを含みfの零点を含まないUの部分区間であるから、Cに含まれる。したがってCは開区間である。a,b∈Cについてf(a)とf(b)の符号が異なるならば、aとbを端点とする閉区間はCに含まれ、§D1.12 系 1.2によりその閉区間にfの零点が存在する。これはCがfの零点を含まないことに反するから、fはC上でf(x)と同じ符号をとる。
t∈I(x)とし、y=φt(x)とおく。f(y)=0ならば、§E13.15 命題 4.3 (2)によりγ(y)={y}であり、y∈γ(x)であるから§E13.15 命題 4.3 (1)によりγ(x)=γ(y)={y}である。このときx=yであり、これはf(x)=0に反する。したがってγ(x)はfの零点を含まない。t1<t2をI(x)の点とすると、§D1.12 系 1.2を[t1,t2]上の連続関数s↦φs(x)に適用することにより、φt1(x)とφt2(x)の間の値はすべてγ(x)に属する。したがってγ(x)はxを含みfの零点を含まないUの部分区間であり、γ(x)⊆Cである。
主張 (2)を示す。t1<t2をI(x)の点とする。s↦φs(x)はI(x)上で微分可能であり、その導関数はs↦f(φs(x))である。§D1.14 定理 3.1により、c∈(t1,t2)が存在してφt2(x)−φt1(x)=f(φc(x))(t2−t1)である。主張 (1)によりφc(x)∈Cであり、f(φc(x))はf(x)と同じ符号をもつ。
主張 (3)を示す。§E13.15 定理 5.3 (2)により、[0,∞)⊆I(x)であり、ω(x)は空でなく、ω(x)⊆Kである。f(x)>0ならばσ=1、f(x)<0ならばσ=−1とおくと、主張 (2)によりs↦σφs(x)はI(x)上で狭義単調増加である。ℓ∈ω(x)と、φtm(x)→ℓを満たす数列tm→∞を取る。s≥0について、tm>sならばσφtm(x)>σφs(x)であるから、m→∞としてσℓ≥σφs(x)である。ε>0に対してmをσℓ−σφtm(x)<εを満たすように取ると、t≥tmについてσφtm(x)≤σφt(x)≤σℓであるから、∣φt(x)−ℓ∣<εである。したがってφt(x)→ℓ(t→∞)である。ω(x)の任意の元はある数列tm′→∞に沿うφtm′(x)の極限であるからℓに等しく、ω(x)={ℓ}である。§E13.15 定理 5.3 (1)によりω(x)は不変集合であるからγ(ℓ)⊆{ℓ}であり、§E13.15 命題 4.3 (2)によりf(ℓ)=0である。ℓ∈ω(x)⊆Kである。▨
定義 1.3.U⊆R2を開集合、F:U→R2を局所 Lipschitz 連続な写像とし、自律系x′=F(x)を考える。状態x=(x1,x2)を点として表す平面を相平面 (phase plane) という。各点x∈Uに速度ベクトルF(x)を対応させたものを、Fの定めるベクトル場 (vector field) という。
- Fの流れ、軌道、平衡点は§E13.15 定義 4.1のn=2の場合であり、相平面にはx∈Uを通る軌道γ(x)をtが増加する向きを付けて描く。
2 平衡点の安定性
定義 2.1.n∈N≥1とし、U⊆Rnを開集合、F:U→Rnを局所 Lipschitz 連続な写像とする。
- Fの平衡点は§E13.15 定義 4.1の平衡点である。
- x∗∈UをFの平衡点とする。x∗の開近傍W⊆Uで、Wに含まれるFの平衡点がx∗だけであるものが存在するとき、x∗をFの孤立平衡点 (isolated equilibrium point) という。
定義 2.2.n∈N≥1とし、U⊆Rnを開集合、F:U→Rnを局所 Lipschitz 連続な写像、φをFの流れ、x∗∈UをFの平衡点とする。Rnのノルム∥⋅∥を一つ固定する。
x∗がLyapunov 安定 (Lyapunov stable) であるとは、任意のε>0に対してδ>0が存在し、∥x−x∗∥<δを満たす任意のx∈Uについて、[0,∞)⊆I(x)であり、すべてのt≥0で∥φt(x)−x∗∥<εが成り立つことをいう。
x∗が漸近安定 (asymptotically stable) であるとは、x∗が Lyapunov 安定であり、さらにρ>0が存在して、∥x−x∗∥<ρを満たす任意のx∈Uについて、[0,∞)⊆I(x)かつφt(x)→x∗(t→∞)が成り立つことをいう。Lyapunov 安定でない平衡点を不安定 (unstable) という。
補題 2.3.n∈N≥1とし、U⊆Rnを開集合、F:U→Rnを局所 Lipschitz 連続な写像、x∗∈UをFの平衡点とする。∥⋅∥1、∥⋅∥2をRnの二つのノルムとする。このとき、x∗が∥⋅∥1について Lyapunov 安定であることは、∥⋅∥2について Lyapunov 安定であるための必要十分条件であり、x∗が∥⋅∥1について漸近安定であることは、∥⋅∥2について漸近安定であるための必要十分条件である。特に、x∗が不安定であるかどうかはノルムの選び方によらない。
証明.φをFの流れとする。§E10.1 命題 2.2により、定数0<c≤Cが存在して、すべてのv∈Rnでc∥v∥1≤∥v∥2≤C∥v∥1が成り立つ。x∗が∥⋅∥1について Lyapunov 安定であるとし、ε>0を取る。ε/Cに対して∥⋅∥1についてのδ>0を取る。x∈Uが∥x−x∗∥2<cδを満たすならば∥x−x∗∥1<δであるから、[0,∞)⊆I(x)であり、すべてのt≥0で∥φt(x)−x∗∥2≤C∥φt(x)−x∗∥1<εである。したがってx∗は∥⋅∥2について Lyapunov 安定である。さらにx∗が∥⋅∥1について漸近安定であり、ρ>0が漸近安定性の定義の条件を∥⋅∥1について満たすとする。x∈Uが∥x−x∗∥2<cρを満たすならば∥x−x∗∥1<ρであるから、[0,∞)⊆I(x)であり、∥φt(x)−x∗∥2≤C∥φt(x)−x∗∥1→0(t→∞)である。したがってx∗は∥⋅∥2について漸近安定である。∥⋅∥1と∥⋅∥2の役割を入れ替えると、逆の含意を得る。不安定性は Lyapunov 安定性の否定である。▨
補題 2.4.n∈N≥1とし、U,V⊆Rnを開集合、F:U→Rn、G:V→Rnを局所 Lipschitz 連続な写像とし、それらの流れをφ、ψとする。x∗∈UをFの平衡点、y∗∈VをGの平衡点とし、Fはx∗でGのy∗と局所位相共役であるとする。このとき、次が成り立つ。
- x∗が Lyapunov 安定ならば、y∗は Lyapunov 安定である。
- x∗が漸近安定ならば、y∗は漸近安定である。
証明. 局所位相共役の定義の開近傍W⊆U、W′⊆Vと同相写像h:W→W′を取り、Rnのノルムを一つ固定する。p∈Rn、r>0についてB(p,r)={x∈Rn:∥x−p∥<r}と書く。
主張 2.4.1.ε>0とする。このとき、ε1>0が存在して、次が成り立つ。x∈Wが[0,∞)⊆IF(x)を満たし、すべてのs≥0でφs(x)∈B(x∗,ε1)であるならば、[0,∞)⊆IG(h(x))であり、すべてのt≥0でψt(h(x))=h(φt(x))∈B(y∗,ε)である。
証明.Wは開集合で、hはx∗で連続でありh(x∗)=y∗であるから、ε1>0が存在してB(x∗,ε1)⊆Wかつh(B(x∗,ε1))⊆B(y∗,ε)である。xを仮定のとおりに取り、t≥0とすると、0とtの間のすべてのsでφs(x)∈B(x∗,ε1)⊆Wである。局所位相共役の定義により、t∈IG(h(x))であり、ψt(h(x))=h(φt(x))∈B(y∗,ε)である。▨
主張 (1)を示す。x∗が Lyapunov 安定であるとし、ε>0を取る。εに対する主張 2.4.1のε1を取り、Lyapunov 安定性の定義によりε1に対するδ1>0を取る。W′は開集合で、h−1はy∗で連続でありh−1(y∗)=x∗であるから、δ>0が存在してB(y∗,δ)⊆W′かつh−1(B(y∗,δ))⊆B(x∗,δ1)である。y∈V、∥y−y∗∥<δとし、x=h−1(y)とおく。∥x−x∗∥<δ1であるから、[0,∞)⊆IF(x)であり、すべてのs≥0でφs(x)∈B(x∗,ε1)である。主張 2.4.1により、[0,∞)⊆IG(y)であり、すべてのt≥0で∥ψt(y)−y∗∥<εである。
主張 (2)を示す。x∗が漸近安定であるとし、ρ>0を漸近安定性の定義の条件を満たすように取る。主張 (1)によりy∗は Lyapunov 安定である。ε=1に対する主張 2.4.1のε1を取り、ε1に対する Lyapunov 安定性の定義のδ1>0を取って、ρ1=min{δ1,ρ}とおく。h−1はy∗で連続であるから、ρ′>0が存在してB(y∗,ρ′)⊆W′かつh−1(B(y∗,ρ′))⊆B(x∗,ρ1)である。y∈V、∥y−y∗∥<ρ′とし、x=h−1(y)とおく。∥x−x∗∥<δ1であるから、主張 2.4.1により[0,∞)⊆IG(y)であり、すべてのt≥0でψt(y)=h(φt(x))である。∥x−x∗∥<ρであるからφt(x)→x∗であり、hはx∗で連続であるからψt(y)→h(x∗)=y∗(t→∞)である。▨
系 2.5.n∈N≥1とし、U,V⊆Rnを開集合、F:U→Rn、G:V→Rnを局所 Lipschitz 連続な写像、h:U→VをFとGの流れの位相共役とする。x∗∈UをFの平衡点とする。このとき、h(x∗)はGの平衡点である。さらに、x∗が Lyapunov 安定であることは、h(x∗)が Lyapunov 安定であるための必要十分条件であり、x∗が漸近安定であることは、h(x∗)が漸近安定であるための必要十分条件である。
証明.FとGの流れをφ、ψとする。§E13.15 命題 7.2 (1)によりh(x∗)はGの平衡点である。x∈U、t∈IF(x)ならば、位相共役の定義によりt∈IG(h(x))かつh(φt(x))=ψt(h(x))であるから、W=U、W′=VとしてFはx∗でGのh(x∗)と局所位相共役である。y∈V、t∈IG(y)とし、x=h−1(y)とおくと、IF(x)=IG(y)であるからt∈IF(x)であり、h(φt(x))=ψt(y)によりh−1(ψt(y))=φt(h−1(y))である。したがって、W=V、W′=Uとして、Gはh(x∗)でFのx∗と局所位相共役である。この二つの局所位相共役に補題 2.4を適用すると、結論を得る。▨
例 2.6.U=R、f(x)=x(1−x)とする。fはR上でC1級であるから、§E4.3 定理 2.1により局所 Lipschitz 連続である。fの平衡点は0と1であり、fは(0,1)上で正、(−∞,0)と(1,∞)上で負である。fの流れをφとし、Rには絶対値をノルムとして用いる。
0<x<1とする。補題 1.2のCは(0,1)であるから、補題 1.2 (1)と補題 1.2 (2)により、t∈I(x)、t≥0についてφt(x)∈[x,1)である。補題 1.2 (3)をK=[x,1]に適用すると、[0,∞)⊆I(x)であり、φt(x)は[x,1]に属する平衡点に収束する。その平衡点は1である。x>1とすると、C=(1,∞)であり、補題 1.2 (1)と補題 1.2 (2)によりt∈I(x)、t≥0についてφt(x)∈(1,x]である。補題 1.2 (3)をK=[1,x]に適用すると、[0,∞)⊆I(x)であり、φt(x)→1(t→∞)である。x=1ならば、§E13.15 命題 4.3 (2)によりすべてのt∈Rでφt(x)=1である。以上により、x>0ならば、すべてのt≥0で∣φt(x)−1∣≤∣x−1∣であり、φt(x)→1(t→∞)である。
ε>0に対してδ=min{ε,1}とおくと、∣x−1∣<δならばx>0であるから、[0,∞)⊆I(x)であり、すべてのt≥0で∣φt(x)−1∣<εである。また∣x−1∣<1ならばφt(x)→1である。したがって1は漸近安定である。一方、任意のδ>0について、0<x<min{δ,1/2}を満たすxはφt(x)→1を満たすから、t≥0が存在して∣φt(x)−0∣>1/2である。したがって0はε=1/2について Lyapunov 安定性の条件を満たさず、不安定である。
3 線形系の安定性と二次元線形系の分類
命題 3.1.n∈N≥1、A∈Rn×nとし、F(x)=Axで定まるF:Rn→Rnの流れをφとする。このとき、次が成り立つ。
- φは完備な流れであり、すべてのt∈R、x∈Rnについてφt(x)=etAxである。
- S∈Rn×nを可逆行列とし、G(y)=SAS−1yでG:Rn→Rnを定める。このとき、h(x)=Sxで定まるh:Rn→RnはFとGの流れの位相共役である。
証明.主張 (1)を示す。Rnのノルムと対応する作用素ノルムについて∥Ax−Ay∥≤∥A∥∥x−y∥であるから、Fは Lipschitz 連続である。x∈Rnとし、z(t)=etAx(t∈R)とおく。§E13.12 定義 1.2によりetAは実行列であり、§E13.12 定理 2.1によりzはR上でC1級で、z′(t)=AetAx=F(z(t))、z(0)=xである。§E13.15 補題 1.2をJ=R、τ=0に適用すると、R⊆I(x)であり、R上でz=φ(⋅,x)である。
主張 (2)を示す。hは連続な逆写像y↦S−1yをもつ同相写像である。Gの流れをψとする。主張 (1)をAとSAS−1に適用すると、FとGの流れはともに完備であり、ψt(y)=etSAS−1yである。§E13.12 命題 4.1を行列tAとSに適用するとetSAS−1=SetAS−1であるから、すべてのt∈R、x∈RnについてIG(h(x))=R=IF(x)であり、ψt(h(x))=SetAx=h(φt(x))である。▨
補題 3.2.d∈N≥1とし、A∈Rd×dの複素固有値の実部がすべて負であるとして、その最大値をαとする。Rdのノルムを一つ固定し、行列には対応する作用素ノルムを用いる。このとき、0<γ<−αを満たす任意のγについて、M≥1が存在して
∥etA∥≤Me−γt(t≥0)(3.2.1) が成り立つ。特に、式 (3.2.1)を満たすM≥1とγ>0が存在する。
証明. 作用素ノルムはRd×d上のノルムであるから、§E13.12 系 6.6をK=R、B=Aに適用して結論を得る。▨
命題 3.3.n∈N≥1、A∈Rn×nとし、Aの相異なる複素固有値全体の集合をσ(A)と書き、F(x)=AxでF:Rn→Rnを定める。Rnのノルムを一つ固定し、行列には対応する作用素ノルムを用いる。このとき、原点はFの平衡点であり、次が成り立つ。
- 原点が Lyapunov 安定であることは、supt≥0∥etA∥<∞であるための必要十分条件である。さらに、supt≥0∥etA∥<∞であることは、すべてのβ∈σ(A)がℜβ≤0を満たし、ℜβ=0を満たすすべてのβ∈σ(A)についてCnの部分空間としてGβ(A)=ker(A−βI)が成り立つための必要十分条件である。
- 原点が漸近安定であることは、すべてのβ∈σ(A)がℜβ<0を満たすための必要十分条件である。このとき、すべてのx∈RnについてetAx→0(t→∞)である。
証明.命題 3.1 (1)により、Fの流れφは完備であり、すべてのt∈R、x∈Rnについてφt(x)=etAxである。
主張 (1)を示す。C=supt≥0∥etA∥<∞とする。e0A=IであるからC≥1である。ε>0に対してδ=ε/Cとおくと、∥x∥<δならば、すべてのt≥0で∥φt(x)∥≤C∥x∥<εである。したがって原点は Lyapunov 安定である。逆に、原点が Lyapunov 安定であるとし、ε=1に対するδ>0を取る。x=0についてy=δx/(2∥x∥)とおくと、∥y∥<δであるから、すべてのt≥0で∥etAy∥<1であり、∥etAx∥<2∥x∥/δである。したがってsupt≥0∥etA∥≤2/δである。作用素ノルムはRn×n上のノルムであるから、§E13.12 補題 6.4 (2)をK=R、B=Aに適用すると、残りの同値を得る。
主張 (2)を示す。すべてのβ∈σ(A)がℜβ<0を満たすとする。補題 3.2により、M≥1とγ>0が存在して、すべてのt≥0で∥etA∥≤Me−γtである。したがってsupt≥0∥etA∥≤Mであり、主張 (1)により原点は Lyapunov 安定である。さらに、すべてのx∈Rnで∥etAx∥≤Me−γt∥x∥→0(t→∞)であるから、原点は漸近安定である。
逆に、原点が漸近安定であるとし、ρ>0を漸近安定性の定義の条件を満たすように取る。β∈σ(A)がℜβ≥0を満たすと仮定する。βに対するAの固有ベクトルv∈Cnを取り、v=p+iq(p,q∈Rn)と書く。§E13.12 補題 6.2 (3)によりetAv=etβvであり、etAは実行列であるから、etApとetAqはそれぞれetβvの実部と虚部である。vj=0となる成分vjを取ると、すべてのt≥0で∣etβvj∣=etℜβ∣vj∣≥∣vj∣>0である。etAp→0かつetAq→0(t→∞)ならば、§E10.1 命題 2.2により両者の各成分が0に収束し、etβvj→0となって前の不等式に反する。したがって、w=pまたはw=qのいずれかについて、etAwは0に収束しない。このときw=0であり、w′=ρw/(2∥w∥)とおくと、∥w′∥<ρであるが、etAw′=ρetAw/(2∥w∥)は0に収束しない。これはρの取り方に反する。したがって、すべてのβ∈σ(A)はℜβ<0を満たす。▨
定義 3.4.A∈R2×2とし、線形系x′=Axの平衡点である原点を考える。Aの複素固有値を重複度を込めてλ1,λ2とする。λ1,λ2が相異なる実数で同じ符号をもつとき、原点を結節点 (node) という。実数λ=0についてA=λIであるとき、原点を星形結節点 (star node) という。λ1=λ2が0でない実数であり、A=λ1Iであるとき、原点を非正則結節点 (improper node) という。λ1,λ2が実数でλ1λ2<0を満たすとき、原点を鞍点 (saddle point) という。λ1,λ2が実数でなく、その実部が0でないとき、原点を渦状点 (spiral point) という。λ1,λ2が実数でなく、その実部が0であるとき、原点を中心 (center of a planar linear system) という。
定理 3.5.A∈R2×2とし、F(x)=Axで定まるF:R2→R2の流れをφとし、∥⋅∥をR2のノルムとする。このとき、原点はFの平衡点であり、次が成り立つ。
- 原点が結節点であるとき、Aの二つの固有値がともに負ならば原点は漸近安定であり、ともに正ならば原点は不安定である。
- 原点が鞍点であるとき、原点は不安定である。さらに、Aの負の固有値をλsとすると、x∈R2について、Ax=λsxならばφt(x)→0(t→∞)であり、Ax=λsxならば∥φt(x)∥→∞(t→∞)である。
- 原点が星形結節点または非正則結節点であるとき、Aの重複固有値が負ならば原点は漸近安定であり、正ならば原点は不安定である。
- 原点が渦状点であるとき、Aの固有値の実部が負ならば原点は漸近安定であり、正ならば原点は不安定である。原点が中心であるとき、原点は Lyapunov 安定であるが漸近安定でない。
- Aが固有値0をもつとき、原点はFの孤立平衡点でない。Aの固有値が0と実数λ=0であるとき、λ<0ならば原点は Lyapunov 安定であるが漸近安定でなく、λ>0ならば原点は不安定である。Aの固有値が0だけであるとき、原点が Lyapunov 安定であることは、A=0であるための必要十分条件であり、このとき原点は漸近安定でない。
証明.命題 3.1 (1)により、すべてのt∈R、x∈R2についてφt(x)=etAxである。Aの相異なる複素固有値全体の集合をσ(A)とする。
主張 3.5.1. 次の三つが成り立つ。σ(A)のすべての元の実部が負ならば、原点は漸近安定である。σ(A)が実部0の元をもつならば、原点は漸近安定でない。σ(A)が実部が正の元をもつならば、原点は不安定である。
主張 3.5.2.Aが相異なる二つの複素固有値β1,β2をもつならば、j=1,2についてC2の部分空間としてGβj(A)=ker(A−βjI)である。
証明.§E13.12 定義 6.1によりC2=Gβ1(A)⊕Gβ2(A)であり、各Gβj(A)はβjに対する固有ベクトルを含むから零でない。したがって各Gβj(A)は一次元であり、零でない部分空間ker(A−βjI)⊆Gβj(A)と一致する。▨
主張 3.5.1により、主張 (1)、主張 (3)、主張 (4)の渦状点についての主張、および主張 (2)の不安定性が従う。原点が中心であるとき、ω=0が存在してσ(A)={iω,−iω}であるから、主張 3.5.2と命題 3.3 (1)により原点は Lyapunov 安定であり、主張 3.5.1により漸近安定でない。これで主張 (4)は示された。
主張 (2)の後半を示す。Aの正の固有値をλuとする。A−λsIとA−λuIは実行列で可逆でないから、実の固有ベクトルvs,vu∈R2∖{0}が存在してAvs=λsvs、Avu=λuvuである。vu=cvsを満たすc∈Rが存在するならばλuvu=Avu=cλsvs=λsvuとなってλu=λsに反するから、vs,vuはR2の基底である。x=avs+bvu(a,b∈R)と書くと、(A−λsI)x=b(λu−λs)vuであるから、Ax=λsxであることとb=0であることは同値である。§E13.12 補題 6.2 (3)により
φt(x)=aeλstvs+beλutvu(t∈R)である。b=0ならば∥φt(x)∥=∣a∣eλst∥vs∥→0(t→∞)である。N(avs+bvu)=max{∣a∣,∣b∣}はR2のノルムであるから、§E10.1 命題 2.2により定数K>0が存在して、すべてのy∈R2でN(y)≤K∥y∥である。b=0ならば∥φt(x)∥≥K−1N(φt(x))≥K−1∣b∣eλut→∞(t→∞)である。
主張 (5)を示す。Aは固有値0をもつ実行列であるから、v∈R2∖{0}が存在してAv=0である。原点の任意の開近傍Wについて、十分小さいs>0に対してsv∈Wであり、sv=0はFの平衡点である。したがって原点は孤立平衡点でない。Aの固有値が0とλ=0であるとき、λはAのトレースに等しいから実数である。主張 3.5.2によりG0(A)=kerAであるから、λ<0ならば命題 3.3 (1)により原点は Lyapunov 安定であり、0∈σ(A)であるから主張 3.5.1により漸近安定でない。λ>0ならば主張 3.5.1により原点は不安定である。Aの固有値が0だけであるとき、§E13.12 定義 6.1によりG0(A)=C2であるから、命題 3.3 (1)により、原点が Lyapunov 安定であることはkerA=C2、すなわちA=0と同値である。このとき0∈σ(A)であるから、主張 3.5.1により原点は漸近安定でない。▨
例 3.6 (鞍点).
A=(1221)の特性多項式は(z−1)2−4=(z−3)(z+1)であるから、固有値は3,−1であり、原点は鞍点である。A(1,1)T=3(1,1)T、A(1,−1)T=−(1,−1)Tであるから、x=u(1,1)T+v(1,−1)Tと書くと、x′=Axの流れは
φt(x)=ue3t(1,1)T+ve−t(1,−1)Tである。定理 3.5 (2)により、φt(x)→0(t→∞)となるのはxが直線x1+x2=0上にあるときに限り、それ以外のxについては∥φt(x)∥→∞である。
例 3.7 (安定な渦状点).
A=(−11−1−1),J=(01−10)とおくと、A=−I+Jであり、Aの固有値は−1±iであるから、原点は渦状点である。−IとJは可換であるから、§E13.12 命題 3.2と§E13.12 例 4.4により
etA=e−tetJ=e−t(costsint−sintcost)である。したがって、Euclid ノルム∥⋅∥2について∥etAx∥2=e−t∥x∥2であり、原点は漸近安定である。x=0の軌道上の点etAxは、原点のまわりの偏角がxの偏角からtだけ増え、原点からの距離がe−t倍になった点である。
4 非線形系の線形化
注意 4.1.n∈N≥1とし、U⊆Rnを開集合、F:U→Rnを写像、x∗∈UをF(x∗)=0を満たす点とし、Fはx∗で微分可能であるとする。A=DF(x∗)とおき、x∗+u∈Uを満たすu∈Rnについてr(u)=F(x∗+u)−Auと定める。F(x∗)=0であるから、微分可能性の定義により∥r(u)∥/∥u∥→0(u→0)である。xがx′=F(x)の古典解ならば、u=x−x∗は
u′=F(x∗+u)=Au+r(u),A=DF(x∗),∥u∥∥r(u)∥⟶0(u→0)(4.1.1) を満たす。
定義 4.2.n∈N≥1とし、U⊆Rnを開集合、F:U→Rnを写像、x∗∈UをF(x∗)=0を満たす点とし、Fはx∗で微分可能であるとする。Jacobi 行列DF(x∗)が実部0の複素固有値をもたないとき、x∗をFの双曲型平衡点 (hyperbolic equilibrium point) という。
補題 4.3.n∈N≥1とし、A∈Rn×nは実部が0の複素固有値をもたないとする。∥⋅∥をRnのノルムとし、行列には対応する作用素ノルムを用いる。Ps,Pu,Es,Euを、Aに対する§E13.12 補題 6.7(K=R、d=n)の射影と部分空間とする。K≥1、α>0は、すべてのt≥0について∥etAPs∥≤Ke−αtを満たし、すべてのt≤0について∥etAPu∥≤Keαtを満たすとし、0<β<αとする。
∥x∥s=t≥0supeβt∥etAx∥(x∈Es),∥x∥u=t≤0supe−βt∥etAx∥(x∈Eu)とおき、∥x∥∗=max{∥Psx∥s,∥Pux∥u}(x∈Rn)、p=max{∥Ps∥,∥Pu∥}とおく。このとき、次が成り立つ。
- ∥⋅∥s、∥⋅∥uはそれぞれEs、Euのノルムであり、y∈Esについて∥y∥≤∥y∥s≤K∥y∥、y∈Euについて∥y∥≤∥y∥u≤K∥y∥である。さらに、∥⋅∥∗はRnのノルムであり、すべてのx∈Rnについて21∥x∥≤∥x∥∗≤Kp∥x∥である。
- すべてのt≥0についてetAEs=Es、etAEu=Euであり、y∈Esならば∥etAy∥s≤e−βt∥y∥s、y∈Euならば∥e−tAy∥u≤e−βt∥y∥uである。
証明.主張 (1)を示す。§E13.12 補題 6.7 (4)により、y∈EsならばPsy=yであり、y∈EuならばPuy=yである。y∈Es、t≥0について
eβt∥etAy∥=eβt∥etAPsy∥≤Ke−(α−β)t∥y∥≤K∥y∥であり、t=0での値は∥y∥であるから、∥y∥≤∥y∥s≤K∥y∥である。y∈Eu、t≤0について
e−βt∥etAy∥=e−βt∥etAPuy∥≤Ke(α−β)t∥y∥≤K∥y∥であるから、∥y∥≤∥y∥u≤K∥y∥である。各tについてy↦eβt∥etAy∥とy↦e−βt∥etAy∥は半ノルムであるから、それらの有限な上限∥⋅∥s、∥⋅∥uも半ノルムであり、∥y∥以上であることからノルムである。Ps、Puは線形であるから∥⋅∥∗はRnの半ノルムである。§E13.12 補題 6.7 (2)によりPs+Pu=Iであるから、x∈Rnについて
∥x∥≤∥Psx∥+∥Pux∥≤∥Psx∥s+∥Pux∥u≤2∥x∥∗であり、∥x∥∗≤Kmax{∥Psx∥,∥Pux∥}≤Kp∥x∥である。前の不等式により∥⋅∥∗はノルムである。
主張 (2)を示す。t≥0とする。§E13.12 補題 6.7 (4)によりetAEs=Es、etAEu=Eu、e−tAEu=Euであり、§E13.12 定理 2.1によりs∈RについてesAetA=e(s+t)Aである。y∈Esについて
∥etAy∥s=s≥0supeβs∥e(s+t)Ay∥=e−βts≥tsupeβs∥esAy∥≤e−βt∥y∥sであり、y∈Euについて
∥e−tAy∥u=s≤0supe−βs∥e(s−t)Ay∥=e−βts≤−tsupe−βs∥esAy∥≤e−βt∥y∥uである。▨
補題 4.4.Eを有限次元実ベクトル空間とし、E=Eu⊕EsをEu={0}を満たす直和分解とする。EuとEsにノルム∥⋅∥を与え、zu∈Eu、zs∈Esについて∥zu+zs∥=max{∥zu∥,∥zs∥}によりEのノルムを定める。線形写像B:E→EはB(Eu)⊆Eu、B(Es)⊆Esを満たし、Bu=B∣Eu:Eu→Euは可逆であるとし、Bs=B∣Es:Es→Esとおく。実数a,bは0≤a<b、b>1を満たし、
∥Bsv∥≤a∥v∥(v∈Es),∥Bu−1v∥≤b−1∥v∥(v∈Eu)(4.4.1) が成り立つとする。N⊆Eを0の開近傍、G:N→EをG(z)=Bz+g(z)の形の写像とし、g(0)=0かつ∥g(z)∥/∥z∥→0(z→0)であるとする。このとき、r>0が存在して、{z∈E:∥z∥<r}⊆Nであり、任意のδ>0に対して、∥z0∥<δを満たすz0∈Eと整数k≥1で、zj=G(zj−1)(1≤j≤k)により定まるz1,…,zkが、0≤j≤k−1で∥zj∥<rを満たし、∥zk∥≥rを満たすものが存在する。
証明.η>0をa+η<b−ηかつμ=b−η>1を満たすように取る。g(0)=0と∥g(z)∥/∥z∥→0により、r>0が存在して、∥z∥<rならばz∈Nかつ∥g(z)∥≤η∥z∥である。z∈Eをz=zu+zs(zu∈Eu、zs∈Es)と分解し、錐C={z∈E:∥zs∥≤∥zu∥}を考える。ノルムの定義により∥g(z)u∥≤∥g(z)∥、∥g(z)s∥≤∥g(z)∥である。
主張 4.4.1.z∈C、∥z∥<rならば、G(z)∈Cであり、∥G(z)∥≥μ∥z∥である。
証明.∥z∥=∥zu∥であり、BはEuとEsを不変にするからG(z)u=Buzu+g(z)u、G(z)s=Bszs+g(z)sである。式 (4.4.1)から∥zu∥=∥Bu−1Buzu∥≤b−1∥Buzu∥であるから、
∥G(z)u∥∥G(z)s∥≥∥Buzu∥−∥g(z)u∥≥(b−η)∥zu∥,≤∥Bszs∥+∥g(z)s∥≤a∥zs∥+η∥zu∥≤(a+η)∥zu∥である。a+η<b−ηであるからG(z)∈Cであり、∥G(z)∥=∥G(z)u∥≥μ∥z∥である。▨
δ>0とし、z0∈Euを0<∥z0∥<min{δ,r}を満たすように取る。z0∈Cである。j≥1について、z0,…,zj−1が定まってCに属し、ノルムがr未満であり、∥zj−1∥≥μj−1∥z0∥であるならば、主張 4.4.1によりzj=G(zj−1)が定まってCに属し、∥zj∥≥μj∥z0∥である。μ>1であるからμj∥z0∥≥rを満たすjが存在し、したがって∥zk∥≥rとなる最小のk≥1が存在する。このkについて、0≤j≤k−1で∥zj∥<rである。▨
定理 4.5 (Lyapunov の間接法).n∈N≥1とし、U⊆Rnを開集合、F:U→Rnを局所 Lipschitz 連続な写像、φをFの流れ、x∗∈UをFの平衡点とする。x∗の開近傍W⊆Uが存在してF∣WはC1級であるとし、A=DF(x∗)とおく。このとき、次が成り立つ。
- Aの複素固有値の実部がすべて負ならば、x∗は漸近安定である。さらに、その実部の最大値をαとすると、Rnの任意のノルム∥⋅∥と0<γ<−αを満たす任意のγについて、δ>0とM≥1が存在して、∥x−x∗∥<δを満たす任意のx∈Uについて[0,∞)⊆I(x)であり、すべてのt≥0で∥φt(x)−x∗∥≤Me−γt∥x−x∗∥である。
- Aが実部が正の複素固有値を少なくとも一つもつならば、x∗は不安定である。
- x∗が双曲型平衡点ならば、x∗が漸近安定であることは、Aの複素固有値の実部がすべて負であるための必要十分条件であり、そうでないときx∗は不安定である。
証明.主張 (1)を示す。∥⋅∥をRnのノルム、0<γ<−αとし、行列には対応する作用素ノルムを用いる。γ<γ′<−αを満たすγ′を取る。補題 3.2により、M≥1が存在して、すべてのt≥0で∥etA∥≤Me−γ′tである。η>0をMη≤γ′−γを満たすように取る。注意 4.1のrについて、ρ0>0が存在して、閉球B={y∈Rn:∥y−x∗∥≤ρ0}はWに含まれ、∥u∥≤ρ0ならば∥r(u)∥≤η∥u∥である。δ=ρ0/Mとおく。
x∈U、∥x−x∗∥<δとし、u(t)=φt(x)−x∗(t∈I(x))、β+=supI(x)とおく。[0,t]上で∥u∥<ρ0が成り立つt∈[0,β+)全体の上限をτとする。∥u(0)∥<δ≤ρ0とuの連続性によりτ>0である。[0,τ)上でuはC1級であり、式 (4.1.1)によりu′=Au+b、b(t)=r(u(t))を満たし、bは[0,τ)上で連続である。§E13.12 系 2.5を区間[0,τ)とt0=0に適用すると、t∈[0,τ)について
u(t)=etAu(0)+∫0te(t−s)Ar(u(s))dsである。t∈[0,τ)で∥u(t)∥<ρ0であるから、w(t)=eγ′t∥u(t)∥は
w(t)≤M∥u(0)∥+∫0tMηw(s)ds(0≤t<τ)を満たす。0<t1<τを取り、§E13.7 系 1.2を区間[0,t1]、定数M∥u(0)∥とMη、零関数に適用すると、0≤t≤t1でw(t)≤M∥u(0)∥eMηtである。t1は任意であるから、t∈[0,τ)について
∥u(t)∥≤M∥u(0)∥e−(γ′−Mη)t≤Me−γt∥u(0)∥(4.5.1) である。τ<β+と仮定すると、uの連続性と式 (4.5.1)により∥u(τ)∥≤M∥u(0)∥<Mδ=ρ0であり、τを含む開区間上で∥u∥<ρ0であるから、τより大きいt∈[0,β+)で[0,t]上∥u∥<ρ0が成り立つ。これはτが上限であることに反するから、τ=β+である。したがって、すべてのt∈I(x)、t≥0についてφt(x)∈Bである。BはUに含まれるコンパクト集合であるから、§E13.15 補題 3.1 (1)により[0,∞)⊆I(x)であり、式 (4.5.1)はすべてのt≥0で成り立つ。ε>0について、∥x−x∗∥<min{δ,ε/M}ならば、すべてのt≥0で∥φt(x)−x∗∥<εである。また∥x−x∗∥<δならばφt(x)→x∗(t→∞)である。したがってx∗は漸近安定である。
主張 (2)を示す。Aの複素固有値で実部が正のもの全体にわたる実部の最小値をλ+とし、c=λ+/2、A′=A−cIとおく。μ∈CがAの固有値であることとμ−cがA′の固有値であることは同値であるから、A′の固有値の実部は、実部が正のAの固有値に対してはc以上であり、それ以外のAの固有値に対しては−c以下である。特にA′は実部0の固有値をもたず、実部が正の固有値をもつ。Rnのノルム∥⋅∥0を一つ固定し、行列には対応する作用素ノルムを用いる。§E13.12 補題 6.7をK=RとA′に適用してPs,Pu,Es,Euを取り、§E13.12 補題 6.7 (5)によりK≥1とκ>0を取って、t≥0で∥etA′Ps∥0≤Ke−κt、t≤0で∥etA′Pu∥0≤Keκtとする。§E13.12 補題 6.7 (3)により、Rn=Es⊕Euであり、Eu={0}である。補題 4.3を、A′、∥⋅∥0、Kと、同補題のα,βをそれぞれκ,κ/2として適用し、Esのノルム∥⋅∥s、Euのノルム∥⋅∥uとRnのノルム∥⋅∥∗を取る。
cIとA′は可換であるから、§E13.12 命題 3.2により、すべてのt∈RでetA=ectetA′である。B=eAとおく。補題 4.3 (2)をt=1に適用すると、BはEsとEuを不変にし、v∈Esならば
∥Bv∥s=ec∥eA′v∥s≤ec−κ/2∥v∥sであり、v∈Euならば
∥e−Av∥u=e−c∥e−A′v∥u≤e−(c+κ/2)∥v∥uである。§E13.12 定理 2.1によりB∣Euは逆写像e−A∣Euをもつ。a=ec−κ/2、b=ec+κ/2とおくと、c>0かつκ>0であるから0≤a<bかつb>1である。§E13.12 補題 6.7 (4)により、zs∈Es、zu∈EuについてPs(zu+zs)=zs、Pu(zu+zs)=zuであるから、∥zu+zs∥∗=max{∥zu∥u,∥zs∥s}である。したがって、E=Rnに∥⋅∥∗を、EsとEuにそれぞれ∥⋅∥sと∥⋅∥uを与えると、Bとa,bは式 (4.4.1)を満たす。
F∣Wの流れをψ、その最大存在区間をIW(y)とし、W1={y∈W:1∈IW(y)}とおく。y∈Wについて、すべてのs∈IW(y)でψs(y)∈Wであるから、§E13.15 補題 7.3 (2)をF∣W、N=W、Fに適用すると、IW(y)⊆I(y)であり、IW(y)上でψs(y)=φs(y)である。§E13.15 定理 2.1 (2)によりW1は開集合であり、§E13.15 系 2.2 (4)によりIW(x∗)=Rであるからx∗∈W1である。N0={z∈Rn:x∗+z∈W1}とおき、G:N0→RnをG(z)=ψ1(x∗+z)−x∗で定める。§E13.15 系 2.2 (3)によりGはC1級であり、§E13.15 系 2.2 (4)によりG(0)=0、DG(0)=eAである。g(z)=G(z)−Bzとおくと、g(0)=0であり、微分可能性の定義により∥g(z)∥∗/∥z∥∗→0(z→0)である。
主張 4.5.1.k≥1とし、z0,…,zk−1∈N0がzj=G(zj−1)(1≤j≤k)を満たすとする。このとき、y=x∗+z0についてk∈I(y)かつφk(y)=x∗+zkである。
証明.y=x∗+z0とおく。0≤j≤k−1についてj∈I(y)かつφj(y)=x∗+zjであるとすると、x∗+zj∈W1であるから1∈IW(x∗+zj)⊆I(x∗+zj)であり、φ1(x∗+zj)=ψ1(x∗+zj)=x∗+zj+1である。§E13.15 定理 1.3によりI(x∗+zj)=I(y)−jであるからj+1∈I(y)であり、φj+1(y)=φ1(φj(y))=x∗+zj+1である。j=0では仮定は成り立つから、帰納法によりk∈I(y)かつφk(y)=x∗+zkである。▨
x∗が Lyapunov 安定であると仮定する。補題 2.3により、x∗はノルム∥⋅∥∗について Lyapunov 安定である。補題 4.4をE=Rn、B、N0、Gに適用してr>0を取り、Lyapunov 安定性の定義によりε=rに対するδ>0を取る。補題 4.4により、∥z0∥∗<δを満たすz0とk≥1が存在して、0≤j≤k−1で∥zj∥∗<rであり、∥zk∥∗≥rである。{z:∥z∥∗<r}⊆N0であるからz0,…,zk−1∈N0であり、主張 4.5.1によりy=x∗+z0∈Uについてk∈I(y)かつ∥φk(y)−x∗∥∗=∥zk∥∗≥rである。一方、∥y−x∗∥∗<δであるから、Lyapunov 安定性により∥φk(y)−x∗∥∗<rである。この二つは両立しないから、x∗は Lyapunov 安定でない。
主張 (3)を示す。x∗が双曲型平衡点であるとする。Aの複素固有値の実部がすべて負ならば、主張 (1)によりx∗は漸近安定である。そうでないならば、Aは実部0の固有値をもたないから実部が正の固有値をもち、主張 (2)によりx∗は不安定であって、特に漸近安定でない。▨
5 演習
問題 5.1 (零固有値と Jordan 構造). 二つの行列
A0=(0000),N=(0010)について、線形系x′=A0xとx′=Nxの解を求めよ。両行列の固有値はともに0,0であるが、前者の原点は Lyapunov 安定であって漸近安定ではなく、後者の原点は不安定であることを証明せよ。
解答.
A0とNの特性多項式はともにz2である。初期値を(u0,v0)Tとすると、命題 3.1 (1)により、A0に対する解はすべての実数tに対して
u(t)=u0,v(t)=v0である。Euclid ノルムを用いる。任意のε>0に対してδ=εと取れば、原点からの距離はε未満に保たれるため、原点は Lyapunov 安定である。非零の初期値から出る解は原点へ収束しないため、原点は漸近安定ではない。
N2=0であるから、§E13.12 定理 5.1によりetN=I+tNであり、解は
u(t)=u0+tv0,v(t)=v0(t∈R)となる。ε=1と固定する。任意のδ>0に対して0<a<min{δ,1}を取り、(u0,v0)=(0,a)と置くと、初期値のノルムはδ未満であるが、t=2/aではu(t)=2となり、解のノルムは1より大きい。したがって、Nに対する原点は不安定である。▨
問題 5.2 (重複固有値をもつ安定な線形系の軌道).λ<0とし、
J=(λ01λ)と置く。Jに対する線形系u′=λu+v、v′=λvを考える。初期値(u0,v0)に対する解を求め、v0=0ならば十分大きいtでu(t)=0となり、v(t)/u(t)→0が成り立つことを示せ。また、行列λIに対する系と比較し、両系の原点が漸近安定であることと、λIに対する非定数解の軌道が原点を通る直線上にあることを示せ。
解答.
§E13.12 定理 5.1から
u(t)=eλt(u0+tv0),v(t)=eλtv0(t∈R)を得る。v0=0のとき、t>∣u0/v0∣ならばt+u0/v0>0であるからu(t)=0であり、
u(t)v(t)=t+u0/v01⟶0(t→∞)となる。v0=0の場合はv(t)=0であり、解は固有直線v=0上にある。
行列λIに対する解は
(u(t),v(t))=eλt(u0,v0)(t∈R)である。非零の初期値に対して、軌道は初期値と原点を結ぶ直線上にあり、方向は一定である。一方、Jに対するv0=0の解では、比v(t)/u(t)が変化する。この解は、原点へ近づくにつれて固有直線v=0の方向へ近づく。λIに対する原点は星形結節点であり、Jに対する原点は非正則結節点である。両行列の重複固有値λは負であるから、定理 3.5 (3)により両系の原点は漸近安定である。▨
問題 5.3 (Jacobi 行列による局所安定性の判定). 非線形系
x′y′=−x+xy,=−2y+x2について、原点が平衡点であることを確認し、Jacobi 行列の固有値から原点の局所安定性を判定せよ。また、原点以外の平衡点をすべて求め、原点がすべての初期値から出る解を吸引するわけではないことを示せ。
解答.
右辺をF(x,y)と書くとF(0,0)=0であるから、原点は平衡点である。FはR2上でC1級であり、各点を中心とする閉球の上でDFは有界であるから、§E4.3 定理 2.1によりFは局所 Lipschitz 連続である。Jacobi 行列は
DF(x,y)=(−1+y2xx−2)であり、
DF(0,0)=(−100−2)の固有値は−1,−2である。二つの固有値の実部がともに負であるから、定理 4.5 (1)により原点は漸近安定である。
平衡点はx(−1+y)=0と−2y+x2=0の共通解である。x=0ならば第二式からy=0である。x=0ならば第一式からy=1であり、第二式からx2=2、すなわちx=±2である。したがって原点以外の平衡点は(2,1)と(−2,1)だけである。これらの点から出る解は定数解であり、原点へ収束しない。ゆえに原点はすべての初期値から出る解を吸引するわけではない。▨
問題 5.4.問題 5.3の系の平衡点(2,1)と(−2,1)が不安定であることを証明せよ。また、系
x′=x+y2,y′=−y3について、原点がただ一つの平衡点であり、双曲型平衡点ではないが不安定であることを証明せよ。
解答.
問題 5.3の解答のとおり、その系の右辺Fは局所 Lipschitz 連続であり、R2上でC1級であって、
DF(±2,1)=(0±22±2−2)(複号同順)である。その特性多項式はz(z+2)−4=z2+2z−4であるから、固有値は−1±5であり、−1+5>0である。定理 4.5 (2)により、(2,1)と(−2,1)はいずれも不安定である。
後半の系の右辺をH(x,y)=(x+y2,−y3)と書く。HはR2上の多項式写像であるからC1級であり、各点を中心とする閉球の上でDHは有界であるから、§E4.3 定理 2.1により局所 Lipschitz 連続である。H(x,y)=0ならばy3=0からy=0であり、x+y2=0からx=0である。したがって原点がただ一つの平衡点である。
DH(x,y)=(102y−3y2),DH(0,0)=(1000)の固有値は1と0である。固有値0の実部は0であるから、原点は双曲型平衡点ではなく、定理 4.5 (3)の仮定は満たされない。一方、DH(0,0)は実部が正の固有値1をもつから、定理 4.5 (2)により原点は不安定である。▨
問題 5.5 (非双曲型平衡点で線形化だけでは安定性が定まらないことの証明).r2=x2+y2と書き、次の三つの系を考える。
-
x′=−y,y′=x,
-
x′=−y−xr2,y′=x−yr2,
-
x′=−y+xr2,y′=x+yr2.
三つの系が原点で同じ Jacobi 行列をもち、その固有値の実部が零であることを示せ。そのうえで、(I)の系の原点は Lyapunov 安定であるが漸近安定ではなく、(II)の系の原点は漸近安定であり、(III)の系の原点は不安定であることを証明せよ。
解答.
三つの系の右辺はR2上の多項式写像であり、局所 Lipschitz 連続である。いずれの系でも原点は平衡点であり、xr2、yr2の偏導関数は原点で0であるから、原点における Jacobi 行列は
A=(01−10)である。固有値は±iであり、実部は零である。Euclid ノルムを用いる。
各系の流れをφとし、z∈R2∖{0}についてs(t)=∥φt(z)∥22(t∈I(z))とおく。φt(z)=0となるtが存在するならば、§E13.15 命題 4.3 (2)と§E13.15 命題 4.3 (1)によりγ(z)={0}となってz=0に反するから、I(z)上でs>0である。φt(z)=(x,y)と書くと、s′=2xx′+2yy′である。
(I)の系ではs′=2x(−y)+2yx=0であるから、§D1.14 定理 3.1によりsはI(z)上で定数であり、前方軌道は半径∥z∥2の円周に含まれる。§E13.15 補題 3.1 (1)により[0,∞)⊆I(z)であり、すべてのt≥0で∥φt(z)∥2=∥z∥2である。任意のε>0に対してδ=εと取ると Lyapunov 安定性の条件が成り立つ。z=0から出る解は原点へ収束しないから、原点は漸近安定でない。
(II)の系では
s′=2x(−y−xs)+2y(x−ys)=−2s2であるから、(1/s)′=−s′/s2=2であり、t∈I(z)について1/s(t)=1/s(0)+2tである。したがってt∈I(z)、t≥0についてs(t)≤s(0)であり、前方軌道は半径∥z∥2の閉円板に含まれる。§E13.15 補題 3.1 (1)により[0,∞)⊆I(z)であり、t≥0について
s(t)=1+2s(0)ts(0)≤s(0),s(t)→0(t→∞)である。任意のε>0に対してδ=εと取ると Lyapunov 安定性の条件が成り立ち、すべてのzについてφt(z)→0であるから、原点は漸近安定である。
(III)の系では
s′=2x(−y+xs)+2y(x+ys)=2s2であるから、(1/s)′=−2であり、t∈I(z)について1/s(t)=1/s(0)−2tである。ε=1とし、δ>0を任意に取る。0<∥z∥2<min{δ,1}を満たすzを取り、t1=(1/s(0)−1)/2とおくと、s(0)<1であるからt1>0である。t1∈/I(z)ならば[0,∞)⊆I(z)である。t1∈I(z)ならば1/s(t1)=1であり、∥φt1(z)∥2=1である。いずれの場合も、zについて Lyapunov 安定性の条件はε=1に対して成り立たない。δは任意であったから、原点は不安定である。▨
問題 5.6.
A=(−11−4−1),S=(1002)とし、F(x)=Axで定まるF:R2→R2の流れをφとする。SAS−1を求め、命題 3.1 (2)と系 2.5を用いて、原点がFの漸近安定な平衡点であることを証明せよ。また、Euclid ノルムを∥⋅∥2とし、p=(1,−1)について、関数t↦∥φt(p)∥2が0を含むある開区間で狭義単調増加であることを示せ。
解答.
S−1=(1001/2)であるから
SAS−1=(1002)(−11−2−1/2)=(−12−2−1)である。B=SAS−1、J=(01−10)とおくとB=−I+2Jであり、−Iと2Jは可換であるから、§E13.12 命題 3.2と§E13.12 例 4.4(ω=2)によりetB=e−te2tJはe−tと角2tの回転行列の積である。命題 3.1 (1)によりG(y)=Byの流れはψt(y)=etByであり、すべてのt∈Rで∥ψt(y)∥2=e−t∥y∥2である。任意のε>0に対してδ=εと取ると Lyapunov 安定性の条件が成り立ち、すべてのyについてψt(y)→0であるから、原点はGの漸近安定な平衡点である。命題 3.1 (2)によりh(x)=SxはFとGの流れの位相共役であり、h(0)=0であるから、系 2.5により原点はFの漸近安定な平衡点である。
q(t)=∥φt(p)∥22とおく。φt(p)=(x1,x2)と書くと、qはR上でC1級であり、
q′(t)=2⟨φt(p),Aφt(p)⟩=2(−x12−3x1x2−x22)である。t=0では(x1,x2)=(1,−1)であるからq′(0)=2(−1+3−1)=2>0である。q′は連続であるから、0を含む開区間J0が存在してJ0上でq′>0であり、§D1.14 定理 3.1によりqはJ0上で狭義単調増加である。∥φt(p)∥2=q(t)であるから、t↦∥φt(p)∥2もJ0上で狭義単調増加である。一方、Sφt(p)=h(φt(p))=ψt(Sp)であるから、∥Sφt(p)∥2=e−t∥Sp∥2は狭義単調減少である。▨