フロケ理論

周期的な駆動力を受ける振り子や交流電流で駆動される振動回路など、力が周期的なシステムでは、システム全体の挙動は必ずしも完全に周期的であるとは限りません。例えば、ブランコで押されている子供を考えてみましょう。動きは規則的で周期的な押圧によって駆動されていますが、ブランコは前後に振動しながら徐々に高い位置まで上昇します。これにより、根底にある周期性と成長が組み合わさった状態が生まれます。

フロケ理論は、このようなシステムを解析する方法を提供します。その本質的な洞察はブランコの例と同様です。つまり、解は2つの部分に分解できます。周期成分(繰り返し動作を反映)と指数因子(成長、減衰、または中立安定性を反映)です。この分解により、時間周期システムにおける長期的な挙動と安定性の解析が可能になります。

正式には、フロケ理論は、次の形式の 周期線形微分方程式の解のクラスに関連する常微分方程式の一分野である。

×˙t×{\displaystyle {\dot {x}}=A(t)x,}

は周期を持つ周期関数であり、解の安定性の状態を定義します。 ×Rn{\displaystyle x\in {R^{n}}}tRn×n{\displaystyle \displaystyle A(t)\in {R^{n\times n}}}T{\displaystyle T}

フロケ理論の主要定理であるフロケの定理は、ガストン・フロケ (1883)によって提唱され、この共通線形系の各基本行列解の標準形を与える。この定理は、 による座標変換を与え、周期系を定数の実係数を持つ従来の線形系に変換する。この定理は、凝縮系物理学における結晶など、周期ポテンシャルを持つ物理系に適用され、ブロッホの定理として知られる。 y質問1t×{\displaystyle \displaystyle y=Q^{-1}(t)x}質問t+2T質問t{\displaystyle \displaystyle Q(t+2T)=Q(t)}

線型微分方程式の解はベクトル空間を形成することに注意してください。行列は、列が解集合の基底を形成する場合、基本行列解と呼ばれます。すべての列が線型独立な解であり、が恒等解となるような行列が存在する場合、行列は主基本行列解と呼ばれます。主基本行列は、 を用いて基本行列から構築できます。初期条件 を持つ線型微分方程式の解はであり、は任意の基本行列解です。 ϕt{\displaystyle \phi \,(t)}Φt{\displaystyle \Phi (t)}t0{\displaystyle t_{0}}Φt0{\displaystyle \Phi (t_{0})}Φtϕtϕ1t0{\displaystyle \Phi (t)=\phi \,(t){\phi \,}^{-1}(t_{0})}×0×0{\displaystyle x(0)=x_{0}}×tϕtϕ10×0{\displaystyle x(t)=\phi \,(t){\phi \,}^{-1}(0)x_{0}}ϕt{\displaystyle \phi \,(t)}

フロケの定理

を線型一階微分方程式とする。ここで は長さ の列ベクトルであり、周期行列( のすべての実数に対して となる)である。をこの微分方程式の基本行列解とする。すると、すべての に対して、 ×˙t×{\displaystyle {\dot {x}}=A(t)x}×t{\displaystyle x(t)}n{\displaystyle n}t{\displaystyle A(t)}n×n{\displaystyle n\times n}T{\displaystyle T}A(t+T)=A(t){\displaystyle A(t+T)=A(t)}t{\displaystyle t}ϕ(t){\displaystyle \phi \,(t)}tR{\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}

eTB=ϕ1(0)ϕ(T),{\displaystyle e^{TB}=\phi ^{-1}(0)\phi (T),}

周期行列関数(周期)が存在し、 T{\displaystyle T}tP(t){\displaystyle t\mapsto P(t)}

ϕ(t)=P(t)etB for all tR.{\displaystyle \phi (t)=P(t)e^{tB}{\text{ for all }}t\in \mathbb {R} .}

この表現は、基本行列解のフロケ正規形と呼ばれます。 ϕ(t){\displaystyle \phi \,(t)}

さらに、が の任意の値に対して実数行列である場合、を満たす 実数行列も少なくとも1つ存在する。 ϕ(t){\displaystyle \phi (t)}t{\displaystyle t}R{\displaystyle R}

ϕ1(0)ϕ(2T)=(ϕ1(0)ϕ(T))2=e2TR.{\displaystyle \phi ^{-1}(0)\phi (2T)=(\phi ^{-1}(0)\phi (T))^{2}=e^{2TR}.}

の任意の選択に対して、周期 を持つ周期行列関数が存在し、 R{\displaystyle R}2T{\displaystyle 2T}tQ(t){\displaystyle t\mapsto Q(t)}

ϕ(t)=Q(t)etR for all tR.{\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)etR{\displaystyle \phi \,(t)=Q(t)e^{tR}}y=Q1(t)x{\displaystyle y=Q^{-1}(t)x}y˙=Ry{\displaystyle {\dot {y}}=Ry}Q(t){\displaystyle Q(t)}y(t){\displaystyle y(t)}x(t){\displaystyle x(t)}R{\displaystyle R}

の固有値、システムの特性乗数と呼ばれます。行列は一意ではありませんが、 の固有値はの選択に対して同じです。これらは、(線形)ポアンカレ写像の固有値でもあります。フロケ指数(特性指数と呼ばれることもあります)は、がシステムの特性乗数となるような複素数です。 (ただし は整数)であるため、フロケ指数は一意ではないことに注意してください。フロケ指数の実部は、リアプノフ指数と呼ばれます。すべてのリアプノフ指数が負の場合、零解は漸近的に安定であり、リアプノフ指数が非正の場合、リアプノフ安定であり、それ以外の場合は不安定です。 eTB{\displaystyle e^{TB}}B{\displaystyle B}eTB{\displaystyle e^{TB}}B{\displaystyle B}x(t)x(t+T){\displaystyle x(t)\to x(t+T)}μ{\displaystyle \mu }eμT{\displaystyle e^{\mu T}}e(μ+2πikT)T=eμT{\displaystyle e^{(\mu +{\frac {2\pi ik}{T}})T}=e^{\mu T}}k{\displaystyle k}

参考文献

  1. ^ arXiv:1211.5623