logo

力学におけるリュービルの定理の証明 📂動力学

力学におけるリュービルの定理の証明

定理

ユークリッド空間 Rn\mathbb{R}^{n} と関数 f:RnRnf : \mathbb{R}^{n} \to \mathbb{R}^{n} について次のようなベクトル場微分方程式で与えられているとする。 x˙=f(x) \dot{x} = f(x) このシステムのフロー ϕt()\phi_t ( \cdot ) と領域 D0RnD_{0} \subset \mathbb{R}^{n} について Dt:=ϕt(D0)D_{t} := \phi_{t} \left( D_{0} \right) をフローに従って時間 tt が経過して移動した領域、そのボリュームを V(t)V(Dt)V(t) \equiv V \left( D_{t} \right) として示す。もし f=0\nabla \cdot f = 0 ならすべての D0RnD_{0} \subset \mathbb{R}^{n}tRt \in \mathbb{R} について次が成り立つ。 V(Dt)=V(D0) V \left( D_{t} \right) = V \left( D_{0} \right)

数式の説明

f\nabla \cdot fベクトル場の発散で、ベクトル場が広がったり集まったりする様子を示す。

VVベクトル場で与えられた領域のボリュームを意味する。リウヴィルの定理は言葉より数式の方が簡単だが、簡単に言えば、発散がどこでも 00 ならばフローが移動した領域のボリュームが変わらないということである。

証明 1

戦略: ベクトル微積分学を少し使えばよい。次のようにベクトル関数と各ベクトルの軸を表せば、解析学と線形代数の基本的なツールによる計算だけで直接演繹される。 f:=(f1,,fn)x:=(x1,,xn) f := \left( f_{1} , \cdots , f_{n} \right) \\ x := \left( x_{1} , \cdots , x_{n} \right)


Part 1. t0=0t_{0} = 0 のとき dVdtt=t0=D0fdx\displaystyle \left.{{ d V } \over { d t }}\right|_{t = t_{0}} = \int_{D_{0}} \nabla \cdot f dx

ボリュームの定義: uRn\textbf{u} \in \mathbb{R}^{n} がベクトル関数 f:RnRn\textbf{f} : \mathbb{R}^{n} \to \mathbb{R}^{n} によって f(u)=(f1(u),,fn(u))\textbf{f} \left( \textbf{u} \right) = \left( f_{1} (\textbf{u}) , \cdots , f_{n} (\textbf{u}) \right) として変換されるとき、DD のボリュームは次のようになる。 V(D)=Df(u)udu1du2dun V(D) = \int_{D} \left| {{ \partial \textbf{f} (\textbf{u}) } \over { \partial \textbf{u} }} \right| d u_{1} d u_{2} \cdots d u_{n}

ボリュームの定義に従って V(t)=D0detϕt(x)xdx V(t) = \int_{D_{0}} \det {{ \partial \phi_{t} (x) } \over { \partial x }} dx フロー ϕt(x)\phi_{t} (x)t=0t=0 でテイラー展開すると ϕt(x)=x+f(x)t+O(t2) \phi_{t} (x) = x + f(x) t + O \left( t^{2} \right) xx偏微分すると n×nn\times n 単位行列 EE に対して ϕt(x)x=E+fxt+O(t2) {{ \partial \phi_{t} (x) } \over { \partial x }} = E + {{ \partial f } \over { \partial x }} t + O \left( t^{2} \right) のように示すことができる。すると、その行列式detϕt(x)x=det[E+fxt]+O(t2)=1+tr[fx]t+O(t2)=1+(df1dx1++dfndxn)t+O(t2)=1+tf+O(t2) \begin{align*} \det {{ \partial \phi_{t} (x) } \over { \partial x }} =& \det \left[ E + {{ \partial f } \over { \partial x }} t \right] + O \left( t^{2} \right) \\ =& 1 + \text{tr} \left[ {{ \partial f } \over { \partial x }} \right] t + O \left( t^{2} \right) \\ =& 1 + \left( {{ d f_{1} } \over { d x_{1} }} + \cdots + {{ d f_{n} } \over { d x_{n} }} \right) t + O \left( t^{2} \right) \\ =& 1 + t \nabla \cdot f + \mathcal{O} \left( t^{2} \right) \end{align*} ここで tr\text{tr} はトレースTraceで、行列の対角成分の合計を意味する。両辺に D0dx\int_{D_{0}} \cdot dx を取ると V(t)=D0detϕt(x)xdx=D01dx+D0tfdx+O(t2)=V(D0)+tD0fdx+O(t2) \begin{align*} V(t) =& \int_{D_{0}} \det {{ \partial \phi_{t} (x) } \over { \partial x }} dx \\ =& \int_{D_{0}} 1 dx + \int_{D_{0}} t \nabla \cdot f dx + O \left( t^{2} \right) \\ =& V \left( D_{0} \right) + t \int_{D_{0}} \nabla \cdot f dx + O \left( t^{2} \right) \end{align*} V(D0)=V(0)V \left( D_{0} \right) = V(0) を左辺に移すと V(t)V(0)=tD0fdx+O(t2) V(t) - V(0) = t \int_{D_{0}} \nabla \cdot f dx + O \left( t^{2} \right) 両辺を (t0)(t - 0) で割ると V(t)V(0)t0=D0fdx+O(t1) {{ V(t) - V(0) } \over { t-0 }} = \int_{D_{0}} \nabla \cdot f dx + O \left( t^{1} \right) テイラー近似を t=0t = 0 近傍で用いたので、t0t \to 0 のとき O(t1)0O \left( t^{1} \right) \to 0 であるから、 dVdtt=0=D0fdx \left.{{ d V } \over { d t }}\right|_{t = 0} = \int_{D_{0}} \nabla \cdot f dx


Part 2. t=t0t = t_{0} への拡張

y:=ϕt0(x)y := \phi_{t_{0}} (x) と仮定すると、Part 1.の議論を同じように進めて次を得る。 dVdtt=t0=D0fdy \left.{{ d V } \over { d t }}\right|_{t = t_{0}} = \int_{D_{0}} \nabla \cdot f dy ある定数 cRc \in \mathbb{R} に対して f=c\nabla \cdot f = c なら、 dVdtt=t0=D0fdy=D0cdy=cD0dy \left.{{ d V } \over { d t }}\right|_{t = t_{0}} = \int_{D_{0}} \nabla \cdot f dy = \int_{D_{0}} c dy = c \int_{D_{0}} dy yyxx を時間 t0t_{0} だけフローに従って変換されたものであるから、D0dy\int_{D_{0}} dyt0t_{0} 時点のボリュームとなる。これは任意の t0Rt_{0} \in \mathbb{R} に対しても同様であるから、 V=cV V’ = c V 上記の微分方程式は自明な解 V(t)=ectV(0)V(t) = e^{ct} V (0) を持つ。


Part 3.

上記 Part 2. の最後の式で c=0c = 0 なら V(t)=ectV(0)=V(0)V(t) = e^{ct} V (0) = V(0) である。

参考文献


  1. Wiggins. (2003). Introduction to Applied Nonlinear Dynamical Systems and Chaos Second Edition(2nd Edition): 99~100. ↩︎