変数分離のソースを表示

{{Differential equations}}
'''変数分離'''（へんすうぶんり、{{lang-en-short|Separation of variables}}）は、[[常微分方程式]]や[[偏微分方程式]]を解くための手法。方程式を変形することにより、2つあるいはそれ以上の変数が式の右辺・左辺に分かれるようにすること。

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

== 常微分方程式 ==
次の形に書かれる[[常微分方程式]]を考える。
{{Indent|<math> \frac{d}{dx} f(x) = g(x)\,h(f(x)) </math>}}
あるいは ''y'' = ''f'' (''x'' ) と書くことにより、もっと簡単に
{{Indent|<math> \frac{dy}{dx} = g(x)h(y)    \qquad\qquad (1) </math>}}
ここで、''h'' (''y'' ) &ne; 0 のとき、両辺を ''h'' (''y'' ) で割って
{{Indent|<math> \frac{1}{h(y)} \frac{dy}{dx} = g(x) </math>}}
となる。この両辺を ''x'' で積分すると
{{Indent|<math> \int \frac{1}{h(y)} \frac{dy}{dx} dx = \int g(x)\,dx + C 
   \qquad\qquad (2) </math>}}
で、[[置換積分]]の法則により
{{Indent|<math> \int \frac{1}{h(y)} dy = \int g(x)\,dx + C </math>}}
となる。

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

（注意：両辺の積分に対し
{{Indent|<math> \int \frac{1}{h(y)} dy + C_1 = \int g(x)\,dx + C_2 </math>}}
のように積分定数をそれぞれ書く必要はない。これは ''C'' = ''C''<sub>2</sub> - ''C''<sub>1</sub> として定数を一つにまとめることが出来るからである。）

=== 例1 ===
常微分方程式
{{Indent|<math> \frac{df(x)}{dx} = f(x)\,(1-f(x)) </math>}}
は、より簡単に
{{Indent|<math> \frac{dy}{dx} = y(1-y) </math>}}
と書けるが、ここで ''g'' (''x'' ) = 1, ''h'' (''y'' ) = ''y'' (1-''y'' ) とすれば、この微分方程式は(1)式の形になる。よってこの微分方程式は変数分離が可能である。

上記の説明により、''dy'' と ''dx'' を分けて扱うことができる。すなわち両辺に ''dx'' をかける。それから両辺を ''y'' (1-''y'' ) でわると
{{Indent|<math> \frac{dy}{y(1-y)} = dx </math>}}
となる。これで ''x'' と ''y'' を分離することができた。つまり、''x'' は右辺のみにあり、''y'' は左辺のみにある状態になった。

両辺を積分して
{{Indent|<math> \int \frac{dy}{y(1-y)} = \int dx </math>}}
となる。これを[[部分分数分解]]して
{{Indent|<math> \int \left( \frac{1}{y} + \frac{1}{1-y} \right) dy = \int dx </math>}}
そして積分を計算すると
{{Indent|<math> \log{y}-\log(1-y)=x+C</math>}}
ここで ''C'' は積分定数である。多少の計算により、''y'' について解くことができて
{{Indent|<math> y = \frac{1}{1 + Be^{-x}} </math>}}
となる。''B'' は任意の定数である。この解を ''x'' で微分すれば、この解が正しいことを確かめることができる。その結果はもともとの微分方程式と一致するはずだ。

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

=== 例2 ===

変数分離を用いて解ける2階非線形常微分方程式の例<ref name="diffequ-2005ngi">{{cite book|和書|author=長島 隆廣 |title=常微分方程式80余例とその厳密解|publisher= 近代文芸社|year=2005|isbn= 4-7733-7282-6}} 国立国会図書館蔵書, 請求記号：MA117-H55（東京　本館書庫）</ref>。

::<math>x\frac{\;d^2 y\;}{dx^2}+ \bigl(1+P(y)\bigr) \frac{\;dy\;}{dx}=0.</math>
             <!-- ★ d^2y/dx^2 +(1+P(x))(dy/dx)=0, (80余例) p.88 (コメント) --> 
この微分方程式は，このまま両辺を ''x'' で積分し，部分積分法を適用して整理すると，変数分離を用いて解くことができる。一般解は，

::<math>x=C_2 \exp{}\Biggl( \int \frac{\;dy\;}{\;C_1 - \displaystyle \int P(y)\,dy\;} \Biggr)</math>
             <!-- ★ (80余例) p.88 (コメント) --> 

と表示される<ref name="diffequ-2005ngi"/>。ここに，''P''(''y'') は既知関数であり，''C''<sub>1</sub>, ''C''<sub>2</sub> は積分定数である。
ただし，''C''<sub>2</sub> &ne; 0 とする。求積法で解ける微分方程式は，変数分離を用いることが多い<ref name="diffequ-2018rmap">{{Cite web|和書|author=長島 隆廣|title=常微分方程式80余例と求積法による解法|date=2018年12月|publisher= researchmap |format=PDF|url=https://researchmap.jp/multidatabases/multidatabase_contents/detail/263160/16f8fddfba5ab789f6475ac2962bfd31?frame_id=539358|accessdate=2020-06-29}}</ref>。

== 偏微分方程式 ==
''n'' 変数関数
{{Indent|<math> F(x_1, x_2, \dots , x_n) </math>}}
についての[[偏微分方程式]]を解くにあたって、その解の形を
{{Indent|<math> F = F_1(x_1)\,F_2(x_2) \cdots F_n(x_n) </math>}}
あるいは
{{Indent|<math> F = f_1(x_1) + f_2(x_2) + \cdots + f_n(x_n) </math>}}
のように仮定すると、偏微分方程式がいくつかの[[常微分方程式]]になる場合がある。多くの場合、個々の変数に対して、微分方程式からは決定できない'''分離定数'''が現れることになる。

=== 例1 ===
未知関数 ''F'' (''x'', ''y'', ''z'' ) と、それが満たす偏微分方程式
{{Indent|<math> \frac{\partial F}{\partial x} + \frac{\partial F}{\partial y}  
         + \frac{\partial F}{\partial z}  =  0    \qquad\qquad (1) </math>}}
を考える。関数 ''F'' (''x'', ''y'', ''z'' ) が
{{Indent|<math> F(x,y,z) = X(x) + Y(y) + Z(z)    \qquad\qquad (2) </math>}}
の形に書けると仮定すると、(1)式は
{{Indent|<math> \frac{dX}{dx} + \frac{dY}{dy} + \frac{dZ}{dz}  =  0 </math>}}
となる。なぜなら &part;''F'' /&part;''x'' = ''dX'' /''dx'' などが成り立つからである。

いま、''X''<nowiki>'</nowiki> (''x'' ) は ''x'' のみに依存し、''Y''<nowiki>'</nowiki> (''y'' ) は ''y'' のみに依存し、そして''Z''<nowiki>'</nowiki> (''z'' ) についても同様である。また、微分方程式 (1) は任意の ''x'', ''y'', ''z'' について成り立つ。これより、それぞれの項が定数になることがわかる。すなわち
{{Indent|<math>
  \frac{dX}{dx} = c_1, \quad 
  \frac{dY}{dy} = c_2, \quad 
  \frac{dZ}{dz} = c_3  \qquad\qquad (3)
</math>}}
となる。定数 ''c''<sub>1</sub>, ''c''<sub>2</sub>, ''c''<sub>3</sub> は
{{Indent|<math> c_1 + c_2 + c_3 = 0    \qquad\qquad (4) </math>}}
を満たす。(3) 式は3つの微分方程式のセットである。この場合、それぞれの微分方程式は単に積分するだけで解を得ることができて、答えは
{{Indent|<math> F(x,y,z) = c_1 x + c_2 y + c_3 z + c_4    \qquad\qquad (5) </math>}}
となる。積分定数 ''c''<sub>4</sub> は初期条件によって定まる。

=== 例2 ===
以下の偏微分方程式を考える。；
{{Indent|<math>
  \nabla^2 v + \lambda v
    =  \frac{\partial^2 v}{\partial x^2}
         + \frac{\partial^2 v}{\partial y^2} + \lambda v
    =  0
</math>}}
まず解の形を
{{Indent|<math> v = X(x) Y(y) </math>}}
とおく。これ以外の解は、このような解の[[線形結合]]になっていると考える。

これを微分方程式に代入すると
{{Indent|<math> X''(x) Y(y) + X(x) Y''(y) + \lambda X(x) Y(y) = 0 </math>}}
となる。この両辺を ''X'' (''x'' ) で割って
{{Indent|<math> \frac{X''(x) Y(y)}{X(x)} + Y''(y) + \lambda Y(y) = 0 </math>}}
更に ''Y'' (''y'' ) で割って
{{Indent|<math> \frac{X''(x)}{X(x)} + \frac{Y''(y) + \lambda Y(y)}{Y(y)} = 0 </math>}}
すると ''X''<nowiki>''</nowiki> (''x'' )/''X'' (''x'' ) は ''x'' のみの関数で、もう一つの項は ''y'' のみの関数だから、分離定数を用いて
{{Indent|<math> \frac{X''(x)}{X(x)} = -\frac{Y''(y) + \lambda Y(y)}{Y(y)} = k </math>}}
と書ける。これによってニつの2階線型常微分方程式
{{Indent|<math> \frac{X''(x)}{X(x)} = k, \quad X''(x) = kX(x) </math>}}
および
{{Indent|<math> \frac{Y''(y) + \lambda Y(y)}{Y(y)} = -k, \quad 
   Y''(y) + (\lambda+k)Y(y) = 0 </math>}}
が得られ、それぞれ解くことができる。もとの問題が[[境界値問題]]であるなら、その境界条件を用いて解を定めることができる。

== 参考文献 ==
<references />

== 関連項目 ==
*[[求積法]]
*[[微分方程式]]
*[[常微分方程式]]

{{DEFAULTSORT:へんすうふんり}}
[[Category:微分方程式]]
[[Category:常微分方程式]]
[[Category:物理数学]]
[[Category:数学に関する記事]]