定数変化法のソースを表示

[[数学]]における'''係数変化法'''（けいすうへんかほう、{{lang-en-short|''variation of parameters''}}）または'''定数変化法'''（じょうすうへんかほう、ていすうへんかほう、{{lang-en-short|''variation of constants''}}）は[[線型微分方程式|線型非斉次]]な[[常微分方程式]]の一般解法である。[[ジョゼフ＝ルイ・ラグランジュ|ラグランジュ]]の定数変化法と呼ばれることもある。

一階の非斉次線型微分方程式は、かなり労力の少ない[[積分因子]]や[[未定係数法]]を通じて解けるのが普通であるが、それらは推測から来る経験則として利用するもので、しかもすべての非斉次微分方程式に対してうまくいくわけではない。

定数変化法は線型[[偏微分方程式]]にも拡張することができて、具体的に[[熱方程式]]、[[波動方程式]]、[[振動板方程式]]などの線型発展方程式の非斉次問題が解ける。この設定での定数変化法を用いた解法は、むしろ[[デュアメルの原理]]としてよく知られている。この呼称は、非斉次熱方程式の解法として定数変化法を初めて適用した[[ジャン＝マリー・デュアメル]]に因むものであり、一般の定数変化法をデュアメルの原理と呼ぶこともある。

== 解法の説明 ==
階数 ''n'' の非斉次常微分方程式
: <math>\text{(i) }\quad y^{(n)}(x) + \sum_{i=0}^{n-1} a_i(x) y^{(i)}(x) = b(x)</math>
が与えられたとき、''y''<sub>1</sub>, …, ''y''<sub>''n''</sub> を対応する斉次方程式
: <math>\text{(ii) }\quad y^{(n)}(x) + \sum_{i=0}^{n-1} a_i(x) y^{(i)}(x) = 0</math>
の解の基本系とすると、もとの非斉次方程式のひとつの特殊解が
:<math>\text{(iii) }\quad y_p(x) = \sum_{i=1}^{n} c_i(x) y_i(x)</math>
で与えられる。ここで、''c''<sub>''i''</sub>(''x'') は連続函数で方程式
:<math>\text{(iv) }\quad \sum_{i=1}^{n} c_i'(x) y_i^{(j)}(x) = 0\quad (j = 0,\ldots, n-2)</math>
を満足する。(iii) を (i) に代入して (iv) を適用すれば
:<math>\text{(v) }\quad \sum_{i=1}^n c_i'(x) y_i^{(n-1)}(x) = b(x)</math>
を得る。''y''<sub>''i''</sub>(''x'') たちは線型独立だから、条件を満たすにはすべての ''x'' および ''i'' に対して ''c''<sub>''i''</sub>&prime; = 0 でなければならない。従って、''b''(''x'') = 0 の場合には、すべての ''c''<sub>''i''</sub>(''x'') が ''x'' に無関係な定数になる。

この ''n'' 本の線型方程式系は[[クラメルの公式]]を用いて解くことができて、
:<math>c_i'(x) = \frac{W_i(x)}{W(x)}, \quad (i=1,\ldots,n)</math>
が導かれる。ただし、''W''(''x'') は解の基本系の[[ロンスキー行列式]]で、''W''<sub>''i''</sub>(''x'') は基本系のロンスキー行列式の第 ''i''-列を (0, 0, …, ''b''(''x'')) で置き換えたものとする。

ゆえに、非斉次方程式の特殊解は
:<math>\sum_{i=1}^n \left[\int \frac{W_i(x)}{W(x)} dx\right] y_i(x)</math>
と書くことができる。

== 例 ==
=== 特定の二階方程式 ===
方程式
: <math> y''+4y'+4y=\cosh{x}</math>
を解くことを考える。一般解を求めるために、斉次方程式
: <math>y''+4y'+4y=0</math>
を解くと、この固有多項式は
: <math>\lambda^2+4\lambda+4=(\lambda+2)^2</math>
で、固有値 &minus;2 は重根であるから、''u''<sub>1</sub> = ''e''<sup>&minus;2''x''</sup> および ''u''<sub>2</sub> = ''xe''<sup>&minus;2''x''</sup> が基本解となる。これらのロンスキー行列式は
: <math>\begin{vmatrix}
   e^{-2x} & xe^{-2x} \\
 -2e^{-2x} & -e^{-2x}(2x-1)
\end{vmatrix} = e^{-4x}</math>
である。これは 0 でないから、この二つの函数は確かに斉次方程式の一般解を生成する。

従って、''A''(''x'')''u''<sub>1</sub> + ''B''(''x'')''u''<sub>2</sub> が非斉次方程式の一般解となるような ''A''(''x''), ''B''(''x'') を求めればよいが、それには積分
:<math>A(x) = - \int {1\over W}\, u_2(x) b(x)\,dx,\quad B(x) = \int {1 \over W}\, u_1(x)b(x)\,dx</math>
を計算すればよい。結局、
:<math>\begin{cases}
  A(x) = - \int xe^{2x}\cosh{x}\,dx = -{1\over 18}e^x(9(x-1)+e^{2x}(3x-1))+C_1\\
  B(x) = \int e^{2x}\cosh{x}\,dx ={1\over 6}e^{x}(3+e^{2x})+C_2
\end{cases}</math>
が求まる。ただし、''C''<sub>1</sub>, ''C''<sub>2</sub> は積分定数である。

=== 一般の二階方程式 ===
微分方程式
:<math>u''+p(x)u'+q(x)u=f(x)</math>
を解くにあたって、''D'' を[[微分演算子]]として線型微分作用素
:<math>L=D^2+p(x)D+q(x)</math>
を定義すると、''L'' および ''f''(''x'') が既知として、方程式 ''Lu'' = ''f'' を ''u'' に関して解けばよい、ということになる。

定数変化法を用いるために、まずは対応する斉次方程式
: <math>u''+p(x)u'+q(x)u=0</math>
を解かねばならない。この方程式は二階であるから、線型独立な二つの解 ''u''<sub>1</sub>, ''u''<sub>2</sub> が得られれば、定数変化法を適用することができる。

求める微分方程式の一般解 ''u''<sub>''G''</sub> は
:<math>u_G(x)=A(x)u_1(x)+B(x)u_2(x)</math>
の形をしているはずである。ただし、''A''(''x''), ''B''(''x'') は未知で、''u''<sub>1</sub>(''x''), ''u''<sub>2</sub>(''x'') は斉次方程式の解である。''A''(''x'') と ''B''(''x'') がともに定数ならば ''Lu''<sub>''G''</sub> = 0 となるのは明らかである。''A'' = ''A''(''x''), ''B'' = ''B''(''x'') は
:<math>A'(x)u_1(x)+B'(x)u_2(x)=0</math>
となるものと仮定すると、
:<math>u_G'(x) = A(x)u_1'(x)+B(x)u_2'(x)</math>
となり、さらに微分して
:<math>u_G''(x) = A(x)u_1''(x)+B(x)u_2''(x)+A'(x)u_1'(x)+B'(x)u_2'(x)</math>
を得る。従って、''L'' の ''u''<sub>''G''</sub> への作用は
:<math>Lu_G = A(x)Lu_1(x)+B(x)Lu_2(x)+A'(x)u_1'(x)+B'(x)u_2'(x)</math>
と書くことができるが、''u''<sub>1</sub> と ''u''<sub>2</sub> は斉次方程式の解だから
:<math>Lu_G(= f) = A'(x)u_1'(x)+B'(x)u_2'(x)</math>
となる。

以上から連立方程式
:<math>\begin{pmatrix}u_1(x) & u_2(x)\\u_1'(x) & u_2'(x)\end{pmatrix}\begin{pmatrix}A'(x)\\B'(x)\end{pmatrix}=\begin{pmatrix}0\\f\end{pmatrix}</math>
が得られたので、''A''(''x''), ''B''(''x'') を求めるために、これを ''A''&prime;, ''B''&prime; について解くと
:<math>\begin{pmatrix}A'(x) \\B'(x)\end{pmatrix}
={1\over W}\begin{pmatrix}
u_2'(x)  & -u_2(x) \\
-u_1'(x) & u_1(x) \end{pmatrix}
\begin{pmatrix}0\\f\end{pmatrix}</math>
を得る。ただし ''W'' は ''u''<sub>1</sub> と ''u''<sub>2</sub> の[[ロンスキー行列式]]である（''u''<sub>1</sub> と ''u''<sub>2</sub> が線型独立であるという仮定から ''W'' が 0 でないことは分かっている）。ゆえに
:<math>A(x) = - \int {1\over W}\, u_2(x) f(x)\,dx,\quad B(x) = \int {1 \over W}\, u_1(x)f(x)\,dx</math>
を得る。

斉次方程式が比較的容易に解ける限り、この方法で非斉次方程式の一般解の係数を計算することができて、非斉次方程式の完全な一般解を決定することができる。

''A''(''x'') も ''B''(''x'') も任意定数（[[積分定数]]）を除いて定まる点に注意。元々の方程式が二階だったので、積分定数が2個出ることは予期されることである。''A''(''x'') または ''B''(''x'') に定数を加えても、''L'' は[[線型作用素|線型]]だから、''Lu''<sub>''G''</sub>(''x'') の値は変わらない。

== 参考文献 ==
* {{cite book | last1=Coddington | first1=Earl A. | last2=Levinson | first2=Norman | title=Theory of Ordinary Differential Equations | publisher=[[マグロウヒル|McGraw-Hill]] | location=New York | year=1955}}
* {{cite book | title = Elementary Differential Equations and Boundary Value Problems 8th Edition | first1 = W. E. | last1 = Boyce | first2 = R. C. | last2 = DiPrima | publisher = Wiley Interscience | year = 1965}}, pages 186-192, 237-241
* {{cite book | surname = Teschl | given = Gerald |authorlink=Gerald Teschl | title = Ordinary Differential Equations and Dynamical Systems | publisher=[[アメリカ数学会|American Mathematical Society]] | place = [[Providence, Rhode Island|Providence]] | url = http://www.mat.univie.ac.at/~gerald/ftp/book-ode/}}

== 外部リンク ==
*[http://tutorial.math.lamar.edu/classes/de/VariationofParameters.aspx Online Notes / Proof] by Paul Dawkins, [[Lamar University]].
* {{PlanetMath|urlname=VariationOfParameters|title=variation of parameters}}
*[http://www.math-cs.ucmo.edu/~mjms/2007.1/srinivasan10.pdf Motivation of method via celestial mechanics]

{{DEFAULTSORT:ていすうへんかほう}}
[[Category:微分方程式]]
[[Category:常微分方程式]]
[[Category:数学に関する記事]]