積分因子

提供: testwiki
ナビゲーションに移動 検索に移動

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

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

カラテオドリの定理

n 次元多様体 M の領域 U で定義された 1-形式 ψ について テンプレート:Indent が成立するような関数 f, g (ただし f0) が存在するとき、1/fψ の積分因子と呼ぶ[2][3][4]。1-形式 ψ の積分因子の存在に関して、カラテオドリの定理は以下の3命題が同値であることを主張する[4]

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

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

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

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

y+P(x)y=Q(x)(1)

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

M(x)y+M(x)P(x)y=M(x)Q(x)(2)
(M(x)y)=M(x)Q(x)(3)

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

yM(x)=Q(x)M(x)dx

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

y=Q(x)M(x)dxM(x)

が得られる。

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

M(x)y+M(x)y=M(x)y+M(x)P(x)y

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

M(x)=M(x)P(x)(4)

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

M(x)M(x)=P(x)(5)
(lnM(x))=P(x)

したがって

M(x)=exp(P(x)dx)

となる。

脚注

テンプレート:Reflist

外部リンク

テンプレート:Mathanalysis-stub

  1. テンプレート:Cite
  2. 2.0 2.1 引用エラー: 無効な <ref> タグです。「MartinasBrodszky2000」という名前の注釈に対するテキストが指定されていません
  3. 3.0 3.1 引用エラー: 無効な <ref> タグです。「PoglianiBerberan-Santos2000」という名前の注釈に対するテキストが指定されていません
  4. 4.0 4.1 引用エラー: 無効な <ref> タグです。「Boyling1968」という名前の注釈に対するテキストが指定されていません
  5. 引用エラー: 無効な <ref> タグです。「Caratheodory1909」という名前の注釈に対するテキストが指定されていません
  6. 引用エラー: 無効な <ref> タグです。「Bernstein1960」という名前の注釈に対するテキストが指定されていません
  7. テンプレート:Cite book. Section 6.3.