§E13.23Hopf 分岐

最終更新

自律系x′=F(x,μ)x'=F(x,\mu)の平衡点は、その点での線形化の固有値の実部がすべて負であるとき漸近安定である。一方、パラメータμ\muの変化によって固有値が虚軸を横切るときには、平衡点の安定性や、平衡点の近傍にある平衡点と周期解の個数が変わることがある。実固有値が00を通る場合には、非退化条件の下で、§E13.22 分岐 で扱った平衡点の分岐が起こる。これに対して、複素共役な一対の固有値が虚軸を横切る場合には、平衡点の近傍に周期解が現れることがある。この周期解の発生を Hopf 分岐という。たとえば、平面上の系ζ′=(μ+i)ζ−∣ζ∣2ζ\zeta'=(\mu+i)\zeta-\lvert\zeta\rvert^2\zeta(ζ∈C\zeta\in\C)では、μ\muが00を超えると原点の周りに円周状の周期軌道が現れる。本記事では、Hopf 分岐の定理と、標準形の例について解説する。

1 Banach 空間の陰関数定理

補題 1.1.X,YX,YをR\R上の Banach 空間、(P,d)(P,d)を距離空間とし、p0∈Pp_0\in P、x0∈Xx_0\in X、r>0r>0とする。r′>0r'>0についてB(x0,r′)={x∈X:∥x−x0∥<r′}B(x_0,r')=\{x\in X:\lVert x-x_0\rVert<r'\}、B‾(x0,r′)={x∈X:∥x−x0∥≤r′}\overline B(x_0,r')=\{x\in X:\lVert x-x_0\rVert\le r'\}と書く。写像H:P×B(x0,r)→YH:P\times B(x_0,r)\to YはH(p0,x0)=0H(p_0,x_0)=0を満たし、各x∈B(x0,r)x\in B(x_0,r)についてp↦H(p,x)p\mapsto H(p,x)は連続であるとする。L:X→YL:X\to Yは全単射な有界線形作用素で、逆写像L−1:Y→XL^{-1}:Y\to Xも有界線形作用素であるとする。さらに、任意のη>0\eta>0に対してρ∈(0,r]\rho\in(0,r]が存在して、d(p,p0)<ρd(p,p_0)<\rhoかつx,y∈B(x0,ρ)x,y\in B(x_0,\rho)ならば

∥H(p,x)−H(p,y)−L(x−y)∥≤η∥x−y∥\lVert H(p,x)-H(p,y)-L(x-y)\rVert\le\eta\lVert x-y\rVert

が成り立つとする。このとき、ρ1>0\rho_1>0とρ2∈(0,r)\rho_2\in(0,r)が存在して、次が成り立つ。

  1. d(p,p0)<ρ1d(p,p_0)<\rho_1を満たす各p∈Pp\in Pについて、H(p,x)=0H(p,x)=0を満たすx∈B‾(x0,ρ2)x\in\overline B(x_0,\rho_2)がただ一つ存在する。これをx(p)x(p)と書く。
  2. x(p0)=x0x(p_0)=x_0であり、p↦x(p)p\mapsto x(p)は{p∈P:d(p,p0)<ρ1}\{p\in P:d(p,p_0)<\rho_1\}の上で連続である。

証明.η>0\eta>0を∥L−1∥η≤1/2\lVert L^{-1}\rVert\eta\le1/2を満たすように取り、このη\etaに対する仮定のρ∈(0,r]\rho\in(0,r]を取る。ρ2=ρ/2\rho_2=\rho/2とおく。p↦H(p,x0)p\mapsto H(p,x_0)はp0p_0で連続でH(p0,x0)=0H(p_0,x_0)=0であるから、ρ1∈(0,ρ]\rho_1\in(0,\rho]が存在して、d(p,p0)<ρ1d(p,p_0)<\rho_1ならば∥L−1∥∥H(p,x0)∥≤ρ2/2\lVert L^{-1}\rVert\lVert H(p,x_0)\rVert\le\rho_2/2である。P1={p∈P:d(p,p0)<ρ1}P_1=\{p\in P:d(p,p_0)<\rho_1\}とおき、P1P_1にはddの制限を距離として入れる。

p∈P1p\in P_1について、Tp:B‾(x0,ρ2)→XT_p:\overline B(x_0,\rho_2)\to XをTp(x)=x−L−1H(p,x)T_p(x)=x-L^{-1}H(p,x)で定める。B‾(x0,ρ2)⊆B(x0,ρ)\overline B(x_0,\rho_2)\subseteq B(x_0,\rho)であるから、x,y∈B‾(x0,ρ2)x,y\in\overline B(x_0,\rho_2)について

∥Tp(x)−Tp(y)∥=∥L−1(H(p,x)−H(p,y)−L(x−y))∥≤∥L−1∥η∥x−y∥≤12∥x−y∥\lVert T_p(x)-T_p(y)\rVert=\bigl\lVert L^{-1}\bigl(H(p,x)-H(p,y)-L(x-y)\bigr)\bigr\rVert\le\lVert L^{-1}\rVert\eta\lVert x-y\rVert\le\tfrac12\lVert x-y\rVert

である。Tp(x0)−x0=−L−1H(p,x0)T_p(x_0)-x_0=-L^{-1}H(p,x_0)であるから、x∈B‾(x0,ρ2)x\in\overline B(x_0,\rho_2)について∥Tp(x)−x0∥≤∥Tp(x)−Tp(x0)∥+∥L−1H(p,x0)∥≤ρ2/2+ρ2/2=ρ2\lVert T_p(x)-x_0\rVert\le\lVert T_p(x)-T_p(x_0)\rVert+\lVert L^{-1}H(p,x_0)\rVert\le\rho_2/2+\rho_2/2=\rho_2であり、TpT_pはB‾(x0,ρ2)\overline B(x_0,\rho_2)をそれ自身へ写す。B‾(x0,ρ2)\overline B(x_0,\rho_2)は Banach 空間XXの空でない閉集合であるから、§E2.5 定理 3.2により完備である。

主張 (1)を示す。§E13.20 補題 1.3 (1)(完備距離空間B‾(x0,ρ2)\overline B(x_0,\rho_2)、パラメータの集合P1P_1、θ=12\theta=\frac12)により、各p∈P1p\in P_1についてTpT_pはB‾(x0,ρ2)\overline B(x_0,\rho_2)にただ一つの不動点をもつ。L−1L^{-1}は単射であるから、x∈B‾(x0,ρ2)x\in\overline B(x_0,\rho_2)がTpT_pの不動点であることとH(p,x)=0H(p,x)=0であることは同値である。

主張 (2)を示す。x0∈B‾(x0,ρ2)x_0\in\overline B(x_0,\rho_2)はH(p0,x0)=0H(p_0,x_0)=0を満たすから、一意性によりx(p0)=x0x(p_0)=x_0である。q∈P1q\in P_1とする。p∈P1p\in P_1についてTp(x(q))=x(q)−L−1H(p,x(q))T_p(x(q))=x(q)-L^{-1}H(p,x(q))であり、p↦H(p,x(q))p\mapsto H(p,x(q))は連続でL−1L^{-1}は有界であるから、p↦Tp(x(q))p\mapsto T_p(x(q))はqqで連続である。§E13.20 補題 1.3 (2)により、p↦x(p)p\mapsto x(p)はqqで連続である。▨

2 周期関数の空間の線形化作用素と最小周期

補題 2.1.n∈N≥1n\in\NN、ω0>0\omega_0>0とし、A0∈Rn×nA_0\in\R^{n\times n}は±iω0\pm i\omega_0を代数的に単純な固有値にもち、すべてのk∈Z∖{1,−1}k\in\Z\setminus\{1,-1\}についてikω0ik\omega_0を固有値にもたないとする。v,w∈Cnv,w\in\C^nはA0v=iω0vA_0v=i\omega_0v、w∗A0=iω0w∗w^*A_0=i\omega_0w^*、w∗v=1w^*v=1を満たすとする。Rn\R^nには Euclid ノルム∣⋅∣\lvert\cdot\rvertを、行列にはその作用素ノルムを用いる。

C2π={u∈C(R;Rn):u(s+2π)=u(s) (s∈R)},∥u∥∞=sup⁡s∈R∣u(s)∣,C_{2\pi}=\{u\in C(\R;\R^n):u(s+2\pi)=u(s)\ (s\in\R)\},\qquad\lVert u\rVert_\infty=\sup_{s\in\R}\lvert u(s)\rvert,C2π1=C2π∩C1(R;Rn),∥u∥C1=∥u∥∞+∥u′∥∞C^1_{2\pi}=C_{2\pi}\cap C^1(\R;\R^n),\qquad\lVert u\rVert_{C^1}=\lVert u\rVert_\infty+\lVert u'\rVert_\infty

とおき、s∈Rs\in\Rについて

ϕ1(s)=ℜ(eisv),ϕ2(s)=ℑ(eisv),ψ1(s)=ℜ(eisw),ψ2(s)=ℑ(eisw)\phi_1(s)=\Re(e^{is}v),\qquad\phi_2(s)=\Im(e^{is}v),\qquad\psi_1(s)=\Re(e^{is}w),\qquad\psi_2(s)=\Im(e^{is}w)

とおく。L0:C2π1→C2πL_0:C^1_{2\pi}\to C_{2\pi}をL0u=ω0u′−A0uL_0u=\omega_0u'-A_0uで定め、X1={u∈C2π1:∫02πϕjTu ds=0 (j=1,2)}X_1=\{u\in C^1_{2\pi}:\int_0^{2\pi}\phi_j^{\mathsf T}u\,ds=0\ (j=1,2)\}とおく。このとき、次が成り立つ。

  1. ϕ1,ϕ2∈C2π1\phi_1,\phi_2\in C^1_{2\pi}、ϕ1′=−ϕ2\phi_1'=-\phi_2、ϕ2′=ϕ1\phi_2'=\phi_1であり、ϕ1,ϕ2\phi_1,\phi_2はker⁡L0\ker L_0の基底である。
  2. ψ1,ψ2\psi_1,\psi_2は、ω0ψ′=−A0Tψ\omega_0\psi'=-A_0^{\mathsf T}\psiを満たすψ∈C2π1\psi\in C^1_{2\pi}の全体のなすR\R線形空間の基底である。
  3. すべてのB∈Rn×nB\in\R^{n\times n}について、β=w∗Bv\beta=w^*Bvとおくと ∫02πψ1TBϕ1 ds=∫02πψ2TBϕ2 ds=πℜβ,∫02πψ1TBϕ2 ds=−∫02πψ2TBϕ1 ds=πℑβ\int_0^{2\pi}\psi_1^{\mathsf T}B\phi_1\,ds=\int_0^{2\pi}\psi_2^{\mathsf T}B\phi_2\,ds=\pi\Re\beta,\qquad \int_0^{2\pi}\psi_1^{\mathsf T}B\phi_2\,ds=-\int_0^{2\pi}\psi_2^{\mathsf T}B\phi_1\,ds=\pi\Im\beta である。特に、i,j∈{1,2}i,j\in\{1,2\}について∫02πψiTϕj ds=πδij\int_0^{2\pi}\psi_i^{\mathsf T}\phi_j\,ds=\pi\delta_{ij}である。
  4. im⁡L0={h∈C2π:∫02πψjTh ds=0 (j=1,2)}\im L_0=\{h\in C_{2\pi}:\int_0^{2\pi}\psi_j^{\mathsf T}h\,ds=0\ (j=1,2)\}である。X1X_1はC2π1C^1_{2\pi}の閉部分空間で、C2π1=ker⁡L0⊕X1C^1_{2\pi}=\ker L_0\oplus X_1であり、L0∣X1:X1→im⁡L0L_0|_{X_1}:X_1\to\im L_0は全単射である。さらに、あるC>0C>0が存在して、すべてのu∈X1u\in X_1で∥u∥C1≤C∥L0u∥∞\lVert u\rVert_{C^1}\le C\lVert L_0u\rVert_\inftyである。
  5. X=X1×R×RX=X_1\times\R\times\Rにノルム∥(u,a,b)∥X=max⁡{∥u∥C1,∣a∣,∣b∣}\lVert(u,a,b)\rVert_X=\max\{\lVert u\rVert_{C^1},\lvert a\rvert,\lvert b\rvert\}を入れると、XXは Banach 空間である。B∈Rn×nB\in\R^{n\times n}がℜ(w∗Bv)≠0\Re(w^*Bv)\ne0を満たすとし、LB:X→C2πL_B:X\to C_{2\pi}をLB(u,a,b)=L0u+aϕ1′−bBϕ1L_B(u,a,b)=L_0u+a\phi_1'-bB\phi_1で定める。このとき、LBL_Bは全単射な有界線形作用素であり、逆写像LB−1L_B^{-1}も有界線形作用素である。

証明.A=A0/ω0A=A_0/\omega_0とおき、LFu=u′−AuL_Fu=u'-Auとおく。L0=ω0LFL_0=\omega_0L_Fであるから、ker⁡L0=ker⁡LF\ker L_0=\ker L_F、im⁡L0=im⁡LF\im L_0=\im L_Fである。AAの固有値はA0A_0の固有値をω0\omega_0で割ったものであり、代数的重複度も等しいから、±i\pm iはAAの特性多項式の単根であり、AAの固有値でiZi\Zに属するものは±i\pm iの他に無い。したがって§E13.13 命題 4.7 (4)(K=R\mathbb K=\R、T=2πT=2\pi、ω=1\omega=1)により、u′=Auu'=Auの2π2\pi周期の実数値の局所絶対連続解の全体と、ψ′=−ATψ\psi'=-A^{\mathsf T}\psiの2π2\pi周期の実数値の局所絶対連続解の全体は、ともに二次元のR\R線形空間である。

主張 (3)を示す。ϕ=ϕ1+iϕ2\phi=\phi_1+i\phi_2、ψ=ψ1+iψ2\psi=\psi_1+i\psi_2とおくとϕ(s)=eisv\phi(s)=e^{is}v、ψ(s)=eisw\psi(s)=e^{is}wである。B∈Rn×nB\in\R^{n\times n}についてψ∗Bϕ=w∗Bv=β\psi^*B\phi=w^*Bv=\beta、ψTBϕ=e2iswTBv\psi^{\mathsf T}B\phi=e^{2is}w^{\mathsf T}Bvであり、∫02πe2is ds=0\int_0^{2\pi}e^{2is}\,ds=0であるから

∫02πψ∗Bϕ ds=2πβ,∫02πψTBϕ ds=0\int_0^{2\pi}\psi^*B\phi\,ds=2\pi\beta,\qquad\int_0^{2\pi}\psi^{\mathsf T}B\phi\,ds=0

である。aij=∫02πψiTBϕj dsa_{ij}=\int_0^{2\pi}\psi_i^{\mathsf T}B\phi_j\,dsとおくと、ψ∗=ψ1T−iψ2T\psi^*=\psi_1^{\mathsf T}-i\psi_2^{\mathsf T}、ψT=ψ1T+iψ2T\psi^{\mathsf T}=\psi_1^{\mathsf T}+i\psi_2^{\mathsf T}であるから、二つの積分はそれぞれa11+a22+i(a12−a21)a_{11}+a_{22}+i(a_{12}-a_{21})とa11−a22+i(a12+a21)a_{11}-a_{22}+i(a_{12}+a_{21})に等しい。実部と虚部を比べてa11=a22=πℜβa_{11}=a_{22}=\pi\Re\beta、a12=−a21=πℑβa_{12}=-a_{21}=\pi\Im\betaを得る。B=InB=I_nのときβ=w∗v=1\beta=w^*v=1である。

主張 (1)を示す。ω0(eisv)′=iω0eisv=A0eisv\omega_0(e^{is}v)'=i\omega_0e^{is}v=A_0e^{is}vであり、A0A_0とω0\omega_0は実であるから、実部と虚部を取ってω0ϕj′=A0ϕj\omega_0\phi_j'=A_0\phi_j(j=1,2j=1,2)を得る。ϕj\phi_jはC1C^1級で2π2\pi周期であるからϕj∈ker⁡L0\phi_j\in\ker L_0である。ϕ1′+iϕ2′=ieisv=−ϕ2+iϕ1\phi_1'+i\phi_2'=ie^{is}v=-\phi_2+i\phi_1である。c1,c2∈Rc_1,c_2\in\Rがc1ϕ1+c2ϕ2=0c_1\phi_1+c_2\phi_2=0を満たすならば、主張 (3)により0=∫02πψiT(c1ϕ1+c2ϕ2) ds=πci0=\int_0^{2\pi}\psi_i^{\mathsf T}(c_1\phi_1+c_2\phi_2)\,ds=\pi c_i(i=1,2i=1,2)であるから、ϕ1,ϕ2\phi_1,\phi_2は一次独立である。ker⁡L0\ker L_0の元はu′=Auu'=Auの2π2\pi周期の実数値の局所絶対連続解であるから、第一段によりdim⁡Rker⁡L0≤2\dim_\R\ker L_0\le2であり、ϕ1,ϕ2\phi_1,\phi_2はker⁡L0\ker L_0の基底である。

主張 (2)を示す。w∗A0=iω0w∗w^*A_0=i\omega_0w^*の共役転置を取ると、A0A_0は実行列であるからA0Tw=−iω0wA_0^{\mathsf T}w=-i\omega_0wである。したがってω0(eisw)′=iω0eisw=−A0Teisw\omega_0(e^{is}w)'=i\omega_0e^{is}w=-A_0^{\mathsf T}e^{is}wであり、実部と虚部を取ってω0ψj′=−A0Tψj\omega_0\psi_j'=-A_0^{\mathsf T}\psi_jを得る。c1,c2∈Rc_1,c_2\in\Rがc1ψ1+c2ψ2=0c_1\psi_1+c_2\psi_2=0を満たすならば、主張 (3)により0=∫02π(c1ψ1+c2ψ2)Tϕj ds=πcj0=\int_0^{2\pi}(c_1\psi_1+c_2\psi_2)^{\mathsf T}\phi_j\,ds=\pi c_j(j=1,2j=1,2)である。ω0ψ′=−A0Tψ\omega_0\psi'=-A_0^{\mathsf T}\psiはψ′=−ATψ\psi'=-A^{\mathsf T}\psiと同値であるから、第一段によりその解でC2π1C^1_{2\pi}に属するものの全体の次元は22以下であり、ψ1,ψ2\psi_1,\psi_2はその基底である。

主張 (4)を示す。AT=A∗A^{\mathsf T}=A^*であるから、ψ′=−ATψ\psi'=-A^{\mathsf T}\psiはu′=Auu'=Auの随伴系である。§E13.13 系 4.6をK=R\mathbb K=\R、d=nd=n、T=2πT=2\pi、定数係数AAについてLFL_Fに適用する。同系のN=ker⁡LF\mathcal N=\ker L_Fは主張 (1)によりϕ1,ϕ2\phi_1,\phi_2で張られ、N∗\mathcal N^*は主張 (2)によりψ1,ψ2\psi_1,\psi_2で張られる。したがって§E13.13 系 4.6 (3)によりim⁡L0=im⁡LF\im L_0=\im L_Fは主張の集合に等しく、§E13.13 系 4.6 (4)のX1X_1は本補題のX1X_1に一致する。§E13.13 系 4.6 (4)によりC2π1=ker⁡L0⊕X1C^1_{2\pi}=\ker L_0\oplus X_1であり、L0∣X1=ω0LF∣X1L_0|_{X_1}=\omega_0L_F|_{X_1}はX1X_1からim⁡L0\im L_0への全単射であって、CF>0C_F>0が存在して、すべてのu∈X1u\in X_1について∥u∥C1≤CF∥LFu∥∞=(CF/ω0)∥L0u∥∞\lVert u\rVert_{C^1}\le C_F\lVert L_Fu\rVert_\infty=(C_F/\omega_0)\lVert L_0u\rVert_\inftyである。ℓj(u)=∫02πϕjTu ds\ell_j(u)=\int_0^{2\pi}\phi_j^{\mathsf T}u\,dsは Cauchy–Schwarz の不等式により∣ℓj(u)∣≤2π∥ϕj∥∞∥u∥C1\lvert\ell_j(u)\rvert\le2\pi\lVert\phi_j\rVert_\infty\lVert u\rVert_{C^1}を満たす線形写像であるから連続であり、X1=ker⁡ℓ1∩ker⁡ℓ2X_1=\ker\ell_1\cap\ker\ell_2は閉部分空間である。

主張 (5)を示す。β=w∗Bv\beta=w^*Bvとおく。§E13.13 系 4.6 (1)と主張 (4)と§E10.1 命題 1.4 (1)によりX1X_1は Banach 空間である。(um,am,bm)(u_m,a_m,b_m)がXXの Cauchy 列ならば、(um)(u_m)、(am)(a_m)、(bm)(b_m)はそれぞれX1X_1、R\R、R\Rの Cauchy 列であって収束するから、(um,am,bm)(u_m,a_m,b_m)はXXで収束する。§E13.13 系 4.6 (1)によりL0=ω0LFL_0=\omega_0L_Fは有界であり、∥LB(u,a,b)∥∞≤(∥L0∥+∥ϕ2∥∞+∥B∥∥ϕ1∥∞)∥(u,a,b)∥X\lVert L_B(u,a,b)\rVert_\infty\le(\lVert L_0\rVert+\lVert\phi_2\rVert_\infty+\lVert B\rVert\lVert\phi_1\rVert_\infty)\lVert(u,a,b)\rVert_Xであるから、LBL_Bは有界線形作用素である。

h∈C2πh\in C_{2\pi}と(u,a,b)∈X(u,a,b)\in Xについて、ϕ1′=−ϕ2\phi_1'=-\phi_2であるから、LB(u,a,b)=hL_B(u,a,b)=hであることはL0u=h+aϕ2+bBϕ1L_0u=h+a\phi_2+bB\phi_1と同値である。主張 (4)により、これを満たすu∈X1u\in X_1が存在することはh+aϕ2+bBϕ1∈im⁡L0h+a\phi_2+bB\phi_1\in\im L_0と同値であり、そのときuuはただ一つである。hj=∫02πψjTh dsh_j=\int_0^{2\pi}\psi_j^{\mathsf T}h\,dsとおくと、主張 (3)により、h+aϕ2+bBϕ1∈im⁡L0h+a\phi_2+bB\phi_1\in\im L_0は

h1+πbℜβ=0,h2+πa−πbℑβ=0h_1+\pi b\Re\beta=0,\qquad h_2+\pi a-\pi b\Im\beta=0

と同値である。ℜβ≠0\Re\beta\ne0であるから、この連立一次方程式はただ一つの解b=−h1/(πℜβ)b=-h_1/(\pi\Re\beta)、a=−h2/π+bℑβa=-h_2/\pi+b\Im\betaをもつ。したがってLBL_Bは全単射である。∣hj∣≤2π∥ψj∥∞∥h∥∞\lvert h_j\rvert\le2\pi\lVert\psi_j\rVert_\infty\lVert h\rVert_\inftyであるから、hhによらないC1>0C_1>0が存在して、すべてのh∈C2πh\in C_{2\pi}について∣a∣+∣b∣≤C1∥h∥∞\lvert a\rvert+\lvert b\rvert\le C_1\lVert h\rVert_\inftyであり、主張 (4)のCCについて∥u∥C1≤C(1+C1(∥ϕ2∥∞+∥B∥∥ϕ1∥∞))∥h∥∞\lVert u\rVert_{C^1}\le C\bigl(1+C_1(\lVert\phi_2\rVert_\infty+\lVert B\rVert\lVert\phi_1\rVert_\infty)\bigr)\lVert h\rVert_\inftyである。線形な全単射の逆写像は線形であるから、LB−1L_B^{-1}は有界線形作用素である。▨

補題 2.2.n∈N≥1n\in\NNとし、p∈C(R;Rn)p\in C(\R;\R^n)はすべてのs∈Rs\in\Rでp(s+2π)=p(s)p(s+2\pi)=p(s)を満たし、∫02πe−isv∗p(s) ds≠0\int_0^{2\pi}e^{-is}v^*p(s)\,ds\ne0を満たすv∈Cnv\in\C^nが存在するとする。このとき、22以上のすべての整数mmについて、p(s+2π/m)≠p(s)p(s+2\pi/m)\ne p(s)を満たすs∈Rs\in\Rが存在する。特にppは定数でない。

証明.m≥2m\ge2を整数とし、c=2π/mc=2\pi/mがすべてのs∈Rs\in\Rでp(s+c)=p(s)p(s+c)=p(s)を満たすと仮定する。I0=∫0ce−iσv∗p(σ) dσI_0=\int_0^ce^{-i\sigma}v^*p(\sigma)\,d\sigmaとおく。j=0,…,m−1j=0,\ldots,m-1について、s=jc+σs=jc+\sigmaと置換しp(jc+σ)=p(σ)p(jc+\sigma)=p(\sigma)を用いると∫jc(j+1)ce−isv∗p(s) ds=e−ijcI0\int_{jc}^{(j+1)c}e^{-is}v^*p(s)\,ds=e^{-ijc}I_0である。ζ=e−2πi/m\zeta=e^{-2\pi i/m}はζ≠1\zeta\ne1、ζm=1\zeta^m=1を満たすから

∫02πe−isv∗p(s) ds=I0∑j=0m−1ζj=I01−ζm1−ζ=0\int_0^{2\pi}e^{-is}v^*p(s)\,ds=I_0\sum_{j=0}^{m-1}\zeta^j=I_0\frac{1-\zeta^m}{1-\zeta}=0

である。これは∫02πe−isv∗p(s) ds≠0\int_0^{2\pi}e^{-is}v^*p(s)\,ds\ne0に反する。定数関数はすべてのssでp(s+π)=p(s)p(s+\pi)=p(s)を満たすから、m=2m=2の場合によりppは定数でない。▨

3 Hopf 分岐定理

定理 3.1 (Hopf 分岐定理).n∈N≥1n\in\NN、n≥2n\ge2とし、U⊆RnU\subseteq\R^nを開集合、Λ⊆R\Lambda\subseteq\Rを00を含む開区間、F∈C2(U×Λ,Rn)F\in C^2(U\times\Lambda,\R^n)とする。x0∈Ux_0\in UはF(x0,0)=0F(x_0,0)=0を満たし、A0=DxF(x0,0)A_0=D_xF(x_0,0)とω0>0\omega_0>0は次を満たすとする。

  1. ±iω0\pm i\omega_0はA0A_0の代数的に単純な固有値である。
  2. すべてのk∈Z∖{1,−1}k\in\Z\setminus\{1,-1\}について、ikω0ik\omega_0はA0A_0の固有値でない。

Λ0⊆Λ\Lambda_0\subseteq\Lambdaを00を含む開区間とし、ξ:Λ0→U\xi:\Lambda_0\to Uを、ξ(0)=x0\xi(0)=x_0とすべてのμ∈Λ0\mu\in\Lambda_0についてのF(ξ(μ),μ)=0F(\xi(\mu),\mu)=0を満たすC2C^2級の写像とする(条件 (b)のk=0k=0によりA0A_0は可逆であり、§E13.22 定理 3.3 (1)によりこのようなΛ0\Lambda_0とξ\xiが存在する)。A(μ)=DxF(ξ(μ),μ)A(\mu)=D_xF(\xi(\mu),\mu)とおき、§E13.22 補題 4.2をAA、μ0=0\mu_0=0、λ0=iω0\lambda_0=i\omega_0に適用して得る固有値の分枝をλ\lambdaとする。さらに次を仮定する。

  1. ℜλ′(0)≠0\Re\lambda'(0)\ne0である。

v∈Cn∖{0}v\in\C^n\setminus\{0\}をA0v=iω0vA_0v=i\omega_0vを満たすベクトルとし、ϕ1(s)=ℜ(eisv)\phi_1(s)=\Re(e^{is}v)、ϕ2(s)=ℑ(eisv)\phi_2(s)=\Im(e^{is}v)とおく。Rn\R^nには Euclid ノルム∣⋅∣\lvert\cdot\rvertを用い、C2π1C^1_{2\pi}を2π2\pi周期のC1C^1級写像u:R→Rnu:\R\to\R^nの全体にノルム∥u∥C1=sup⁡s∣u(s)∣+sup⁡s∣u′(s)∣\lVert u\rVert_{C^1}=\sup_s\lvert u(s)\rvert+\sup_s\lvert u'(s)\rvertを入れた空間とする。このとき、ε0>0\varepsilon_0>0、ρ>0\rho>0と、連続写像(−ε0,ε0)∋ε↦(μ(ε),ω(ε),uε)∈Λ0×(0,∞)×C2π1(-\varepsilon_0,\varepsilon_0)\ni\varepsilon\mapsto(\mu(\varepsilon),\omega(\varepsilon),u_\varepsilon)\in\Lambda_0\times(0,\infty)\times C^1_{2\pi}が存在して、次が成り立つ。

  1. μ(0)=0\mu(0)=0、ω(0)=ω0\omega(0)=\omega_0、u0=0u_0=0である。すべてのε∈(−ε0,ε0)\varepsilon\in(-\varepsilon_0,\varepsilon_0)について、∣μ(ε)∣≤ρ\lvert\mu(\varepsilon)\rvert\le\rho、∣ω(ε)−ω0∣≤ρ\lvert\omega(\varepsilon)-\omega_0\rvert\le\rho、∥uε−εϕ1∥C1≤ρ∣ε∣\lVert u_\varepsilon-\varepsilon\phi_1\rVert_{C^1}\le\rho\lvert\varepsilon\rvert、∫02πϕjT(uε−εϕ1) ds=0\int_0^{2\pi}\phi_j^{\mathsf T}(u_\varepsilon-\varepsilon\phi_1)\,ds=0(j=1,2j=1,2)であり、すべてのs∈Rs\in\Rでξ(μ(ε))+uε(s)∈U\xi(\mu(\varepsilon))+u_\varepsilon(s)\in Uかつω(ε)uε′(s)=F(ξ(μ(ε))+uε(s),μ(ε))\omega(\varepsilon)u_\varepsilon'(s)=F\bigl(\xi(\mu(\varepsilon))+u_\varepsilon(s),\mu(\varepsilon)\bigr)である。さらに、ε→0\varepsilon\to0(ε≠0\varepsilon\ne0)のとき∥uε−εϕ1∥C1/∣ε∣→0\lVert u_\varepsilon-\varepsilon\phi_1\rVert_{C^1}/\lvert\varepsilon\rvert\to0である。
  2. 0<∣ε∣<ε00<\lvert\varepsilon\rvert<\varepsilon_0とし、xε(t)=ξ(μ(ε))+uε(ω(ε)t)x_\varepsilon(t)=\xi(\mu(\varepsilon))+u_\varepsilon(\omega(\varepsilon)t)(t∈Rt\in\R)とおく。このとき、xεx_\varepsilonはx′=F(x,μ(ε))x'=F(x,\mu(\varepsilon))のR\R上の古典解であり、xε(0)x_\varepsilon(0)はF(⋅,μ(ε))F(\cdot,\mu(\varepsilon))の流れの周期点で、その最小周期は2π/ω(ε)2\pi/\omega(\varepsilon)である。
  3. 0<∣ε∣<ε00<\lvert\varepsilon\rvert<\varepsilon_0とし、(μ,ω,u)∈Λ0×R×C2π1(\mu,\omega,u)\in\Lambda_0\times\R\times C^1_{2\pi}が∣μ∣≤ρ\lvert\mu\rvert\le\rho、∣ω−ω0∣≤ρ\lvert\omega-\omega_0\rvert\le\rho、∥u−εϕ1∥C1≤ρ∣ε∣\lVert u-\varepsilon\phi_1\rVert_{C^1}\le\rho\lvert\varepsilon\rvert、∫02πϕjT(u−εϕ1) ds=0\int_0^{2\pi}\phi_j^{\mathsf T}(u-\varepsilon\phi_1)\,ds=0(j=1,2j=1,2)を満たし、すべてのs∈Rs\in\Rでξ(μ)+u(s)∈U\xi(\mu)+u(s)\in Uかつωu′(s)=F(ξ(μ)+u(s),μ)\omega u'(s)=F(\xi(\mu)+u(s),\mu)を満たすとする。このとき(μ,ω,u)=(μ(ε),ω(ε),uε)(\mu,\omega,u)=(\mu(\varepsilon),\omega(\varepsilon),u_\varepsilon)である。

証明.§E13.22 補題 4.1 (2)により、w∗A0=iω0w∗w^*A_0=i\omega_0w^*とw∗v≠0w^*v\ne0を満たすw∈Cnw\in\C^nが存在する。wwをw∗v‾\overline{w^*v}で割って、w∗v=1w^*v=1とする。DxFD_xFとξ\xiはC1C^1級であるから、A:Λ0→Rn×nA:\Lambda_0\to\R^{n\times n}はC1C^1級である。B1=A′(0)B_1=A'(0)とおく。§E13.22 補題 4.2 (2)によりλ′(0)=w∗B1v\lambda'(0)=w^*B_1vであるから、条件 (c)によりℜ(w∗B1v)≠0\Re(w^*B_1v)\ne0である。補題 2.1をA0A_0、ω0\omega_0、vv、wwに適用し、同補題のC2πC_{2\pi}、∥⋅∥∞\lVert\cdot\rVert_\infty、L0L_0、X1X_1を用いる。補題 2.1 (5)により、X=X1×R×RX=X_1\times\R\times\Rはノルム∥(z,ω,μ)∥X=max⁡{∥z∥C1,∣ω∣,∣μ∣}\lVert(z,\omega,\mu)\rVert_X=\max\{\lVert z\rVert_{C^1},\lvert\omega\rvert,\lvert\mu\rvert\}について Banach 空間であり、L=LB1L=L_{B_1}、すなわちL(z,ω,μ)=L0z+ωϕ1′−μB1ϕ1L(z,\omega,\mu)=L_0z+\omega\phi_1'-\mu B_1\phi_1は全単射な有界線形作用素で、L−1L^{-1}も有界線形作用素である。

W={(y,μ)∈Rn×Λ0:ξ(μ)+y∈U}W=\{(y,\mu)\in\R^n\times\Lambda_0:\xi(\mu)+y\in U\}は開集合であり、G(y,μ)=F(ξ(μ)+y,μ)G(y,\mu)=F(\xi(\mu)+y,\mu)はWW上のC2C^2級の写像で、すべてのμ∈Λ0\mu\in\Lambda_0についてG(0,μ)=0G(0,\mu)=0、DyG(0,μ)=A(μ)D_yG(0,\mu)=A(\mu)である。δ>0\delta>0を[−2δ,2δ]⊆Λ0[-2\delta,2\delta]\subseteq\Lambda_0を満たすように取る。{0}×[−2δ,2δ]\{0\}\times[-2\delta,2\delta]は開集合WWのコンパクト部分集合であるから、{(y,μ):∣y∣≤2r, ∣μ∣≤2δ}⊆W\{(y,\mu):\lvert y\rvert\le2r,\ \lvert\mu\rvert\le2\delta\}\subseteq Wを満たすr>0r>0が存在する。そのようなrrを一つ取る。R1=∥ϕ1∥∞+1R_1=\lVert\phi_1\rVert_\infty+1、ε1=r/R1\varepsilon_1=r/R_1とし、

E={(a,ν)∈Rn×R:∣a∣<2R1, ∣ν∣<2δ},D=(−ε1,ε1)×EE=\{(a,\nu)\in\R^n\times\R:\lvert a\rvert<2R_1,\ \lvert\nu\rvert<2\delta\},\qquad D=(-\varepsilon_1,\varepsilon_1)\times E

とおく。(ε,a,ν)∈D(\varepsilon,a,\nu)\in Dならば∣εa∣<2r\lvert\varepsilon a\rvert<2rかつ∣ν∣<2δ\lvert\nu\rvert<2\deltaであるから(εa,ν)∈W(\varepsilon a,\nu)\in Wであり、τ∈[0,1]\tau\in[0,1]について(τε,a,ν)∈D(\tau\varepsilon,a,\nu)\in Dである。したがって、(ε,a,ν)∈D(\varepsilon,a,\nu)\in Dについて

k(ε,a,ν)=∫01DyG(τεa,ν) a dτk(\varepsilon,a,\nu)=\int_0^1D_yG(\tau\varepsilon a,\nu)\,a\,d\tau

が定まる。

主張 3.1.1.(ε,a,ν)∈D(\varepsilon,a,\nu)\in Dについて、ε≠0\varepsilon\ne0ならばk(ε,a,ν)=ε−1G(εa,ν)k(\varepsilon,a,\nu)=\varepsilon^{-1}G(\varepsilon a,\nu)であり、k(0,a,ν)=A(ν)ak(0,a,\nu)=A(\nu)aである。各ε∈(−ε1,ε1)\varepsilon\in(-\varepsilon_1,\varepsilon_1)について、(a,ν)↦k(ε,a,ν)(a,\nu)\mapsto k(\varepsilon,a,\nu)はEE上でC1C^1級であり、その(a,ν)∈E(a,\nu)\in Eでの導関数は

Θ(ε,a,ν)=(DyG(εa,ν), ∫01∂μDyG(τεa,ν) a dτ)∈Rn×(n+1)\Theta(\varepsilon,a,\nu)=\Bigl(D_yG(\varepsilon a,\nu),\ \int_0^1\partial_\mu D_yG(\tau\varepsilon a,\nu)\,a\,d\tau\Bigr)\in\R^{n\times(n+1)}

である。kkとΘ\ThetaはDD上で連続である。

証明.(ε,(a,ν))∈(−ε1,ε1)×E(\varepsilon,(a,\nu))\in(-\varepsilon_1,\varepsilon_1)\times Eについてg(ε,(a,ν))=G(εa,ν)g(\varepsilon,(a,\nu))=G(\varepsilon a,\nu)とおく。§E4.3 定理 1.1によりggはC2C^2級であり、∂εg(ε,(a,ν))=DyG(εa,ν)a\partial_\varepsilon g(\varepsilon,(a,\nu))=D_yG(\varepsilon a,\nu)a、∂νg(ε,(a,ν))=∂μG(εa,ν)\partial_\nu g(\varepsilon,(a,\nu))=\partial_\mu G(\varepsilon a,\nu)である。すべての(a,ν)∈E(a,\nu)\in Eについてg(0,(a,ν))=G(0,ν)=0g(0,(a,\nu))=G(0,\nu)=0であるから、§E13.22 補題 7.1を、滑らかさの指数22、p=n+1p=n+1、m=nm=n、同補題の区間JJと開集合Λ\Lambdaをそれぞれ(−ε1,ε1)(-\varepsilon_1,\varepsilon_1)とEEとして、ggに適用することができる。∂εg(τε,(a,ν))=DyG(τεa,ν)a\partial_\varepsilon g(\tau\varepsilon,(a,\nu))=D_yG(\tau\varepsilon a,\nu)aであるから、同補題でggから定まる関数はkkに等しく、同補題によりkkはDD上でC1C^1級である。§E13.22 補題 7.1 (1)により、ε≠0\varepsilon\ne0ならばk(ε,a,ν)=ε−1G(εa,ν)k(\varepsilon,a,\nu)=\varepsilon^{-1}G(\varepsilon a,\nu)である。§E13.22 補題 7.1 (2)のi=0i=0、β=0\beta=0の場合によりk(0,a,ν)=∂εg(0,(a,ν))=DyG(0,ν)a=A(ν)ak(0,a,\nu)=\partial_\varepsilon g(0,(a,\nu))=D_yG(0,\nu)a=A(\nu)aである。

kkはDD上でC1C^1級であるから、各ε\varepsilonについて(a,ν)↦k(ε,a,ν)(a,\nu)\mapsto k(\varepsilon,a,\nu)はEE上でC1C^1級であり、その導関数はDD上で連続である。ε≠0\varepsilon\ne0のとき、k(ε,a,ν)=ε−1G(εa,ν)k(\varepsilon,a,\nu)=\varepsilon^{-1}G(\varepsilon a,\nu)と§E4.3 定理 1.1によりaaについての導関数はDyG(εa,ν)D_yG(\varepsilon a,\nu)であり、ε=0\varepsilon=0のとき、k(0,a,ν)=A(ν)ak(0,a,\nu)=A(\nu)aのaaについての導関数はA(ν)=DyG(0,ν)A(\nu)=D_yG(0,\nu)である。§E13.22 補題 7.1 (2)のi=0i=0でβ\betaをν\nuの座標の単位多重指数とする場合により

∂νk(ε,a,ν)=∫01∂ε∂νg(τε,(a,ν)) dτ=∫01Dy∂μG(τεa,ν) a dτ\partial_\nu k(\varepsilon,a,\nu)=\int_0^1\partial_\varepsilon\partial_\nu g(\tau\varepsilon,(a,\nu))\,d\tau=\int_0^1D_y\partial_\mu G(\tau\varepsilon a,\nu)\,a\,d\tau

であり、GGはC2C^2級であるから§E4.4 定理 2.1によりDy∂μG=∂μDyGD_y\partial_\mu G=\partial_\mu D_yGである。したがって導関数はΘ(ε,a,ν)\Theta(\varepsilon,a,\nu)であり、kkとΘ\ThetaはDD上で連続である。▨

x∗=(0,ω0,0)∈Xx_*=(0,\omega_0,0)\in X、rX=min⁡{1,δ,ω0/2}r_X=\min\{1,\delta,\omega_0/2\}、P=(−ε1/2,ε1/2)P=(-\varepsilon_1/2,\varepsilon_1/2)とし、PPにはR\Rの距離を入れる。BX(x∗,r′)B_X(x_*,r')、B‾X(x∗,r′)\overline B_X(x_*,r')をXXの開球・閉球とする。(z,ω,μ)∈BX(x∗,rX)(z,\omega,\mu)\in B_X(x_*,r_X)、ε∈P\varepsilon\in P、s∈Rs\in\Rについて∣ϕ1(s)+z(s)∣<R1\lvert\phi_1(s)+z(s)\rvert<R_1、∣μ∣<δ\lvert\mu\rvert<\deltaであるから、(ε,ϕ1(s)+z(s),μ)(\varepsilon,\phi_1(s)+z(s),\mu)はコンパクト集合D′=[−ε1/2,ε1/2]×{a∈Rn:∣a∣≤R1}×[−δ,δ]⊆DD'=[-\varepsilon_1/2,\varepsilon_1/2]\times\{a\in\R^n:\lvert a\rvert\le R_1\}\times[-\delta,\delta]\subseteq Dに属する。H:P×BX(x∗,rX)→C2πH:P\times B_X(x_*,r_X)\to C_{2\pi}を

H(ε;z,ω,μ)(s)=ω(ϕ1′(s)+z′(s))−k(ε,ϕ1(s)+z(s),μ)(s∈R)H(\varepsilon;z,\omega,\mu)(s)=\omega\bigl(\phi_1'(s)+z'(s)\bigr)-k\bigl(\varepsilon,\phi_1(s)+z(s),\mu\bigr)\qquad(s\in\R)

で定める。右辺はssについて連続で2π2\pi周期である。kkはD′D'上で一様連続であるから、HHは連続である。k(0,a,0)=A0ak(0,a,0)=A_0aであるから、補題 2.1 (1)によりH(0;x∗)=ω0ϕ1′−A0ϕ1=L0ϕ1=0H(0;x_*)=\omega_0\phi_1'-A_0\phi_1=L_0\phi_1=0である。

主張 3.1.2. 任意のη>0\eta>0に対してσ∈(0,rX]\sigma\in(0,r_X]が存在して、∣ε∣<σ\lvert\varepsilon\rvert<\sigmaかつx=(z,ω,μ), x~=(z~,ω~,μ~)∈BX(x∗,σ)x=(z,\omega,\mu),\ \tilde x=(\tilde z,\tilde\omega,\tilde\mu)\in B_X(x_*,\sigma)ならば、∥H(ε;x)−H(ε;x~)−L(x−x~)∥∞≤4η∥x−x~∥X\lVert H(\varepsilon;x)-H(\varepsilon;\tilde x)-L(x-\tilde x)\rVert_\infty\le4\eta\lVert x-\tilde x\rVert_Xである。

証明.η>0\eta>0とする。主張 3.1.1によりΘ\Thetaはコンパクト集合D′D'上で一様連続であるから、σ′>0\sigma'>0が存在して、D′D'の二点の Euclid 距離がσ′\sigma'未満ならば、その二点でのΘ\Thetaの値の差の作用素ノルムはη\eta以下である。σ=min⁡{rX,ε1/2,σ′/3,η}\sigma=\min\{r_X,\varepsilon_1/2,\sigma'/3,\eta\}とおき、ε\varepsilon、xx、x~\tilde xを主張のとおりに取る。p=ϕ1+zp=\phi_1+z、p~=ϕ1+z~\tilde p=\phi_1+\tilde zとおく。L(x−x~)=ω0(z−z~)′−A0(z−z~)+(ω−ω~)ϕ1′−(μ−μ~)B1ϕ1L(x-\tilde x)=\omega_0(z-\tilde z)'-A_0(z-\tilde z)+(\omega-\tilde\omega)\phi_1'-(\mu-\tilde\mu)B_1\phi_1であるから、H(ε;x)−H(ε;x~)−L(x−x~)=I1−I2H(\varepsilon;x)-H(\varepsilon;\tilde x)-L(x-\tilde x)=I_1-I_2であり、ここで

I1=ωp′−ω~p~′−ω0(z−z~)′−(ω−ω~)ϕ1′=(ω−ω0)(z−z~)′+(ω−ω~)z~′,I_1=\omega p'-\tilde\omega\tilde p'-\omega_0(z-\tilde z)'-(\omega-\tilde\omega)\phi_1'=(\omega-\omega_0)(z-\tilde z)'+(\omega-\tilde\omega)\tilde z',I2(s)=k(ε,p(s),μ)−k(ε,p~(s),μ~)−A0(z(s)−z~(s))−(μ−μ~)B1ϕ1(s)I_2(s)=k(\varepsilon,p(s),\mu)-k(\varepsilon,\tilde p(s),\tilde\mu)-A_0\bigl(z(s)-\tilde z(s)\bigr)-(\mu-\tilde\mu)B_1\phi_1(s)

である。∣ω−ω0∣<σ\lvert\omega-\omega_0\rvert<\sigma、∥z~∥C1<σ\lVert\tilde z\rVert_{C^1}<\sigmaであるから、∥I1∥∞≤2σ∥x−x~∥X≤2η∥x−x~∥X\lVert I_1\rVert_\infty\le2\sigma\lVert x-\tilde x\rVert_X\le2\eta\lVert x-\tilde x\rVert_Xである。

s∈Rs\in\Rを固定し、EE上でfs(a,ν)=k(ε,a,ν)−A0a−νB1ϕ1(s)f_s(a,\nu)=k(\varepsilon,a,\nu)-A_0a-\nu B_1\phi_1(s)とおく。I2(s)=fs(p(s),μ)−fs(p~(s),μ~)I_2(s)=f_s(p(s),\mu)-f_s(\tilde p(s),\tilde\mu)である。Θ(0,ϕ1(s),0)=(A0,B1ϕ1(s))\Theta(0,\phi_1(s),0)=(A_0,B_1\phi_1(s))であるから、主張 3.1.1によりfsf_sはEE上でC1C^1級で、(a,ν)(a,\nu)での導関数はΘ(ε,a,ν)−Θ(0,ϕ1(s),0)\Theta(\varepsilon,a,\nu)-\Theta(0,\phi_1(s),0)である。(p~(s),μ~)(\tilde p(s),\tilde\mu)と(p(s),μ)(p(s),\mu)を結ぶ線分の点は、∣a′∣<σ\lvert a'\rvert<\sigma、∣ν′∣<σ\lvert\nu'\rvert<\sigmaを満たす(ϕ1(s)+a′,ν′)(\phi_1(s)+a',\nu')の形であり、EEに属する。(ε,ϕ1(s)+a′,ν′)(\varepsilon,\phi_1(s)+a',\nu')と(0,ϕ1(s),0)(0,\phi_1(s),0)はともにD′D'に属し、その Euclid 距離は∣ε∣+∣a′∣+∣ν′∣\lvert\varepsilon\rvert+\lvert a'\rvert+\lvert\nu'\rvert以下であり、∣ε∣+∣a′∣+∣ν′∣<3σ≤σ′\lvert\varepsilon\rvert+\lvert a'\rvert+\lvert\nu'\rvert<3\sigma\le\sigma'であるから、線分上でfsf_sの導関数の作用素ノルムはη\eta以下である。§E4.3 定理 2.1(Rn+1\R^{n+1}の Euclid ノルム)により

∣I2(s)∣≤η(∣z(s)−z~(s)∣2+∣μ−μ~∣2)1/2≤2η∥x−x~∥X\lvert I_2(s)\rvert\le\eta\bigl(\lvert z(s)-\tilde z(s)\rvert^2+\lvert\mu-\tilde\mu\rvert^2\bigr)^{1/2}\le2\eta\lVert x-\tilde x\rVert_X

である。ssは任意であるから∥I1−I2∥∞≤4η∥x−x~∥X\lVert I_1-I_2\rVert_\infty\le4\eta\lVert x-\tilde x\rVert_Xである。▨

任意のη>0\eta>0について主張 3.1.2をη/4\eta/4に適用すると、補題 1.1の剰余の条件がη\etaについて成り立つ。したがって、同補題の仮定はXX、Y=C2πY=C_{2\pi}、PP、p0=0p_0=0、x∗x_*、rXr_X、HH、LLについて成り立つ。同補題のρ1\rho_1、ρ2∈(0,rX)\rho_2\in(0,r_X)を取り、ε0=min⁡{ρ1,ε1/2}\varepsilon_0=\min\{\rho_1,\varepsilon_1/2\}、ρ=ρ2\rho=\rho_2とおく。ε∈(−ε0,ε0)\varepsilon\in(-\varepsilon_0,\varepsilon_0)について、H(ε;x)=0H(\varepsilon;x)=0を満たすx∈B‾X(x∗,ρ)x\in\overline B_X(x_*,\rho)をただ一つ取り、x=(zε,ω(ε),μ(ε))x=(z_\varepsilon,\omega(\varepsilon),\mu(\varepsilon))と書く。補題 1.1 (2)によりε↦(zε,ω(ε),μ(ε))\varepsilon\mapsto(z_\varepsilon,\omega(\varepsilon),\mu(\varepsilon))は連続で、(z0,ω(0),μ(0))=x∗(z_0,\omega(0),\mu(0))=x_*である。uε=ε(ϕ1+zε)u_\varepsilon=\varepsilon(\phi_1+z_\varepsilon)とおく。

主張 (1)を示す。ε↦uε\varepsilon\mapsto u_\varepsilonはC2π1C^1_{2\pi}への連続写像であり、μ(0)=0\mu(0)=0、ω(0)=ω0\omega(0)=\omega_0、u0=0u_0=0である。∣μ(ε)∣≤ρ<δ\lvert\mu(\varepsilon)\rvert\le\rho<\deltaであるからμ(ε)∈Λ0\mu(\varepsilon)\in\Lambda_0であり、ω(ε)≥ω0−ρ>ω0/2>0\omega(\varepsilon)\ge\omega_0-\rho>\omega_0/2>0である。∥uε−εϕ1∥C1=∣ε∣∥zε∥C1≤ρ∣ε∣\lVert u_\varepsilon-\varepsilon\phi_1\rVert_{C^1}=\lvert\varepsilon\rvert\lVert z_\varepsilon\rVert_{C^1}\le\rho\lvert\varepsilon\rvertであり、zε∈X1z_\varepsilon\in X_1であるから∫02πϕjT(uε−εϕ1) ds=ε∫02πϕjTzε ds=0\int_0^{2\pi}\phi_j^{\mathsf T}(u_\varepsilon-\varepsilon\phi_1)\,ds=\varepsilon\int_0^{2\pi}\phi_j^{\mathsf T}z_\varepsilon\,ds=0である。∣uε(s)∣<(ε1/2)R1=r/2\lvert u_\varepsilon(s)\rvert<(\varepsilon_1/2)R_1=r/2であるから(uε(s),μ(ε))∈W(u_\varepsilon(s),\mu(\varepsilon))\in W、すなわちξ(μ(ε))+uε(s)∈U\xi(\mu(\varepsilon))+u_\varepsilon(s)\in Uである。ε≠0\varepsilon\ne0のとき、主張 3.1.1によりεk(ε,ϕ1(s)+zε(s),μ(ε))=G(uε(s),μ(ε))\varepsilon k(\varepsilon,\phi_1(s)+z_\varepsilon(s),\mu(\varepsilon))=G(u_\varepsilon(s),\mu(\varepsilon))であるから、H(ε;zε,ω(ε),μ(ε))=0H(\varepsilon;z_\varepsilon,\omega(\varepsilon),\mu(\varepsilon))=0の両辺にε\varepsilonを掛けると

ω(ε)uε′(s)=G(uε(s),μ(ε))=F(ξ(μ(ε))+uε(s),μ(ε))(s∈R)\omega(\varepsilon)u_\varepsilon'(s)=G(u_\varepsilon(s),\mu(\varepsilon))=F\bigl(\xi(\mu(\varepsilon))+u_\varepsilon(s),\mu(\varepsilon)\bigr)\qquad(s\in\R)

を得る。ε=0\varepsilon=0のときはu0=0u_0=0とF(ξ(0),0)=F(x0,0)=0F(\xi(0),0)=F(x_0,0)=0により両辺は00である。ε≠0\varepsilon\ne0について∥uε−εϕ1∥C1/∣ε∣=∥zε∥C1\lVert u_\varepsilon-\varepsilon\phi_1\rVert_{C^1}/\lvert\varepsilon\rvert=\lVert z_\varepsilon\rVert_{C^1}であり、ε→0\varepsilon\to0のときzε→z0=0z_\varepsilon\to z_0=0である。

主張 (3)を示す。ε\varepsilonと(μ,ω,u)(\mu,\omega,u)を主張のとおりに取り、z=(u−εϕ1)/εz=(u-\varepsilon\phi_1)/\varepsilonとおく。z∈C2π1z\in C^1_{2\pi}、∫02πϕjTz ds=0\int_0^{2\pi}\phi_j^{\mathsf T}z\,ds=0であるからz∈X1z\in X_1であり、∥z∥C1≤ρ\lVert z\rVert_{C^1}\le\rho、∣ω−ω0∣≤ρ\lvert\omega-\omega_0\rvert\le\rho、∣μ∣≤ρ\lvert\mu\rvert\le\rhoであるから(z,ω,μ)∈B‾X(x∗,ρ)(z,\omega,\mu)\in\overline B_X(x_*,\rho)である。ρ<rX\rho<r_Xにより、すべてのssで(ε,ϕ1(s)+z(s),μ)∈D(\varepsilon,\phi_1(s)+z(s),\mu)\in Dであり、ξ(μ)+u(s)∈U\xi(\mu)+u(s)\in Uにより(u(s),μ)∈W(u(s),\mu)\in W、G(u(s),μ)=F(ξ(μ)+u(s),μ)G(u(s),\mu)=F(\xi(\mu)+u(s),\mu)である。主張 3.1.1によりk(ε,ϕ1(s)+z(s),μ)=ε−1G(u(s),μ)=ε−1ωu′(s)=ω(ϕ1′(s)+z′(s))k(\varepsilon,\phi_1(s)+z(s),\mu)=\varepsilon^{-1}G(u(s),\mu)=\varepsilon^{-1}\omega u'(s)=\omega(\phi_1'(s)+z'(s))であるから、H(ε;z,ω,μ)=0H(\varepsilon;z,\omega,\mu)=0である。補題 1.1 (1)の一意性により(z,ω,μ)=(zε,ω(ε),μ(ε))(z,\omega,\mu)=(z_\varepsilon,\omega(\varepsilon),\mu(\varepsilon))であり、u=ε(ϕ1+zε)=uεu=\varepsilon(\phi_1+z_\varepsilon)=u_\varepsilonである。

主張 (2)を示す。0<∣ε∣<ε00<\lvert\varepsilon\rvert<\varepsilon_0とし、μ=μ(ε)\mu=\mu(\varepsilon)、ω=ω(ε)\omega=\omega(\varepsilon)と書く。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 連続である。F(⋅,μ)F(\cdot,\mu)の流れをφ\varphiとする。xεx_\varepsilonはC1C^1級であり、主張 (1)によりxε′(t)=ωuε′(ωt)=F(xε(t),μ)x_\varepsilon'(t)=\omega u_\varepsilon'(\omega t)=F(x_\varepsilon(t),\mu)であるから、xεx_\varepsilonはR\R上の古典解である。q=xε(0)q=x_\varepsilon(0)とおくと、§E13.15 補題 1.2(J=RJ=\R、τ=0\tau=0)によりI(q)=RI(q)=\Rであり、すべてのt∈Rt\in\Rでφt(q)=xε(t)\varphi_t(q)=x_\varepsilon(t)である。

e−isv∗=ϕ1(s)T−iϕ2(s)Te^{-is}v^*=\phi_1(s)^{\mathsf T}-i\phi_2(s)^{\mathsf T}であり、zε∈X1z_\varepsilon\in X_1であるから∫02πe−isv∗zε(s) ds=0\int_0^{2\pi}e^{-is}v^*z_\varepsilon(s)\,ds=0である。ϕ1(s)=12(eisv+e−isvˉ)\phi_1(s)=\frac12(e^{is}v+e^{-is}\bar v)と∫02πe−2is ds=0\int_0^{2\pi}e^{-2is}\,ds=0により

∫02πe−isv∗uε(s) ds=ε∫02π12(v∗v+e−2isv∗vˉ) ds=πε v∗v≠0\int_0^{2\pi}e^{-is}v^*u_\varepsilon(s)\,ds=\varepsilon\int_0^{2\pi}\tfrac12\bigl(v^*v+e^{-2is}v^*\bar v\bigr)\,ds=\pi\varepsilon\,v^*v\ne0

である。補題 2.2によりuεu_\varepsilonは定数でなく、したがってxεx_\varepsilonは定数でない。F(q,μ)=0F(q,\mu)=0ならば、§E13.15 命題 4.3 (2)によりγ(q)={q}\gamma(q)=\{q\}であり、xεx_\varepsilonは定数となるから、F(q,μ)≠0F(q,\mu)\ne0である。φ2π/ω(q)=ξ(μ)+uε(2π)=q\varphi_{2\pi/\omega}(q)=\xi(\mu)+u_\varepsilon(2\pi)=qであるから、qqは周期点である。§E13.15 命題 4.3 (3)により、{t∈R:φt(q)=q}=T∗Z\{t\in\R:\varphi_t(q)=q\}=T_*\Zを満たすT∗>0T_*>0が存在し、2π/ω∈T∗Z2\pi/\omega\in T_*\Zであるから、2π/ω=mT∗2\pi/\omega=mT_*を満たす整数m≥1m\ge1が存在する。§E13.15 定理 1.3により、すべてのttでxε(t+T∗)=φt(φT∗(q))=xε(t)x_\varepsilon(t+T_*)=\varphi_t(\varphi_{T_*}(q))=x_\varepsilon(t)であるから、すべてのssでuε(s+2π/m)=uε(s)u_\varepsilon(s+2\pi/m)=u_\varepsilon(s)である。補題 2.2によりm=1m=1であり、qqの最小周期は2π/ω2\pi/\omegaである。▨

4 分岐の向きと安定性

命題 4.1.定理 3.1の仮定と記号の下で、さらにA0A_0の±iω0\pm i\omega_0以外の固有値の実部がすべて負であるとする。このとき、δ>0\delta>0が存在して(−δ,δ)⊆Λ0(-\delta,\delta)\subseteq\Lambda_0であり、0<∣μ∣<δ0<\lvert\mu\rvert<\deltaを満たすすべてのμ\muについて、F(⋅,μ)F(\cdot,\mu)の平衡点ξ(μ)\xi(\mu)は、μℜλ′(0)<0\mu\Re\lambda'(0)<0ならば漸近安定であり、μℜλ′(0)>0\mu\Re\lambda'(0)>0ならば不安定である。

証明.λ\lambdaの定義域をΛλ\Lambda_\lambdaとする。DxFD_xFとξ\xiはC1C^1級であるからAAはC1C^1級であり、§E13.22 補題 4.2によりΛλ⊆Λ0\Lambda_\lambda\subseteq\Lambda_0は00を含む開区間で、λ\lambdaはC1C^1級であり、§E13.22 補題 4.2 (1)によりλ(0)=iω0\lambda(0)=i\omega_0であって、μ∈Λλ\mu\in\Lambda_\lambdaについてλ(μ)\lambda(\mu)はA(μ)A(\mu)の固有値である。a=min⁡{−ℜν:νa=\min\{-\Re\nu:\nuはA0A_0の固有値でν≠±iω0}\nu\ne\pm i\omega_0\}(そのような固有値が無ければa=1a=1)とする。A0A_0の固有値は±iω0\pm i\omega_0と実部が−a-a以下の固有値であるから、直線{z∈C:ℜz=−a/2}\{z\in\C:\Re z=-a/2\}の上にA0A_0の固有値は無く、定理 3.1 条件 (a)により、ℜz>−a/2\Re z>-a/2にあるA0A_0の固有値の個数(代数的重複度を込める)は22である。AAは00で連続であるから、§E13.22 補題 3.2 (4)(T=Λ0T=\Lambda_0、t0=0t_0=0、c=−a/2c=-a/2)により、δ1>0\delta_1>0が存在して、μ∈Λ0\mu\in\Lambda_0、∣μ∣<δ1\lvert\mu\rvert<\delta_1ならば、ℜz>−a/2\Re z>-a/2にあるA(μ)A(\mu)の固有値の個数(代数的重複度を込める)は22であり、A(μ)A(\mu)のその他の固有値の実部は−a/2-a/2より小さい。λ\lambdaは連続でλ(0)=iω0\lambda(0)=i\omega_0であるから、δ2∈(0,δ1]\delta_2\in(0,\delta_1]が存在して、(−δ2,δ2)⊆Λλ(-\delta_2,\delta_2)\subseteq\Lambda_\lambdaであり、∣μ∣<δ2\lvert\mu\rvert<\delta_2ならば∣λ(μ)−iω0∣<min⁡{a/2,ω0}\lvert\lambda(\mu)-i\omega_0\rvert<\min\{a/2,\omega_0\}、したがってℜλ(μ)>−a/2\Re\lambda(\mu)>-a/2かつℑλ(μ)>0\Im\lambda(\mu)>0である。A(μ)A(\mu)は実行列であるからλ(μ)‾\overline{\lambda(\mu)}もA(μ)A(\mu)の固有値であり、ℑλ(μ)>0\Im\lambda(\mu)>0によりλ(μ)‾≠λ(μ)\overline{\lambda(\mu)}\ne\lambda(\mu)である。したがって∣μ∣<δ2\lvert\mu\rvert<\delta_2ならば、ℜz>−a/2\Re z>-a/2にあるA(μ)A(\mu)の固有値はλ(μ)\lambda(\mu)とλ(μ)‾\overline{\lambda(\mu)}であり、A(μ)A(\mu)のその他の固有値の実部は−a/2-a/2より小さい。

ℜλ(0)=0\Re\lambda(0)=0であり、λ\lambdaは00で微分可能であるから、μ→0\mu\to0(μ≠0\mu\ne0)のときℜλ(μ)/μ→ℜλ′(0)\Re\lambda(\mu)/\mu\to\Re\lambda'(0)である。定理 3.1 条件 (c)によりℜλ′(0)≠0\Re\lambda'(0)\ne0であるから、δ∈(0,δ2]\delta\in(0,\delta_2]が存在して、0<∣μ∣<δ0<\lvert\mu\rvert<\deltaならばℜλ(μ)/μ\Re\lambda(\mu)/\muとℜλ′(0)\Re\lambda'(0)の符号は等しく、ℜλ(μ)\Re\lambda(\mu)の符号はμℜλ′(0)\mu\Re\lambda'(0)の符号に等しい。

0<∣μ∣<δ0<\lvert\mu\rvert<\deltaとする。F(ξ(μ),μ)=0F(\xi(\mu),\mu)=0であるから、ξ(μ)\xi(\mu)はF(⋅,μ)F(\cdot,\mu)の平衡点である。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として適用することができる。DxF(ξ(μ),μ)=A(μ)D_xF(\xi(\mu),\mu)=A(\mu)である。μℜλ′(0)<0\mu\Re\lambda'(0)<0ならば、ℜλ(μ)‾=ℜλ(μ)<0\Re\overline{\lambda(\mu)}=\Re\lambda(\mu)<0であり、A(μ)A(\mu)の固有値の実部はすべて負であるから、§E13.16 定理 4.5 (1)によりξ(μ)\xi(\mu)は漸近安定である。μℜλ′(0)>0\mu\Re\lambda'(0)>0ならば、λ(μ)\lambda(\mu)は実部が正の固有値であるから、§E13.16 定理 4.5 (2)によりξ(μ)\xi(\mu)は不安定である。▨

定理 4.2.定理 3.1の仮定と記号の下で、さらにF∈C3(U×Λ,Rn)F\in C^3(U\times\Lambda,\R^n)とし、q,p∈Cnq,p\in\C^nをA0q=iω0qA_0q=i\omega_0q、q∗q=1q^*q=1、p∗A0=iω0p∗p^*A_0=i\omega_0p^*、p∗q=1p^*q=1を満たすベクトルとする。y1,y2,y3∈Cny_1,y_2,y_3\in\C^nについて

b(y1,y2)=∑j,k=1n∂xj∂xkF(x0,0) y1,jy2,k,c(y1,y2,y3)=∑j,k,l=1n∂xj∂xk∂xlF(x0,0) y1,jy2,ky3,lb(y_1,y_2)=\sum_{j,k=1}^n\partial_{x_j}\partial_{x_k}F(x_0,0)\,y_{1,j}y_{2,k},\qquad c(y_1,y_2,y_3)=\sum_{j,k,l=1}^n\partial_{x_j}\partial_{x_k}\partial_{x_l}F(x_0,0)\,y_{1,j}y_{2,k}y_{3,l}

とおき(yi,jy_{i,j}はyiy_iの第jj成分)、

ℓ1=12ω0ℜ(p∗c(q,q,qˉ)−2p∗b(q,A0−1b(q,qˉ))+p∗b(qˉ,(2iω0In−A0)−1b(q,q)))\ell_1=\frac1{2\omega_0}\Re\Bigl(p^*c(q,q,\bar q)-2p^*b\bigl(q,A_0^{-1}b(q,\bar q)\bigr)+p^*b\bigl(\bar q,(2i\omega_0I_n-A_0)^{-1}b(q,q)\bigr)\Bigr)

とおく(定理 3.1 条件 (b)によりA0A_0と2iω0In−A02i\omega_0I_n-A_0は可逆である)。定理 3.1 条件 (a)により、A0q=iω0qA_0q=i\omega_0qを満たすqqの全体とp∗A0=iω0p∗p^*A_0=i\omega_0p^*を満たすppの全体はともに一次元の複素線形空間であるから、qqは絶対値11の複素数κ\kappa倍を除いて定まり、qqをκq\kappa qに替えるとppはκp\kappa pに替わる。このとき右辺の三つの項はいずれも∣κ∣4=1\lvert\kappa\rvert^4=1倍になるから、ℓ1\ell_1はq,pq,pの選び方によらない。ℓ1\ell_1を(x0,0)(x_0,0)における第一 Lyapunov 係数 (first Lyapunov coefficient) という。0<∣ε∣<ε00<\lvert\varepsilon\rvert<\varepsilon_0について、Γε\Gamma_\varepsilonを定理 3.1 (2)のxε(0)x_\varepsilon(0)を通る周期軌道、φε\varphi^\varepsilonをF(⋅,μ(ε))F(\cdot,\mu(\varepsilon))の流れ、Iε(y)I_\varepsilon(y)をyyの最大存在区間とする。このとき、次が成り立つ。

  1. ℓ1≠0\ell_1\ne0ならば、ε0\varepsilon_0を縮めると、0<∣ε∣<ε00<\lvert\varepsilon\rvert<\varepsilon_0のときμ(ε)≠0\mu(\varepsilon)\ne0であり、μ(ε)\mu(\varepsilon)の符号は−ℓ1ℜλ′(0)-\ell_1\Re\lambda'(0)の符号に等しい。
  2. ℓ1<0\ell_1<0であり、A0A_0の±iω0\pm i\omega_0以外の固有値の実部がすべて負であるとする。このとき、ε0\varepsilon_0を縮めると、0<∣ε∣<ε00<\lvert\varepsilon\rvert<\varepsilon_0について次の二つが成り立つ。Γε\Gamma_\varepsilonの任意の開近傍V′V'に対してΓε\Gamma_\varepsilonの開近傍V′′V''が存在して、y∈V′′y\in V''ならば[0,∞)⊆Iε(y)[0,\infty)\subseteq I_\varepsilon(y)かつφtε(y)∈V′\varphi^\varepsilon_t(y)\in V'(t≥0t\ge0)である。Γε\Gamma_\varepsilonの開近傍VVが存在して、y∈Vy\in Vならば[0,∞)⊆Iε(y)[0,\infty)\subseteq I_\varepsilon(y)であり、t→∞t\to\inftyのときdist⁡(φtε(y),Γε)→0\operatorname{dist}(\varphi^\varepsilon_t(y),\Gamma_\varepsilon)\to0である。
  3. ℓ1>0\ell_1>0ならば、ε0\varepsilon_0を縮めると、0<∣ε∣<ε00<\lvert\varepsilon\rvert<\varepsilon_0について、Γε\Gamma_\varepsilonの開近傍V′V'で次を満たすものが存在する。Γε\Gamma_\varepsilonのすべての開近傍V′′V''は、t∈Iε(y)t\in I_\varepsilon(y)、t≥0t\ge0、φtε(y)∉V′\varphi^\varepsilon_t(y)\notin V'を満たすttをもつ点yyを含む。

注意 4.3.定理 4.2の証明は本記事では与えない。A0A_0の±iω0\pm i\omega_0以外の固有値の実部がすべて00でない場合には、パラメータを変数に加えた系(x,μ)′=(F(x,μ),0)(x,\mu)'=(F(x,\mu),0)に(x0,0)(x_0,0)で§E13.21 定理 5.2を適用すると、同定理のEcE^cは三次元であり、得られる中心多様体の各μ\muでの切り口は二次元である。この場合の証明には、その二次元の切り口への制限と、制限した系の三次までの標準形の計算を用いる。定理 4.2 (2)だけがA0A_0の±iω0\pm i\omega_0以外の固有値の実部がすべて負であることを仮定し、定理 4.2 (1)と定理 4.2 (3)は定理 3.1の仮定とF∈C3(U×Λ,Rn)F\in C^3(U\times\Lambda,\R^n)だけを仮定するので、±iω0\pm i\omega_0以外にA0A_0の純虚数の固有値がある場合も含む。この場合には中心多様体の切り口は二次元でなく、二つの項の証明には、定理 3.1の証明の方程式H(ε;z,ω,μ)=0H(\varepsilon;z,\omega,\mu)=0の解(zε,ω(ε),μ(ε))(z_\varepsilon,\omega(\varepsilon),\mu(\varepsilon))のε\varepsilonについての二次までの展開と、周期軌道Γε\Gamma_\varepsilonの Floquet 乗数のうち11の近くにある自明でないもののε\varepsilonについての展開を用いる。

5 例

例 5.1.σ∈{1,−1}\sigma\in\{1,-1\}とし、U=R2U=\R^2、Λ=R\Lambda=\R、x=(x1,x2)x=(x_1,x_2)について

F(x,μ)=(μx1−x2−σx1(x12+x22)x1+μx2−σx2(x12+x22))F(x,\mu)=\begin{pmatrix}\mu x_1-x_2-\sigma x_1(x_1^2+x_2^2)\\x_1+\mu x_2-\sigma x_2(x_1^2+x_2^2)\end{pmatrix}

とする。R2\R^2をC\Cと同一視してζ=x1+ix2\zeta=x_1+ix_2と書くと、x′=F(x,μ)x'=F(x,\mu)はζ′=(μ+i)ζ−σ∣ζ∣2ζ\zeta'=(\mu+i)\zeta-\sigma\lvert\zeta\rvert^2\zetaである。

F(0,μ)=0F(0,\mu)=0であるからΛ0=R\Lambda_0=\R、ξ=0\xi=0とすることができ、A(μ)=(μ−11μ)A(\mu)=\begin{pmatrix}\mu&-1\\1&\mu\end{pmatrix}の固有値はμ±i\mu\pm iである。ω0=1\omega_0=1について定理 3.1 条件 (a)と定理 3.1 条件 (b)が成り立ち、§E13.22 補題 4.2 (3)により固有値の分枝はλ(μ)=μ+i\lambda(\mu)=\mu+iであるから、ℜλ′(0)=1\Re\lambda'(0)=1である。v=(1,−i)Tv=(1,-i)^{\mathsf T}はA0v=ivA_0v=ivを満たし、ϕ1(s)=(cos⁡s,sin⁡s)T\phi_1(s)=(\cos s,\sin s)^{\mathsf T}である。ε∈R\varepsilon\in\Rについて、ζ\zetaの表示でu(s)=εeisu(s)=\varepsilon e^{is}とおくと、u=εϕ1u=\varepsilon\phi_1であり、

(σε2+i)εeis−σε2⋅εeis=iεeis=u′(s)(\sigma\varepsilon^2+i)\varepsilon e^{is}-\sigma\varepsilon^2\cdot\varepsilon e^{is}=i\varepsilon e^{is}=u'(s)

であるから、(μ,ω,u)=(σε2,1,εϕ1)(\mu,\omega,u)=(\sigma\varepsilon^2,1,\varepsilon\phi_1)はωu′=F(u,μ)\omega u'=F(u,\mu)を満たす。u−εϕ1=0u-\varepsilon\phi_1=0であるから、0<∣ε∣<ε00<\lvert\varepsilon\rvert<\varepsilon_0、ε2≤ρ\varepsilon^2\le\rhoならば、定理 3.1 (3)によりμ(ε)=σε2\mu(\varepsilon)=\sigma\varepsilon^2、ω(ε)=1\omega(\varepsilon)=1、uε=εϕ1u_\varepsilon=\varepsilon\phi_1である。

ζ0∈C∖{0}\zeta_0\in\C\setminus\{0\}を初期値とする最大解をζ\zeta、その最大存在区間をIIとし、R(t)=∣ζ(t)∣2R(t)=\lvert\zeta(t)\rvert^2とおくと、R′=2ℜ(ζˉζ′)=2R(μ−σR)R'=2\Re(\bar\zeta\zeta')=2R(\mu-\sigma R)である。ζ=0\zeta=0は解であるから、一意性によりII上でR>0R>0である。h=R/∣ζ0∣h=\sqrt R/\lvert\zeta_0\rvertはh′=(μ−σR)hh'=(\mu-\sigma R)hを満たすから、η(t)=eith(t)ζ0\eta(t)=e^{it}h(t)\zeta_0は∣η∣2=R\lvert\eta\rvert^2=Rによりη′=iη+(μ−σR)η=(μ+i)η−σ∣η∣2η\eta'=i\eta+(\mu-\sigma R)\eta=(\mu+i)\eta-\sigma\lvert\eta\rvert^2\etaを満たし、η(0)=ζ0\eta(0)=\zeta_0である。一意性によりII上でζ(t)=eith(t)ζ0\zeta(t)=e^{it}h(t)\zeta_0である。RRは一次元の自律系R′=2R(μ−σR)R'=2R(\mu-\sigma R)の解であり、R′(t1)=0R'(t_1)=0を満たすt1t_1があればR(t1)R(t_1)は平衡点で、一意性によりRRは定数である。したがってRRは定数であるか狭義単調である。周期解の上でRRは周期的であるから定数であり、R>0R>0によりR=σμR=\sigma\muである。

σ=1\sigma=1とする。μ≤0\mu\le0ならば、R>0R>0でR′<0R'<0であるから、原点以外に周期解は存在しない。μ>0\mu>0ならば、円周Γμ={∣ζ∣=μ}\Gamma_\mu=\{\lvert\zeta\rvert=\sqrt\mu\}は唯一の周期軌道である。R=μR=\muはR′=2R(μ−R)R'=2R(\mu-R)の定常解であるから、R(0)<μR(0)<\muならばII上でR<μR<\muかつR′>0R'>0であり、R(0)>μR(0)>\muならばII上でR>μR>\muかつR′<0R'<0である。いずれの場合もt≥0t\ge0でR(t)R(t)はR(0)R(0)とμ\muの間にあるから、§E13.15 補題 3.1 (1)により[0,∞)⊆I[0,\infty)\subseteq Iであり、単調で有界なRRの極限R∞>0R_\infty>0は2R∞(μ−R∞)=02R_\infty(\mu-R_\infty)=0を満たすからR∞=μR_\infty=\muである。したがってt→∞t\to\inftyのときdist⁡(ζ(t),Γμ)=∣R(t)−μ∣→0\operatorname{dist}(\zeta(t),\Gamma_\mu)=\lvert\sqrt{R(t)}-\sqrt\mu\rvert\to0である。さらにθ∈R\theta\in\Rについてζ(2πj+θ)=eiθh(2πj+θ)ζ0→eiθμ ζ0/∣ζ0∣\zeta(2\pi j+\theta)=e^{i\theta}h(2\pi j+\theta)\zeta_0\to e^{i\theta}\sqrt\mu\,\zeta_0/\lvert\zeta_0\rvert(j→∞j\to\infty)であるから、ζ0∉Γμ\zeta_0\notin\Gamma_\muのω\omega極限集合はΓμ\Gamma_\muに等しく、Γμ\Gamma_\muは極限閉軌道である。

σ=−1\sigma=-1とする。μ≥0\mu\ge0ならばR>0R>0でR′>0R'>0であるから、原点以外に周期解は存在しない。μ<0\mu<0ならばΓμ={∣ζ∣=−μ}\Gamma_\mu=\{\lvert\zeta\rvert=\sqrt{-\mu}\}は唯一の周期軌道である。0<R(0)<−μ0<R(0)<-\muならば、t≥0t\ge0でRRは狭義単調減少で0<R<−μ0<R<-\muにとどまり、同じ議論によりR(t)→0R(t)\to0(t→∞t\to\infty)である。したがってΓμ\Gamma_\muのすべての開近傍は、前方軌道が原点へ近づく点を含む。

q=p=12(1,−i)Tq=p=\frac1{\sqrt2}(1,-i)^{\mathsf T}はA0q=iqA_0q=iq、q∗q=1q^*q=1、p∗A0=ip∗p^*A_0=ip^*、p∗q=1p^*q=1を満たす。b=0b=0であり、c(y1,y2,y3)=−2σ(y1(y2Ty3)+y2(y1Ty3)+y3(y1Ty2))c(y_1,y_2,y_3)=-2\sigma\bigl(y_1(y_2^{\mathsf T}y_3)+y_2(y_1^{\mathsf T}y_3)+y_3(y_1^{\mathsf T}y_2)\bigr)であるから、qTq=0q^{\mathsf T}q=0、qTqˉ=1q^{\mathsf T}\bar q=1によりc(q,q,qˉ)=−4σqc(q,q,\bar q)=-4\sigma q、ℓ1=−2σ\ell_1=-2\sigmaである。μ(ε)=σε2\mu(\varepsilon)=\sigma\varepsilon^2の符号と、σ=1\sigma=1の周期軌道の吸引とσ=−1\sigma=-1の周期軌道の不安定性は、定理 4.2の主張と一致する。

例 5.2.U=R2U=\R^2、Λ=R\Lambda=\R、x=(x1,x2)x=(x_1,x_2)について

F(x,μ)=(μx1−x2+x12x1+μx2+x12)F(x,\mu)=\begin{pmatrix}\mu x_1-x_2+x_1^2\\x_1+\mu x_2+x_1^2\end{pmatrix}

とする。このとき、FFはC∞C^\infty級でF(0,μ)=0F(0,\mu)=0であるから、Λ0=R\Lambda_0=\R、ξ=0\xi=0とすることができ、A(μ)=(μ−11μ)A(\mu)=\begin{pmatrix}\mu&-1\\1&\mu\end{pmatrix}は例 5.1のA(μ)A(\mu)と同じである。したがってω0=1\omega_0=1について定理 3.1 条件 (a)、定理 3.1 条件 (b)、定理 3.1 条件 (c)が成り立ち、ℜλ′(0)=1\Re\lambda'(0)=1である。q=p=12(1,−i)Tq=p=\frac1{\sqrt2}(1,-i)^{\mathsf T}は定理 4.2のq,pq,pの条件を満たす。

F(⋅,0)F(\cdot,0)の二階偏導関数のうち00でないものは∂x1∂x1F(0,0)=(2,2)T\partial_{x_1}\partial_{x_1}F(0,0)=(2,2)^{\mathsf T}だけであり、三階偏導関数はすべて00である。したがって、e1=(1,0)Te_1=(1,0)^{\mathsf T}とおくと、y,z∈C2y,z\in\C^2についてb(y,z)=2(e1Ty)(e1Tz)(1,1)Tb(y,z)=2(e_1^{\mathsf T}y)(e_1^{\mathsf T}z)(1,1)^{\mathsf T}であり、c=0c=0である。e1Tq=e1Tqˉ=1/2e_1^{\mathsf T}q=e_1^{\mathsf T}\bar q=1/\sqrt2、p∗(1,1)T=(1+i)/2p^*(1,1)^{\mathsf T}=(1+i)/\sqrt2であり、

A0−1=(01−10),(2iI2−A0)−1=−13(2i−112i)A_0^{-1}=\begin{pmatrix}0&1\\-1&0\end{pmatrix},\qquad(2iI_2-A_0)^{-1}=-\frac13\begin{pmatrix}2i&-1\\1&2i\end{pmatrix}

である。b(q,qˉ)=(1,1)Tb(q,\bar q)=(1,1)^{\mathsf T}、A0−1(1,1)T=(1,−1)TA_0^{-1}(1,1)^{\mathsf T}=(1,-1)^{\mathsf T}、b(q,(1,−1)T)=2(1,1)Tb\bigl(q,(1,-1)^{\mathsf T}\bigr)=\sqrt2(1,1)^{\mathsf T}であるから、p∗b(q,A0−1b(q,qˉ))=1+ip^*b\bigl(q,A_0^{-1}b(q,\bar q)\bigr)=1+iである。b(q,q)=(1,1)Tb(q,q)=(1,1)^{\mathsf T}、(2iI2−A0)−1(1,1)T=13(1−2i,−1−2i)T(2iI_2-A_0)^{-1}(1,1)^{\mathsf T}=\frac13(1-2i,-1-2i)^{\mathsf T}、b(qˉ,13(1−2i,−1−2i)T)=23(1−2i)(1,1)Tb\bigl(\bar q,\frac13(1-2i,-1-2i)^{\mathsf T}\bigr)=\frac{\sqrt2}3(1-2i)(1,1)^{\mathsf T}であるから、p∗b(qˉ,(2iI2−A0)−1b(q,q))=13(1−2i)(1+i)=1−i3p^*b\bigl(\bar q,(2iI_2-A_0)^{-1}b(q,q)\bigr)=\frac13(1-2i)(1+i)=1-\frac i3である。したがって

ℓ1=12ℜ(0−2(1+i)+1−i3)=−12\ell_1=\frac12\Re\Bigl(0-2(1+i)+1-\frac i3\Bigr)=-\frac12

である。

ℓ1<0\ell_1<0、ℜλ′(0)=1\Re\lambda'(0)=1であり、A0A_0の固有値は±i\pm iだけである。本記事で証明を与えていない定理 4.2を認めると、定理 4.2 (1)と定理 4.2 (2)により、ε0\varepsilon_0を縮めると、0<∣ε∣<ε00<\lvert\varepsilon\rvert<\varepsilon_0についてμ(ε)>0\mu(\varepsilon)>0であり、周期軌道Γε\Gamma_\varepsilonは定理 4.2 (2)の二つの性質をもつ。一方、命題 4.1により、δ>0\delta>0が存在して、0<μ<δ0<\mu<\deltaならば原点はF(⋅,μ)F(\cdot,\mu)の不安定な平衡点であり、−δ<μ<0-\delta<\mu<0ならば漸近安定な平衡点である。したがって、定理 4.2を認めると、この例の周期軌道は、原点が不安定である側のμ\muに現れる。

例 5.3.σ∈{1,−1}\sigma\in\{1,-1\}とし、U=R2U=\R^2、Λ=R\Lambda=\Rとする。R2\R^2をC\Cと同一視し、ζ=x1+ix2\zeta=x_1+ix_2についてFFをζ′=(μ2+i)ζ+σ∣ζ∣2ζ\zeta'=(\mu^2+i)\zeta+\sigma\lvert\zeta\rvert^2\zetaがx′=F(x,μ)x'=F(x,\mu)を表すように定める。このとき、F(0,μ)=0F(0,\mu)=0であり、A(μ)=DxF(0,μ)=(μ2−11μ2)A(\mu)=D_xF(0,\mu)=\begin{pmatrix}\mu^2&-1\\1&\mu^2\end{pmatrix}の固有値はμ2±i\mu^2\pm iである。ω0=1\omega_0=1について定理 3.1 条件 (a)と定理 3.1 条件 (b)が成り立つが、固有値の分枝はλ(μ)=μ2+i\lambda(\mu)=\mu^2+iであり、ℜλ′(0)=0\Re\lambda'(0)=0であるから定理 3.1 条件 (c)は成り立たない。R=∣ζ∣2R=\lvert\zeta\rvert^2はR′=2R(μ2+σR)R'=2R(\mu^2+\sigma R)を満たす。

σ=1\sigma=1とする。すべてのμ∈R\mu\in\Rについて、ζ(0)≠0\zeta(0)\ne0の解のRRは狭義単調増加であり、周期解は存在しない。したがって、定理 3.1 (2)の結論は、μ(ε)\mu(\varepsilon)をどのように選んでも成り立たない。

σ=−1\sigma=-1とする。μ≠0\mu\ne0について円周{∣ζ∣=∣μ∣}\{\lvert\zeta\rvert=\lvert\mu\rvert\}は周期軌道である。ε∈R\varepsilon\in\Rについてμ(ε)=ε\mu(\varepsilon)=\varepsilon、ω(ε)=1\omega(\varepsilon)=1、uε=εϕ1u_\varepsilon=\varepsilon\phi_1(ϕ1(s)=(cos⁡s,sin⁡s)T\phi_1(s)=(\cos s,\sin s)^{\mathsf T})とおくと、ζ\zetaの表示で(ε2+i)εeis−ε2⋅εeis=iεeis(\varepsilon^2+i)\varepsilon e^{is}-\varepsilon^2\cdot\varepsilon e^{is}=i\varepsilon e^{is}であるから、ω(ε)uε′=F(uε,μ(ε))\omega(\varepsilon)u_\varepsilon'=F(u_\varepsilon,\mu(\varepsilon))であり、ε≠0\varepsilon\ne0についてxε(t)=uε(t)x_\varepsilon(t)=u_\varepsilon(t)はx′=F(x,ε)x'=F(x,\varepsilon)の最小周期2π2\piの周期解である。この族は定理 3.1 (2)の結論を満たす。定理 3.1 条件 (a)と定理 3.1 条件 (b)の下で、定理 3.1 条件 (c)は周期解の族が存在するための十分条件であるが、必要条件ではない。この例では、族の存否は三次の項の符号σ\sigmaで決まる。

6 演習

問題 6.1.定理 3.1のε0\varepsilon_0、ρ\rhoと(μ(ε),ω(ε),uε)(\mu(\varepsilon),\omega(\varepsilon),u_\varepsilon)について、0<∣ε∣<ε00<\lvert\varepsilon\rvert<\varepsilon_0ならばμ(−ε)=μ(ε)\mu(-\varepsilon)=\mu(\varepsilon)、ω(−ε)=ω(ε)\omega(-\varepsilon)=\omega(\varepsilon)であり、すべてのs∈Rs\in\Rでu−ε(s)=uε(s+π)u_{-\varepsilon}(s)=u_\varepsilon(s+\pi)であることを証明せよ。さらに、定理 3.1 (2)のx−ε(0)x_{-\varepsilon}(0)とxε(0)x_\varepsilon(0)がF(⋅,μ(ε))F(\cdot,\mu(\varepsilon))の同じ周期軌道に属することを証明せよ。

解答.

0<∣ε∣<ε00<\lvert\varepsilon\rvert<\varepsilon_0とし、u~(s)=uε(s+π)\tilde u(s)=u_\varepsilon(s+\pi)とおく。ei(s+π)=−eise^{i(s+\pi)}=-e^{is}であるからϕj(s+π)=−ϕj(s)\phi_j(s+\pi)=-\phi_j(s)(j=1,2j=1,2)であり、u~(s)−(−ε)ϕ1(s)=(uε−εϕ1)(s+π)\tilde u(s)-(-\varepsilon)\phi_1(s)=(u_\varepsilon-\varepsilon\phi_1)(s+\pi)である。したがってu~∈C2π1\tilde u\in C^1_{2\pi}、∥u~−(−ε)ϕ1∥C1=∥uε−εϕ1∥C1≤ρ∣−ε∣\lVert\tilde u-(-\varepsilon)\phi_1\rVert_{C^1}=\lVert u_\varepsilon-\varepsilon\phi_1\rVert_{C^1}\le\rho\lvert-\varepsilon\rvertであり、被積分関数の周期性により

∫02πϕj(s)T(u~(s)+εϕ1(s)) ds=−∫02πϕj(s+π)T(uε−εϕ1)(s+π) ds=−∫02πϕjT(uε−εϕ1) ds=0\int_0^{2\pi}\phi_j(s)^{\mathsf T}\bigl(\tilde u(s)+\varepsilon\phi_1(s)\bigr)\,ds=-\int_0^{2\pi}\phi_j(s+\pi)^{\mathsf T}(u_\varepsilon-\varepsilon\phi_1)(s+\pi)\,ds=-\int_0^{2\pi}\phi_j^{\mathsf T}(u_\varepsilon-\varepsilon\phi_1)\,ds=0

である。定理 3.1 (1)の等式をs+πs+\piで用いると、すべてのssでξ(μ(ε))+u~(s)∈U\xi(\mu(\varepsilon))+\tilde u(s)\in Uかつω(ε)u~′(s)=F(ξ(μ(ε))+u~(s),μ(ε))\omega(\varepsilon)\tilde u'(s)=F(\xi(\mu(\varepsilon))+\tilde u(s),\mu(\varepsilon))である。∣μ(ε)∣≤ρ\lvert\mu(\varepsilon)\rvert\le\rho、∣ω(ε)−ω0∣≤ρ\lvert\omega(\varepsilon)-\omega_0\rvert\le\rhoであるから、(μ(ε),ω(ε),u~)(\mu(\varepsilon),\omega(\varepsilon),\tilde u)は−ε-\varepsilonについて定理 3.1 (3)の仮定を満たし、(μ(ε),ω(ε),u~)=(μ(−ε),ω(−ε),u−ε)(\mu(\varepsilon),\omega(\varepsilon),\tilde u)=(\mu(-\varepsilon),\omega(-\varepsilon),u_{-\varepsilon})である。

φ\varphiをF(⋅,μ(ε))F(\cdot,\mu(\varepsilon))の流れとする。前段により、すべてのttでx−ε(t)=ξ(μ(ε))+uε(ω(ε)t+π)=xε(t+π/ω(ε))x_{-\varepsilon}(t)=\xi(\mu(\varepsilon))+u_\varepsilon(\omega(\varepsilon)t+\pi)=x_\varepsilon\bigl(t+\pi/\omega(\varepsilon)\bigr)である。xεx_\varepsilonはR\R上の古典解であるから、§E13.15 補題 1.2によりxε(t)=φt(xε(0))x_\varepsilon(t)=\varphi_t(x_\varepsilon(0))であり、x−ε(0)=φπ/ω(ε)(xε(0))∈γ(xε(0))x_{-\varepsilon}(0)=\varphi_{\pi/\omega(\varepsilon)}(x_\varepsilon(0))\in\gamma(x_\varepsilon(0))である。§E13.15 命題 4.3 (1)によりγ(x−ε(0))=γ(xε(0))\gamma(x_{-\varepsilon}(0))=\gamma(x_\varepsilon(0))である。▨

問題 6.2.J=(0−110)J=\begin{pmatrix}0&-1\\1&0\end{pmatrix}、A0=(J002J)∈R4×4A_0=\begin{pmatrix}J&0\\0&2J\end{pmatrix}\in\R^{4\times4}、ω0=1\omega_0=1とする。±i\pm iがA0A_0の代数的に単純な固有値であり、2i2iがA0A_0の固有値であることを確かめよ。v=(1,−i,0,0)Tv=(1,-i,0,0)^{\mathsf T}とし、C2πC_{2\pi}、C2π1C^1_{2\pi}、ϕ1,ϕ2\phi_1,\phi_2、L0L_0、X1X_1を補題 2.1と同じ式で定めるとき、ker⁡L0\ker L_0が四次元であること、および、すべてのB∈R4×4B\in\R^{4\times4}についてX1×R×R∋(u,a,b)↦L0u+aϕ1′−bBϕ1∈C2πX_1\times\R\times\R\ni(u,a,b)\mapsto L_0u+a\phi_1'-bB\phi_1\in C_{2\pi}が単射でないことを証明せよ。

解答.

det⁡(zI2−J)=z2+1\det(zI_2-J)=z^2+1、det⁡(zI2−2J)=z2+4\det(zI_2-2J)=z^2+4であるから、A0A_0の特性多項式は(z2+1)(z2+4)(z^2+1)(z^2+4)であり、根±i\pm i、±2i\pm2iはすべて単根である。したがって±i\pm iは代数的に単純な固有値であり、2i2iは固有値である。

ϕ1(s)=ℜ(eisv)=(cos⁡s,sin⁡s,0,0)T\phi_1(s)=\Re(e^{is}v)=(\cos s,\sin s,0,0)^{\mathsf T}、ϕ2(s)=(sin⁡s,−cos⁡s,0,0)T\phi_2(s)=(\sin s,-\cos s,0,0)^{\mathsf T}である。e1=ϕ1e_1=\phi_1、e2=ϕ2e_2=\phi_2、e3(s)=(0,0,cos⁡2s,sin⁡2s)Te_3(s)=(0,0,\cos2s,\sin2s)^{\mathsf T}、e4(s)=(0,0,sin⁡2s,−cos⁡2s)Te_4(s)=(0,0,\sin2s,-\cos2s)^{\mathsf T}とおくと

e1′=(−sin⁡s,cos⁡s,0,0)T=A0e1,e2′=(cos⁡s,sin⁡s,0,0)T=A0e2,e3′=(0,0,−2sin⁡2s,2cos⁡2s)T=A0e3,e4′=(0,0,2cos⁡2s,2sin⁡2s)T=A0e4\begin{aligned} e_1'&=(-\sin s,\cos s,0,0)^{\mathsf T}=A_0e_1,&\quad e_2'&=(\cos s,\sin s,0,0)^{\mathsf T}=A_0e_2,\\ e_3'&=(0,0,-2\sin2s,2\cos2s)^{\mathsf T}=A_0e_3,&\quad e_4'&=(0,0,2\cos2s,2\sin2s)^{\mathsf T}=A_0e_4 \end{aligned}

であるから、e1,…,e4∈ker⁡L0e_1,\ldots,e_4\in\ker L_0である。e1(0)=(1,0,0,0)Te_1(0)=(1,0,0,0)^{\mathsf T}、e2(0)=(0,−1,0,0)Te_2(0)=(0,-1,0,0)^{\mathsf T}、e3(0)=(0,0,1,0)Te_3(0)=(0,0,1,0)^{\mathsf T}、e4(0)=(0,0,0,−1)Te_4(0)=(0,0,0,-1)^{\mathsf T}は一次独立であるから、e1,…,e4e_1,\ldots,e_4は一次独立である。§E13.13 命題 4.7をA0A_0、T=2πT=2\pi、ω=1\omega=1に適用すると、ker⁡(A0−inI)≠{0}\ker(A_0-inI)\ne\{0\}となるnnは±1,±2\pm1,\pm2であり、各核は§E13.22 補題 4.1 (1)により一次元であるから、§E13.13 命題 4.7 (1)と§E13.13 命題 4.7 (3)によりu′=A0uu'=A_0uの2π2\pi周期の実数値の局所絶対連続解の全体は四次元である。ker⁡L0\ker L_0はその部分空間であり、四つの一次独立な元を含むから、ker⁡L0\ker L_0は四次元である。

e3e_3とϕ1,ϕ2\phi_1,\phi_2は各ssで直交するから∫02πϕjTe3 ds=0\int_0^{2\pi}\phi_j^{\mathsf T}e_3\,ds=0であり、e3∈X1e_3\in X_1である。B∈R4×4B\in\R^{4\times4}について、(e3,0,0)≠0(e_3,0,0)\ne0はL0e3+0⋅ϕ1′−0⋅Bϕ1=0L_0e_3+0\cdot\phi_1'-0\cdot B\phi_1=0を満たすから、写像は単射でない。補題 2.1の仮定のうち、k∈Z∖{1,−1}k\in\Z\setminus\{1,-1\}についてikω0ik\omega_0が固有値でないという条件を外すと、補題 2.1 (5)の結論が成り立たないA0A_0が存在する。▨

参考文献

  1. Lawrence Perko, Differential Equations and Dynamical Systems, 3rd ed., Texts in Applied Mathematics 7, Springer, New York, 2001.Hopf 分岐定理の仮定と第一 Lyapunov 係数による分岐の向きと安定性を参考にした。
  2. Gerald Teschl, Ordinary Differential Equations and Dynamical Systems, American Mathematical Society, Providence, Rhode Island, 2012.周期解の分岐と標準形の例の扱いを参考にした。
  3. Yuri A. Kuznetsov, Elements of Applied Bifurcation Theory, 3rd ed., Applied Mathematical Sciences 112, Springer, New York, 2004.第一 Lyapunov 係数の n 次元の公式を参考にした。

前提記事

15 本の記事・単元を表示