積分因子

微分方程式
ナビエ–ストークス方程式。障害物の周囲の気流のシミュレーションに用いられる。
ナビエ–ストークス方程式。障害物の周囲の気流のシミュレーションに用いられる。
分類
タイプ
変数のタイプにより
  • 独立変数と従属変数(英語版)
特徴
過程との関係
  • 差分 (離散類似)
  • 確率
    • 確率偏(英語版)
  • 遅延(英語版)
一般的な話題

積分因子 (せきぶんいんし、: integrating factor) とは微分方程式の解法に用いられる関数である。常微分方程式の解法で最もよく用いられ、積分因子を掛けることにより不完全微分から完全微分(積分するとスカラー場を与える)を得ることができる。特に熱力学の分野で用いられ、そこではエントロピーを完全微分にするために温度が積分因子となる。

2変数の方程式の場合には積分因子は必ず存在する[1]

カラテオドリの定理

n {\displaystyle n} 次元多様体 M {\displaystyle M} の領域 U {\displaystyle U} で定義された 1-形式 ψ {\displaystyle \psi } について

ψ = f d g {\displaystyle \psi =f\,dg}

が成立するような関数 f {\displaystyle f} , g {\displaystyle g} (ただし f 0 {\displaystyle f\neq 0} ) が存在するとき、 1 / f {\displaystyle 1/f} ψ {\displaystyle \psi } の積分因子と呼ぶ[2][3][4]。1-形式 ψ {\displaystyle \psi } の積分因子の存在に関して、カラテオドリの定理は以下の3命題が同値であることを主張する[4]

  1. M {\displaystyle M} の任意の点 x {\displaystyle x} について、点 x {\displaystyle x} の近傍 V {\displaystyle V} が存在し、 V {\displaystyle V} 内の任意の x {\displaystyle x} の近傍 W {\displaystyle W} に、区分的 C {\displaystyle C^{\infty }} -級曲線 γ {\displaystyle \gamma } であって γ ˙ {\displaystyle {\dot {\gamma }}} が定義されるすべての点で ψ { γ ˙ ( t ) } = 0 {\displaystyle \psi \{{\dot {\gamma }}(t)\}=0} を満たすような曲線 γ {\displaystyle \gamma } によって点 x {\displaystyle x} と結ばれることのないような点 y W {\displaystyle y\in W} が存在する。
  2. ψ d ψ = 0 {\displaystyle \psi \wedge d\psi =0}
  3. 任意の点 x M {\displaystyle x\in M} にある近傍 V {\displaystyle V} が存在し、 ψ {\displaystyle \psi } V {\displaystyle V} への制限は積分因子を持つ。

この定理は1909年にコンスタンティン・カラテオドリによって熱力学第二法則の定式化のために導入された[5][6] ψ {\displaystyle \psi } 閉形式(すなわち d ψ = 0 {\displaystyle d\psi =0} が成立するもの)であればポアンカレの補題により積分因子 1 {\displaystyle 1} が存在する[7]。また n 2 {\displaystyle n\leq 2} のときカラテオドリの定理によりすべての 1-形式について積分因子の存在が保証される(これはJohann Friedrich Pfaffによって最初に示された)[2][3]

積分因子を用いた常微分方程式の解法の例

次のような1階の線形常微分方程式を考える。

y + P ( x ) y = Q ( x ) ( 1 ) {\displaystyle y'+P(x)y=Q(x)\quad \quad \quad (1)}

この方程式に対し、適当な積分因子 M ( x ) {\displaystyle M(x)} を (1) 式の両辺に掛け、左辺に積の微分の公式を適用できるようにすると、

M ( x ) y + M ( x ) P ( x ) y = M ( x ) Q ( x ) ( 2 ) {\displaystyle M(x)y'+M(x)P(x)y=M(x)Q(x)\quad \quad \quad (2)}
( M ( x ) y ) = M ( x ) Q ( x ) ( 3 ) {\displaystyle (M(x)y)'=M(x)Q(x)\quad \quad \quad (3)}

となるから、(3) 式を積分して

y M ( x ) = Q ( x ) M ( x ) d x {\displaystyle yM(x)=\int Q(x)M(x)\,dx}

となり、これよりもとの微分方程式の解として

y = Q ( x ) M ( x ) d x M ( x ) {\displaystyle y={\frac {\int Q(x)M(x)\,dx}{M(x)}}\,}

が得られる。

次に積分因子 M ( x ) {\displaystyle M(x)} を具体的に求める。(2), (3) 式それぞれの左辺が等しくなるように M ( x ) {\displaystyle M(x)} をとっていることから、

M ( x ) y + M ( x ) y = M ( x ) y + M ( x ) P ( x ) y {\displaystyle M'(x)y+M(x)y'=M(x)y'+M(x)P(x)y\quad \quad \quad }

となり、 M ( x ) {\displaystyle M(x)} がつぎの微分方程式

M ( x ) = M ( x ) P ( x ) ( 4 ) {\displaystyle M'(x)=M(x)P(x)\quad \quad \quad (4)\,}

を満たすことがわかる。この式を変形すると、

M ( x ) M ( x ) = P ( x ) ( 5 ) {\displaystyle {\frac {M'(x)}{M(x)}}=P(x)\quad \quad \quad (5)\,}
( ln M ( x ) ) = P ( x ) {\displaystyle (\ln M(x))'=P(x)}

したがって

M ( x ) = exp ( P ( x ) d x ) {\displaystyle M(x)=\exp \left(\int P(x)\,dx\right)}

となる。

脚注

  1. ^ 和達三樹; 十河清; 出口哲生『ゼロからの熱力学と統計力学』岩波書店、2005年、26頁。ISBN 4-00-006700-1。 
  2. ^ a b Martinás, Katalin; Brodszky, Ildikó (2000). “Thermodynamics of Gyula Farkas - A new (old) approach to entropy”. Periodica Polytechnica: Chemical Engineering 44 (1): 17-27. https://pp.bme.hu/ch/article/download/286/179/. 
  3. ^ a b Pogliani, Lionello; Berberan-Santos, Mario N. (2000). Journal of Mathematical Chemistry 28 (1/3): 313–324. doi:10.1023/A:1018834326958. ISSN 02599791. 
  4. ^ a b Boyling, J. B. (1968). “Carathéodory's principle and the existence of global integrating factors”. Communications in Mathematical Physics 10 (1): 52–68. doi:10.1007/BF01654133. ISSN 0010-3616. https://projecteuclid.org/download/pdf_1/euclid.cmp/1103840983. 
  5. ^ Carathéodory, C. (1909). “Untersuchungen über die Grundlagen der Thermodynamik”. Mathematische Annalen 67 (3): 355–386. doi:10.1007/BF01450409. ISSN 0025-5831. ; 英訳: http://neo-classical-physics.info/uploads/3/0/6/5/3065888/caratheodory_-_thermodynamics.pdf
  6. ^ Bernstein, B. (1960). “Proof of Carathéodory's Local Theorem and Its Global Application to Thermostatics”. Journal of Mathematical Physics 1 (3): 222–224. doi:10.1063/1.1703655. ISSN 0022-2488. 
  7. ^ Nakahara, Mikio (2003). Geometry, Topology and Physics. CRC Press . Section 6.3.

外部リンク

  • Joakim Munkhammar. "Integrating Factor". mathworld.wolfram.com (英語).
スタブアイコン

この項目は、解析学に関連した書きかけの項目です。この項目を加筆・訂正などしてくださる協力者を求めています(プロジェクト:数学/Portal:数学)。

  • 表示
  • 編集