変数分離

微分方程式
ナビエ–ストークス微分方程式、障害物のまわりの気流をシミュレートするのに用いられる。
範囲 .mw-parser-output .hlist ul,.mw-parser-output .hlist ol{padding-left:0}.mw-parser-output .hlist li,.mw-parser-output .hlist dd,.mw-parser-output .hlist dt{margin-right:0;display:inline-block}.mw-parser-output .hlist dt:after,.mw-parser-output .hlist dd:after,.mw-parser-output .hlist li:after{white-space:normal}.mw-parser-output .hlist li:after,.mw-parser-output .hlist dd:after{content:" · ";font-weight:bold}.mw-parser-output .hlist dt:after{content:": "}.mw-parser-output .hlist-pipe dd:after,.mw-parser-output .hlist-pipe li:after{content:" | ";font-weight:normal}.mw-parser-output .hlist-hyphen dd:after,.mw-parser-output .hlist-hyphen li:after{content:" - ";font-weight:normal}.mw-parser-output .hlist-comma dd:after,.mw-parser-output .hlist-comma li:after{content:"、 ";font-weight:normal}.mw-parser-output .hlist dd:last-child:after,.mw-parser-output .hlist dt:last-child:after,.mw-parser-output .hlist li:last-child:after{content:none}.mw-parser-output .hlist dd dd:first-child:before,.mw-parser-output .hlist dd dt:first-child:before,.mw-parser-output .hlist dd li:first-child:before,.mw-parser-output .hlist dt dd:first-child:before,.mw-parser-output .hlist dt dt:first-child:before,.mw-parser-output .hlist dt li:first-child:before,.mw-parser-output .hlist li dd:first-child:before,.mw-parser-output .hlist li dt:first-child:before,.mw-parser-output .hlist li li:first-child:before{content:" (";font-weight:normal}.mw-parser-output .hlist dd dd:last-child:after,.mw-parser-output .hlist dd dt:last-child:after,.mw-parser-output .hlist dd li:last-child:after,.mw-parser-output .hlist dt dd:last-child:after,.mw-parser-output .hlist dt dt:last-child:after,.mw-parser-output .hlist dt li:last-child:after,.mw-parser-output .hlist li dd:last-child:after,.mw-parser-output .hlist li dt:last-child:after,.mw-parser-output .hlist li li:last-child:after{content:") ";font-weight:normal}.mw-parser-output .hlist ol{counter-reset:listitem}.mw-parser-output .hlist ol>li{counter-increment:listitem}.mw-parser-output .hlist ol>li:before{content:" "counter(listitem)" ";white-space:nowrap}.mw-parser-output .hlist dd ol>li:first-child:before,.mw-parser-output .hlist dt ol>li:first-child:before,.mw-parser-output .hlist li ol>li:first-child:before{content:" ("counter(listitem)" "} 自然科学 工学 天文学 物理学 化学 生物学 地質学 応用数学 連続体力学 カオス理論 力学系 社会科学 経済学 個体群動態論
分類
タイプ 常 偏 微分代数（英語版） 積分微分 分数階 線型 非線型（英語版） 変数のタイプにより 独立変数と従属変数（英語版） 自律（英語版） 複素（英語版） Coupled/Decoupled 完全（英語版） 斉次（英語版）/非斉次（英語版） 特徴 階数（英語版） 作用素 記法
過程との関係 差分 (離散類似) 確率 遅延（英語版）
解
一般的な話題 ピカール・リンデレフの定理（英語版） ロンスキアン Phase portrait（英語版） 相空間 リャプノフ/漸近（英語版）/指数安定性（英語版） 収束率（英語版） 級数（英語版）/積分解 数値積分 ディラックのデルタ関数
解法（英語版） Inspection 特性曲線法 オイラー 有限差分 (クランク・ニコルソン（英語版）) 有限要素 有限体積 ガリョルキン（英語版） 積分因子 積分変換 摂動論 ルンゲ・クッタ 変数分離 未定係数（英語版）
人物 アイザック・ニュートン レオンハルト・オイラー エミール・ピカール ユゼフ・マリア・ハーネー＝ウロンスキー エルンスト・リンデレフ（英語版） ルドルフ・リプシッツ（英語版） オーギュスタン＝ルイ・コーシー ジョン・クランク（英語版） フィリス・ニコルソン（英語版） カール・ダーフィト・トルメ・ルンゲ マルチン・ヴィルヘルム・クッタ（英語版）
表 話 編 歴

変数分離（へんすうぶんり、separation of variables）は、常微分方程式や偏微分方程式を解くための手法。方程式を変形することにより、2つあるいはそれ以上の変数が式の右辺・左辺に分かれるようにすること。

常微分方程式に対して用いるときと、偏微分方程式に対して用いるときは、そのやり方がかなり異なっているが、それぞれの変数に依存する部分を両辺に分けるという点では共通している。

目次

• 
  1 常微分方程式
  
  
  
  • 
    1.1 例１
    
  
  
  • 
    1.2 例２
    
  
  
  
  


• 
  2 偏微分方程式
  
  
  
  • 
    2.1 例1
    
  
  
  • 
    2.2 例2
    
  
  
  
  


• 
  3 関連項目
  


• 
  4 参考文献
  

常微分方程式

次の形に書かれる常微分方程式を考える。

ddxf(x)=g(x)h(f(x)){displaystyle {frac {d}{dx}}f(x)=g(x),h(f(x))}

あるいは y = f (x ) と書くことにより、もっと簡単に

dydx=g(x)h(y)(1){displaystyle {frac {dy}{dx}}=g(x)h(y)qquad qquad (1)}

ここで、h (y ) ≠ 0 のとき、両辺を h (y ) で割って

1h(y)dydx=g(x){displaystyle {frac {1}{h(y)}}{frac {dy}{dx}}=g(x)}

となる。この両辺を x で積分すると

∫1h(y)dydxdx=∫g(x)dx+C(2){displaystyle int {frac {1}{h(y)}}{frac {dy}{dx}}dx=int g(x),dx+Cqquad qquad (2)}

で、置換積分の法則により

∫1h(y)dy=∫g(x)dx+C{displaystyle int {frac {1}{h(y)}}dy=int g(x),dx+C}

となる。

この両辺の積分を実行すれば、微分方程式の解が求まる。この手続きは実際のところ、導関数 dy /dx を分数とみなして分母を払うのと同じことである。そうすることによって解くのがもっと簡単になる。具体的なやり方は以下の例で示す。

（注意：両辺の積分に対し

∫1h(y)dy+C1=∫g(x)dx+C2{displaystyle int {frac {1}{h(y)}}dy+C_{1}=int g(x),dx+C_{2}}

のように積分定数をそれぞれ書く必要はない。これは C = C₂ - C₁ として定数を一つにまとめることが出来るからである。）

例１

常微分方程式

df(x)dx=f(x)(1−f(x)){displaystyle {frac {df(x)}{dx}}=f(x),(1-f(x))}

は、より簡単に

dydx=y(1−y){displaystyle {frac {dy}{dx}}=y(1-y)}

と書けるが、ここで g (x ) = 1, h (y ) = y (1-y ) とすれば、この微分方程式は(1)式の形になる。よってこの微分方程式は変数分離が可能である。

上記の説明により、dy と dx を分けて扱うことができる。すなわち両辺に dx をかける。それから両辺を y (1-y ) でわると

dyy(1−y)=dx{displaystyle {frac {dy}{y(1-y)}}=dx}

となる。これで x と y を分離することができた。つまり、x は右辺のみにあり、y は左辺のみにある状態になった。

両辺を積分して

∫dyy(1−y)=∫dx{displaystyle int {frac {dy}{y(1-y)}}=int dx}

となる。これを部分分数分解して

∫(1y+11−y)dy=∫dx{displaystyle int left({frac {1}{y}}+{frac {1}{1-y}}right)dy=int dx}

そして積分を計算すると

log⁡y−log⁡(1−y)=x+C{displaystyle log {y}-log(1-y)=x+C}

ここで C は積分定数である。多少の計算により、y について解くことができて

y=11+Be−x{displaystyle y={frac {1}{1+Be^{-x}}}}

となる。B は任意の定数である。この解を x で微分すれば、この解が正しいことを確かめることができる。その結果はもともとの微分方程式と一致するはずだ。

ところで、両辺を y (1-y ) で割るにあたって、y (x ) = 0 や y (x ) = 1 が微分方程式の解になるかどうかを検討する必要がある。そのような解は特異解となりうる。

例２

変数分離を用いて解ける２階非線型常微分方程式の例^[1]。

xd2ydx2+(1+P(y))dydx=0.{displaystyle x{frac {;d^{2}y;}{dx^{2}}}+{bigl (}1+P(y){bigr )}{frac {;dy;}{dx}}=0.}

この微分方程式は，このまま両辺を x で積分し，部分積分法を適用して整理すると，変数分離を用いて解くことができる。一般解は，

x=C2exp⁡(∫dyC1−∫P(y)dy){displaystyle x=C_{2}exp {}{Biggl (}int {frac {;dy;}{;C_{1}-displaystyle int P(y),dy;}}{Biggr )}}

と表示される^[1]。ここに，P(y) は既知関数であり，C₁, C₂ は積分定数である。
ただし，C₂ ≠ 0 とする。求積法で解ける微分方程式は，変数分離を用いることが多い。

偏微分方程式

n 変数関数

F(x1,x2,…,xn){displaystyle F(x_{1},x_{2},dots ,x_{n})}

についての偏微分方程式を解くにあたって、その解の形を

F=F1(x1)F2(x2)⋯Fn(xn){displaystyle F=F_{1}(x_{1}),F_{2}(x_{2})cdots F_{n}(x_{n})}

あるいは

F=f1(x1)+f2(x2)+⋯+fn(xn){displaystyle F=f_{1}(x_{1})+f_{2}(x_{2})+cdots +f_{n}(x_{n})}

のように仮定すると、偏微分方程式がいくつかの常微分方程式になる場合がある。多くの場合、個々の変数に対して、微分方程式からは決定できない分離定数が現れることになる。

例1

未知関数 F (x, y, z ) と、それが満たす偏微分方程式

∂F∂x+∂F∂y+∂F∂z=0(1){displaystyle {frac {partial F}{partial x}}+{frac {partial F}{partial y}}+{frac {partial F}{partial z}}=0qquad qquad (1)}

を考える。関数 F (x, y, z ) が

F(x,y,z)=X(x)+Y(y)+Z(z)(2){displaystyle F(x,y,z)=X(x)+Y(y)+Z(z)qquad qquad (2)}

の形に書けると仮定すると、(1)式は

dXdx+dYdy+dZdz=0{displaystyle {frac {dX}{dx}}+{frac {dY}{dy}}+{frac {dZ}{dz}}=0}

となる。なぜなら ∂F /∂x = dX /dx などが成り立つからである。

いま、X' (x ) は x のみに依存し、Y' (y ) は y のみに依存し、そしてZ' (z ) についても同様である。また、微分方程式(1)は任意の x, y, z について成り立つ。これより、それぞれの項が定数になることがわかる。すなわち

dXdx=c1,dYdy=c2,dZdz=c3(3){displaystyle {frac {dX}{dx}}=c_{1},quad {frac {dY}{dy}}=c_{2},quad {frac {dZ}{dz}}=c_{3}qquad qquad (3)}

となる。定数 c₁, c₂, c₃ は

c1+c2+c3=0(4){displaystyle c_{1}+c_{2}+c_{3}=0qquad qquad (4)}

を満たす。(3)式は3つの微分方程式のセットである。この場合、それぞれの微分方程式は単に積分するだけで解を得ることができて、答えは

F(x,y,z)=c1x+c2y+c3z+c4(5){displaystyle F(x,y,z)=c_{1}x+c_{2}y+c_{3}z+c_{4}qquad qquad (5)}

となる。積分定数 c₄ は初期条件によって定まる。

例2

次の微分方程式を考える。

∇2v+λv=∂2v∂x2+∂2v∂y2+λv=0{displaystyle nabla ^{2}v+lambda v={frac {partial ^{2}v}{partial x^{2}}}+{frac {partial ^{2}v}{partial y^{2}}}+lambda v=0}

まず解の形を

v=X(x)Y(y){displaystyle v=X(x)Y(y)}

とおく。これ以外の解は、このような解の線形結合になっていると考える。

これを微分方程式に代入すると

X″(x)Y(y)+X(x)Y″(y)+λX(x)Y(y)=0{displaystyle X''(x)Y(y)+X(x)Y''(y)+lambda X(x)Y(y)=0}

となる。この両辺を X (x ) でわって

X″(x)Y(y)X(x)+Y″(y)+λY(y)=0{displaystyle {frac {X''(x)Y(y)}{X(x)}}+Y''(y)+lambda Y(y)=0}

さらに Y (y ) でわって

X″(x)X(x)+Y″(y)+λY(y)Y(y)=0{displaystyle {frac {X''(x)}{X(x)}}+{frac {Y''(y)+lambda Y(y)}{Y(y)}}=0}

すると X'' (x )/X (x ) は x のみの関数で、もう一つの項は y のみの関数だから、分離定数を用いて

X″(x)X(x)=−Y″(y)+λY(y)Y(y)=k{displaystyle {frac {X''(x)}{X(x)}}=-{frac {Y''(y)+lambda Y(y)}{Y(y)}}=k}

と書ける。これによって2つの常微分方程式

X″(x)X(x)=k,X″(x)=kX(x){displaystyle {frac {X''(x)}{X(x)}}=k,quad X''(x)=kX(x)}

および

Y″(y)+λY(y)Y(y)=−k,Y″(y)+(λ+k)Y(y)=0{displaystyle {frac {Y''(y)+lambda Y(y)}{Y(y)}}=-k,quad Y''(y)+(lambda +k)Y(y)=0}

が得られ、それぞれ解くことができる。もとの問題が境界値問題であるなら、その境界条件を用いて解を定めることができる。

関連項目

• 求積法


• 微分方程式


• 常微分方程式

参考文献

1. ^ ^ab長島 隆廣 『常微分方程式80余例とその厳密解』 近代文芸社、2005年 ISBN 4-7733-7282-6. 国立国会図書館蔵書, 請求記号：MA117-H55（東京　本館書庫）
   

This page is only for reference, If you need detailed information, please check here