§E13.22分岐

最終更新

パラメータμ\muを含む自律系x′=F(x,μ)x'=F(x,\mu)について、μ\muを固定するとき、F(x0,μ)=0F(x_0,\mu)=0を満たす点x0x_0を平衡点といい、定数関数x(t)=x0x(t)=x_0はx′=F(x,μ)x'=F(x,\mu)の解である。FFがC1C^1級で、μ=μ0\mu=\mu_0の平衡点x0x_0におけるDxF(x0,μ0)D_xF(x_0,\mu_0)が可逆であれば、μ0\mu_0の近くの各μ\muについて、x0x_0の近くの平衡点はμ\muのC1C^1級関数としてただ一つに定まる。DxF(x0,μ0)D_xF(x_0,\mu_0)が双曲型であれば、その平衡点におけるDxFD_xFの実部が負の固有値の個数と実部が正の固有値の個数も保たれる。DxFD_xFが固有値00をもつ平衡点については、この持続は成り立つとは限らない。たとえばx′=μ−x2x'=\mu-x^2は、μ>0\mu>0では平衡点を二個、μ<0\mu<0では零個もち、μ=0\mu=0で平衡点の個数が変わる。点の近くで平衡点の集合がパラメータの連続関数のグラフとして表されないとき、その点を分岐点という。本記事では、平衡点の持続を証明し、DxFD_xFの固有値00が代数的に単純な場合に、平衡点の方程式を一変数の方程式へ帰着して、サドルノード分岐、トランスクリティカル分岐、ピッチフォーク分岐と安定性の交代を証明する。

1 分岐点

定義 1.1.n∈N≥1n\in\NNとし、U⊆RnU\subseteq\R^nを開集合、Λ⊆R\Lambda\subseteq\Rを開区間、F:U×Λ→RnF:U\times\Lambda\to\R^nを連続写像とし、(x0,μ0)∈U×Λ(x_0,\mu_0)\in U\times\LambdaがF(x0,μ0)=0F(x_0,\mu_0)=0を満たすとする。μ0\mu_0を含む開区間Λ0⊆Λ\Lambda_0\subseteq\Lambda、x0x_0の開近傍B⊆UB\subseteq Uと連続写像ξ:Λ0→B\xi:\Lambda_0\to Bの組であって、(y,μ)∈B×Λ0(y,\mu)\in B\times\Lambda_0についてF(y,μ)=0F(y,\mu)=0であることとy=ξ(μ)y=\xi(\mu)であることが同値になるものが存在しないとき、(x0,μ0)(x_0,\mu_0)を族x′=F(x,μ)x'=F(x,\mu)の分岐点 (bifurcation point) という。(x0,μ0)(x_0,\mu_0)が分岐点であるとき、族x′=F(x,μ)x'=F(x,\mu)の平衡点は(x0,μ0)(x_0,\mu_0)で分岐 (bifurcation) するという。

注意 1.2.定義 1.1の分岐点は平衡点の集合だけで定まる変種であり、パラメータの変化による相図の位相的な型の変化で分岐を定める定式化とは一致しない。たとえばn=2n=2、U=R2U=\R^2、Λ=R\Lambda=\R、F(x,μ)=A(μ)xF(x,\mu)=A(\mu)x、A(μ)=(μ−11μ)A(\mu)=\begin{pmatrix}\mu&-1\\1&\mu\end{pmatrix}とする。このとき、det⁡A(μ)=μ2+1>0\det A(\mu)=\mu^2+1>0であるから、各μ\muについて平衡点は原点だけであり、Λ0=R\Lambda_0=\R、B=R2B=\R^2、ξ≡0\xi\equiv0の組により(0,0)(0,0)は定義 1.1の意味の分岐点でない。一方、A(μ)A(\mu)の固有値はμ±i\mu\pm iであるから、原点はμ<0\mu<0で漸近安定、μ>0\mu>0で不安定である(F(⋅,μ)F(\cdot,\mu)はC1C^1級であり、§E13.4 補題 4.2により局所 Lipschitz 連続であるから、§E13.16 定理 4.5 (1)と§E13.16 定理 4.5 (2)をW=R2W=\R^2として適用する)。

2 CkC^k級の陰関数定理

補題 2.1.n,m,k∈N≥1n,m,k\in\NNとし、U⊆Rn×RmU\subseteq\R^n\times\R^mを開集合、F:U→RmF:U\to\R^mをCkC^k級の写像とする。(a,b)∈U(a,b)\in UがF(a,b)=0F(a,b)=0を満たし、DyF(a,b)D_yF(a,b)が可逆であるとする。このとき、aaの開近傍A⊆RnA\subseteq\R^n、bbの開近傍B⊆RmB\subseteq\R^mとCkC^k級の写像g:A→Bg:A\to Bが存在して、次が成り立つ。

  1. (x,y)∈A×B(x,y)\in A\times Bについて、F(x,y)=0F(x,y)=0であることとy=g(x)y=g(x)であることは同値である。特にg(a)=bg(a)=bである。
  2. x∈Ax\in AについてDyF(x,g(x))D_yF(x,g(x))は可逆であり、Dg(x)=−DyF(x,g(x))−1DxF(x,g(x))Dg(x)=-D_yF(x,g(x))^{-1}D_xF(x,g(x))である。

証明.§E4.8 定理 2.1により、AA、BBとC1C^1級の写像g:A→Bg:A\to Bで主張 (1)と主張 (2)を満たすものが存在する。k=1k=1のときはこれで主張を得る。

C0C^0級は連続を意味する。§E4.4 定理 1.7により、j≥0j\ge0について、開集合の上の写像がCj+1C^{j+1}級であることは、その一階偏導関数がすべて存在してCjC^j級であることと同値である。この同値と、∂i(φψ)=φ ∂iψ+ψ ∂iφ\partial_i(\varphi\psi)=\varphi\,\partial_i\psi+\psi\,\partial_i\varphi、∂i(1/φ)=−∂iφ/φ2\partial_i(1/\varphi)=-\partial_i\varphi/\varphi^2から、jjについての帰納法により、CjC^j級の実数値関数の和と積、および零にならないCjC^j級の実数値関数の逆数はCjC^j級である。

0≤j≤k−10\le j\le k-1について、UUの上の任意のCjC^j級の実数値関数Φ\Phiについてx↦Φ(x,g(x))x\mapsto\Phi(x,g(x))がAAの上でCjC^j級であることを、(Cj)(\mathrm C_j)と書く。j=0j=0のとき、(C0)(\mathrm C_0)はggの連続性と§E4.1 命題 4.1 (4)から従い、ggは上で得たとおりC1C^1級である。1≤j≤k−11\le j\le k-1とし、(Cj−1)(\mathrm C_{j-1})が成り立ち、ggがCjC^j級であると仮定する。Φ\PhiをUUの上のCjC^j級の実数値関数とする。連鎖律(§E4.3 定理 1.1)により

∂xi(Φ(x,g(x)))=(∂xiΦ)(x,g(x))+∑l=1m(∂ylΦ)(x,g(x)) ∂xigl(x)\partial_{x_i}\bigl(\Phi(x,g(x))\bigr)=(\partial_{x_i}\Phi)(x,g(x))+\sum_{l=1}^m(\partial_{y_l}\Phi)(x,g(x))\,\partial_{x_i}g_l(x)

である。∂xiΦ\partial_{x_i}\Phi、∂ylΦ\partial_{y_l}\PhiはCj−1C^{j-1}級であるから、(Cj−1)(\mathrm C_{j-1})により右辺の(∂xiΦ)(x,g(x))(\partial_{x_i}\Phi)(x,g(x))、(∂ylΦ)(x,g(x))(\partial_{y_l}\Phi)(x,g(x))はCj−1C^{j-1}級であり、ggはCjC^j級であるから∂xigl\partial_{x_i}g_lはCj−1C^{j-1}級である。したがってx↦Φ(x,g(x))x\mapsto\Phi(x,g(x))の一階偏導関数はCj−1C^{j-1}級であり、x↦Φ(x,g(x))x\mapsto\Phi(x,g(x))はCjC^j級である。これで(Cj)(\mathrm C_j)を得る。DxFD_xF、DyFD_yFの成分はUUの上のCk−1C^{k-1}級、したがってCjC^j級の関数であるから、(Cj)(\mathrm C_j)によりx↦DxF(x,g(x))x\mapsto D_xF(x,g(x))、x↦DyF(x,g(x))x\mapsto D_yF(x,g(x))の成分はCjC^j級である。m=1m=1のときDyF(x,g(x))−1=1/DyF(x,g(x))D_yF(x,g(x))^{-1}=1/D_yF(x,g(x))であり、m≥2m\ge2のとき§D3.4 定理 7.4によりDyF(x,g(x))−1D_yF(x,g(x))^{-1}は余因子行列を行列式で割ったものである。余因子と行列式は成分の多項式であり、行列式はAAの上で00でないから、DyF(x,g(x))−1D_yF(x,g(x))^{-1}の成分はCjC^j級である。主張 (2)によりggの一階偏導関数はCjC^j級であり、ggはCj+1C^{j+1}級である。jjについての帰納法により、0≤j≤k−10\le j\le k-1について(Cj)(\mathrm C_j)が成り立ち、ggはCj+1C^{j+1}級である。j=k−1j=k-1とすれば、ggはCkC^k級である。▨

3 平衡点の持続

定義 3.1.n∈N≥1n\in\NNとし、A∈Cn×nA\in\C^{n\times n}の特性多項式をχA(z)=det⁡(zIn−A)\chi_A(z)=\det(zI_n-A)とする。λ∈C\lambda\in\CがAAの固有値であることはχA(λ)=0\chi_A(\lambda)=0と同値である(§E3.20 命題 3.1)。固有値λ\lambdaのχA\chi_Aの根としての重複度を、λ\lambdaの代数的重複度 (algebraic multiplicity) という。代数的重複度が11である固有値を、代数的に単純 (algebraically simple) であるという。§D3.12 定理 1.3により、AAは代数的重複度を込めてちょうどnn個の固有値をもつ。実行列は複素行列とみなす。

補題 3.2.n,p∈N≥1n,p\in\NNとし、T⊆RpT\subseteq\R^p、t0∈Tt_0\in Tとし、A:T→Rn×nA:T\to\R^{n\times n}はt0t_0で連続であるとする。このとき、次が成り立つ。

  1. Γ⊆C\Gamma\subseteq\Cを辺が座標軸に平行な閉長方形とし、∂Γ\partial\Gammaの上にA(t0)A(t_0)の固有値は無いとする。このとき、δ>0\delta>0が存在して、t∈Tt\in T、∥t−t0∥<δ\lVert t-t_0\rVert<\deltaならば、∂Γ\partial\Gammaの上にA(t)A(t)の固有値は無く、Γ\Gammaの内部にあるA(t)A(t)の固有値の個数(代数的重複度を込める)はA(t0)A(t_0)のそれに等しい。
  2. r∈N≥1r\in\NNとし、Γ1,…,Γr⊆C\Gamma_1,\ldots,\Gamma_r\subseteq\Cを辺が座標軸に平行な互いに素な閉長方形とする。A(t0)A(t_0)の固有値はすべてΓ1,…,Γr\Gamma_1,\ldots,\Gamma_rのいずれかの内部にあるとする。このとき、δ>0\delta>0が存在して、t∈Tt\in T、∥t−t0∥<δ\lVert t-t_0\rVert<\deltaならば、A(t)A(t)の固有値はすべてΓ1,…,Γr\Gamma_1,\ldots,\Gamma_rのいずれかの内部にあり、各jjについてΓj\Gamma_jの内部にあるA(t)A(t)の固有値の個数(代数的重複度を込める)はA(t0)A(t_0)のそれに等しい。
  3. Γ⊆C\Gamma\subseteq\Cを辺が座標軸に平行で実軸について対称な閉長方形とし、t∈Tt\in Tについて、Γ\Gammaの内部にあるA(t)A(t)の固有値の個数(代数的重複度を込める)は11であるとする。このとき、その固有値は実数で代数的に単純である。特に、∂Γ\partial\Gammaの上にA(t0)A(t_0)の固有値は無く、Γ\Gammaの内部にあるA(t0)A(t_0)の固有値の個数が11であるならば、主張 (1)のδ\deltaについて、t∈Tt\in T、∥t−t0∥<δ\lVert t-t_0\rVert<\deltaのとき、Γ\Gammaの内部にあるA(t)A(t)の固有値はちょうど一個であり、それは実数で代数的に単純である。
  4. c∈Rc\in\Rとし、直線{z∈C:ℜz=c}\{z\in\C:\Re z=c\}の上にA(t0)A(t_0)の固有値は無いとする。このとき、δ>0\delta>0が存在して、t∈Tt\in T、∥t−t0∥<δ\lVert t-t_0\rVert<\deltaならば、直線{z∈C:ℜz=c}\{z\in\C:\Re z=c\}の上にA(t)A(t)の固有値は無く、ℜz<c\Re z<cにあるA(t)A(t)の固有値の個数とℜz>c\Re z>cにあるA(t)A(t)の固有値の個数(ともに代数的重複度を込める)は、それぞれA(t0)A(t_0)のそれに等しい。

証明.主張 (1)を示す。f=χA(t0)f=\chi_{A(t_0)}、gt=χA(t)−fg_t=\chi_{A(t)}-fとおく。B∈Rn×nB\in\R^{n\times n}についてχB(z)=zn+∑j=0n−1cj(B)zj\chi_B(z)=z^n+\sum_{j=0}^{n-1}c_j(B)z^jと書くと、cj(B)c_j(B)はBBの成分の多項式であるからcj(A(t))→cj(A(t0))c_j(A(t))\to c_j(A(t_0))(t→t0t\to t_0、t∈Tt\in T)である。∂Γ\partial\Gammaはコンパクトでffは∂Γ\partial\Gammaの上で零点をもたない連続関数であるから、m=min⁡z∈∂Γ∣f(z)∣>0m=\min_{z\in\partial\Gamma}\lvert f(z)\rvert>0である。K=max⁡z∈∂Γmax⁡{1,∣z∣}nK=\max_{z\in\partial\Gamma}\max\{1,\lvert z\rvert\}^nとおくと、z∈∂Γz\in\partial\Gammaについて∣gt(z)∣≤K∑j=0n−1∣cj(A(t))−cj(A(t0))∣\lvert g_t(z)\rvert\le K\sum_{j=0}^{n-1}\lvert c_j(A(t))-c_j(A(t_0))\rvertである。したがってδ>0\delta>0が存在して、∥t−t0∥<δ\lVert t-t_0\rVert<\deltaならば∂Γ\partial\Gammaの上で∣gt∣<m≤∣f∣\lvert g_t\rvert<m\le\lvert f\rvertである。このとき∂Γ\partial\Gammaの上で∣f+gt∣≥∣f∣−∣gt∣>0\lvert f+g_t\rvert\ge\lvert f\rvert-\lvert g_t\rvert>0であるから、∂Γ\partial\Gammaの上にχA(t)=f+gt\chi_{A(t)}=f+g_tの零点は無い。§E5.14 系 3.3(D=CD=\C)により、ffとf+gtf+g_tはΓ\Gammaの内部に位数込みで同数の零点をもつ。零点の位数は固有値の代数的重複度である。

主張 (2)を示す。Γj\Gamma_jは互いに素であり、A(t0)A(t_0)の固有値はすべていずれかのΓj\Gamma_jの内部にあるから、∂Γj\partial\Gamma_jの上にA(t0)A(t_0)の固有値は無く、Γj\Gamma_jの内部にあるA(t0)A(t_0)の固有値の個数のjjについての和はnnである。各Γj\Gamma_jに主張 (1)を適用してδj>0\delta_j>0を取り、δ=min⁡jδj\delta=\min_j\delta_jとおく。∥t−t0∥<δ\lVert t-t_0\rVert<\deltaならば、Γj\Gamma_jの内部にあるA(t)A(t)の固有値の個数はA(t0)A(t_0)のそれに等しく、その和はnnである。Γj\Gamma_jの内部は互いに素であり、A(t)A(t)は代数的重複度を込めてちょうどnn個の固有値をもつから、A(t)A(t)の固有値はすべていずれかのΓj\Gamma_jの内部にある。

主張 (3)を示す。Γ\Gammaの内部にあるA(t)A(t)の固有値をλ\lambdaとすると、その代数的重複度は11である。A(t)A(t)は実行列であるからχA(t)\chi_{A(t)}は実係数であり、χA(t)(λˉ)=χA(t)(λ)‾=0\chi_{A(t)}(\bar\lambda)=\overline{\chi_{A(t)}(\lambda)}=0である。Γ\Gammaは実軸について対称であるからλˉ\bar\lambdaもΓ\Gammaの内部にある。λ∉R\lambda\notin\Rならばλˉ≠λ\bar\lambda\ne\lambdaはΓ\Gammaの内部にある第二の固有値であり、Γ\Gammaの内部の固有値が一個であることに反する。したがってλ∈R\lambda\in\Rである。後半は主張 (1)と前半から従う。

主張 (4)を示す。η=min⁡{∣ℜλ−c∣:λ\eta=\min\{\lvert\Re\lambda-c\rvert:\lambdaはA(t0)A(t_0)の固有値}>0\}>0、ρ=1+∣c∣+η+max⁡{∣λ∣:λ\rho=1+\lvert c\rvert+\eta+\max\{\lvert\lambda\rvert:\lambdaはA(t0)A(t_0)の固有値}\}とし、閉長方形

Lc={z∈C:−ρ≤ℜz≤c−η/2, ∣ℑz∣≤ρ},Rc={z∈C:c+η/2≤ℜz≤ρ, ∣ℑz∣≤ρ}L_c=\{z\in\C:-\rho\le\Re z\le c-\eta/2,\ \lvert\Im z\rvert\le\rho\},\qquad R_c=\{z\in\C:c+\eta/2\le\Re z\le\rho,\ \lvert\Im z\rvert\le\rho\}

を取る。−ρ<c−η/2<c+η/2<ρ-\rho<c-\eta/2<c+\eta/2<\rhoであるから、LcL_cとRcR_cは互いに素である。A(t0)A(t_0)の固有値λ\lambdaは∣ℜλ∣≤∣λ∣<ρ\lvert\Re\lambda\rvert\le\lvert\lambda\rvert<\rho、∣ℑλ∣<ρ\lvert\Im\lambda\rvert<\rhoを満たし、ℜλ≤c−η\Re\lambda\le c-\etaまたはℜλ≥c+η\Re\lambda\ge c+\etaであるから、ℜλ<c\Re\lambda<cならばLcL_cの内部にあり、ℜλ>c\Re\lambda>cならばRcR_cの内部にある。主張 (2)(r=2r=2、Γ1=Lc\Gamma_1=L_c、Γ2=Rc\Gamma_2=R_c)により、δ>0\delta>0が存在して、t∈Tt\in T、∥t−t0∥<δ\lVert t-t_0\rVert<\deltaならば、A(t)A(t)の固有値はすべてLcL_cまたはRcR_cの内部にあり、LcL_cの内部にあるA(t)A(t)の固有値の個数とRcR_cの内部にあるA(t)A(t)の固有値の個数は、それぞれA(t0)A(t_0)のそれに等しい。LcL_cの内部の点の実部はccより小さく、RcR_cの内部の点の実部はccより大きいから、直線{z∈C:ℜz=c}\{z\in\C:\Re z=c\}の上にA(t)A(t)の固有値は無く、ℜz<c\Re z<cにあるA(t)A(t)の固有値はLcL_cの内部にあるものであり、ℜz>c\Re z>cにあるA(t)A(t)の固有値はRcR_cの内部にあるものである。t=t0t=t_0についても同じであるから、主張を得る。▨

定理 3.3.n,k∈N≥1n,k\in\NNとし、U⊆RnU\subseteq\R^nを開集合、Λ⊆R\Lambda\subseteq\Rを開区間、F∈Ck(U×Λ,Rn)F\in C^k(U\times\Lambda,\R^n)とする。(x0,μ0)∈U×Λ(x_0,\mu_0)\in U\times\LambdaがF(x0,μ0)=0F(x_0,\mu_0)=0を満たし、DxF(x0,μ0)D_xF(x_0,\mu_0)が可逆であるとする。このとき、次が成り立つ。

  1. μ0\mu_0を含む開区間Λ0⊆Λ\Lambda_0\subseteq\Lambda、x0x_0の開近傍B⊆UB\subseteq UとCkC^k級の写像ξ:Λ0→B\xi:\Lambda_0\to Bが存在して、(y,μ)∈B×Λ0(y,\mu)\in B\times\Lambda_0についてF(y,μ)=0F(y,\mu)=0であることとy=ξ(μ)y=\xi(\mu)であることは同値である。ξ(μ0)=x0\xi(\mu_0)=x_0であり、μ∈Λ0\mu\in\Lambda_0についてDxF(ξ(μ),μ)D_xF(\xi(\mu),\mu)は可逆で、ξ′(μ)=−DxF(ξ(μ),μ)−1DμF(ξ(μ),μ)\xi'(\mu)=-D_xF(\xi(\mu),\mu)^{-1}D_\mu F(\xi(\mu),\mu)である。
  2. 主張 (1)のΛ0\Lambda_0とξ\xiについて、DxF(x0,μ0)D_xF(x_0,\mu_0)が双曲型行列であるとする。このとき、Λ0\Lambda_0を縮めれば、すべてのμ∈Λ0\mu\in\Lambda_0についてDxF(ξ(μ),μ)D_xF(\xi(\mu),\mu)は双曲型行列であり、実部が負の固有値の個数と実部が正の固有値の個数(ともに代数的重複度を込める)はμ\muによらない。したがって、DxF(ξ(μ),μ)D_xF(\xi(\mu),\mu)で不変な部分空間Es,EuE^s,E^uであって、Rn=Es⊕Eu\R^n=E^s\oplus E^uを満たし、DxF(ξ(μ),μ)D_xF(\xi(\mu),\mu)のEsE^sへの制限の固有値の実部がすべて負、EuE^uへの制限の固有値の実部がすべて正であるもの(§E13.12 補題 6.7 (3)により存在する)について、dim⁡Es\dim E^sとdim⁡Eu\dim E^uはμ\muによらない。
  3. (x0,μ0)(x_0,\mu_0)は族x′=F(x,μ)x'=F(x,\mu)の分岐点でない。

証明.主張 (1)を示す。U~={(μ,y)∈R×Rn:(y,μ)∈U×Λ}\tilde U=\{(\mu,y)\in\R\times\R^n:(y,\mu)\in U\times\Lambda\}は開集合であり、F~(μ,y)=F(y,μ)\tilde F(\mu,y)=F(y,\mu)はU~\tilde Uの上のCkC^k級の写像で、F~(μ0,x0)=0\tilde F(\mu_0,x_0)=0、DyF~(μ0,x0)=DxF(x0,μ0)D_y\tilde F(\mu_0,x_0)=D_xF(x_0,\mu_0)は可逆である。補題 2.1(n=1n=1、m=nm=n)により、μ0\mu_0の開近傍AA、x0x_0の開近傍BBとCkC^k級の写像ξ:A→B\xi:A\to Bが存在して、補題 2.1 (1)と補題 2.1 (2)が成り立つ。A∩ΛA\cap\Lambdaに含まれるμ0\mu_0を含む開区間をΛ0\Lambda_0とし、ξ\xiをΛ0\Lambda_0へ制限すれば、主張を得る。

主張 (2)を示す。Aμ=DxF(ξ(μ),μ)A_\mu=D_xF(\xi(\mu),\mu)とおく。DxFD_xFは連続でξ\xiは連続であるから、μ↦Aμ\mu\mapsto A_\muは連続である。Aμ0A_{\mu_0}の固有値のうち実部が負のものの個数をnsn_s、実部が正のものの個数をnun_uとする(代数的重複度を込める)。Aμ0A_{\mu_0}は双曲型であるから、直線{z∈C:ℜz=0}\{z\in\C:\Re z=0\}の上にAμ0A_{\mu_0}の固有値は無い。補題 3.2 (4)(T=Λ0T=\Lambda_0、t0=μ0t_0=\mu_0、c=0c=0)により、Λ0\Lambda_0を縮めれば、すべてのμ∈Λ0\mu\in\Lambda_0について、実部が00のAμA_\muの固有値は無く、実部が負の固有値の個数はnsn_s、実部が正の固有値の個数はnun_uである。したがってAμA_\muは双曲型であり、これらの個数はμ\muによらない。

μ∈Λ0\mu\in\Lambda_0とし、Es,EuE^s,E^uを、AμA_\muで不変でRn=Es⊕Eu\R^n=E^s\oplus E^uを満たし、Aμ∣EsA_\mu|_{E^s}の固有値の実部がすべて負、Aμ∣EuA_\mu|_{E^u}の固有値の実部がすべて正である部分空間とする。そのようなEs,EuE^s,E^uは、§E13.12 補題 6.7(K=R\mathbb K=\R、d=nd=n)をAμA_\muに適用したものとして存在する(§E13.12 補題 6.7 (3))。χAμ=χAμ∣EsχAμ∣Eu\chi_{A_\mu}=\chi_{A_\mu|_{E^s}}\chi_{A_\mu|_{E^u}}であり、χAμ∣Es\chi_{A_\mu|_{E^s}}の根の実部はすべて負、χAμ∣Eu\chi_{A_\mu|_{E^u}}の根の実部はすべて正であるから、dim⁡Es=ns\dim E^s=n_s、dim⁡Eu=nu\dim E^u=n_uである。

主張 (3)を示す。主張 (1)のΛ0\Lambda_0、BB、ξ\xiについて、Λ0\Lambda_0はμ0\mu_0を含む開区間、BBはx0x_0の開近傍、ξ:Λ0→B\xi:\Lambda_0\to BはCkC^k級であるから連続であり、(y,μ)∈B×Λ0(y,\mu)\in B\times\Lambda_0についてF(y,μ)=0F(y,\mu)=0であることとy=ξ(μ)y=\xi(\mu)であることは同値である。FFは連続であるから、定義 1.1により(x0,μ0)(x_0,\mu_0)は分岐点でない。▨

4 単純固有値の摂動

補題 4.1.n∈N≥1n\in\NNとし、A∈Cn×nA\in\C^{n\times n}の固有値λ0\lambda_0が代数的に単純であるとする。このとき、次が成り立つ。

  1. ker⁡(A−λ0In)\ker(A-\lambda_0I_n)は一次元である。
  2. w∗(A−λ0In)=0w^*(A-\lambda_0I_n)=0を満たすw∈Cn∖{0}w\in\C^n\setminus\{0\}が存在し、そのようなすべてのwwとすべてのv0∈ker⁡(A−λ0In)∖{0}v_0\in\ker(A-\lambda_0I_n)\setminus\{0\}についてw∗v0≠0w^*v_0\ne0である。
  3. wwを主張 (2)のものとするとき、im⁡(A−λ0In)={y∈Cn:w∗y=0}\operatorname{im}(A-\lambda_0I_n)=\{y\in\C^n:w^*y=0\}である。
  4. AAが実行列でλ0∈R\lambda_0\in\Rならば、主張 (1)の核の生成元v0v_0と主張 (2)のwwを、実ベクトルに取ることができる。

証明.λ0\lambda_0は固有値であるから、v0∈ker⁡(A−λ0In)∖{0}v_0\in\ker(A-\lambda_0I_n)\setminus\{0\}が存在する。b1=v0b_1=v_0を含むCn\C^nの基底b1,…,bnb_1,\ldots,b_nを取り、P=(b1,…,bn)P=(b_1,\ldots,b_n)とすると、

T=P−1AP=(λ0β∗0A′)(β∈Cn−1, A′∈C(n−1)×(n−1))T=P^{-1}AP=\begin{pmatrix}\lambda_0&\beta^*\\0&A'\end{pmatrix}\qquad(\beta\in\C^{n-1},\ A'\in\C^{(n-1)\times(n-1)})

である。§E3.20 定理 2.2によりχA=χT=(z−λ0)χA′\chi_A=\chi_T=(z-\lambda_0)\chi_{A'}である。λ0\lambda_0はχA\chi_Aの単根であるからχA′(λ0)≠0\chi_{A'}(\lambda_0)\ne0であり、A′−λ0In−1A'-\lambda_0I_{n-1}は可逆である。

主張 (1)を示す。(a,c)∈C×Cn−1(a,c)\in\C\times\C^{n-1}について(T−λ0In)(a,c)=(β∗c,(A′−λ0In−1)c)(T-\lambda_0I_n)(a,c)=(\beta^*c,(A'-\lambda_0I_{n-1})c)であるから、これが00であることからc=0c=0を得る。したがってker⁡(T−λ0In)=Ce1\ker(T-\lambda_0I_n)=\C e_1であり、ker⁡(A−λ0In)=Pker⁡(T−λ0In)=Cv0\ker(A-\lambda_0I_n)=P\ker(T-\lambda_0I_n)=\C v_0である。

主張 (2)を示す。行ベクトルu∗=(α,c∗)u^*=(\alpha,c^*)についてu∗(T−λ0In)=(0,αβ∗+c∗(A′−λ0In−1))u^*(T-\lambda_0I_n)=(0,\alpha\beta^*+c^*(A'-\lambda_0I_{n-1}))であるから、これが00であることはc∗=−αβ∗(A′−λ0In−1)−1c^*=-\alpha\beta^*(A'-\lambda_0I_{n-1})^{-1}と同値である。したがってu≠0u\ne0ならばα≠0\alpha\ne0であり、このようなu∗u^*が存在する。w∗=u∗P−1w^*=u^*P^{-1}とおくと、w∗(A−λ0In)=u∗(T−λ0In)P−1w^*(A-\lambda_0I_n)=u^*(T-\lambda_0I_n)P^{-1}であるから、w≠0w\ne0がw∗(A−λ0In)=0w^*(A-\lambda_0I_n)=0を満たすことと、u∗=w∗P≠0u^*=w^*P\ne0がu∗(T−λ0In)=0u^*(T-\lambda_0I_n)=0を満たすことは同値である。このときw∗v0=u∗P−1Pe1=α≠0w^*v_0=u^*P^{-1}Pe_1=\alpha\ne0である。

主張 (3)を示す。主張 (1)によりrank⁡(A−λ0In)=n−1\operatorname{rank}(A-\lambda_0I_n)=n-1である。w∗(A−λ0In)=0w^*(A-\lambda_0I_n)=0であるからim⁡(A−λ0In)⊆{y:w∗y=0}\operatorname{im}(A-\lambda_0I_n)\subseteq\{y:w^*y=0\}であり、右辺はw≠0w\ne0によりn−1n-1次元であるから、両者は等しい。

主張 (4)を示す。AAとλ0\lambda_0が実のとき、実行列A−λ0InA-\lambda_0I_nは特異であるから、実ベクトルv0≠0v_0\ne0で(A−λ0In)v0=0(A-\lambda_0I_n)v_0=0を満たすものと、実ベクトルw≠0w\ne0でwT(A−λ0In)=0w^{\mathsf T}(A-\lambda_0I_n)=0を満たすものが存在する。これらは複素の条件も満たす。▨

補題 4.2.n,k∈N≥1n,k\in\NNとし、Λ⊆R\Lambda\subseteq\Rを開区間、A∈Ck(Λ,Rn×n)A\in C^k(\Lambda,\R^{n\times n})、μ0∈Λ\mu_0\in\Lambdaとする。λ0∈C\lambda_0\in\CをA(μ0)A(\mu_0)の代数的に単純な固有値とし、v0∈Cn∖{0}v_0\in\C^n\setminus\{0\}はA(μ0)v0=λ0v0A(\mu_0)v_0=\lambda_0v_0を、w∈Cn∖{0}w\in\C^n\setminus\{0\}はw∗A(μ0)=λ0w∗w^*A(\mu_0)=\lambda_0w^*を満たすとする。このときw∗v0≠0w^*v_0\ne0であり、μ0\mu_0を含む開区間Λ0⊆Λ\Lambda_0\subseteq\LambdaとCkC^k級の写像λ:Λ0→C\lambda:\Lambda_0\to\C、v:Λ0→Cnv:\Lambda_0\to\C^nが存在して、次が成り立つ。

  1. λ(μ0)=λ0\lambda(\mu_0)=\lambda_0、v(μ0)=v0v(\mu_0)=v_0であり、μ∈Λ0\mu\in\Lambda_0についてA(μ)v(μ)=λ(μ)v(μ)A(\mu)v(\mu)=\lambda(\mu)v(\mu)である。
  2. λ′(μ0)=w∗A′(μ0)v0w∗v0\lambda'(\mu_0)=\dfrac{w^*A'(\mu_0)v_0}{w^*v_0}である。
  3. r>0r>0が存在して、Λ0\Lambda_0を縮めれば、すべてのμ∈Λ0\mu\in\Lambda_0について、開正方形Q={z∈C:∣ℜ(z−λ0)∣<r, ∣ℑ(z−λ0)∣<r}Q=\{z\in\C:\lvert\Re(z-\lambda_0)\rvert<r,\ \lvert\Im(z-\lambda_0)\rvert<r\}の中にあるA(μ)A(\mu)の固有値はλ(μ)\lambda(\mu)だけであり、λ(μ)\lambda(\mu)は代数的に単純である。
  4. Λ0\Lambda_0を主張 (3)のとおりに縮めたものとし、λ0∈R\lambda_0\in\Rとする。このとき、すべてのμ∈Λ0\mu\in\Lambda_0についてλ(μ)∈R\lambda(\mu)\in\Rである。さらにv0∈Rnv_0\in\R^nならば、Λ0\Lambda_0を縮めればv(μ)∈Rnv(\mu)\in\R^nである。

証明.補題 4.1 (2)によりw∗v0≠0w^*v_0\ne0である。λ0∈R\lambda_0\in\Rかつv0∈Rnv_0\in\R^nの場合には、補題 4.1 (1)と補題 4.1 (4)により、左零空間は一次元であって実ベクトルwRw_Rを含むから、wRTv0≠0w_R^{\mathsf T}v_0\ne0であり、w~=wR/(wRTv0)\tilde w=w_R/(w_R^{\mathsf T}v_0)とおく。それ以外の場合はw~=w/w∗v0‾\tilde w=w/\overline{w^*v_0}とおく。どちらの場合もw~∗A(μ0)=λ0w~∗\tilde w^*A(\mu_0)=\lambda_0\tilde w^*、w~∗v0=1\tilde w^*v_0=1である。

G:Λ×Cn×C→Cn×CG:\Lambda\times\C^n\times\C\to\C^n\times\CをG(μ,v,λ)=(A(μ)v−λv, w~∗v−1)G(\mu,v,\lambda)=(A(\mu)v-\lambda v,\ \tilde w^*v-1)で定める。Cn×C\C^n\times\CをR2n+2\R^{2n+2}とみなすと、GGはCkC^k級であり、G(μ0,v0,λ0)=0G(\mu_0,v_0,\lambda_0)=0である。(v,λ)(v,\lambda)についての微分はC\C線形写像L(h,ν)=((A(μ0)−λ0In)h−νv0, w~∗h)L(h,\nu)=((A(\mu_0)-\lambda_0I_n)h-\nu v_0,\ \tilde w^*h)である。L(h,ν)=0L(h,\nu)=0とする。w~∗(A(μ0)−λ0In)=0\tilde w^*(A(\mu_0)-\lambda_0I_n)=0であるから、第一成分にw~∗\tilde w^*を掛けて−νw~∗v0=−ν=0-\nu\tilde w^*v_0=-\nu=0を得る。したがって(A(μ0)−λ0In)h=0(A(\mu_0)-\lambda_0I_n)h=0であり、補題 4.1 (1)によりh=cv0h=cv_0(c∈Cc\in\C)である。0=w~∗h=c0=\tilde w^*h=cであるからh=0h=0である。LLは単射であるから、有限次元のR\R線形空間の間の同型である。補題 2.1により、μ0\mu_0を含む開区間Λ0\Lambda_0とCkC^k級の写像(v,λ):Λ0→Cn×C(v,\lambda):\Lambda_0\to\C^n\times\Cが存在してG(μ,v(μ),λ(μ))=0G(\mu,v(\mu),\lambda(\mu))=0、(v(μ0),λ(μ0))=(v0,λ0)(v(\mu_0),\lambda(\mu_0))=(v_0,\lambda_0)である。したがって主張 (1)が成り立つ。

主張 (2)を示す。A(μ)v(μ)=λ(μ)v(μ)A(\mu)v(\mu)=\lambda(\mu)v(\mu)をμ0\mu_0で微分するとA′(μ0)v0+A(μ0)v′(μ0)=λ′(μ0)v0+λ0v′(μ0)A'(\mu_0)v_0+A(\mu_0)v'(\mu_0)=\lambda'(\mu_0)v_0+\lambda_0v'(\mu_0)である。w∗w^*を掛け、w∗A(μ0)=λ0w∗w^*A(\mu_0)=\lambda_0w^*を用いるとw∗A′(μ0)v0=λ′(μ0)w∗v0w^*A'(\mu_0)v_0=\lambda'(\mu_0)w^*v_0を得る。

主張 (3)を示す。A(μ0)A(\mu_0)のλ0\lambda_0以外の固有値とλ0\lambda_0との距離の最小値をddとし(他の固有値が無ければd=1d=1)、r=d/2r=d/2とする。閉正方形Qˉ={z:∣ℜ(z−λ0)∣≤r, ∣ℑ(z−λ0)∣≤r}\bar Q=\{z:\lvert\Re(z-\lambda_0)\rvert\le r,\ \lvert\Im(z-\lambda_0)\rvert\le r\}はλ0\lambda_0以外の固有値を含まず、λ0\lambda_0はQQの内部にある。補題 3.2 (1)(T=ΛT=\Lambda、Γ=Qˉ\Gamma=\bar Q)により、Λ0\Lambda_0を縮めれば、すべてのμ∈Λ0\mu\in\Lambda_0についてQQの中にあるA(μ)A(\mu)の固有値は代数的重複度を込めてちょうど11個である。λ\lambdaは連続でλ(μ0)=λ0\lambda(\mu_0)=\lambda_0であるから、Λ0\Lambda_0を縮めればλ(μ)∈Q\lambda(\mu)\in Qであり、λ(μ)\lambda(\mu)はA(μ)A(\mu)の固有値であるから、これがQQの中の唯一の固有値で、代数的に単純である。

主張 (4)を示す。λ0∈R\lambda_0\in\Rとし、μ∈Λ0\mu\in\Lambda_0とする。Qˉ\bar Qは実軸について対称であり、主張 (3)により、その内部QQにあるA(μ)A(\mu)の固有値は代数的重複度を込めてλ(μ)\lambda(\mu)の一個だけであるから、補題 3.2 (3)(T=ΛT=\Lambda、Γ=Qˉ\Gamma=\bar Q)によりλ(μ)∈R\lambda(\mu)\in\Rである。v0∈Rnv_0\in\R^nとする。このときw~\tilde wは実ベクトルであり、v(μ)‾\overline{v(\mu)}はA(μ)v(μ)‾=λ(μ)v(μ)‾A(\mu)\overline{v(\mu)}=\lambda(\mu)\overline{v(\mu)}とw~Tv(μ)‾=1\tilde w^{\mathsf T}\overline{v(\mu)}=1を満たす。A(μ)−λ(μ)InA(\mu)-\lambda(\mu)I_nの核は補題 4.1 (1)により一次元であるから、v(μ)‾=cv(μ)\overline{v(\mu)}=cv(\mu)と書くと、w~T\tilde w^{\mathsf T}を掛けてc=1c=1を得る。したがってv(μ)∈Rnv(\mu)\in\R^nである。▨

5 Lyapunov–Schmidt 還元

定理 5.1.n,k∈N≥1n,k\in\NNとし、U⊆RnU\subseteq\R^nを開集合、Λ⊆R\Lambda\subseteq\Rを開区間、F∈Ck(U×Λ,Rn)F\in C^k(U\times\Lambda,\R^n)とする。(x0,μ0)∈U×Λ(x_0,\mu_0)\in U\times\LambdaがF(x0,μ0)=0F(x_0,\mu_0)=0を満たすとし、A0=DxF(x0,μ0)A_0=D_xF(x_0,\mu_0)、N=ker⁡A0N=\ker A_0、m=dim⁡Nm=\dim N、R=im⁡A0R=\operatorname{im}A_0とおく。Rn=N⊕N′=R⊕R′\R^n=N\oplus N'=R\oplus R'を満たす部分空間N′,R′N',R'を取り、Q:Rn→RnQ:\R^n\to\R^nを像がR′R'、核がRRである射影、P=In−QP=I_n-Qとする。このとき、0∈N0\in Nの開近傍ZZ、μ0\mu_0を含む開区間Λ0⊆Λ\Lambda_0\subseteq\Lambda、0∈N′0\in N'の開近傍YYとCkC^k級の写像W:Z×Λ0→YW:Z\times\Lambda_0\to Yが存在して、x0+Z+Y⊆Ux_0+Z+Y\subseteq Uであり、b(z,μ)=QF(x0+z+W(z,μ),μ)b(z,\mu)=QF(x_0+z+W(z,\mu),\mu)((z,μ)∈Z×Λ0(z,\mu)\in Z\times\Lambda_0)とおくと、次が成り立つ。

  1. W(0,μ0)=0W(0,\mu_0)=0、DzW(0,μ0)=0D_zW(0,\mu_0)=0、b(0,μ0)=0b(0,\mu_0)=0、Dzb(0,μ0)=0D_zb(0,\mu_0)=0である。
  2. V=x0+Z+YV=x_0+Z+Yはx0x_0の開近傍であり、(x,μ)∈V×Λ0(x,\mu)\in V\times\Lambda_0について、F(x,μ)=0F(x,\mu)=0であることと、z∈Zz\in Zが存在してx=x0+z+W(z,μ)x=x_0+z+W(z,\mu)かつb(z,μ)=0b(z,\mu)=0であることは同値である。このときzzはx−x0x-x_0のNN成分である。
  3. (z,μ)∈Z×Λ0(z,\mu)\in Z\times\Lambda_0についてPF(x0+z+W(z,μ),μ)=0PF(x_0+z+W(z,\mu),\mu)=0であり、したがってF(x0+z+W(z,μ),μ)=b(z,μ)∈R′F(x_0+z+W(z,\mu),\mu)=b(z,\mu)\in R'である。

証明.n=mn=mのときA0=0A_0=0、R={0}R=\{0\}、N′={0}N'=\{0\}、Q=InQ=I_nであり、ZZをx0+Z⊆Ux_0+Z\subseteq Uを満たす開近傍、Y={0}Y=\{0\}、W=0W=0とすれば、b(z,μ)=F(x0+z,μ)b(z,\mu)=F(x_0+z,\mu)であって、主張は明らかである。以下m<nm<nとする。

A0A_0のN′N'への制限A0∣N′:N′→RA_0|_{N'}:N'\to Rは、N′∩ker⁡A0=N′∩N={0}N'\cap\ker A_0=N'\cap N=\{0\}により単射であり、dim⁡N′=n−m=dim⁡R\dim N'=n-m=\dim Rであるから同型である。NN、N′N'、RRを基底によってRm\R^m、Rn−m\R^{n-m}、Rn−m\R^{n-m}と同一視する。H(z,μ,y)=PF(x0+z+y,μ)∈RH(z,\mu,y)=PF(x_0+z+y,\mu)\in R(z∈Nz\in N、y∈N′y\in N')は(z,μ,y)(z,\mu,y)の開集合の上のCkC^k級の写像であり、H(0,μ0,0)=PF(x0,μ0)=0H(0,\mu_0,0)=PF(x_0,\mu_0)=0である。η∈N′\eta\in N'についてDyH(0,μ0,0)η=PA0η=A0ηD_yH(0,\mu_0,0)\eta=PA_0\eta=A_0\etaである(A0η∈RA_0\eta\in RでありPPはRRの上で恒等写像である)から、DyH(0,μ0,0)=A0∣N′D_yH(0,\mu_0,0)=A_0|_{N'}は可逆である。補題 2.1(変数(z,μ)∈N×R(z,\mu)\in N\times\R、未知数y∈N′y\in N')により、(0,μ0)(0,\mu_0)の開近傍と00の開近傍YYとCkC^k級の写像WWが存在して、その開近傍の中でH(z,μ,y)=0H(z,\mu,y)=0であることとy=W(z,μ)y=W(z,\mu)であることは同値である。開近傍を、Z×Λ0Z\times\Lambda_0の形の部分集合へ縮める。Z×Λ0×YZ\times\Lambda_0\times YはHHの定義域に含まれるから、x0+Z+Y⊆Ux_0+Z+Y\subseteq Uである。WWは連続でW(0,μ0)=0W(0,\mu_0)=0であるから、ZZとΛ0\Lambda_0を縮めてW(Z×Λ0)⊆YW(Z\times\Lambda_0)\subseteq Yとする。これらの縮小の後も、(z,μ,y)∈Z×Λ0×Y(z,\mu,y)\in Z\times\Lambda_0\times YについてH(z,μ,y)=0H(z,\mu,y)=0であることとy=W(z,μ)y=W(z,\mu)であることは同値である。

主張 (3)はH(z,μ,W(z,μ))=0H(z,\mu,W(z,\mu))=0から従う。

主張 (1)を示す。H(0,μ0,0)=0H(0,\mu_0,0)=0であるからW(0,μ0)=0W(0,\mu_0)=0であり、b(0,μ0)=QF(x0,μ0)=0b(0,\mu_0)=QF(x_0,\mu_0)=0である。H(z,μ0,W(z,μ0))=0H(z,\mu_0,W(z,\mu_0))=0をzzについて00で微分すると、h∈Nh\in NについてPA0(h+DzW(0,μ0)h)=0PA_0(h+D_zW(0,\mu_0)h)=0である。A0h=0A_0h=0でありA0DzW(0,μ0)h∈RA_0D_zW(0,\mu_0)h\in RであるからA0DzW(0,μ0)h=0A_0D_zW(0,\mu_0)h=0であり、DzW(0,μ0)h∈N′∩N={0}D_zW(0,\mu_0)h\in N'\cap N=\{0\}である。またQA0=0QA_0=0(im⁡A0=R=ker⁡Q\operatorname{im}A_0=R=\ker Q)であるから、Dzb(0,μ0)h=QA0(h+DzW(0,μ0)h)=0D_zb(0,\mu_0)h=QA_0(h+D_zW(0,\mu_0)h)=0である。

主張 (2)を示す。N⊕N′→RnN\oplus N'\to\R^nは線形同型であるから、VVは開集合であり、x∈Vx\in Vはx=x0+z+yx=x_0+z+y(z∈Zz\in Z、y∈Yy\in Y)とただ一通りに書かれ、zzはx−x0x-x_0のNN成分である。F(x,μ)=0F(x,\mu)=0ならばH(z,μ,y)=PF(x,μ)=0H(z,\mu,y)=PF(x,\mu)=0であるからy=W(z,μ)y=W(z,\mu)であり、b(z,μ)=QF(x,μ)=0b(z,\mu)=QF(x,\mu)=0である。逆にx=x0+z+W(z,μ)x=x_0+z+W(z,\mu)かつb(z,μ)=0b(z,\mu)=0ならば、F(x,μ)=PF(x,μ)+QF(x,μ)=0+b(z,μ)=0F(x,\mu)=PF(x,\mu)+QF(x,\mu)=0+b(z,\mu)=0である。▨

定義 5.2.n∈N≥1n\in\NNとし、A∈Rn×nA\in\R^{n\times n}の固有値00が代数的に単純であるとする。v,w∈Rnv,w\in\R^nがAv=0Av=0、wTA=0w^{\mathsf T}A=0、wTv=1w^{\mathsf T}v=1を満たすとき、組(v,w)(v,w)をAAの零固有対 (zero eigenpair) という。零固有対は存在する(補題 4.1 (4)と補題 4.1 (2))。

U⊆RnU\subseteq\R^nを開集合、Λ⊆R\Lambda\subseteq\Rを開区間、F∈C2(U×Λ,Rn)F\in C^2(U\times\Lambda,\R^n)、(x0,μ0)∈U×Λ(x_0,\mu_0)\in U\times\Lambdaとする。ξ1,ξ2∈Rn\xi_1,\xi_2\in\R^nについて、Dx2F(x0,μ0)[ξ1,ξ2]=∂σ1∂σ2F(x0+σ1ξ1+σ2ξ2,μ0)∣σ1=σ2=0D_x^2F(x_0,\mu_0)[\xi_1,\xi_2]=\partial_{\sigma_1}\partial_{\sigma_2}F(x_0+\sigma_1\xi_1+\sigma_2\xi_2,\mu_0)\big|_{\sigma_1=\sigma_2=0}と書く。v∈Rnv\in\R^nについてDμDxF(x0,μ0)[v]=ddμ(DxF(x0,μ)v)∣μ=μ0D_\mu D_xF(x_0,\mu_0)[v]=\dfrac{d}{d\mu}\bigl(D_xF(x_0,\mu)v\bigr)\Big|_{\mu=\mu_0}と書く。F∈C3F\in C^3のとき、Dx3F(x0,μ0)[v,v,v]=d3dσ3F(x0+σv,μ0)∣σ=0D_x^3F(x_0,\mu_0)[v,v,v]=\dfrac{d^3}{d\sigma^3}F(x_0+\sigma v,\mu_0)\Big|_{\sigma=0}と書く。

系 5.3.n,k∈N≥1n,k\in\NNとし、U⊆RnU\subseteq\R^nを開集合、Λ⊆R\Lambda\subseteq\Rを開区間、F∈Ck(U×Λ,Rn)F\in C^k(U\times\Lambda,\R^n)とし、F(x0,μ0)=0F(x_0,\mu_0)=0を満たす(x0,μ0)∈U×Λ(x_0,\mu_0)\in U\times\Lambdaについて、A0=DxF(x0,μ0)A_0=D_xF(x_0,\mu_0)の固有値00が代数的に単純であるとする。(v,w)(v,w)をA0A_0の零固有対とする。このとき、ε0>0\varepsilon_0>0、μ0\mu_0を含む開区間Λ0⊆Λ\Lambda_0\subseteq\Lambda、x0x_0の開近傍V⊆UV\subseteq Uと、CkC^k級の写像W:(−ε0,ε0)×Λ0→ker⁡wTW:(-\varepsilon_0,\varepsilon_0)\times\Lambda_0\to\ker w^{\mathsf T}が存在して、x(s,μ)=x0+sv+W(s,μ)x(s,\mu)=x_0+sv+W(s,\mu)、g(s,μ)=wTF(x(s,μ),μ)g(s,\mu)=w^{\mathsf T}F(x(s,\mu),\mu)とおくと、次が成り立つ。

  1. (s,μ)∈(−ε0,ε0)×Λ0(s,\mu)\in(-\varepsilon_0,\varepsilon_0)\times\Lambda_0についてx(s,μ)∈Vx(s,\mu)\in Vである。(y,μ)∈V×Λ0(y,\mu)\in V\times\Lambda_0について、F(y,μ)=0F(y,\mu)=0であることと、s∈(−ε0,ε0)s\in(-\varepsilon_0,\varepsilon_0)が存在してy=x(s,μ)y=x(s,\mu)かつg(s,μ)=0g(s,\mu)=0であることは同値である。このときs=wT(y−x0)s=w^{\mathsf T}(y-x_0)である。
  2. F(x(s,μ),μ)=g(s,μ)vF(x(s,\mu),\mu)=g(s,\mu)v((s,μ)∈(−ε0,ε0)×Λ0(s,\mu)\in(-\varepsilon_0,\varepsilon_0)\times\Lambda_0)である。
  3. W(0,μ0)=0W(0,\mu_0)=0、∂sW(0,μ0)=0\partial_sW(0,\mu_0)=0、g(0,μ0)=0g(0,\mu_0)=0、∂sg(0,μ0)=0\partial_sg(0,\mu_0)=0、∂μg(0,μ0)=wTDμF(x0,μ0)\partial_\mu g(0,\mu_0)=w^{\mathsf T}D_\mu F(x_0,\mu_0)である。
  4. k≥2k\ge2ならば、∂s2g(0,μ0)=wTDx2F(x0,μ0)[v,v]\partial_s^2g(0,\mu_0)=w^{\mathsf T}D_x^2F(x_0,\mu_0)[v,v]である。さらに∂s2W(0,μ0)\partial_s^2W(0,\mu_0)は、ξ∈ker⁡wT\xi\in\ker w^{\mathsf T}かつA0ξ=−Dx2F(x0,μ0)[v,v]+(wTDx2F(x0,μ0)[v,v])vA_0\xi=-D_x^2F(x_0,\mu_0)[v,v]+\bigl(w^{\mathsf T}D_x^2F(x_0,\mu_0)[v,v]\bigr)vを満たす唯一のξ\xiである。
  5. k≥3k\ge3ならば、ξ=∂s2W(0,μ0)\xi=\partial_s^2W(0,\mu_0)について∂s3g(0,μ0)=wT(Dx3F(x0,μ0)[v,v,v]+3Dx2F(x0,μ0)[v,ξ])\partial_s^3g(0,\mu_0)=w^{\mathsf T}\bigl(D_x^3F(x_0,\mu_0)[v,v,v]+3D_x^2F(x_0,\mu_0)[v,\xi]\bigr)である。
  6. すべてのμ∈Λ\mu\in\LambdaについてF(x0,μ)=0F(x_0,\mu)=0ならば、μ∈Λ0\mu\in\Lambda_0についてW(0,μ)=0W(0,\mu)=0、g(0,μ)=0g(0,\mu)=0であり、k≥2k\ge2ならば∂μg(0,μ0)=0\partial_\mu g(0,\mu_0)=0、∂s∂μg(0,μ0)=wTDμDxF(x0,μ0)[v]\partial_s\partial_\mu g(0,\mu_0)=w^{\mathsf T}D_\mu D_xF(x_0,\mu_0)[v]である。

証明.N=Rv=ker⁡A0N=\R v=\ker A_0である(補題 4.1 (1)の実の核はC\C上の核に含まれ、vvを含む)。wTA0=0w^{\mathsf T}A_0=0であるからR=im⁡A0⊆ker⁡wTR=\operatorname{im}A_0\subseteq\ker w^{\mathsf T}であり、dim⁡R=n−1=dim⁡ker⁡wT\dim R=n-1=\dim\ker w^{\mathsf T}によりR=ker⁡wTR=\ker w^{\mathsf T}である。wTv=1w^{\mathsf T}v=1によりRn=Rv⊕ker⁡wT\R^n=\R v\oplus\ker w^{\mathsf T}であるから、定理 5.1をN′=ker⁡wTN'=\ker w^{\mathsf T}、R′=RvR'=\R vについて適用する。このときQy=(wTy)vQy=(w^{\mathsf T}y)v、Py=y−(wTy)vPy=y-(w^{\mathsf T}y)vである。Z={sv:∣s∣<ε0}Z=\{sv:\lvert s\rvert<\varepsilon_0\}となるようにε0\varepsilon_0を取り、z=svz=svとssを同一視してW(s,μ)W(s,\mu)と書く。b(sv,μ)=(wTF(x(s,μ),μ))v=g(s,μ)vb(sv,\mu)=(w^{\mathsf T}F(x(s,\mu),\mu))v=g(s,\mu)vである。定理 5.1 (2)から主張 (1)が従い(wT(y−x0)=s+wTW(s,μ)=sw^{\mathsf T}(y-x_0)=s+w^{\mathsf T}W(s,\mu)=sによる)、定理 5.1 (3)から主張 (2)が従う。定理 5.1 (1)からW(0,μ0)=0W(0,\mu_0)=0、∂sW(0,μ0)=0\partial_sW(0,\mu_0)=0、g(0,μ0)=0g(0,\mu_0)=0、∂sg(0,μ0)=0\partial_sg(0,\mu_0)=0を得る。

g(s,μ)=wTF(x(s,μ),μ)g(s,\mu)=w^{\mathsf T}F(x(s,\mu),\mu)を微分し、wTA0=0w^{\mathsf T}A_0=0と∂sx(0,μ0)=v+∂sW(0,μ0)=v\partial_sx(0,\mu_0)=v+\partial_sW(0,\mu_0)=vを用いる。∂μg=wT(DxF ∂μx+DμF)\partial_\mu g=w^{\mathsf T}(D_xF\,\partial_\mu x+D_\mu F)に(0,μ0)(0,\mu_0)を代入して∂μg(0,μ0)=wT(A0∂μW(0,μ0)+DμF(x0,μ0))=wTDμF(x0,μ0)\partial_\mu g(0,\mu_0)=w^{\mathsf T}(A_0\partial_\mu W(0,\mu_0)+D_\mu F(x_0,\mu_0))=w^{\mathsf T}D_\mu F(x_0,\mu_0)を得る。これで主張 (3)は示された。

k≥2k\ge2とする。∂s2g=wT(Dx2F[∂sx,∂sx]+DxF ∂s2x)\partial_s^2g=w^{\mathsf T}(D_x^2F[\partial_sx,\partial_sx]+D_xF\,\partial_s^2x)に(0,μ0)(0,\mu_0)を代入して∂s2g(0,μ0)=wTDx2F(x0,μ0)[v,v]\partial_s^2g(0,\mu_0)=w^{\mathsf T}D_x^2F(x_0,\mu_0)[v,v]を得る。主張 (2)をssについて二回微分して(0,μ0)(0,\mu_0)を代入すると、Dx2F(x0,μ0)[v,v]+A0∂s2W(0,μ0)=∂s2g(0,μ0)vD_x^2F(x_0,\mu_0)[v,v]+A_0\partial_s^2W(0,\mu_0)=\partial_s^2g(0,\mu_0)vである(∂sx(0,μ0)=v\partial_sx(0,\mu_0)=v、∂s2x=∂s2W\partial_s^2x=\partial_s^2W)。wTW=0w^{\mathsf T}W=0であるから∂s2W(0,μ0)∈ker⁡wT\partial_s^2W(0,\mu_0)\in\ker w^{\mathsf T}であり、上の式はξ\xiについての方程式を与える。A0A_0をker⁡wT\ker w^{\mathsf T}へ制限した写像は、ker⁡A0∩ker⁡wT=Rv∩ker⁡wT={0}\ker A_0\cap\ker w^{\mathsf T}=\R v\cap\ker w^{\mathsf T}=\{0\}により単射であるから、ξ\xiは唯一である。これで主張 (4)は示された。

k≥3k\ge3とする。∂s3g=wT(Dx3F[∂sx,∂sx,∂sx]+3Dx2F[∂sx,∂s2x]+DxF ∂s3x)\partial_s^3g=w^{\mathsf T}(D_x^3F[\partial_sx,\partial_sx,\partial_sx]+3D_x^2F[\partial_sx,\partial_s^2x]+D_xF\,\partial_s^3x)に(0,μ0)(0,\mu_0)を代入すると、wTA0=0w^{\mathsf T}A_0=0により主張 (5)を得る。

主張 (6)を示す。μ∈Λ0\mu\in\Lambda_0についてF(x0,μ)=0F(x_0,\mu)=0でありx0∈Vx_0\in Vであるから、主張 (1)によりs∈(−ε0,ε0)s\in(-\varepsilon_0,\varepsilon_0)が存在してx0=x(s,μ)x_0=x(s,\mu)、g(s,μ)=0g(s,\mu)=0である。s=wT(x0−x0)=0s=w^{\mathsf T}(x_0-x_0)=0であるから、x0=x0+W(0,μ)x_0=x_0+W(0,\mu)でありW(0,μ)=0W(0,\mu)=0、g(0,μ)=0g(0,\mu)=0である。したがって∂μW(0,μ0)=0\partial_\mu W(0,\mu_0)=0かつ∂μg(0,μ0)=0\partial_\mu g(0,\mu_0)=0である。k≥2k\ge2のとき、g=wTF(x(s,μ),μ)g=w^{\mathsf T}F(x(s,\mu),\mu)をμ\mu、ssの順に微分すると

∂s∂μg=wT(Dx2F[∂sx,∂μx]+DμDxF[∂sx]+DxF ∂s∂μx)\partial_s\partial_\mu g=w^{\mathsf T}\bigl(D_x^2F[\partial_sx,\partial_\mu x]+D_\mu D_xF[\partial_sx]+D_xF\,\partial_s\partial_\mu x\bigr)

である。(0,μ0)(0,\mu_0)では∂sx=v\partial_sx=v、∂μx=∂μW(0,μ0)=0\partial_\mu x=\partial_\mu W(0,\mu_0)=0であるから、∂s∂μg(0,μ0)=wTDμDxF(x0,μ0)[v]\partial_s\partial_\mu g(0,\mu_0)=w^{\mathsf T}D_\mu D_xF(x_0,\mu_0)[v]である。▨

6 サドルノード分岐

補題 6.1.ε1>0\varepsilon_1>0とし、μ^∈C2((−ε1,ε1))\hat\mu\in C^2((-\varepsilon_1,\varepsilon_1))がμ^(0)=μ0\hat\mu(0)=\mu_0、μ^′(0)=0\hat\mu'(0)=0、μ^′′(0)≠0\hat\mu''(0)\ne0を満たすとする。σ=sgn⁡μ^′′(0)\sigma=\operatorname{sgn}\hat\mu''(0)とおく。このとき、ε∈(0,ε1]\varepsilon\in(0,\varepsilon_1]が存在して、次が成り立つ。

  1. 0<∣s∣<ε0<\lvert s\rvert<\varepsilonならばσsμ^′(s)>0\sigma s\hat\mu'(s)>0である。
  2. ε′∈(0,ε]\varepsilon'\in(0,\varepsilon]とする。このとき、δ>0\delta>0が存在して、μ∈(μ0−δ,μ0+δ)\mu\in(\mu_0-\delta,\mu_0+\delta)について、{s∈(−ε′,ε′):μ^(s)=μ}\{s\in(-\varepsilon',\varepsilon'):\hat\mu(s)=\mu\}の元の個数は、σ(μ−μ0)>0\sigma(\mu-\mu_0)>0のとき22(一方は負、他方は正)、μ=μ0\mu=\mu_0のとき11(s=0s=0)、σ(μ−μ0)<0\sigma(\mu-\mu_0)<0のとき00である。

証明.主張 (1)を示す。μ^′′\hat\mu''は連続であるから、ε∈(0,ε1)\varepsilon\in(0,\varepsilon_1)を、σμ^′′>0\sigma\hat\mu''>0(∣s∣<ε\lvert s\rvert<\varepsilon)となるように取る。μ^′(s)=∫0sμ^′′(t) dt\hat\mu'(s)=\int_0^s\hat\mu''(t)\,dtであるから、0<s<ε0<s<\varepsilonならばσμ^′(s)>0\sigma\hat\mu'(s)>0、−ε<s<0-\varepsilon<s<0ならばσμ^′(s)<0\sigma\hat\mu'(s)<0である。

主張 (2)を示す。ε′∈(0,ε]\varepsilon'\in(0,\varepsilon]とする。主張 (1)により、h(s)=σ(μ^(s)−μ0)h(s)=\sigma(\hat\mu(s)-\mu_0)は[0,ε′)[0,\varepsilon')の上で狭義単調増加、(−ε′,0](-\varepsilon',0]の上で狭義単調減少であり、h(0)=0h(0)=0である。したがって0<∣s∣<ε′0<\lvert s\rvert<\varepsilon'ならばh(s)>0h(s)>0である。L+=sup⁡0≤s<ε′h(s)>0L_+=\sup_{0\le s<\varepsilon'}h(s)>0、L−=sup⁡−ε′<s≤0h(s)>0L_-=\sup_{-\varepsilon'<s\le0}h(s)>0(+∞+\inftyであってもよい)とし、δ=min⁡{L+,L−,1}\delta=\min\{L_+,L_-,1\}とおく。0<t<δ0<t<\deltaについて、中間値の定理と狭義単調性により、h(s)=th(s)=tを満たすs∈(0,ε′)s\in(0,\varepsilon')がちょうど一つ、s∈(−ε′,0)s\in(-\varepsilon',0)がちょうど一つ存在する。t=0t=0のときh(s)=0h(s)=0を満たすssは00だけであり、t<0t<0のときh(s)=th(s)=tを満たすssは無い。t=σ(μ−μ0)t=\sigma(\mu-\mu_0)とおけば、主張を得る。▨

定理 6.2.n∈N≥1n\in\NNとし、U⊆RnU\subseteq\R^nを開集合、Λ⊆R\Lambda\subseteq\Rを開区間、F∈C2(U×Λ,Rn)F\in C^2(U\times\Lambda,\R^n)とし、(x0,μ0)∈U×Λ(x_0,\mu_0)\in U\times\LambdaがF(x0,μ0)=0F(x_0,\mu_0)=0を満たすとする。A0=DxF(x0,μ0)A_0=D_xF(x_0,\mu_0)の固有値00が代数的に単純であり、(v,w)(v,w)をA0A_0の零固有対とする。

θμ=wTDμF(x0,μ0)≠0,θxx=wTDx2F(x0,μ0)[v,v]≠0\theta_\mu=w^{\mathsf T}D_\mu F(x_0,\mu_0)\ne0,\qquad \theta_{xx}=w^{\mathsf T}D_x^2F(x_0,\mu_0)[v,v]\ne0

であるとする。x(s,μ)x(s,\mu)、ggを系 5.3(k=2k=2)のものとする。このとき、ε>0\varepsilon>0、x0x_0の開近傍V⊆UV\subseteq U、δ>0\delta>0と、C2C^2級の写像x^:(−ε,ε)→V\hat x:(-\varepsilon,\varepsilon)\to V、μ^:(−ε,ε)→Λ\hat\mu:(-\varepsilon,\varepsilon)\to\Lambdaが存在して、次が成り立つ。

  1. x^(s)=x(s,μ^(s))\hat x(s)=x(s,\hat\mu(s))かつg(s,μ^(s))=0g(s,\hat\mu(s))=0(∣s∣<ε\lvert s\rvert<\varepsilon)であり、x^(0)=x0\hat x(0)=x_0、x^′(0)=v\hat x'(0)=v、μ^(0)=μ0\hat\mu(0)=\mu_0、μ^′(0)=0\hat\mu'(0)=0、μ^′′(0)=−θxx/θμ≠0\hat\mu''(0)=-\theta_{xx}/\theta_\mu\ne0である。
  2. (y,μ)∈V×(μ0−δ,μ0+δ)(y,\mu)\in V\times(\mu_0-\delta,\mu_0+\delta)について、F(y,μ)=0F(y,\mu)=0であることと、s∈(−ε,ε)s\in(-\varepsilon,\varepsilon)が存在して(y,μ)=(x^(s),μ^(s))(y,\mu)=(\hat x(s),\hat\mu(s))であることは同値である。
  3. σ=sgn⁡μ^′′(0)\sigma=\operatorname{sgn}\hat\mu''(0)とする。μ∈(μ0−δ,μ0+δ)\mu\in(\mu_0-\delta,\mu_0+\delta)について、VVの中にあるF(⋅,μ)F(\cdot,\mu)の平衡点の個数は、σ(μ−μ0)>0\sigma(\mu-\mu_0)>0のとき22、μ=μ0\mu=\mu_0のとき11、σ(μ−μ0)<0\sigma(\mu-\mu_0)<0のとき00である。
  4. (x0,μ0)(x_0,\mu_0)は族x′=F(x,μ)x'=F(x,\mu)の分岐点である。

証明.系 5.3(k=2k=2)のε0\varepsilon_0、Λ0\Lambda_0、V0V_0、WW、x(s,μ)x(s,\mu)、ggを取る。系 5.3 (3)と系 5.3 (4)により、g(0,μ0)=0g(0,\mu_0)=0、∂sg(0,μ0)=0\partial_sg(0,\mu_0)=0、∂μg(0,μ0)=θμ≠0\partial_\mu g(0,\mu_0)=\theta_\mu\ne0、∂s2g(0,μ0)=θxx\partial_s^2g(0,\mu_0)=\theta_{xx}である。補題 2.1(パラメータss、未知数μ\mu)により、ε1∈(0,ε0]\varepsilon_1\in(0,\varepsilon_0]、μ0\mu_0を含む開区間M⊆Λ0M\subseteq\Lambda_0とC2C^2級の写像μ^:(−ε1,ε1)→M\hat\mu:(-\varepsilon_1,\varepsilon_1)\to Mが存在して、(s,μ)∈(−ε1,ε1)×M(s,\mu)\in(-\varepsilon_1,\varepsilon_1)\times Mについてg(s,μ)=0g(s,\mu)=0であることとμ=μ^(s)\mu=\hat\mu(s)であることは同値である。μ^(0)=μ0\hat\mu(0)=\mu_0であり、補題 2.1 (2)によりμ^′(s)=−∂sg(s,μ^(s))/∂μg(s,μ^(s))\hat\mu'(s)=-\partial_sg(s,\hat\mu(s))/\partial_\mu g(s,\hat\mu(s))であるからμ^′(0)=0\hat\mu'(0)=0である。g(s,μ^(s))=0g(s,\hat\mu(s))=0を二回微分してs=0s=0を代入すると、∂s2g+2∂s∂μg μ^′+∂μ2g μ^′2+∂μg μ^′′=0\partial_s^2g+2\partial_s\partial_\mu g\,\hat\mu'+\partial_\mu^2g\,\hat\mu'^2+\partial_\mu g\,\hat\mu''=0であり、μ^′(0)=0\hat\mu'(0)=0によりμ^′′(0)=−θxx/θμ≠0\hat\mu''(0)=-\theta_{xx}/\theta_\mu\ne0を得る。

補題 6.1により、ε∈(0,ε1]\varepsilon\in(0,\varepsilon_1]を、その主張が成り立つように取り、δ1>0\delta_1>0を補題 6.1 (2)(ε′=ε\varepsilon'=\varepsilon)のδ\deltaとする。δ>0\delta>0を、δ≤δ1\delta\le\delta_1かつ(μ0−δ,μ0+δ)⊆M(\mu_0-\delta,\mu_0+\delta)\subseteq Mとなるように取り、V={y∈V0:∣wT(y−x0)∣<ε}V=\{y\in V_0:\lvert w^{\mathsf T}(y-x_0)\rvert<\varepsilon\}、x^(s)=x(s,μ^(s))\hat x(s)=x(s,\hat\mu(s))とおく。x^\hat xはC2C^2級であり、wT(x^(s)−x0)=sw^{\mathsf T}(\hat x(s)-x_0)=sであるからx^(s)∈V\hat x(s)\in Vである。x^(0)=x0+W(0,μ0)=x0\hat x(0)=x_0+W(0,\mu_0)=x_0、x^′(0)=v+∂sW(0,μ0)+∂μW(0,μ0)μ^′(0)=v\hat x'(0)=v+\partial_sW(0,\mu_0)+\partial_\mu W(0,\mu_0)\hat\mu'(0)=vである。これで主張 (1)は示された。

主張 (2)を示す。(y,μ)∈V×(μ0−δ,μ0+δ)(y,\mu)\in V\times(\mu_0-\delta,\mu_0+\delta)がF(y,μ)=0F(y,\mu)=0を満たすとする。系 5.3 (1)により、s=wT(y−x0)∈(−ε,ε)s=w^{\mathsf T}(y-x_0)\in(-\varepsilon,\varepsilon)についてy=x(s,μ)y=x(s,\mu)、g(s,μ)=0g(s,\mu)=0であり、μ∈M\mu\in Mであるからμ=μ^(s)\mu=\hat\mu(s)である。したがって(y,μ)=(x^(s),μ^(s))(y,\mu)=(\hat x(s),\hat\mu(s))である。逆に、(y,μ)=(x^(s),μ^(s))(y,\mu)=(\hat x(s),\hat\mu(s))ならば、系 5.3 (2)によりF(y,μ)=g(s,μ^(s))v=0F(y,\mu)=g(s,\hat\mu(s))v=0である。

主張 (3)を示す。主張 (2)により、VVの中のF(⋅,μ)F(\cdot,\mu)の平衡点はx^(s)\hat x(s)(s∈(−ε,ε)s\in(-\varepsilon,\varepsilon)、μ^(s)=μ\hat\mu(s)=\mu)である。wTx^(s)=wTx0+sw^{\mathsf T}\hat x(s)=w^{\mathsf T}x_0+sによりs↦x^(s)s\mapsto\hat x(s)は単射であるから、平衡点の個数は{s∈(−ε,ε):μ^(s)=μ}\{s\in(-\varepsilon,\varepsilon):\hat\mu(s)=\mu\}の元の個数に等しく、補題 6.1 (2)により主張を得る。

主張 (4)を示す。μ0\mu_0を含む開区間Λ′⊆Λ\Lambda'\subseteq\Lambda、x0x_0の開近傍B′⊆UB'\subseteq Uと連続写像ξ:Λ′→B′\xi:\Lambda'\to B'の組であって、(y,μ)∈B′×Λ′(y,\mu)\in B'\times\Lambda'についてF(y,μ)=0F(y,\mu)=0であることとy=ξ(μ)y=\xi(\mu)であることが同値になるものが存在すると仮定する。F(x0,μ0)=0F(x_0,\mu_0)=0であるからξ(μ0)=x0\xi(\mu_0)=x_0である。ξ\xiは連続でVVはx0x_0の開近傍であるから、μ∈Λ′\mu\in\Lambda'を、∣μ−μ0∣<δ\lvert\mu-\mu_0\rvert<\delta、σ(μ−μ0)<0\sigma(\mu-\mu_0)<0、ξ(μ)∈V\xi(\mu)\in Vを満たすように取ることができる。このときξ(μ)\xi(\mu)はVVの中にあるF(⋅,μ)F(\cdot,\mu)の平衡点であり、主張 (3)の、σ(μ−μ0)<0\sigma(\mu-\mu_0)<0のときVVの中の平衡点は零個であることと両立しない。したがって定義 1.1の組は存在せず、(x0,μ0)(x_0,\mu_0)は分岐点である。▨

7 自明な枝をもつ場合の還元

補題 7.1.k,p,m∈N≥1k,p,m\in\NNとし、J⊆RJ\subseteq\Rを00を含む開区間、Λ⊆Rp\Lambda\subseteq\R^pを開集合、g∈Ck(J×Λ,Rm)g\in C^k(J\times\Lambda,\R^m)とし、すべてのμ∈Λ\mu\in\Lambdaについてg(0,μ)=0g(0,\mu)=0であるとする。G(s,μ)=∫01∂sg(τs,μ) dτG(s,\mu)=\displaystyle\int_0^1\partial_sg(\tau s,\mu)\,d\tau((s,μ)∈J×Λ(s,\mu)\in J\times\Lambda。積分は成分ごとにとる)とおく。多重指数β=(β1,…,βp)∈Z≥0p\beta=(\beta_1,\ldots,\beta_p)\in\Z_{\ge0}^pについて∣β∣=β1+⋯+βp\lvert\beta\rvert=\beta_1+\cdots+\beta_p、∂μβ=∂μ1β1⋯∂μpβp\partial_\mu^\beta=\partial_{\mu_1}^{\beta_1}\cdots\partial_{\mu_p}^{\beta_p}と書く。このとき、GGはCk−1C^{k-1}級(k=1k=1のときは連続)であり、次が成り立つ。

  1. g(s,μ)=sG(s,μ)g(s,\mu)=sG(s,\mu)である。
  2. i∈Z≥0i\in\Z_{\ge0}、β∈Z≥0p\beta\in\Z_{\ge0}^p、i+∣β∣≤k−1i+\lvert\beta\rvert\le k-1について、∂si∂μβG(s,μ)=∫01τi ∂si+1∂μβg(τs,μ) dτ\partial_s^i\partial_\mu^\beta G(s,\mu)=\displaystyle\int_0^1\tau^i\,\partial_s^{i+1}\partial_\mu^\beta g(\tau s,\mu)\,d\tauである。特に∂si∂μβG(0,μ)=1i+1∂si+1∂μβg(0,μ)\partial_s^i\partial_\mu^\beta G(0,\mu)=\dfrac{1}{i+1}\partial_s^{i+1}\partial_\mu^\beta g(0,\mu)である。

証明.主張 (1)を示す。(s,μ)∈J×Λ(s,\mu)\in J\times\Lambdaを固定する。JJは00とssを含む区間であるから、{τ∈R:τs∈J}\{\tau\in\R:\tau s\in J\}は[0,1][0,1]を含む開区間であり、その上でggの各成分glg_lについてφl(τ)=gl(τs,μ)\varphi_l(\tau)=g_l(\tau s,\mu)はC1C^1級である。§E4.5 補題 1.2(r=0r=0)によりφl(1)=φl(0)+∫01φl′(τ) dτ\varphi_l(1)=\varphi_l(0)+\int_0^1\varphi_l'(\tau)\,d\tauである。φl(0)=gl(0,μ)=0\varphi_l(0)=g_l(0,\mu)=0、φl′(τ)=s ∂sgl(τs,μ)\varphi_l'(\tau)=s\,\partial_sg_l(\tau s,\mu)であるから、g(s,μ)=sG(s,μ)g(s,\mu)=sG(s,\mu)である。

主張 (2)を示す。i+∣β∣≤k−1i+\lvert\beta\rvert\le k-1についてhi,β(τ,s,μ)=τi∂si+1∂μβg(τs,μ)h_{i,\beta}(\tau,s,\mu)=\tau^i\partial_s^{i+1}\partial_\mu^\beta g(\tau s,\mu)とおく。ggはCkC^k級であるから、hi,βh_{i,\beta}は[0,1]×J×Λ[0,1]\times J\times\Lambdaの上で連続であり、J×ΛJ\times\Lambdaの任意のコンパクト部分集合KKについて[0,1]×K[0,1]\times Kの上で一様連続である。J×ΛJ\times\Lambdaの各点はコンパクトな近傍をもつから、∫01hi,β dτ\int_0^1h_{i,\beta}\,d\tauはJ×ΛJ\times\Lambdaの上で連続である。e1,…,epe_1,\ldots,e_pをZ≥0p\Z_{\ge0}^pの標準的な単位ベクトルとする。i+∣β∣+1≤k−1i+\lvert\beta\rvert+1\le k-1のとき、hi,βh_{i,\beta}のssについての偏導関数はhi+1,βh_{i+1,\beta}、μq\mu_qについての偏導関数は(§E4.4 定理 2.1によりggのkk階までの偏導関数の順序を入れ替えて)hi,β+eqh_{i,\beta+e_q}であり、いずれも連続である。§D1.26 補題 2.1を、各成分について、τ∈[0,1]\tau\in[0,1]を積分変数、ssまたはμq\mu_qを一つのパラメータ(他の変数は固定し、パラメータの区間には、固定した点におけるパラメータの値を含む有界な開区間であって、その閉包の上でパラメータを動かしても点がJ×ΛJ\times\Lambdaにとどまるものを取る)として適用すると、∫01hi,β dτ\int_0^1h_{i,\beta}\,d\tauのssについての偏導関数は∫01hi+1,β dτ\int_0^1h_{i+1,\beta}\,d\tau、μq\mu_qについての偏導関数は∫01hi,β+eq dτ\int_0^1h_{i,\beta+e_q}\,d\tauである。G=∫01h0,0 dτG=\int_0^1h_{0,0}\,d\tauから出発して、偏微分の回数についての帰納法により、GGのk−1k-1階までのすべての反復偏導関数は、偏微分の順序によらず、ssについての回数iiとμ\muについての回数を表すβ\betaだけで定まる∫01hi,β dτ\int_0^1h_{i,\beta}\,d\tauに等しく、連続である。したがって主張 (2)の積分表示が成り立ち、k≥2k\ge2ならば§E4.4 定理 1.7によりGGはCk−1C^{k-1}級である。s=0s=0を代入すると、∫01τi dτ=1/(i+1)\int_0^1\tau^i\,d\tau=1/(i+1)から特別な場合を得る。▨

補題 7.2.n,k∈N≥1n,k\in\NN、k≥2k\ge2とし、U⊆RnU\subseteq\R^nを開集合、Λ⊆R\Lambda\subseteq\Rを開区間、F∈Ck(U×Λ,Rn)F\in C^k(U\times\Lambda,\R^n)とし、x0∈Ux_0\in Uがすべてのμ∈Λ\mu\in\LambdaについてF(x0,μ)=0F(x_0,\mu)=0を満たし、μ0∈Λ\mu_0\in\LambdaについてA0=DxF(x0,μ0)A_0=D_xF(x_0,\mu_0)の固有値00が代数的に単純であるとする。(v,w)(v,w)をA0A_0の零固有対とし、θxμ=wTDμDxF(x0,μ0)[v]≠0\theta_{x\mu}=w^{\mathsf T}D_\mu D_xF(x_0,\mu_0)[v]\ne0であるとする。ε0\varepsilon_0、Λ0\Lambda_0、V0V_0、WW、x(s,μ)x(s,\mu)、ggを系 5.3(同じkk)のε0\varepsilon_0、Λ0\Lambda_0、VV、WW、x(s,μ)x(s,\mu)、ggとし、G(s,μ)=∫01∂sg(τs,μ) dτG(s,\mu)=\displaystyle\int_0^1\partial_sg(\tau s,\mu)\,d\tau((s,μ)∈(−ε0,ε0)×Λ0(s,\mu)\in(-\varepsilon_0,\varepsilon_0)\times\Lambda_0)とおく。このとき、GGはCk−1C^{k-1}級であり、ε1∈(0,ε0]\varepsilon_1\in(0,\varepsilon_0]、μ0\mu_0を含む開区間M⊆Λ0M\subseteq\Lambda_0とCk−1C^{k-1}級の写像μ^:(−ε1,ε1)→M\hat\mu:(-\varepsilon_1,\varepsilon_1)\to Mが存在して、次が成り立つ。

  1. g(s,μ)=sG(s,μ)g(s,\mu)=sG(s,\mu)であり、G(0,μ0)=0G(0,\mu_0)=0、∂μG(0,μ0)=θxμ\partial_\mu G(0,\mu_0)=\theta_{x\mu}、∂sG(0,μ0)=12wTDx2F(x0,μ0)[v,v]\partial_sG(0,\mu_0)=\tfrac12w^{\mathsf T}D_x^2F(x_0,\mu_0)[v,v]である。k≥3k\ge3ならば∂s2G(0,μ0)=13∂s3g(0,μ0)\partial_s^2G(0,\mu_0)=\tfrac13\partial_s^3g(0,\mu_0)である。
  2. (s,μ)∈(−ε1,ε1)×M(s,\mu)\in(-\varepsilon_1,\varepsilon_1)\times Mについて、G(s,μ)=0G(s,\mu)=0であることとμ=μ^(s)\mu=\hat\mu(s)であることは同値である。μ^(0)=μ0\hat\mu(0)=\mu_0、μ^′(0)=−∂sG(0,μ0)/θxμ\hat\mu'(0)=-\partial_sG(0,\mu_0)/\theta_{x\mu}である。
  3. ε∈(0,ε1]\varepsilon\in(0,\varepsilon_1]とし、V={y∈V0:∣wT(y−x0)∣<ε}V=\{y\in V_0:\lvert w^{\mathsf T}(y-x_0)\rvert<\varepsilon\}、x^(s)=x(s,μ^(s))\hat x(s)=x(s,\hat\mu(s))(∣s∣<ε\lvert s\rvert<\varepsilon)とおく。このとき、x^\hat xは(−ε,ε)(-\varepsilon,\varepsilon)からVVへのCk−1C^{k-1}級の単射であり、x^(0)=x0\hat x(0)=x_0、x^′(0)=v\hat x'(0)=v、g(s,μ^(s))=0g(s,\hat\mu(s))=0であり、s≠0s\ne0ならばx^(s)≠x0\hat x(s)\ne x_0である。さらに、(y,μ)∈V×M(y,\mu)\in V\times Mについて、F(y,μ)=0F(y,\mu)=0であることと、y=x0y=x_0であるか、またはs∈(−ε,ε)s\in(-\varepsilon,\varepsilon)が存在して(y,μ)=(x^(s),μ^(s))(y,\mu)=(\hat x(s),\hat\mu(s))であることは同値である。
  4. ε′∈(0,ε0]\varepsilon'\in(0,\varepsilon_0]とし、μ~:(−ε′,ε′)→Λ0\tilde\mu:(-\varepsilon',\varepsilon')\to\Lambda_0を、すべてのs∈(−ε′,ε′)s\in(-\varepsilon',\varepsilon')についてg(s,μ~(s))=0g(s,\tilde\mu(s))=0を満たすC1C^1級の写像とする。このとき、∣s∣<ε′\lvert s\rvert<\varepsilon'についてG(s,μ~(s))=0G(s,\tilde\mu(s))=0であり、∂sg(s,μ~(s))=s ∂sG(s,μ~(s))=−s ∂μG(s,μ~(s)) μ~′(s)\partial_sg(s,\tilde\mu(s))=s\,\partial_sG(s,\tilde\mu(s))=-s\,\partial_\mu G(s,\tilde\mu(s))\,\tilde\mu'(s)である。
  5. (x0,μ0)(x_0,\mu_0)は族x′=F(x,μ)x'=F(x,\mu)の分岐点である。

証明.主張 (1)を示す。系 5.3 (6)により、μ∈Λ0\mu\in\Lambda_0についてW(0,μ)=0W(0,\mu)=0、g(0,μ)=0g(0,\mu)=0であり、∂s∂μg(0,μ0)=θxμ\partial_s\partial_\mu g(0,\mu_0)=\theta_{x\mu}である。補題 7.1(p=m=1p=m=1、J=(−ε0,ε0)J=(-\varepsilon_0,\varepsilon_0)、パラメータの集合Λ0\Lambda_0)により、GGはCk−1C^{k-1}級でg=sGg=sGであり、G(0,μ0)=∂sg(0,μ0)G(0,\mu_0)=\partial_sg(0,\mu_0)、∂μG(0,μ0)=∂s∂μg(0,μ0)\partial_\mu G(0,\mu_0)=\partial_s\partial_\mu g(0,\mu_0)、∂sG(0,μ0)=12∂s2g(0,μ0)\partial_sG(0,\mu_0)=\tfrac12\partial_s^2g(0,\mu_0)、k≥3k\ge3ならば∂s2G(0,μ0)=13∂s3g(0,μ0)\partial_s^2G(0,\mu_0)=\tfrac13\partial_s^3g(0,\mu_0)である。系 5.3 (3)により∂sg(0,μ0)=0\partial_sg(0,\mu_0)=0であり、系 5.3 (4)により∂s2g(0,μ0)=wTDx2F(x0,μ0)[v,v]\partial_s^2g(0,\mu_0)=w^{\mathsf T}D_x^2F(x_0,\mu_0)[v,v]である。

主張 (2)を示す。GGはCk−1C^{k-1}級(k−1≥1k-1\ge1)で、G(0,μ0)=0G(0,\mu_0)=0、∂μG(0,μ0)=θxμ≠0\partial_\mu G(0,\mu_0)=\theta_{x\mu}\ne0である。補題 2.1(Ck−1C^{k-1}級、パラメータss、未知数μ\mu)により、00の開近傍SS、μ0\mu_0の開近傍B⊆Λ0B\subseteq\Lambda_0とCk−1C^{k-1}級の写像μ^:S→B\hat\mu:S\to Bが存在して、(s,μ)∈S×B(s,\mu)\in S\times BについてG(s,μ)=0G(s,\mu)=0であることとμ=μ^(s)\mu=\hat\mu(s)であることは同値であり、μ^(0)=μ0\hat\mu(0)=\mu_0である。M⊆BM\subseteq Bをμ0\mu_0を含む開区間とし、μ^\hat\muの連続性によりε1∈(0,ε0]\varepsilon_1\in(0,\varepsilon_0]を、(−ε1,ε1)⊆S(-\varepsilon_1,\varepsilon_1)\subseteq Sかつμ^((−ε1,ε1))⊆M\hat\mu((-\varepsilon_1,\varepsilon_1))\subseteq Mとなるように取る。(−ε1,ε1)×M⊆S×B(-\varepsilon_1,\varepsilon_1)\times M\subseteq S\times Bであるから、主張の同値が成り立つ。補題 2.1 (2)によりμ^′(0)=−∂sG(0,μ0)/∂μG(0,μ0)=−∂sG(0,μ0)/θxμ\hat\mu'(0)=-\partial_sG(0,\mu_0)/\partial_\mu G(0,\mu_0)=-\partial_sG(0,\mu_0)/\theta_{x\mu}である。

主張 (3)を示す。∣s∣<ε\lvert s\rvert<\varepsilonについてμ^(s)∈M⊆Λ0\hat\mu(s)\in M\subseteq\Lambda_0であるから、系 5.3 (1)によりx^(s)=x(s,μ^(s))∈V0\hat x(s)=x(s,\hat\mu(s))\in V_0であり、xxはCkC^k級、μ^\hat\muはCk−1C^{k-1}級であるから、x^\hat xはCk−1C^{k-1}級である。WWはker⁡wT\ker w^{\mathsf T}に値をとりwTv=1w^{\mathsf T}v=1であるから、wT(x^(s)−x0)=sw^{\mathsf T}(\hat x(s)-x_0)=sである。したがってx^(s)∈V\hat x(s)\in Vであり、x^\hat xは単射であり、s≠0s\ne0ならばx^(s)≠x0\hat x(s)\ne x_0である。x^(0)=x0+W(0,μ0)=x0\hat x(0)=x_0+W(0,\mu_0)=x_0である。W(0,μ)=0W(0,\mu)=0(μ∈Λ0\mu\in\Lambda_0)から∂μW(0,μ0)=0\partial_\mu W(0,\mu_0)=0であり、系 5.3 (3)により∂sW(0,μ0)=0\partial_sW(0,\mu_0)=0であるから、x^′(0)=v+∂sW(0,μ0)+∂μW(0,μ0)μ^′(0)=v\hat x'(0)=v+\partial_sW(0,\mu_0)+\partial_\mu W(0,\mu_0)\hat\mu'(0)=vである。主張 (2)によりG(s,μ^(s))=0G(s,\hat\mu(s))=0であり、g(s,μ^(s))=sG(s,μ^(s))=0g(s,\hat\mu(s))=sG(s,\hat\mu(s))=0である。

(y,μ)∈V×M(y,\mu)\in V\times MがF(y,μ)=0F(y,\mu)=0を満たすとする。系 5.3 (1)により、s=wT(y−x0)∈(−ε,ε)s=w^{\mathsf T}(y-x_0)\in(-\varepsilon,\varepsilon)についてy=x(s,μ)y=x(s,\mu)、g(s,μ)=0g(s,\mu)=0である。s=0s=0ならばy=x0+W(0,μ)=x0y=x_0+W(0,\mu)=x_0である。s≠0s\ne0ならばG(s,μ)=g(s,μ)/s=0G(s,\mu)=g(s,\mu)/s=0であり、(s,μ)∈(−ε1,ε1)×M(s,\mu)\in(-\varepsilon_1,\varepsilon_1)\times Mであるから、主張 (2)によりμ=μ^(s)\mu=\hat\mu(s)、y=x(s,μ^(s))=x^(s)y=x(s,\hat\mu(s))=\hat x(s)である。逆に、仮定によりF(x0,μ)=0F(x_0,\mu)=0であり、(y,μ)=(x^(s),μ^(s))(y,\mu)=(\hat x(s),\hat\mu(s))ならば、系 5.3 (2)によりF(y,μ)=g(s,μ^(s))v=0F(y,\mu)=g(s,\hat\mu(s))v=0である。

主張 (4)を示す。s≠0s\ne0ならばG(s,μ~(s))=g(s,μ~(s))/s=0G(s,\tilde\mu(s))=g(s,\tilde\mu(s))/s=0であり、GGとμ~\tilde\muは連続であるから、s=0s=0でもG(0,μ~(0))=0G(0,\tilde\mu(0))=0である。GGはC1C^1級であるから、G(s,μ~(s))=0G(s,\tilde\mu(s))=0をssについて微分して∂sG(s,μ~(s))+∂μG(s,μ~(s))μ~′(s)=0\partial_sG(s,\tilde\mu(s))+\partial_\mu G(s,\tilde\mu(s))\tilde\mu'(s)=0を得る。g=sGg=sGをssについて微分すると∂sg=G+s ∂sG\partial_sg=G+s\,\partial_sGであるから、∂sg(s,μ~(s))=s ∂sG(s,μ~(s))=−s ∂μG(s,μ~(s)) μ~′(s)\partial_sg(s,\tilde\mu(s))=s\,\partial_sG(s,\tilde\mu(s))=-s\,\partial_\mu G(s,\tilde\mu(s))\,\tilde\mu'(s)である。

主張 (5)を示す。μ0\mu_0を含む開区間Λ′⊆Λ\Lambda'\subseteq\Lambda、x0x_0の開近傍B′⊆UB'\subseteq Uと連続写像ξ:Λ′→B′\xi:\Lambda'\to B'の組であって、(y,μ)∈B′×Λ′(y,\mu)\in B'\times\Lambda'についてF(y,μ)=0F(y,\mu)=0であることとy=ξ(μ)y=\xi(\mu)であることが同値になるものが存在すると仮定する。x0∈B′x_0\in B'であり、すべてのμ∈Λ′\mu\in\Lambda'についてF(x0,μ)=0F(x_0,\mu)=0であるから、ξ(μ)=x0\xi(\mu)=x_0(μ∈Λ′\mu\in\Lambda')である。主張 (3)をε=ε1\varepsilon=\varepsilon_1について適用すると、x^\hat xとμ^\hat\muは連続でx^(0)=x0\hat x(0)=x_0、μ^(0)=μ0\hat\mu(0)=\mu_0であるから、s∈(−ε1,ε1)∖{0}s\in(-\varepsilon_1,\varepsilon_1)\setminus\{0\}をx^(s)∈B′\hat x(s)\in B'、μ^(s)∈Λ′\hat\mu(s)\in\Lambda'となるように取ることができる。系 5.3 (2)によりF(x^(s),μ^(s))=g(s,μ^(s))v=0F(\hat x(s),\hat\mu(s))=g(s,\hat\mu(s))v=0であるから、x^(s)=ξ(μ^(s))=x0\hat x(s)=\xi(\hat\mu(s))=x_0である。一方、s≠0s\ne0であるからx^(s)≠x0\hat x(s)\ne x_0である。この二つは両立しないから、定義 1.1の組は存在せず、(x0,μ0)(x_0,\mu_0)は分岐点である。▨

8 トランスクリティカル分岐

定理 8.1.n∈N≥1n\in\NNとし、U⊆RnU\subseteq\R^nを開集合、Λ⊆R\Lambda\subseteq\Rを開区間、F∈C2(U×Λ,Rn)F\in C^2(U\times\Lambda,\R^n)とし、x0∈Ux_0\in Uがすべてのμ∈Λ\mu\in\LambdaについてF(x0,μ)=0F(x_0,\mu)=0を満たし、μ0∈Λ\mu_0\in\LambdaについてA0=DxF(x0,μ0)A_0=D_xF(x_0,\mu_0)の固有値00が代数的に単純であるとする。(v,w)(v,w)をA0A_0の零固有対とする。

θxμ=wTDμDxF(x0,μ0)[v]≠0,θxx=wTDx2F(x0,μ0)[v,v]≠0\theta_{x\mu}=w^{\mathsf T}D_\mu D_xF(x_0,\mu_0)[v]\ne0,\qquad \theta_{xx}=w^{\mathsf T}D_x^2F(x_0,\mu_0)[v,v]\ne0

であるとする。x(s,μ)x(s,\mu)、ggを系 5.3(k=2k=2)のものとする。このとき、ε>0\varepsilon>0、x0x_0の開近傍V⊆UV\subseteq U、δ>0\delta>0と、C1C^1級の写像x^:(−ε,ε)→V\hat x:(-\varepsilon,\varepsilon)\to V、μ^:(−ε,ε)→Λ\hat\mu:(-\varepsilon,\varepsilon)\to\Lambdaが存在して、次が成り立つ。

  1. x^(s)=x(s,μ^(s))\hat x(s)=x(s,\hat\mu(s))かつg(s,μ^(s))=0g(s,\hat\mu(s))=0(∣s∣<ε\lvert s\rvert<\varepsilon)であり、x^(0)=x0\hat x(0)=x_0、x^′(0)=v\hat x'(0)=v、μ^(0)=μ0\hat\mu(0)=\mu_0、μ^′(0)=−θxx/(2θxμ)≠0\hat\mu'(0)=-\theta_{xx}/(2\theta_{x\mu})\ne0であり、s≠0s\ne0ならばx^(s)≠x0\hat x(s)\ne x_0である。
  2. (y,μ)∈V×(μ0−δ,μ0+δ)(y,\mu)\in V\times(\mu_0-\delta,\mu_0+\delta)について、F(y,μ)=0F(y,\mu)=0であることと、y=x0y=x_0であるか、またはs∈(−ε,ε)s\in(-\varepsilon,\varepsilon)が存在して(y,μ)=(x^(s),μ^(s))(y,\mu)=(\hat x(s),\hat\mu(s))であることは同値である。
  3. μ∈(μ0−δ,μ0+δ)\mu\in(\mu_0-\delta,\mu_0+\delta)について、μ≠μ0\mu\ne\mu_0ならば、VVの中にあるF(⋅,μ)F(\cdot,\mu)の平衡点はちょうど二個であり、そのうち一つはx0x_0、他の一つはμ^(s)=μ\hat\mu(s)=\muを満たすs≠0s\ne0についてx^(s)\hat x(s)である。μ=μ0\mu=\mu_0ならば、VVの中にある平衡点はx0x_0だけである。
  4. (x0,μ0)(x_0,\mu_0)は族x′=F(x,μ)x'=F(x,\mu)の分岐点である。

証明.補題 7.2(k=2k=2)のε0\varepsilon_0、Λ0\Lambda_0、V0V_0、x(s,μ)x(s,\mu)、gg、GG、ε1\varepsilon_1、MM、μ^\hat\muを取る。補題 7.2 (1)と補題 7.2 (2)により、μ^\hat\muはC1C^1級であり、μ^(0)=μ0\hat\mu(0)=\mu_0、μ^′(0)=−∂sG(0,μ0)/θxμ=−θxx/(2θxμ)≠0\hat\mu'(0)=-\partial_sG(0,\mu_0)/\theta_{x\mu}=-\theta_{xx}/(2\theta_{x\mu})\ne0である。μ^′\hat\mu'は連続であるから、ε∈(0,ε1]\varepsilon\in(0,\varepsilon_1]を、μ^′\hat\mu'が(−ε,ε)(-\varepsilon,\varepsilon)で零にならないように取る。このときμ^\hat\muは(−ε,ε)(-\varepsilon,\varepsilon)の上で狭義単調であり、μ^((−ε,ε))\hat\mu((-\varepsilon,\varepsilon))はμ0\mu_0を含む開区間である。δ>0\delta>0を、(μ0−δ,μ0+δ)⊆M∩μ^((−ε,ε))(\mu_0-\delta,\mu_0+\delta)\subseteq M\cap\hat\mu((-\varepsilon,\varepsilon))となるように取り、VV、x^\hat xを補題 7.2 (3)のもの(このε\varepsilonについて)とする。同項により主張 (1)が成り立つ。

(μ0−δ,μ0+δ)⊆M(\mu_0-\delta,\mu_0+\delta)\subseteq Mであるから、補題 7.2 (3)の同値から主張 (2)が従う。

主張 (3)を示す。μ∈(μ0−δ,μ0+δ)\mu\in(\mu_0-\delta,\mu_0+\delta)はμ^((−ε,ε))\hat\mu((-\varepsilon,\varepsilon))に属し、μ^\hat\muは狭義単調であるから、μ^(s)=μ\hat\mu(s)=\muを満たすs∈(−ε,ε)s\in(-\varepsilon,\varepsilon)はただ一つである。μ≠μ0=μ^(0)\mu\ne\mu_0=\hat\mu(0)ならばこのssは00でなく、x^(s)≠x0\hat x(s)\ne x_0であるから、平衡点は主張 (2)によりちょうど二個である。μ=μ0\mu=\mu_0ならばs=0s=0であり、x^(0)=x0\hat x(0)=x_0である。

本定理の仮定は補題 7.2(k=2k=2)の仮定を含むから、補題 7.2 (5)により主張 (4)が成り立つ。▨

9 ピッチフォーク分岐

定理 9.1.n∈N≥1n\in\NNとし、U⊆RnU\subseteq\R^nを開集合、Λ⊆R\Lambda\subseteq\Rを開区間、F∈C3(U×Λ,Rn)F\in C^3(U\times\Lambda,\R^n)とし、x0∈Ux_0\in Uがすべてのμ∈Λ\mu\in\LambdaについてF(x0,μ)=0F(x_0,\mu)=0を満たし、μ0∈Λ\mu_0\in\LambdaについてA0=DxF(x0,μ0)A_0=D_xF(x_0,\mu_0)の固有値00が代数的に単純であるとする。(v,w)(v,w)をA0A_0の零固有対とし、wTDx2F(x0,μ0)[v,v]=0w^{\mathsf T}D_x^2F(x_0,\mu_0)[v,v]=0であるとする。ξ∈ker⁡wT\xi\in\ker w^{\mathsf T}をA0ξ=−Dx2F(x0,μ0)[v,v]A_0\xi=-D_x^2F(x_0,\mu_0)[v,v]を満たす唯一の元とし(系 5.3 (4))、

θxμ=wTDμDxF(x0,μ0)[v],θxxx=wT(Dx3F(x0,μ0)[v,v,v]+3Dx2F(x0,μ0)[v,ξ])\theta_{x\mu}=w^{\mathsf T}D_\mu D_xF(x_0,\mu_0)[v],\qquad \theta_{xxx}=w^{\mathsf T}\bigl(D_x^3F(x_0,\mu_0)[v,v,v]+3D_x^2F(x_0,\mu_0)[v,\xi]\bigr)

とおく。θxμ≠0\theta_{x\mu}\ne0かつθxxx≠0\theta_{xxx}\ne0であるとする。x(s,μ)x(s,\mu)、ggを系 5.3(k=3k=3)のものとする。このとき、ε>0\varepsilon>0、x0x_0の開近傍V⊆UV\subseteq U、δ>0\delta>0と、C2C^2級の写像x^:(−ε,ε)→V\hat x:(-\varepsilon,\varepsilon)\to V、μ^:(−ε,ε)→Λ\hat\mu:(-\varepsilon,\varepsilon)\to\Lambdaが存在して、次が成り立つ。

  1. x^(s)=x(s,μ^(s))\hat x(s)=x(s,\hat\mu(s))かつg(s,μ^(s))=0g(s,\hat\mu(s))=0(∣s∣<ε\lvert s\rvert<\varepsilon)であり、x^(0)=x0\hat x(0)=x_0、x^′(0)=v\hat x'(0)=v、μ^(0)=μ0\hat\mu(0)=\mu_0、μ^′(0)=0\hat\mu'(0)=0、μ^′′(0)=−θxxx/(3θxμ)\hat\mu''(0)=-\theta_{xxx}/(3\theta_{x\mu})であり、s≠0s\ne0ならばx^(s)≠x0\hat x(s)\ne x_0である。μ^(s)=μ0+c2s2+o(s2)\hat\mu(s)=\mu_0+c_2s^2+o(s^2)(s→0s\to0)、c2=−θxxx/(6θxμ)≠0c_2=-\theta_{xxx}/(6\theta_{x\mu})\ne0である。
  2. (y,μ)∈V×(μ0−δ,μ0+δ)(y,\mu)\in V\times(\mu_0-\delta,\mu_0+\delta)について、F(y,μ)=0F(y,\mu)=0であることと、y=x0y=x_0であるか、またはs∈(−ε,ε)s\in(-\varepsilon,\varepsilon)が存在して(y,μ)=(x^(s),μ^(s))(y,\mu)=(\hat x(s),\hat\mu(s))であることは同値である。
  3. τ=sgn⁡c2\tau=\operatorname{sgn}c_2とする。μ∈(μ0−δ,μ0+δ)\mu\in(\mu_0-\delta,\mu_0+\delta)について、VVの中にあるF(⋅,μ)F(\cdot,\mu)の平衡点の個数は、τ(μ−μ0)>0\tau(\mu-\mu_0)>0のとき33(x0x_0と、s<0s<0とs>0s>0にそれぞれ対応するx^(s)\hat x(s)の二個)、τ(μ−μ0)≤0\tau(\mu-\mu_0)\le0のとき11(x0x_0)である。
  4. (x0,μ0)(x_0,\mu_0)は族x′=F(x,μ)x'=F(x,\mu)の分岐点である。

証明.補題 7.2(k=3k=3)のε0\varepsilon_0、Λ0\Lambda_0、V0V_0、WW、x(s,μ)x(s,\mu)、gg、GG、ε1\varepsilon_1、MM、μ^\hat\muを取る。wTDx2F(x0,μ0)[v,v]=0w^{\mathsf T}D_x^2F(x_0,\mu_0)[v,v]=0であるから、系 5.3 (4)によりξ=∂s2W(0,μ0)\xi=\partial_s^2W(0,\mu_0)であり、系 5.3 (5)により∂s3g(0,μ0)=θxxx\partial_s^3g(0,\mu_0)=\theta_{xxx}である。補題 7.2 (1)により、GGはC2C^2級であり、∂μG(0,μ0)=θxμ\partial_\mu G(0,\mu_0)=\theta_{x\mu}、∂sG(0,μ0)=0\partial_sG(0,\mu_0)=0、∂s2G(0,μ0)=θxxx/3\partial_s^2G(0,\mu_0)=\theta_{xxx}/3である。補題 7.2 (2)により、μ^\hat\muはC2C^2級であり、μ^(0)=μ0\hat\mu(0)=\mu_0、μ^′(0)=0\hat\mu'(0)=0、∣s∣<ε1\lvert s\rvert<\varepsilon_1についてG(s,μ^(s))=0G(s,\hat\mu(s))=0である。G(s,μ^(s))=0G(s,\hat\mu(s))=0を二回微分してs=0s=0を代入すると、∂s2G+2∂s∂μG μ^′+∂μ2G μ^′2+∂μG μ^′′=0\partial_s^2G+2\partial_s\partial_\mu G\,\hat\mu'+\partial_\mu^2G\,\hat\mu'^2+\partial_\mu G\,\hat\mu''=0であり、μ^′(0)=0\hat\mu'(0)=0によりμ^′′(0)=−∂s2G(0,μ0)/∂μG(0,μ0)=−θxxx/(3θxμ)\hat\mu''(0)=-\partial_s^2G(0,\mu_0)/\partial_\mu G(0,\mu_0)=-\theta_{xxx}/(3\theta_{x\mu})である。μ^\hat\muはC2C^2級であるから、Taylor の定理によりμ^(s)=μ0+12μ^′′(0)s2+o(s2)\hat\mu(s)=\mu_0+\tfrac12\hat\mu''(0)s^2+o(s^2)である。

補題 6.1により、ε∈(0,ε1]\varepsilon\in(0,\varepsilon_1]を、その主張が成り立つように取り(σ=τ\sigma=\tau)、δ1>0\delta_1>0を補題 6.1 (2)(ε′=ε\varepsilon'=\varepsilon)のδ\deltaとする。δ>0\delta>0を、δ≤δ1\delta\le\delta_1、(μ0−δ,μ0+δ)⊆M(\mu_0-\delta,\mu_0+\delta)\subseteq Mとなるように取り、VV、x^\hat xを補題 7.2 (3)のもの(このε\varepsilonについて)とする。同項により、x^\hat xはC2C^2級であり、x^(s)=x(s,μ^(s))\hat x(s)=x(s,\hat\mu(s))、g(s,μ^(s))=0g(s,\hat\mu(s))=0、x^(0)=x0\hat x(0)=x_0、x^′(0)=v\hat x'(0)=vであり、s≠0s\ne0ならばx^(s)≠x0\hat x(s)\ne x_0である。これで主張 (1)は示された。

(μ0−δ,μ0+δ)⊆M(\mu_0-\delta,\mu_0+\delta)\subseteq Mであるから、補題 7.2 (3)の同値から主張 (2)が従う。

主張 (3)を示す。補題 7.2 (3)によりx^\hat xは単射で、s≠0s\ne0ならばx^(s)≠x0\hat x(s)\ne x_0であるから、VVの中の平衡点の個数は、1+#{s∈(−ε,ε)∖{0}:μ^(s)=μ}1+\#\{s\in(-\varepsilon,\varepsilon)\setminus\{0\}:\hat\mu(s)=\mu\}である。補題 6.1 (2)により、この個数はτ(μ−μ0)>0\tau(\mu-\mu_0)>0のとき1+2=31+2=3、τ(μ−μ0)≤0\tau(\mu-\mu_0)\le0のとき11である。

本定理の仮定は補題 7.2(k=3k=3)の仮定を含むから、補題 7.2 (5)により主張 (4)が成り立つ。▨

10 安定性の交代

補題 10.1.n∈N≥1n\in\NNとし、U⊆RnU\subseteq\R^nを開集合、Λ⊆R\Lambda\subseteq\Rを開区間、F∈C2(U×Λ,Rn)F\in C^2(U\times\Lambda,\R^n)とし、F(x0,μ0)=0F(x_0,\mu_0)=0を満たす(x0,μ0)∈U×Λ(x_0,\mu_0)\in U\times\Lambdaについて、A0=DxF(x0,μ0)A_0=D_xF(x_0,\mu_0)の固有値00は代数的に単純であり、00以外の固有値の実部はすべて負であるとする。(v,w)(v,w)をA0A_0の零固有対とし、ε0\varepsilon_0、Λ0\Lambda_0、WW、x(s,μ)x(s,\mu)、ggを系 5.3(k=2k=2)のものとして、A(s,μ)=DxF(x(s,μ),μ)A(s,\mu)=D_xF(x(s,\mu),\mu)とおく。このとき、ε∈(0,ε0]\varepsilon\in(0,\varepsilon_0]、μ0\mu_0を含む開区間Λ1⊆Λ0\Lambda_1\subseteq\Lambda_0とr>0r>0が存在して、(s,μ)∈(−ε,ε)×Λ1(s,\mu)\in(-\varepsilon,\varepsilon)\times\Lambda_1について、次が成り立つ。

  1. A(s,μ)A(s,\mu)の固有値のうち、開正方形{z∈C:∣ℜz∣<r, ∣ℑz∣<r}\{z\in\C:\lvert\Re z\rvert<r,\ \lvert\Im z\rvert<r\}の中にあるものはちょうど一つ(これをλ(s,μ)\lambda(s,\mu)とする)であり、λ(s,μ)\lambda(s,\mu)は実数で代数的に単純である。他のn−1n-1個の固有値(代数的重複度を込める)の実部は−r-rより小さい。
  2. λ(s,μ)\lambda(s,\mu)の符号は∂sg(s,μ)\partial_sg(s,\mu)の符号に等しい。すなわち、λ(s,μ)\lambda(s,\mu)が正、零、負であることは、それぞれ∂sg(s,μ)\partial_sg(s,\mu)が正、零、負であることと同値である。

証明.a=min⁡{−ℜλ:λa=\min\{-\Re\lambda:\lambdaはA0A_0の固有値でλ≠0}>0\lambda\ne0\}>0(他の固有値が無ければa=1a=1)、r=min⁡{a,1}/2r=\min\{a,1\}/2とする。A0A_0の固有値は00と、実部が−a-a以下のもの(−a≤−2r<−r-a\le-2r<-r)であるから、直線{z∈C:ℜz=−r}\{z\in\C:\Re z=-r\}と{z∈C:ℜz=r}\{z\in\C:\Re z=r\}の上にA0A_0の固有値は無く、ℜz>−r\Re z>-rにあるA0A_0の固有値の個数は11、ℜz>r\Re z>rにあるA0A_0の固有値の個数は00である(代数的重複度を込める)。A(s,μ)A(s,\mu)は(0,μ0)(0,\mu_0)で連続であるから、補題 3.2 (4)(T=(−ε0,ε0)×Λ0T=(-\varepsilon_0,\varepsilon_0)\times\Lambda_0、t0=(0,μ0)t_0=(0,\mu_0))をc=−rc=-rとc=rc=rに適用し、得られる二つのδ\deltaの小さい方をδ∗\delta_*とする。ε∈(0,ε0]\varepsilon\in(0,\varepsilon_0]とμ0\mu_0を含む開区間Λ1⊆Λ0\Lambda_1\subseteq\Lambda_0を、(s,μ)∈(−ε,ε)×Λ1(s,\mu)\in(-\varepsilon,\varepsilon)\times\Lambda_1ならば∥(s,μ)−(0,μ0)∥<δ∗\lVert(s,\mu)-(0,\mu_0)\rVert<\delta_*となるように取る。このとき、(s,μ)∈(−ε,ε)×Λ1(s,\mu)\in(-\varepsilon,\varepsilon)\times\Lambda_1について、A(s,μ)A(s,\mu)の固有値のうち、実部が−r-rより大きいものは代数的重複度を込めてちょうど11個であり、実部がrr以上のものは無く、他のn−1n-1個の実部は−r-rより小さい。実部が−r-rより大きい固有値をλ(s,μ)\lambda(s,\mu)とする。A(s,μ)A(s,\mu)は実行列であるから、λ(s,μ)‾\overline{\lambda(s,\mu)}もA(s,μ)A(s,\mu)の固有値であり、その実部はλ(s,μ)\lambda(s,\mu)の実部に等しい。実部が−r-rより大きい固有値は一個であるからλ(s,μ)‾=λ(s,μ)\overline{\lambda(s,\mu)}=\lambda(s,\mu)であり、λ(s,μ)\lambda(s,\mu)は実数で代数的に単純である。−r<λ(s,μ)<r-r<\lambda(s,\mu)<rであるから、λ(s,μ)\lambda(s,\mu)は開正方形{z∈C:∣ℜz∣<r, ∣ℑz∣<r}\{z\in\C:\lvert\Re z\rvert<r,\ \lvert\Im z\rvert<r\}の中にあり、他の固有値の実部は−r-rより小さいから、この開正方形の中にある他の固有値は無い。これで主張 (1)は示された。

主張 (2)を示す。

主張 10.1.1.d∈Z≥0d\in\Z_{\ge0}とし、実係数の多項式p(z)=∏j=1d(z−λj)p(z)=\prod_{j=1}^d(z-\lambda_j)の根λ1,…,λd\lambda_1,\ldots,\lambda_dの実部がすべて負であるとする。このとき、λ1⋯λd\lambda_1\cdots\lambda_dは実数であり、その符号は(−1)d(-1)^dである。

証明.p(zˉ)‾=p(z)\overline{p(\bar z)}=p(z)であるから、実でない根はその共役と同じ重複度をもち、共役な二根の積λλˉ=∣λ∣2\lambda\bar\lambda=\lvert\lambda\rvert^2は正である。実の根は負であり、その個数(重複度を込める)はddから偶数を引いたものである。▨

n=1n=1のとき、ker⁡wT={0}\ker w^{\mathsf T}=\{0\}であるからW=0W=0、x(s,μ)=x0+svx(s,\mu)=x_0+svであり、wv=1wv=1により∂sg(s,μ)=w DxF(x0+sv,μ) v=A(s,μ)\partial_sg(s,\mu)=w\,D_xF(x_0+sv,\mu)\,v=A(s,\mu)である。A(s,μ)A(s,\mu)の固有値はA(s,μ)A(s,\mu)自身であるから、λ(s,μ)=∂sg(s,μ)\lambda(s,\mu)=\partial_sg(s,\mu)である。

n≥2n\ge2とする。e2,…,ene_2,\ldots,e_nをker⁡wT\ker w^{\mathsf T}の基底とし、E=(e2,…,en)∈Rn×(n−1)E=(e_2,\ldots,e_n)\in\R^{n\times(n-1)}とおく。wTv=1w^{\mathsf T}v=1によりv∉ker⁡wTv\notin\ker w^{\mathsf T}であるから、nn次正方行列(v,E)(v,E)は可逆である。wTA0=0w^{\mathsf T}A_0=0であるからA0EA_0Eの列はker⁡wT\ker w^{\mathsf T}に属し、A0E=EBA_0E=EBを満たすB∈R(n−1)×(n−1)B\in\R^{(n-1)\times(n-1)}が存在する。A0v=0A_0v=0であるからA0(v,E)=(v,E)(000B)A_0(v,E)=(v,E)\begin{pmatrix}0&0\\0&B\end{pmatrix}であり、§E3.20 定理 2.2によりχA0(z)=z χB(z)\chi_{A_0}(z)=z\,\chi_B(z)である。00はχA0\chi_{A_0}の単根であるから、χB\chi_Bの根はA0A_0の00以外の固有値(代数的重複度を込める)であり、その実部はすべて負である。det⁡B=(−1)n−1χB(0)\det B=(-1)^{n-1}\chi_B(0)はχB\chi_Bの根の積であるから、主張 10.1.1によりdet⁡B\det Bの符号は(−1)n−1(-1)^{n-1}である。

系 5.3 (2)をssについて微分すると、u(s,μ)=v+∂sW(s,μ)u(s,\mu)=v+\partial_sW(s,\mu)についてA(s,μ)u(s,μ)=∂sg(s,μ)vA(s,\mu)u(s,\mu)=\partial_sg(s,\mu)vである。WWはker⁡wT\ker w^{\mathsf T}に値をとるから∂sW(s,μ)∈ker⁡wT\partial_sW(s,\mu)\in\ker w^{\mathsf T}であり、det⁡(u(s,μ),E)=det⁡(v,E)\det(u(s,\mu),E)=\det(v,E)である。したがって

det⁡A(s,μ) det⁡(v,E)=det⁡(A(s,μ)u(s,μ),A(s,μ)E)=∂sg(s,μ) det⁡(v,A(s,μ)E)\det A(s,\mu)\,\det(v,E)=\det\bigl(A(s,\mu)u(s,\mu),A(s,\mu)E\bigr)=\partial_sg(s,\mu)\,\det\bigl(v,A(s,\mu)E\bigr)

である。β(s,μ)=det⁡(v,A(s,μ)E)/det⁡(v,E)\beta(s,\mu)=\det(v,A(s,\mu)E)/\det(v,E)は連続であり、(v,A0E)=(v,E)(100B)(v,A_0E)=(v,E)\begin{pmatrix}1&0\\0&B\end{pmatrix}によりβ(0,μ0)=det⁡B\beta(0,\mu_0)=\det Bである。ε\varepsilonとΛ1\Lambda_1を縮めて、(s,μ)∈(−ε,ε)×Λ1(s,\mu)\in(-\varepsilon,\varepsilon)\times\Lambda_1についてβ(s,μ)\beta(s,\mu)の符号が(−1)n−1(-1)^{n-1}であるとする。一方、det⁡A(s,μ)\det A(s,\mu)はA(s,μ)A(s,\mu)の固有値の積であり、主張 (1)によりdet⁡A(s,μ)=λ(s,μ) π(s,μ)\det A(s,\mu)=\lambda(s,\mu)\,\pi(s,\mu)と書かれる。ここでπ(s,μ)\pi(s,\mu)は実係数の多項式χA(s,μ)(z)/(z−λ(s,μ))\chi_{A(s,\mu)}(z)/(z-\lambda(s,\mu))の根の積であり、その根の実部はすべて−r-rより小さいから、主張 10.1.1によりπ(s,μ)\pi(s,\mu)の符号は(−1)n−1(-1)^{n-1}である。det⁡A(s,μ)=∂sg(s,μ) β(s,μ)\det A(s,\mu)=\partial_sg(s,\mu)\,\beta(s,\mu)であるから、λ(s,μ) π(s,μ)=∂sg(s,μ) β(s,μ)\lambda(s,\mu)\,\pi(s,\mu)=\partial_sg(s,\mu)\,\beta(s,\mu)であり、π(s,μ)\pi(s,\mu)とβ(s,μ)\beta(s,\mu)の符号は等しく零でないから、λ(s,μ)\lambda(s,\mu)と∂sg(s,μ)\partial_sg(s,\mu)の符号は一致する。▨

系 10.2.n∈N≥1n\in\NNとし、U⊆RnU\subseteq\R^nを開集合、Λ⊆R\Lambda\subseteq\Rを開区間、F∈C2(U×Λ,Rn)F\in C^2(U\times\Lambda,\R^n)とし、F(x0,μ0)=0F(x_0,\mu_0)=0を満たす(x0,μ0)∈U×Λ(x_0,\mu_0)\in U\times\Lambdaについて、A0=DxF(x0,μ0)A_0=D_xF(x_0,\mu_0)の固有値00は代数的に単純であり、00以外の固有値の実部はすべて負であるとする。(v,w)(v,w)をA0A_0の零固有対とし、ε0\varepsilon_0、Λ0\Lambda_0、x(s,μ)x(s,\mu)、ggを系 5.3(k=2k=2)のものとする。ただし主張 (5)では、定理 9.1の仮定によりF∈C3(U×Λ,Rn)F\in C^3(U\times\Lambda,\R^n)であるから、ε0\varepsilon_0、Λ0\Lambda_0、x(s,μ)x(s,\mu)、ggを系 5.3(k=3k=3)のものとする。これらはk=2k=2の場合の結論も満たす。ε1\varepsilon_1、Λ1\Lambda_1を補題 10.1のε\varepsilon、Λ1\Lambda_1とする。このとき、次が成り立つ。

  1. (s,μ)∈(−ε1,ε1)×Λ1(s,\mu)\in(-\varepsilon_1,\varepsilon_1)\times\Lambda_1がg(s,μ)=0g(s,\mu)=0を満たすとき、x(s,μ)x(s,\mu)はF(⋅,μ)F(\cdot,\mu)の平衡点である。∂sg(s,μ)<0\partial_sg(s,\mu)<0ならばこの平衡点は漸近安定であり、∂sg(s,μ)>0\partial_sg(s,\mu)>0ならば不安定である。
  2. 定理 6.2の仮定が成り立つとし、同定理のx^\hat x、μ^\hat\mu、ε\varepsilon、VVを取り、σ=sgn⁡μ^′′(0)\sigma=\operatorname{sgn}\hat\mu''(0)とする。このとき、ε′∈(0,min⁡{ε,ε1}]\varepsilon'\in(0,\min\{\varepsilon,\varepsilon_1\}]が存在して、0<∣s∣<ε′0<\lvert s\rvert<\varepsilon'について、x^(s)\hat x(s)はF(⋅,μ^(s))F(\cdot,\hat\mu(s))の平衡点であり、sθxx<0s\theta_{xx}<0ならば漸近安定、sθxx>0s\theta_{xx}>0ならば不安定である。さらに、δ′>0\delta'>0が存在して、0<σ(μ−μ0)<δ′0<\sigma(\mu-\mu_0)<\delta'を満たすすべてのμ\muについて、VVの中にあるF(⋅,μ)F(\cdot,\mu)の二個の平衡点の一方は漸近安定であり、他方は不安定である。
  3. すべてのμ∈Λ\mu\in\LambdaについてF(x0,μ)=0F(x_0,\mu)=0であり、θxμ=wTDμDxF(x0,μ0)[v]≠0\theta_{x\mu}=w^{\mathsf T}D_\mu D_xF(x_0,\mu_0)[v]\ne0であるとする。このとき、δ0>0\delta_0>0が存在して、0<∣μ−μ0∣<δ00<\lvert\mu-\mu_0\rvert<\delta_0を満たすすべてのμ\muについて、x0x_0は、θxμ(μ−μ0)<0\theta_{x\mu}(\mu-\mu_0)<0ならばF(⋅,μ)F(\cdot,\mu)の漸近安定な平衡点、θxμ(μ−μ0)>0\theta_{x\mu}(\mu-\mu_0)>0ならば不安定な平衡点である。
  4. 定理 8.1の仮定が成り立つとし、同定理のx^\hat x、μ^\hat\mu、ε\varepsilon、δ\delta、VVを取る。このとき、δ′∈(0,δ]\delta'\in(0,\delta]が存在して、0<∣μ−μ0∣<δ′0<\lvert\mu-\mu_0\rvert<\delta'を満たすすべてのμ\muについて、x0x_0は、θxμ(μ−μ0)<0\theta_{x\mu}(\mu-\mu_0)<0ならばF(⋅,μ)F(\cdot,\mu)の漸近安定な平衡点、θxμ(μ−μ0)>0\theta_{x\mu}(\mu-\mu_0)>0ならば不安定な平衡点である。VVの中にあるF(⋅,μ)F(\cdot,\mu)の他方の平衡点x^(s)\hat x(s)(s∈(−ε,ε)s\in(-\varepsilon,\varepsilon)、μ^(s)=μ\hat\mu(s)=\mu)は、θxμ(μ−μ0)>0\theta_{x\mu}(\mu-\mu_0)>0ならば漸近安定、θxμ(μ−μ0)<0\theta_{x\mu}(\mu-\mu_0)<0ならば不安定である。
  5. 定理 9.1の仮定が成り立つとし、同定理のx^\hat x、μ^\hat\mu、ε\varepsilon、δ\delta、VV、τ\tauを取る。このとき、ε′∈(0,min⁡{ε,ε1}]\varepsilon'\in(0,\min\{\varepsilon,\varepsilon_1\}]とδ′∈(0,δ]\delta'\in(0,\delta]が存在して、次が成り立つ。0<∣μ−μ0∣<δ′0<\lvert\mu-\mu_0\rvert<\delta'について、x0x_0は、θxμ(μ−μ0)<0\theta_{x\mu}(\mu-\mu_0)<0ならばF(⋅,μ)F(\cdot,\mu)の漸近安定な平衡点、θxμ(μ−μ0)>0\theta_{x\mu}(\mu-\mu_0)>0ならば不安定な平衡点である。0<∣s∣<ε′0<\lvert s\rvert<\varepsilon'について、x^(s)\hat x(s)は、θxxx<0\theta_{xxx}<0ならばF(⋅,μ^(s))F(\cdot,\hat\mu(s))の漸近安定な平衡点、θxxx>0\theta_{xxx}>0ならば不安定な平衡点である。さらに、0<τ(μ−μ0)<δ′0<\tau(\mu-\mu_0)<\delta'を満たすすべてのμ\muについて、VVの中にあるF(⋅,μ)F(\cdot,\mu)の平衡点は、x0x_0と、0<∣s∣<ε′0<\lvert s\rvert<\varepsilon'かつμ^(s)=μ\hat\mu(s)=\muを満たす二つのssについてのx^(s)\hat x(s)であり、θxxx>0\theta_{xxx}>0ならばx0x_0は漸近安定、二つのx^(s)\hat x(s)は不安定であり、θxxx<0\theta_{xxx}<0ならばx0x_0は不安定、二つのx^(s)\hat x(s)は漸近安定である。

証明.主張 (1)を示す。rrを補題 10.1のrrとする。系 5.3 (2)によりF(x(s,μ),μ)=g(s,μ)v=0F(x(s,\mu),\mu)=g(s,\mu)v=0である。DxF(x(s,μ),μ)=A(s,μ)D_xF(x(s,\mu),\mu)=A(s,\mu)の固有値は、補題 10.1 (1)により、実数λ(s,μ)\lambda(s,\mu)と実部が−r-rより小さい固有値である。F(⋅,μ)F(\cdot,\mu)はUU上のC1C^1級の写像であるから、§E13.4 補題 4.2(Ω=R×U\Omega=\R\times U、f(t,y)=F(y,μ)f(t,y)=F(y,\mu))により局所 Lipschitz 連続であり、§E13.16 定理 4.5をW=UW=Uとして適用することができる。∂sg(s,μ)<0\partial_sg(s,\mu)<0ならば、補題 10.1 (2)によりλ(s,μ)<0\lambda(s,\mu)<0であり、すべての固有値の実部が負であるから、§E13.16 定理 4.5 (1)によりx(s,μ)x(s,\mu)は漸近安定である。∂sg(s,μ)>0\partial_sg(s,\mu)>0ならばλ(s,μ)>0\lambda(s,\mu)>0は実部が正の固有値であるから、§E13.16 定理 4.5 (2)によりx(s,μ)x(s,\mu)は不安定である。

主張 (2)を示す。定理 6.2 (1)により、x^(s)=x(s,μ^(s))\hat x(s)=x(s,\hat\mu(s))、g(s,μ^(s))=0g(s,\hat\mu(s))=0、μ^(0)=μ0\hat\mu(0)=\mu_0、μ^′′(0)=−θxx/θμ≠0\hat\mu''(0)=-\theta_{xx}/\theta_\mu\ne0であるから、σ=−sgn⁡(θxx)sgn⁡(θμ)\sigma=-\operatorname{sgn}(\theta_{xx})\operatorname{sgn}(\theta_\mu)である。系 5.3 (3)により∂μg(0,μ0)=θμ≠0\partial_\mu g(0,\mu_0)=\theta_\mu\ne0であり、∂μg\partial_\mu gとμ^\hat\muは連続である。補題 6.1をμ^\hat\mu(同補題のε1\varepsilon_1を同定理のε\varepsilonとする)に適用し、同補題のε\varepsilonをε2\varepsilon_2とする。ε′∈(0,min⁡{ε2,ε1}]\varepsilon'\in(0,\min\{\varepsilon_2,\varepsilon_1\}]を、∣s∣<ε′\lvert s\rvert<\varepsilon'についてμ^(s)∈Λ1\hat\mu(s)\in\Lambda_1かつsgn⁡∂μg(s,μ^(s))=sgn⁡θμ\operatorname{sgn}\partial_\mu g(s,\hat\mu(s))=\operatorname{sgn}\theta_\muとなるように取る。補題 6.1 (1)により、0<∣s∣<ε′0<\lvert s\rvert<\varepsilon'についてσsμ^′(s)>0\sigma s\hat\mu'(s)>0である。g(s,μ^(s))=0g(s,\hat\mu(s))=0を微分して∂sg(s,μ^(s))=−∂μg(s,μ^(s))μ^′(s)\partial_sg(s,\hat\mu(s))=-\partial_\mu g(s,\hat\mu(s))\hat\mu'(s)であるから、0<∣s∣<ε′0<\lvert s\rvert<\varepsilon'について

sgn⁡∂sg(s,μ^(s))=−sgn⁡(θμ) σsgn⁡(s)=sgn⁡(s)sgn⁡(θxx)\operatorname{sgn}\partial_sg(s,\hat\mu(s))=-\operatorname{sgn}(\theta_\mu)\,\sigma\operatorname{sgn}(s)=\operatorname{sgn}(s)\operatorname{sgn}(\theta_{xx})

である。主張 (1)により、0<∣s∣<ε′0<\lvert s\rvert<\varepsilon'についての主張を得る。ε′≤ε2\varepsilon'\le\varepsilon_2であるから、δ′′>0\delta''>0を補題 6.1 (2)(このε′\varepsilon')のδ\deltaとし、δ\deltaを定理 6.2のδ\deltaとしてδ′=min⁡{δ,δ′′}\delta'=\min\{\delta,\delta''\}とおく。0<σ(μ−μ0)<δ′0<\sigma(\mu-\mu_0)<\delta'とする。補題 6.1 (2)により、μ^(s)=μ\hat\mu(s)=\muを満たすs∈(−ε′,ε′)s\in(-\varepsilon',\varepsilon')は、負のs−s_-と正のs+s_+の二つである。定理 6.2 (2)によりx^(s±)\hat x(s_\pm)はVVの中にあるF(⋅,μ)F(\cdot,\mu)の平衡点であり、系 5.3 (1)によりwT(x^(s±)−x0)=s±w^{\mathsf T}(\hat x(s_\pm)-x_0)=s_\pmであるからx^(s−)≠x^(s+)\hat x(s_-)\ne\hat x(s_+)である。定理 6.2 (3)により、VVの中にあるF(⋅,μ)F(\cdot,\mu)の平衡点はx^(s−)\hat x(s_-)とx^(s+)\hat x(s_+)だけである。s−θxxs_-\theta_{xx}とs+θxxs_+\theta_{xx}の符号は反対であるから、一方は漸近安定であり、他方は不安定である。

主張 (3)を示す。h(μ)=∂sg(0,μ)h(\mu)=\partial_sg(0,\mu)(μ∈Λ0\mu\in\Lambda_0)とおく。ggはC2C^2級であるからhhはC1C^1級であり、系 5.3 (3)によりh(μ0)=0h(\mu_0)=0である。§E4.4 定理 2.1と系 5.3 (6)によりh′(μ0)=∂s∂μg(0,μ0)=θxμ≠0h'(\mu_0)=\partial_s\partial_\mu g(0,\mu_0)=\theta_{x\mu}\ne0である。したがってδ0>0\delta_0>0が存在して、0<∣μ−μ0∣<δ00<\lvert\mu-\mu_0\rvert<\delta_0ならばμ∈Λ1\mu\in\Lambda_1かつsgn⁡h(μ)=sgn⁡(θxμ(μ−μ0))\operatorname{sgn}h(\mu)=\operatorname{sgn}(\theta_{x\mu}(\mu-\mu_0))である。このようなμ\muについて、系 5.3 (6)によりx(0,μ)=x0+W(0,μ)=x0x(0,\mu)=x_0+W(0,\mu)=x_0、g(0,μ)=0g(0,\mu)=0であるから、主張 (1)(s=0s=0)により、x0x_0はθxμ(μ−μ0)<0\theta_{x\mu}(\mu-\mu_0)<0ならば漸近安定、θxμ(μ−μ0)>0\theta_{x\mu}(\mu-\mu_0)>0ならば不安定である。

主張 (4)を示す。定理 8.1の仮定は主張 (3)の仮定を含むから、同項のδ0\delta_0について、0<∣μ−μ0∣<δ00<\lvert\mu-\mu_0\rvert<\delta_0ならばx0x_0についての主張が成り立つ。他方の平衡点について、GGを補題 7.2(k=2k=2)の関数とする。定理 8.1 (1)により、μ^\hat\muはC1C^1級であり、∣s∣<ε\lvert s\rvert<\varepsilonについてg(s,μ^(s))=0g(s,\hat\mu(s))=0であるから、補題 7.2 (4)により∂sg(s,μ^(s))=−s ∂μG(s,μ^(s)) μ^′(s)\partial_sg(s,\hat\mu(s))=-s\,\partial_\mu G(s,\hat\mu(s))\,\hat\mu'(s)である。補題 7.2 (1)により∂μG(0,μ0)=θxμ\partial_\mu G(0,\mu_0)=\theta_{x\mu}であり、定理 8.1 (1)によりμ^′(0)=−θxx/(2θxμ)≠0\hat\mu'(0)=-\theta_{xx}/(2\theta_{x\mu})\ne0である。ε′∈(0,min⁡{ε,ε1}]\varepsilon'\in(0,\min\{\varepsilon,\varepsilon_1\}]を、∣s∣<ε′\lvert s\rvert<\varepsilon'についてμ^(s)∈Λ1\hat\mu(s)\in\Lambda_1、sgn⁡∂μG(s,μ^(s))=sgn⁡θxμ\operatorname{sgn}\partial_\mu G(s,\hat\mu(s))=\operatorname{sgn}\theta_{x\mu}、sgn⁡μ^′(s)=sgn⁡μ^′(0)\operatorname{sgn}\hat\mu'(s)=\operatorname{sgn}\hat\mu'(0)となるように取る。このときμ^\hat\muは(−ε′,ε′)(-\varepsilon',\varepsilon')の上で狭義単調であり、sgn⁡(μ^(s)−μ0)=sgn⁡(s)sgn⁡μ^′(0)\operatorname{sgn}(\hat\mu(s)-\mu_0)=\operatorname{sgn}(s)\operatorname{sgn}\hat\mu'(0)であるから、0<∣s∣<ε′0<\lvert s\rvert<\varepsilon'について

sgn⁡∂sg(s,μ^(s))=−sgn⁡(s)sgn⁡(θxμ)sgn⁡μ^′(0)=−sgn⁡(θxμ(μ^(s)−μ0))\operatorname{sgn}\partial_sg(s,\hat\mu(s))=-\operatorname{sgn}(s)\operatorname{sgn}(\theta_{x\mu})\operatorname{sgn}\hat\mu'(0)=-\operatorname{sgn}\bigl(\theta_{x\mu}(\hat\mu(s)-\mu_0)\bigr)

である。μ^((−ε′,ε′))\hat\mu((-\varepsilon',\varepsilon'))はμ0\mu_0を含む開区間であるから、δ′∈(0,min⁡{δ,δ0}]\delta'\in(0,\min\{\delta,\delta_0\}]を、(μ0−δ′,μ0+δ′)⊆μ^((−ε′,ε′))(\mu_0-\delta',\mu_0+\delta')\subseteq\hat\mu((-\varepsilon',\varepsilon'))となるように取る。0<∣μ−μ0∣<δ′0<\lvert\mu-\mu_0\rvert<\delta'とすると、μ^(s′)=μ\hat\mu(s')=\muを満たすs′∈(−ε′,ε′)s'\in(-\varepsilon',\varepsilon')が存在し、μ≠μ0=μ^(0)\mu\ne\mu_0=\hat\mu(0)によりs′≠0s'\ne0である。定理 8.1 (1)と定理 8.1 (2)によりx^(s′)\hat x(s')はVVの中にあるx0x_0と異なるF(⋅,μ)F(\cdot,\mu)の平衡点であるから、定理 8.1 (3)により、他方の平衡点はx^(s′)\hat x(s')である。主張 (1)により、x^(s′)\hat x(s')はθxμ(μ−μ0)>0\theta_{x\mu}(\mu-\mu_0)>0ならば漸近安定、θxμ(μ−μ0)<0\theta_{x\mu}(\mu-\mu_0)<0ならば不安定である。

主張 (5)を示す。定理 9.1の仮定は主張 (3)の仮定を含むから、同項のδ0\delta_0について、0<∣μ−μ0∣<δ00<\lvert\mu-\mu_0\rvert<\delta_0ならばx0x_0についての主張が成り立つ。

x^(s)\hat x(s)について、GGを補題 7.2(k=3k=3)の関数とする。定理 9.1 (1)により、μ^\hat\muはC2C^2級であり、∣s∣<ε\lvert s\rvert<\varepsilonについてg(s,μ^(s))=0g(s,\hat\mu(s))=0であり、μ^′(0)=0\hat\mu'(0)=0、μ^′′(0)=−θxxx/(3θxμ)=2c2\hat\mu''(0)=-\theta_{xxx}/(3\theta_{x\mu})=2c_2であるから、τ=sgn⁡μ^′′(0)=−sgn⁡(θxxx)sgn⁡(θxμ)\tau=\operatorname{sgn}\hat\mu''(0)=-\operatorname{sgn}(\theta_{xxx})\operatorname{sgn}(\theta_{x\mu})である。補題 7.2 (4)により∂sg(s,μ^(s))=−s ∂μG(s,μ^(s)) μ^′(s)\partial_sg(s,\hat\mu(s))=-s\,\partial_\mu G(s,\hat\mu(s))\,\hat\mu'(s)であり、補題 7.2 (1)により∂μG(0,μ0)=θxμ\partial_\mu G(0,\mu_0)=\theta_{x\mu}である。補題 6.1をμ^\hat\mu(同補題のε1\varepsilon_1を同定理のε\varepsilonとし、σ=τ\sigma=\tau)に適用し、同補題のε\varepsilonをε2\varepsilon_2とする。ε′∈(0,min⁡{ε2,ε1}]\varepsilon'\in(0,\min\{\varepsilon_2,\varepsilon_1\}]を、∣s∣<ε′\lvert s\rvert<\varepsilon'についてμ^(s)∈Λ1\hat\mu(s)\in\Lambda_1かつsgn⁡∂μG(s,μ^(s))=sgn⁡θxμ\operatorname{sgn}\partial_\mu G(s,\hat\mu(s))=\operatorname{sgn}\theta_{x\mu}となるように取る。補題 6.1 (1)により、0<∣s∣<ε′0<\lvert s\rvert<\varepsilon'についてτsμ^′(s)>0\tau s\hat\mu'(s)>0であるから、

sgn⁡∂sg(s,μ^(s))=−sgn⁡(s)sgn⁡(θxμ) τsgn⁡(s)=−τsgn⁡(θxμ)=sgn⁡(θxxx)\operatorname{sgn}\partial_sg(s,\hat\mu(s))=-\operatorname{sgn}(s)\operatorname{sgn}(\theta_{x\mu})\,\tau\operatorname{sgn}(s)=-\tau\operatorname{sgn}(\theta_{x\mu})=\operatorname{sgn}(\theta_{xxx})

であり、主張 (1)により、x^(s)\hat x(s)はθxxx<0\theta_{xxx}<0ならば漸近安定、θxxx>0\theta_{xxx}>0ならば不安定である。

ε′≤ε2\varepsilon'\le\varepsilon_2であるから、δ′′>0\delta''>0を補題 6.1 (2)(このε′\varepsilon')のδ\deltaとし、δ′=min⁡{δ,δ0,δ′′}\delta'=\min\{\delta,\delta_0,\delta''\}とおく。0<τ(μ−μ0)<δ′0<\tau(\mu-\mu_0)<\delta'とする。補題 6.1 (2)により、μ^(s)=μ\hat\mu(s)=\muを満たすs∈(−ε′,ε′)s\in(-\varepsilon',\varepsilon')は負と正の二つであり、系 5.3 (1)によりwT(x^(s)−x0)=sw^{\mathsf T}(\hat x(s)-x_0)=sであるから、二つのx^(s)\hat x(s)は互いに異なり、x0x_0とも異なる。定理 9.1 (2)によりこれらはVVの中にあるF(⋅,μ)F(\cdot,\mu)の平衡点であり、定理 9.1 (3)により、VVの中にある平衡点はx0x_0とこれら二つだけである。sgn⁡(θxμ(μ−μ0))=τsgn⁡(θxμ)=−sgn⁡θxxx\operatorname{sgn}(\theta_{x\mu}(\mu-\mu_0))=\tau\operatorname{sgn}(\theta_{x\mu})=-\operatorname{sgn}\theta_{xxx}であるから、θxxx>0\theta_{xxx}>0ならばx0x_0は漸近安定、二つのx^(s)\hat x(s)は不安定であり、θxxx<0\theta_{xxx}<0ならばx0x_0は不安定、二つのx^(s)\hat x(s)は漸近安定である。▨

注意 10.3.系 10.2 (3)の仮定の下で、A(μ)=DxF(x0,μ)A(\mu)=D_xF(x_0,\mu)(μ∈Λ\mu\in\Lambda)とおく。FFはC2C^2級であるからA∈C1(Λ,Rn×n)A\in C^1(\Lambda,\R^{n\times n})であり、A(μ0)=A0A(\mu_0)=A_0の固有値00は代数的に単純で、A0v=0A_0v=0、wTA0=0w^{\mathsf T}A_0=0、wTv=1w^{\mathsf T}v=1である。補題 4.2(k=1k=1、λ0=0\lambda_0=0、v0=vv_0=v、ww)のλ\lambdaについて、補題 4.2 (4)によりλ(μ)\lambda(\mu)は実数であり、補題 4.2 (2)とA′(μ0)v=DμDxF(x0,μ0)[v]A'(\mu_0)v=D_\mu D_xF(x_0,\mu_0)[v]によりλ′(μ0)=wTA′(μ0)v=θxμ\lambda'(\mu_0)=w^{\mathsf T}A'(\mu_0)v=\theta_{x\mu}である。したがって、横断条件θxμ≠0\theta_{x\mu}\ne0は、自明な平衡点x0x_0における線形化の代数的に単純な実固有値λ(μ)\lambda(\mu)のμ0\mu_0での微分係数が零でないという条件である。代数的に単純な非実固有値λ(μ)\lambda(\mu)がμ=μ0\mu=\mu_0で虚軸の上にある場合に、これに対応する条件はℜλ′(μ0)≠0\Re\lambda'(\mu_0)\ne0である(問題 12.2)。

11 例

例 11.1.n=1n=1、U=Λ=RU=\Lambda=\R、(x0,μ0)=(0,0)(x_0,\mu_0)=(0,0)とする。このとき、以下のいずれでもDxF(0,0)=0D_xF(0,0)=0であり、00は代数的に単純な固有値で、零固有対は(v,w)=(1,1)(v,w)=(1,1)である。

  1. F(x,μ)=μ−x2F(x,\mu)=\mu-x^2のとき、θμ=1\theta_\mu=1、θxx=−2\theta_{xx}=-2であり、定理 6.2の仮定が成り立つ。平衡点の集合はμ=x2\mu=x^2であり、x^(s)=s\hat x(s)=s、μ^(s)=s2\hat\mu(s)=s^2、μ^′′(0)=2=−θxx/θμ\hat\mu''(0)=2=-\theta_{xx}/\theta_\muである。μ>0\mu>0に平衡点が二個(±μ\pm\sqrt\mu)、μ<0\mu<0に零個ある。定理 6.2 (4)により(0,0)(0,0)は分岐点である。DxF(x,μ)=−2xD_xF(x,\mu)=-2xであるから、μ\sqrt\muは漸近安定、−μ-\sqrt\muは不安定であり、系 10.2 (2)の符号sθxx<0s\theta_{xx}<0と一致する。
  2. F(x,μ)=μx−x2F(x,\mu)=\mu x-x^2のとき、F(0,μ)=0F(0,\mu)=0、θxμ=1\theta_{x\mu}=1、θxx=−2\theta_{xx}=-2であり、定理 8.1の仮定が成り立つ。平衡点はx=0x=0とx=μx=\muであり、x^(s)=s\hat x(s)=s、μ^(s)=s\hat\mu(s)=s、μ^′(0)=1=−θxx/(2θxμ)\hat\mu'(0)=1=-\theta_{xx}/(2\theta_{x\mu})である。定理 8.1 (4)により(0,0)(0,0)は分岐点である。DxF(x,μ)=μ−2xD_xF(x,\mu)=\mu-2xであるから、x=0x=0での値はμ\mu、x=μx=\muでの値は−μ-\muであり、μ\muが00を通過するとき二つの平衡点の安定性が入れ替わる。
  3. F(x,μ)=μx−x3F(x,\mu)=\mu x-x^3のとき、F(0,μ)=0F(0,\mu)=0、Dx2F(0,0)=0D_x^2F(0,0)=0であるからξ=0\xi=0、θxμ=1\theta_{x\mu}=1、θxxx=−6\theta_{xxx}=-6、c2=1c_2=1であり、定理 9.1の仮定が成り立つ。平衡点はx=0x=0と、μ>0\mu>0のときx=±μx=\pm\sqrt\muである。定理 9.1 (4)により(0,0)(0,0)は分岐点である。DxF(x,μ)=μ−3x2D_xF(x,\mu)=\mu-3x^2であるから、±μ\pm\sqrt\muでの値は−2μ<0-2\mu<0であり漸近安定、x=0x=0はμ>0\mu>0で不安定、μ<0\mu<0で漸近安定である。F(x,μ)=μx+x3F(x,\mu)=\mu x+x^3のときはθxxx=6\theta_{xxx}=6、c2=−1c_2=-1であり、x=±−μx=\pm\sqrt{-\mu}(μ<0\mu<0)でのDxFD_xFの値は−2μ>0-2\mu>0であるから不安定、x=0x=0はμ<0\mu<0で漸近安定、μ>0\mu>0で不安定である。系 10.2 (5)の符号と一致する。

例 11.2.n=1n=1、U=Λ=RU=\Lambda=\R、(x0,μ0)=(0,0)(x_0,\mu_0)=(0,0)とする。このとき、以下のいずれでもF(0,μ)=0F(0,\mu)=0(μ∈R\mu\in\R)、DxF(0,0)=0D_xF(0,0)=0であり、00は代数的に単純な固有値で、零固有対は(v,w)=(1,1)(v,w)=(1,1)である。

  1. F(x,μ)=μx−x3F(x,\mu)=\mu x-x^3のとき、θxμ=1≠0\theta_{x\mu}=1\ne0であるが、θxx=Dx2F(0,0)=0\theta_{xx}=D_x^2F(0,0)=0である。μ>0\mu>0の平衡点は00、±μ\pm\sqrt\muの三個である。00の任意の開近傍VVと任意のδ>0\delta>0について、0<μ<δ0<\mu<\deltaかつ±μ∈V\pm\sqrt\mu\in Vを満たすμ\muが存在するから、定理 8.1 (3)の結論(μ≠μ0\mu\ne\mu_0で平衡点はちょうど二個)はθxx≠0\theta_{xx}\ne0を外すと成り立たない。
  2. F(x,μ)=−x2F(x,\mu)=-x^2のとき、θxx=−2≠0\theta_{xx}=-2\ne0であるが、θxμ=0\theta_{x\mu}=0である。すべてのμ\muについて平衡点は00だけであるから、定理 8.1 (3)の結論はθxμ≠0\theta_{x\mu}\ne0を外すと成り立たない。
  3. F(x,μ)=μ2x−x2F(x,\mu)=\mu^2x-x^2のとき、θxx=−2≠0\theta_{xx}=-2\ne0であるが、θxμ=∂μ(μ2)∣μ=0=0\theta_{x\mu}=\partial_\mu(\mu^2)\big|_{\mu=0}=0である。μ≠0\mu\ne0の平衡点は00とμ2\mu^2の二個であり、平衡点の個数については定理 8.1 (3)と同じ結論が成り立つ。しかしDxF(x,μ)=μ2−2xD_xF(x,\mu)=\mu^2-2xであるから、DxF(0,μ)=μ2>0D_xF(0,\mu)=\mu^2>0、DxF(μ2,μ)=−μ2<0D_xF(\mu^2,\mu)=-\mu^2<0であり、μ\muの符号によらず00は不安定、μ2\mu^2は漸近安定である(F(⋅,μ)F(\cdot,\mu)はR\R上のC1C^1級の写像であり、§E13.4 補題 4.2により局所 Lipschitz 連続であるから、§E13.16 定理 4.5 (2)と§E13.16 定理 4.5 (1)をW=RW=\Rとして適用する)。したがって、系 10.2 (4)の安定性の交代はθxμ≠0\theta_{x\mu}\ne0を外すと成り立たない。

例 11.3.n=1n=1、U=Λ=RU=\Lambda=\R、F(x,μ)=μ−x3F(x,\mu)=\mu-x^3、(x0,μ0)=(0,0)(x_0,\mu_0)=(0,0)とする。このとき、DxF(0,0)=0D_xF(0,0)=0であり、零固有対は(1,1)(1,1)で、θμ=1≠0\theta_\mu=1\ne0であるが、θxx=Dx2F(0,0)=0\theta_{xx}=D_x^2F(0,0)=0である。平衡点はx=μ1/3x=\mu^{1/3}であり、すべてのμ∈R\mu\in\Rについてちょうど一個である。したがって、定理 6.2 (3)の結論(一方の側に二個、他方の側に零個)はθxx≠0\theta_{xx}\ne0を外すと成り立たない。また、μ↦μ1/3\mu\mapsto\mu^{1/3}はμ=0\mu=0で微分可能でなく、定理 3.3 (1)のC1C^1級の枝は、DxF(0,0)D_xF(0,0)の可逆性を外すと存在しない。一方、ξ(μ)=μ1/3\xi(\mu)=\mu^{1/3}は連続であり、(y,μ)∈R×R(y,\mu)\in\R\times\RについてF(y,μ)=0F(y,\mu)=0であることとy=ξ(μ)y=\xi(\mu)であることは同値であるから、(0,0)(0,0)は分岐点でない。したがって、定理 3.3 (3)の逆は成り立たず、DxF(x0,μ0)D_xF(x_0,\mu_0)が可逆でなくても(x0,μ0)(x_0,\mu_0)が分岐点であるとは限らない。μ≠0\mu\ne0ではDxF(μ1/3,μ)=−3μ2/3<0D_xF(\mu^{1/3},\mu)=-3\mu^{2/3}<0であるから平衡点は漸近安定であり(§E13.4 補題 4.2によりF(⋅,μ)F(\cdot,\mu)は局所 Lipschitz 連続であるから、§E13.16 定理 4.5 (1)をW=RW=\Rとして適用する)、μ=0\mu=0でも、x′=−x3x'=-x^3の原点は漸近安定である(§E13.17 例 7.5 (1))。

例 11.4.n=1n=1、U=Λ=RU=\Lambda=\R、ε>0\varepsilon>0とし、F(x,μ)=μx−x3+εF(x,\mu)=\mu x-x^3+\varepsilonとする。c=(ε/2)1/3c=(\varepsilon/2)^{1/3}とおく。このとき、x=0x=0は平衡点でなく、F(x,μ)=0F(x,\mu)=0はx≠0x\ne0かつμ=ρ(x)=x2−ε/x\mu=\rho(x)=x^2-\varepsilon/xと同値である。ρ′(x)=(2x3+ε)/x2\rho'(x)=(2x^3+\varepsilon)/x^2である。

x>0x>0ではρ′>0\rho'>0であり、ρ\rhoは(0,∞)(0,\infty)からR\Rへの全単射である。したがって各μ∈R\mu\in\Rについてx>0x>0の平衡点がちょうど一個ある。DxF(x,ρ(x))=ρ(x)−3x2=−2x2−ε/x<0D_xF(x,\rho(x))=\rho(x)-3x^2=-2x^2-\varepsilon/x<0であるから、これらは漸近安定である。

x<0x<0では、ρ′(x)<0\rho'(x)<0(x<−cx<-c)、ρ′(x)>0\rho'(x)>0(−c<x<0-c<x<0)であり、ρ(x)→+∞\rho(x)\to+\infty(x→−∞x\to-\inftyまたはx→0−x\to0^-)、ρ(−c)=c2+ε/c=3c2\rho(-c)=c^2+\varepsilon/c=3c^2である(ε=2c3\varepsilon=2c^3)。したがってx<0x<0の平衡点は、μ<3c2\mu<3c^2で零個、μ=3c2\mu=3c^2で一個、μ>3c2\mu>3c^2で二個である。(x1,μ1)=(−c,3c2)(x_1,\mu_1)=(-c,3c^2)ではDxF(x1,μ1)=μ1−3x12=0D_xF(x_1,\mu_1)=\mu_1-3x_1^2=0であり、θμ=DμF=x1=−c≠0\theta_\mu=D_\mu F=x_1=-c\ne0、θxx=Dx2F=−6x1=6c≠0\theta_{xx}=D_x^2F=-6x_1=6c\ne0であるから、定理 6.2の仮定が成り立ち、μ^′′(0)=−θxx/θμ=6>0\hat\mu''(0)=-\theta_{xx}/\theta_\mu=6>0である。ε=0\varepsilon=0のときFFは例 11.1 (3)のμx−x3\mu x-x^3に一致し、F(−x,μ)=−F(x,μ)F(-x,\mu)=-F(x,\mu)が成り立つ。ε>0\varepsilon>0のときこの等式は成り立たず、平衡点の集合は、すべてのμ\muについて平衡点を一個もつx>0x>0の部分と、(−c,3c2)(-c,3c^2)でサドルノード分岐の仮定を満たすx<0x<0の部分からなる。全体の平衡点の個数は、μ<3c2\mu<3c^2で一個、μ=3c2\mu=3c^2で二個、μ>3c2\mu>3c^2で三個である。

12 演習

問題 12.1.a∈Ra\in\Rとし、U=R2U=\R^2、Λ=R\Lambda=\R、F(x,y,μ)=(μx−x3+axy, −y+x2)F(x,y,\mu)=(\mu x-x^3+axy,\ -y+x^2)とする。(x0,μ0)=(0,0)(x_0,\mu_0)=(0,0)で定理 9.1の仮定を確かめ、a≠1a\ne1のとき、分岐する平衡点の枝のμ−μ0\mu-\mu_0のs2s^2の係数c2c_2をaaで表し、枝上の平衡点が漸近安定になるためのaaの条件を求めよ。

解答.

F(0,0,μ)=0F(0,0,\mu)=0である。DxF(x,y,μ)=(μ−3x2+ayax2x−1)D_xF(x,y,\mu)=\begin{pmatrix}\mu-3x^2+ay&ax\\2x&-1\end{pmatrix}であるから、A0=(000−1)A_0=\begin{pmatrix}0&0\\0&-1\end{pmatrix}であり、固有値は00(代数的に単純)と−1-1である。v=w=(1,0)Tv=w=(1,0)^{\mathsf T}は零固有対である。

二階微分はDx2F(0)[v,v]=(∂x2F1,∂x2F2)=(−6x,2)∣0=(0,2)D_x^2F(0)[v,v]=(\partial_x^2F_1,\partial_x^2F_2)=(-6x,2)|_0=(0,2)であり、wTDx2F(0)[v,v]=0w^{\mathsf T}D_x^2F(0)[v,v]=0である。ξ=(0,η)∈ker⁡wT\xi=(0,\eta)\in\ker w^{\mathsf T}についてA0ξ=(0,−η)A_0\xi=(0,-\eta)であり、A0ξ=−(0,2)A_0\xi=-(0,2)からξ=(0,2)\xi=(0,2)を得る。DμDxF(0)[v]=∂μ(μ−3x2+ay, 2x)∣0=(1,0)D_\mu D_xF(0)[v]=\partial_\mu(\mu-3x^2+ay,\,2x)|_0=(1,0)であるからθxμ=1≠0\theta_{x\mu}=1\ne0である。三階微分はDx3F(0)[v,v,v]=(−6,0)D_x^3F(0)[v,v,v]=(-6,0)であり、Dx2F(0)[v,ξ]D_x^2F(0)[v,\xi]の第一成分は∂x∂yF1 vxξy=a⋅1⋅2=2a\partial_x\partial_yF_1\,v_x\xi_y=a\cdot1\cdot2=2a、第二成分は∂x∂yF2=0\partial_x\partial_yF_2=0により00である。したがってθxxx=−6+3⋅2a=6(a−1)\theta_{xxx}=-6+3\cdot2a=6(a-1)であり、a≠1a\ne1ならばθxxx≠0\theta_{xxx}\ne0である。定理 9.1の仮定が成り立ち、c2=−θxxx/(6θxμ)=1−ac_2=-\theta_{xxx}/(6\theta_{x\mu})=1-aである。

A0A_0の零でない固有値−1-1の実部は負であるから、系 10.2 (5)により、分岐する枝上の平衡点はθxxx<0\theta_{xxx}<0、すなわちa<1a<1のとき漸近安定であり、a>1a>1のとき不安定である。a<1a<1のとき枝はμ>μ0\mu>\mu_0の側に存在する(c2>0c_2>0)。a>1a>1のとき枝はμ<μ0\mu<\mu_0の側に存在する(c2<0c_2<0)。▨

問題 12.2.A(μ)=(μ−11μ)A(\mu)=\begin{pmatrix}\mu&-1\\1&\mu\end{pmatrix}(μ∈R\mu\in\R)とする。μ0=0\mu_0=0におけるA(0)A(0)の固有値λ0=i\lambda_0=iについて、補題 4.2のλ′(0)\lambda'(0)を同補題の公式で求め、A(μ)A(\mu)の固有値を直接求めて結果を確かめよ。また、ℜλ′(0)>0\Re\lambda'(0)>0であることから、∣μ∣\lvert\mu\rvertが小さいときℜλ(μ)\Re\lambda(\mu)の符号がμ\muの符号に等しいことを証明せよ。

解答.

A(0)=(0−110)A(0)=\begin{pmatrix}0&-1\\1&0\end{pmatrix}の特性多項式はz2+1z^2+1であり、λ0=i\lambda_0=iは単根であるから代数的に単純である。v0=(1,−i)Tv_0=(1,-i)^{\mathsf T}はA(0)v0=(i,1)T=iv0A(0)v_0=(i,1)^{\mathsf T}=iv_0を満たす。行ベクトルw∗=(1,i)w^*=(1,i)はw∗A(0)=(i,−1)=i w∗w^*A(0)=(i,-1)=i\,w^*を満たし、w=(1,−i)Tw=(1,-i)^{\mathsf T}である。w∗v0=1+1=2≠0w^*v_0=1+1=2\ne0である。A′(μ)=I2A'(\mu)=I_2であるから、λ′(0)=w∗A′(0)v0/(w∗v0)=w∗v0/(w∗v0)=1\lambda'(0)=w^*A'(0)v_0/(w^*v_0)=w^*v_0/(w^*v_0)=1である。

A(μ)A(\mu)の特性多項式は(z−μ)2+1(z-\mu)^2+1であり、固有値はμ±i\mu\pm iである。λ0=i\lambda_0=iから連続な枝はλ(μ)=μ+i\lambda(\mu)=\mu+iであり、λ′(0)=1\lambda'(0)=1は公式の値と一致する。ℜλ′(0)=1>0\Re\lambda'(0)=1>0である。λ\lambdaはC1C^1級でℜλ(0)=0\Re\lambda(0)=0であるから、ℜλ(μ)=ℜλ′(0) μ+o(μ)\Re\lambda(\mu)=\Re\lambda'(0)\,\mu+o(\mu)であり、ℜλ′(0)>0\Re\lambda'(0)>0ならば、∣μ∣\lvert\mu\rvertが小さいときℜλ(μ)\Re\lambda(\mu)の符号はμ\muの符号に等しい。この例ではℜλ(μ)=μ\Re\lambda(\mu)=\muである。▨

問題 12.3.U=Λ=RU=\Lambda=\R、F(x,μ)=μ2−x2F(x,\mu)=\mu^2-x^2、(x0,μ0)=(0,0)(x_0,\mu_0)=(0,0)とする。(0,0)(0,0)で定理 6.2の仮定のうち成り立たないものを挙げよ。次に、F~(y,μ)=F(y+μ,μ)\tilde F(y,\mu)=F(y+\mu,\mu)とおき、族y′=F~(y,μ)y'=\tilde F(y,\mu)が(0,0)(0,0)で定理 8.1の仮定を満たすことを確かめ、(0,0)(0,0)が族x′=F(x,μ)x'=F(x,\mu)の分岐点であることを、族y′=F~(y,μ)y'=\tilde F(y,\mu)についての定理 8.1 (4)から導け。さらに、∣μ∣\lvert\mu\rvertが小さくμ≠0\mu\ne0のとき、F(⋅,μ)F(\cdot,\mu)の平衡点x=μx=\mu、x=−μx=-\muの安定性を系 10.2 (4)から求め、DxFD_xFの値と照合せよ。

解答.

DxF(x,μ)=−2xD_xF(x,\mu)=-2xであるからDxF(0,0)=0D_xF(0,0)=0であり、00は代数的に単純な固有値で、零固有対は(1,1)(1,1)である。θxx=Dx2F(0,0)=−2≠0\theta_{xx}=D_x^2F(0,0)=-2\ne0であるが、θμ=DμF(0,0)=2μ∣μ=0=0\theta_\mu=D_\mu F(0,0)=2\mu|_{\mu=0}=0である。したがって、成り立たない仮定はθμ≠0\theta_\mu\ne0である。

F~(y,μ)=μ2−(y+μ)2=−y2−2μy\tilde F(y,\mu)=\mu^2-(y+\mu)^2=-y^2-2\mu yであるから、すべてのμ\muについてF~(0,μ)=0\tilde F(0,\mu)=0である。DyF~(y,μ)=−2y−2μD_y\tilde F(y,\mu)=-2y-2\muであるからDyF~(0,0)=0D_y\tilde F(0,0)=0であり、零固有対は(1,1)(1,1)である。θxμ=ddμDyF~(0,μ)∣μ=0=−2≠0\theta_{x\mu}=\dfrac{d}{d\mu}D_y\tilde F(0,\mu)\Big|_{\mu=0}=-2\ne0、θxx=Dy2F~(0,0)=−2≠0\theta_{xx}=D_y^2\tilde F(0,0)=-2\ne0であるから、定理 8.1の仮定が成り立ち、定理 8.1 (4)により(0,0)(0,0)は族y′=F~(y,μ)y'=\tilde F(y,\mu)の分岐点である。

(0,0)(0,0)が族x′=F(x,μ)x'=F(x,\mu)の分岐点でないと仮定し、定義 1.1の組(Λ′,B′,ξ)(\Lambda',B',\xi)を取る。F(0,0)=0F(0,0)=0であるからξ(0)=0\xi(0)=0である。r>0r>0を(−2r,2r)⊆B′(-2r,2r)\subseteq B'となるように取り、B′′=(−r,r)B''=(-r,r)とおく。ξ\xiの連続性により、00を含む開区間Λ′′⊆Λ′∩(−r,r)\Lambda''\subseteq\Lambda'\cap(-r,r)を、μ∈Λ′′\mu\in\Lambda''についてξ(μ)−μ∈B′′\xi(\mu)-\mu\in B''となるように取る。ξ~(μ)=ξ(μ)−μ\tilde\xi(\mu)=\xi(\mu)-\muはΛ′′\Lambda''からB′′B''への連続写像である。(y,μ)∈B′′×Λ′′(y,\mu)\in B''\times\Lambda''ならばy+μ∈(−2r,2r)⊆B′y+\mu\in(-2r,2r)\subseteq B'であるから、F~(y,μ)=0\tilde F(y,\mu)=0であることはF(y+μ,μ)=0F(y+\mu,\mu)=0、すなわちy+μ=ξ(μ)y+\mu=\xi(\mu)であることと同値であり、y=ξ~(μ)y=\tilde\xi(\mu)であることと同値である。したがって(Λ′′,B′′,ξ~)(\Lambda'',B'',\tilde\xi)は族y′=F~(y,μ)y'=\tilde F(y,\mu)についての定義 1.1の組であり、(0,0)(0,0)がこの族の分岐点であることと両立しない。よって(0,0)(0,0)は族x′=F(x,μ)x'=F(x,\mu)の分岐点である。

μ\muを固定すると、xxがx′=F(x,μ)x'=F(x,\mu)の解であることと、y=x−μy=x-\muがy′=F~(y,μ)y'=\tilde F(y,\mu)の解であることは同値である。F~(⋅,μ)\tilde F(\cdot,\mu)の流れφ~\tilde\varphiとF(⋅,μ)F(\cdot,\mu)の流れφ\varphiはφt(x)=φ~t(x−μ)+μ\varphi_t(x)=\tilde\varphi_t(x-\mu)+\muを満たし、∣φt(x)−x∗∣=∣φ~t(x−μ)−(x∗−μ)∣\lvert\varphi_t(x)-x_*\rvert=\lvert\tilde\varphi_t(x-\mu)-(x_*-\mu)\rvertであるから、§E13.16 定義 2.2により、x∗x_*がF(⋅,μ)F(\cdot,\mu)の漸近安定な平衡点、または不安定な平衡点であることは、x∗−μx_*-\muがF~(⋅,μ)\tilde F(\cdot,\mu)のそれであることと同値である。F~(⋅,μ)\tilde F(\cdot,\mu)の平衡点はy=0y=0とy=−2μy=-2\muであり、それぞれx=μx=\muとx=−μx=-\muに対応する。n=1n=1であるからA0A_0に00以外の固有値は無く、系 10.2 (4)を族y′=F~(y,μ)y'=\tilde F(y,\mu)に適用すると、θxμμ=−2μ\theta_{x\mu}\mu=-2\muにより、∣μ∣\lvert\mu\rvertが小さくμ≠0\mu\ne0のとき、y=0y=0はμ>0\mu>0で漸近安定、μ<0\mu<0で不安定であり、他方の平衡点y=−2μy=-2\muはμ<0\mu<0で漸近安定、μ>0\mu>0で不安定である。したがって、x=μx=\muはμ>0\mu>0で漸近安定、μ<0\mu<0で不安定であり、x=−μx=-\muはμ<0\mu<0で漸近安定、μ>0\mu>0で不安定である。DxF(μ,μ)=−2μD_xF(\mu,\mu)=-2\mu、DxF(−μ,μ)=2μD_xF(-\mu,\mu)=2\muであるから、§E13.16 定理 4.5 (1)と§E13.16 定理 4.5 (2)(W=RW=\R。F(⋅,μ)F(\cdot,\mu)は§E13.4 補題 4.2により局所 Lipschitz 連続である)から得られる結論と一致する。▨

参考文献

  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.双曲型平衡点の線形化と平衡点の持続の扱いを参考にした。

前提記事

16 本の記事・単元を表示