フロケ理論

数学のフロケ理論（フロケりろん、英: Floquet theory）とは、次の形の線型微分方程式の解のクラスに関する常微分方程式理論の一分野である。 x ˙ = A ( t ) x , {\displaystyle {\dot {x}}=A(t)x,\,} ここで A ( t ) {\displaystyle \displaystyle A(t)} は区分的連続な周期 T {\displaystyle T} の周期関数である。 フロケ理論における主定理であるフロケの定理（Floquet's theorem）は、Gaston Floquet (1883) によるもので、この共通の線型系の各基本解行列に対する標準形（英語版）を与えるものである。それはまた、 Q ( t + 2 T ) = Q ( t ) {\displaystyle \displaystyle Q(t+2T)=Q(t)} を満たすような座標変換 y = Q − 1 ( t ) x {\displaystyle \displaystyle y=Q^{-1}(t)x} を与え、これは周期系を典型的な実係数の線型系へと変換する。 固体物理学において、（3次元へと一般化された）同様の結果はブロッホの定理として知られている。 線型微分方程式の解はベクトル空間を構成することに注意されたい。ある行列 ϕ ( t ) {\displaystyle \phi \,(t)} が基本解行列（fundamental matrix solution）であるとは、その全ての列が線型独立な解であることを言う。ある行列 Φ ( t ) {\displaystyle \Phi (t)} が主基本解行列（principal fundamental matrix solution）であるとは、その全ての列が線型独立な解で、 Φ ( t 0 ) {\displaystyle \Phi (t_{0})} が単位行列となるようなある t 0 {\displaystyle t_{0}} が存在することを言う。主基本行列は、 Φ ( t ) = ϕ ( t ) ϕ − 1 ( t 0 ) {\displaystyle \Phi (t)=\phi \,(t){\phi \,}^{-1}(t_{0})} を使うことで基本行列から構成することが出来る。初期条件が x ( 0 ) = x 0 {\displaystyle x(0)=x_{0}} であるような線型微分方程式の解は、 x ( t ) = ϕ ( t ) ϕ − 1 ( 0 ) x 0 {\displaystyle x(t)=\phi \,(t){\phi \,}^{-1}(0)x_{0}} である。ここで ϕ ( t ) {\displaystyle \phi \,(t)} は任意の基本行列である。 x ˙ = A ( t ) x {\displaystyle {\dot {x}}=A(t)x} を一階の線型微分方程式とし、 x ( t ) {\displaystyle x(t)} は長さ n {\displaystyle n} の列ベクトルとし、 A ( t ) {\displaystyle A(t)} は周期 T {\displaystyle T} の n × n {\displaystyle n\times n} 周期行列とする（すなわち、すべての実数値 t {\displaystyle t} に対して A ( t + T ) = A ( t ) {\displaystyle A(t+T)=A(t)} が成立する）。 ϕ ( t ) {\displaystyle \phi \,(t)} をこの微分方程式のある基本解行列とする。このとき、全ての t ∈ R {\displaystyle t\in \mathbb {R} } に対して ϕ ( t + T ) = ϕ ( t ) ϕ − 1 ( 0 ) ϕ ( T )   {\displaystyle \phi (t+T)=\phi (t)\phi ^{-1}(0)\phi (T)\ } が成立する。ここで ϕ − 1 ( 0 ) ϕ ( T )   {\displaystyle \phi ^{-1}(0)\phi (T)\ } はモノドロミー行列として知られるものである。さらに（複素数値でもあり得る）各行列 B {\displaystyle B} で e T B = ϕ − 1 ( 0 ) ϕ ( T ) ,   {\displaystyle e^{TB}=\phi ^{-1}(0)\phi (T),\ } を満たすようなものに対し、周期 T {\displaystyle T} の周期行列関数 t ↦ P ( t ) {\displaystyle t\mapsto P(t)} で ϕ ( t ) = P ( t ) e t B  for all  t ∈ R   {\displaystyle \phi (t)=P(t)e^{tB}{\text{ for all }}t\in \mathbb {R} \ } を満たすようなものが存在する。また、ある実行列 R {\displaystyle R} と実周期（ 2 T {\displaystyle 2T} -周期）行列函数 t ↦ Q ( t ) {\displaystyle t\mapsto Q(t)} で ϕ ( t ) = Q ( t ) e t R  for all  t ∈ R   {\displaystyle \phi (t)=Q(t)e^{tR}{\text{ for all }}t\in \mathbb {R} \ } を満たすようなものが存在する。以上の議論において、 B {\displaystyle B} 、 P {\displaystyle P} 、 Q {\displaystyle Q} および R {\displaystyle R} は n × n {\displaystyle n\times n} 行列である。 写像 ϕ ( t ) = Q ( t ) e t R {\displaystyle \phi \,(t)=Q(t)e^{tR}} は時間依存の座標変換（ y = Q − 1 ( t ) x {\displaystyle y=Q^{-1}(t)x} ）を導き、その変換の下で元の方程式系は実定数を係数とする線型系 y ˙ = R y {\displaystyle {\dot {y}}=Ry} になる。 Q ( t ) {\displaystyle Q(t)} は連続かつ周期的であるため、有界である。したがって y ( t ) {\displaystyle y(t)} および x ( t ) {\displaystyle x(t)} に対するゼロ解の安定性は R {\displaystyle R} の固有値によって決定される。 表現 ϕ ( t ) = P ( t ) e t B {\displaystyle \phi \,(t)=P(t)e^{tB}} は、基本行列 ϕ ( t ) {\displaystyle \phi \,(t)} に対する「フロケ正規形（Floquet normal form）」と呼ばれる。 e T B {\displaystyle e^{TB}} の固有値は、その方程式系の特性乗数と呼ばれる。それらはまた（線型の）ポアンカレ写像 x ( t ) → x ( t + T ) {\displaystyle x(t)\to x(t+T)} の固有値でもある。フロケ指数（Floquet exponent、しばしば特性指数とも呼ばれる）とは、 e μ T {\displaystyle e^{\mu T}} がその方程式系の特性乗数となるような複素数 μ {\displaystyle \mu } のことを言う。ここで、整数 k {\displaystyle k} に対して e ( μ + 2 π i k T ) T = e μ T {\displaystyle e^{(\mu +{\frac {2\pi ik}{T}})T}=e^{\mu T}} が成立することより、フロケ指数は一意的ではないことに注意されたい。フロケ指数の実部は、リアプノフ指数（Lyapunov exponents）と呼ばれる。すべてのリアプノフ指数が負であればゼロ解は漸近安定となり、非正であればリアプノフ安定となり、それ以外の場合では不安定となる。 フロケ理論は力学系の研究において非常に重要となる。 フロケ理論は、周期的な重力場における調和振動子として月の運動を近似するヒルの方程式（ジョージ・ウィリアム・ヒルによって考えられた）において、安定性を示すものである。 強レーザー場におけるボンドソフト化（英語版）およびボンドハード化（英語版）は、フロケの定理により得られる解によって表現される。 マシュー方程式とは、楕円柱に対する波動方程式に関連するものである。 与えられた a ∈ R , q ∈ C {\displaystyle a\in \mathbb {R} ,q\in \mathbb {C} } に対するマシュー方程式とは、次のようなものである。 d 2 y d w 2 + ( a − 2 q cos ⁡ 2 w ) y = 0. {\displaystyle {\frac {d^{2}y}{dw^{2}}}+(a-2q\cos 2w)y=0.} マシュー方程式は、周期係数を持つ二階の線型微分方程式である。 マシュー函数に関する最も有用な結果は、フロケの定理 [1, 2] である。その定理によれば、任意のペア (a, q) に対するマシュー方程式の解は次の形式で表現される。 y ( w ) = F ν ( w ) = e i w ν P ( w ) {\displaystyle y(w)=F_{\nu }(w)=e^{iw\nu }P(w)\,} または y ( w ) = F ν ( − w ) = e − i w ν P ( − w ) . {\displaystyle y(w)=F_{\nu }(-w)=e^{-iw\nu }P(-w).\,} ここで ν {\displaystyle \nu } は a および q に依存する定数であり、P(.) は w について π {\displaystyle \pi } -周期的である。 この定数 ν {\displaystyle \nu } は特性指数（characteristic exponent）と呼ばれる。 ν {\displaystyle \nu } が整数であるなら、 F ν ( w ) {\displaystyle F_{\nu }(w)} および F ν ( − w ) {\displaystyle F_{\nu }(-w)} は線型独立な解である。さらに、 y ( w + k π ) = e i ν k π y ( w )  or  y ( w + k π ) = e − i ν k π y ( w ) , {\displaystyle y(w+k\pi )=e^{i\nu k\pi }y(w){\text{ or }}y(w+k\pi )=e^{-i\nu k\pi }y(w),\,} がそれぞれ解 F ν ( w ) {\displaystyle F_{\nu }(w)} or F ν ( − w ) {\displaystyle F_{\nu }(-w)} に対して成立する。 ペア (a, q) は | cosh ⁡ ( i ν π ) | < 1 {\displaystyle |\cosh(i\nu \pi )|<1} を満たすようなもので、したがって解 y ( w ) {\displaystyle y(w)} は実軸上で有界であるものと仮定する。マシュー方程式の一般解は、（ q ∈ R {\displaystyle q\in \mathbb {R} } および ν {\displaystyle \nu } は非整数）次のような形式で記述される。 y ( w ) = c 1 e i w ν P ( w ) + c 2 e − i w ν P ( − w ) , {\displaystyle y(w)=c_{1}e^{iw\nu }P(w)+c_{2}e^{-iw\nu }P(-w),\,} ここで c 1 {\displaystyle c_{1}} および c 2 {\displaystyle c_{2}} は任意定数である。 分数次あるいは整数次のすべての有界な解は、振動数が増加するにつれて振幅が減少するような調和振動の無限級数として表現される。 マシュー方程式の、他の非常に重要な性質として、直交性（orthogonality）が挙げられる [3]： a ( ν + 2 p , q ) {\displaystyle a(\nu +2p,q)} および a ( ν + 2 s , q ) {\displaystyle a(\nu +2s,q)} が cos ⁡ ( π ν ) − y ( π = 0 ) = 0 , {\displaystyle \cos(\pi \nu )-y(\pi =0)=0,\,} の単根であるなら、 ∫ 0 π F ν + 2 p ( w ) F ν + 2 s ( − w ) d w = 0 , p ≠ s , {\displaystyle \int _{0}^{\pi }F_{\nu +2p}(w)F_{\nu +2s}(-w)\,dw=0,\qquad p\neq s,} が成立する。すなわち、 ⟨ F ν + 2 p ( w ) , F ν + 2 s ( w ) ⟩ = 0 , p ≠ s , {\displaystyle \langle F_{\nu +2p}(w),F_{\nu +2s}(w)\rangle =0,\qquad p\neq s,} が成立する。ここで <·,·> は 0 から π に対して定義される内積である。 C. Chicone. Ordinary Differential Equations with Applications. Springer-Verlag, New York 1999. Ekeland, Ivar (1990). “One”. Convexity methods in Hamiltonian mechanics. Ergebnisse der Mathematik und ihrer Grenzgebiete (3) [Results in Mathematics and Related Areas (3)]. 19. Berlin: Springer-Verlag. pp. x+247. ISBN 3-540-50613-6. MR1051888  Floquet, Gaston (1883), “Sur les équations différentielles linéaires à coefficients périodiques”, Annales de l'École Normale Supérieure 12: 47–88, http://archive.numdam.org/ARCHIVE/ASENS/ASENS_1883_2_12_/ASENS_1883_2_12__47_0/ASENS_1883_2_12__47_0.pdf  Krasnosel'skii, M.A. (1968), The Operator of Translation along the Trajectories of Differential Equations, Providence: American Mathematical Society , Translation of Mathematical Monographs, 19, 294p. W. Magnus, S. Winkler. Hill's Equation, Dover-Phoenix Editions, ISBN 0-486-49565-5. N.W. McLachlan, Theory and Application of Mathieu Functions, New York: Dover, 1964. Teschl, Gerald (2012). Ordinary Differential Equations and Dynamical Systems. Providence: American Mathematical Society. ISBN 978-0-8218-8328-0. http://www.mat.univie.ac.at/~gerald/ftp/book-ode/  M.S.P. Eastham, "The Spectral Theory of Periodic Differential Equations", Texts in Mathematics, Scottish Academic Press, Edinburgh, 1973. ISBN 978-0-7011-1936-2.