variation of constants

ウィキペディア

定数変化法

出典: フリー百科事典『ウィキペディア（Wikipedia）』 (2017/01/22 08:45 UTC 版)

数学における係数変化法（けいすうへんかほう、英: variation of parameters）または定数変化法（じょうすうへんかほう、ていすうへんかほう、英: variation of constants）は線型非斉次な常微分方程式の一般解法である。ラグランジュの定数変化法と呼ばれることもある。

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

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

解法の説明

階数 n の非斉次常微分方程式

${\text{(i) }}\quad y^{(n)}(x)+\sum _{i=0}^{n-1}a_{i}(x)y^{(i)}(x)=b(x)$

が与えられたとき、y₁, …, y_n を対応する斉次方程式

${\text{(ii) }}\quad y^{(n)}(x)+\sum _{i=0}^{n-1}a_{i}(x)y^{(i)}(x)=0$

の解の基本系とすると、もとの非斉次方程式のひとつの特殊解が

${\text{(iii) }}\quad y_{p}(x)=\sum _{i=1}^{n}c_{i}(x)y_{i}(x)$

で与えられる。ここで、c_i(x) は連続函数で方程式

${\text{(iv) }}\quad \sum _{i=1}^{n}c_{i}'(x)y_{i}^{(j)}(x)=0\quad (j=0,\ldots ,n-2)$

を満足する。(iii) を (i) に代入して (iv) を適用すれば

${\text{(v) }}\quad \sum _{i=1}^{n}c_{i}'(x)y_{i}^{(n-1)}(x)=b(x)$

を得る。y_i(x) たちは線型独立だから、条件を満たすにはすべての x および i に対して c_i′ = 0 でなければならない。従って、b(x) = 0 の場合には、すべての c_i(x) が x に無関係な定数になる。

この n 本の線型方程式系はクラメルの公式を用いて解くことができて、

$c_{i}'(x)={\frac {W_{i}(x)}{W(x)}},\quad (i=1,\ldots ,n)$

が導かれる。ただし、W(x) は解の基本系のロンスキー行列式で、W_i(x) は基本系のロンスキー行列式の第 i-列を (0, 0, …, b(x)) で置き換えたものとする。

ゆえに、非斉次方程式の特殊解は

$\sum _{i=1}^{n}\left[\int {\frac {W_{i}(x)}{W(x)}}dx\right]y_{i}(x)$

と書くことができる。

例

特定の二階方程式

方程式

$y''+4y'+4y=\cosh {x}$

を解くことを考える。一般解を求めるために、斉次方程式

$y''+4y'+4y=0$

を解くと、この固有多項式は

$\lambda ^{2}+4\lambda +4=(\lambda +2)^{2}$

で、固有値 −2 は重根であるから、u₁ = e^−2x および u₂ = xe^−2x が基本解となる。これらのロンスキー行列式は

${\begin{vmatrix}e^{-2x}&xe^{-2x}\\-2e^{-2x}&-e^{-2x}(2x-1)\end{vmatrix}}=e^{-4x}$

である。これは 0 でないから、この二つの函数は確かに斉次方程式の一般解を生成する。

従って、A(x)u₁ + B(x)u₂ が非斉次方程式の一般解となるような A(x), B(x) を求めればよいが、それには積分

$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$

を計算すればよい。結局、

${\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}}$

が求まる。ただし、C₁, C₂ は積分定数である。

一般の二階方程式

微分方程式

$u''+p(x)u'+q(x)u=f(x)$

を解くにあたって、D を微分演算子として線型微分作用素

$L=D^{2}+p(x)D+q(x)$

を定義すると、L および f(x) が既知として、方程式 Lu = f を u に関して解けばよい、ということになる。

定数変化法を用いるために、まずは対応する斉次方程式

$u''+p(x)u'+q(x)u=0$

を解かねばならない。この方程式は二階であるから、線型独立な二つの解 u₁, u₂ が得られれば、定数変化法を適用することができる。

求める微分方程式の一般解 u_G は

$u_{G}(x)=A(x)u_{1}(x)+B(x)u_{2}(x)$

の形をしているはずである。ただし、A(x), B(x) は未知で、u₁(x), u₂(x) は斉次方程式の解である。A(x) と B(x) がともに定数ならば Lu_G = 0 となるのは明らかである。A = A(x), B = B(x) は

$A'(x)u_{1}(x)+B'(x)u_{2}(x)=0$

となるものと仮定すると、

$u_{G}'(x)=A(x)u_{1}'(x)+B(x)u_{2}'(x)$

となり、さらに微分して

$u_{G}''(x)=A(x)u_{1}''(x)+B(x)u_{2}''(x)+A'(x)u_{1}'(x)+B'(x)u_{2}'(x)$

を得る。従って、L の u_G への作用は

$Lu_{G}=A(x)Lu_{1}(x)+B(x)Lu_{2}(x)+A'(x)u_{1}'(x)+B'(x)u_{2}'(x)$

と書くことができるが、u₁ と u₂ は斉次方程式の解だから

$Lu_{G}(=f)=A'(x)u_{1}'(x)+B'(x)u_{2}'(x)$

となる。

以上から連立方程式

${\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}}$

が得られたので、A(x), B(x) を求めるために、これを A′, B′ について解くと

${\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}}$

を得る。ただし W は u₁ と u₂ のロンスキー行列式である（u₁ と u₂ が線型独立であるという仮定から W が 0 でないことは分かっている）。ゆえに

$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$

を得る。

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

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

参考文献

• Coddington, Earl A.; Levinson, Norman (1955). Theory of Ordinary Differential Equations. New York: McGraw-Hill. 
• Boyce, W. E.; DiPrima, R. C. (1965). Elementary Differential Equations and Boundary Value Problems 8th Edition. Wiley Interscience. , pages 186-192, 237-241
• Teschl, Gerald. Ordinary Differential Equations and Dynamical Systems. Providence: American Mathematical Society. http://www.mat.univie.ac.at/~gerald/ftp/book-ode/. 

外部リンク

• Online Notes / Proof by Paul Dawkins, Lamar University.
• variation of parameters - PlanetMath.（英語）
• Motivation of method via celestial mechanics

Weblio日本語例文用例辞書

「variation of constants」の例文・使い方・用例・文例

• Microsoftがβ版をランチするのは「NetShow streaming server」で動画や音声をオンデマンドで提供する。
• 《主に米国で用いられる》 ＝ 《主に英国で用いられる》 an admiral of the fleet 海軍元帥.
• 篏入的 r 音 《英音の India office /ndiərfɪs/の /r/の音》.
• ＝《口語》 These kind of stamps are rare. この種の[こういう]切手は珍しい.
• (英国の)運輸省. the Ministry of Education(, Science and Culture) (日本の)文部省.
• は of の誤植です.
• を off と誤植する.
• あいまい母音 《about, sofa などの /ə/》.
• 副詞的小詞 《on, in, out, over, off など》.
• 迂言的属格 《語尾変化によらず前置詞によって示す属格; たとえば Caesar's の代わりの of Caesar など》.
• çon of garlic [humor]. それにはガーリック[ユーモア]がちょっぴり必要だ.
• 《主に米国で用いられる》 ＝ 《主に英国で用いられる》 the Speaker of the House of Commons 下院議長.
• 《主に米国で用いられる》 ＝ 《主に英国で用いられる》 the Committee of Ways and Means 歳入委員会.
• 初めて読んだ英文小説は“The Vicar of Wakefield”
• （違法罪―a sin of commission―に対する）怠惰罪
• 『each』、『every』、『either』、『neither』、『none』が分配的、つまり集団の中の1つのものを指すのに対し、『which of the men』の『which』は分離的である
• 『hot off the press（最新情報）』は『hot（最新の）』の拡張感覚を示している
• 『Each made a list of the books that had influenced him』における制限節は、リストに載った本を制限節で定義された特定の本だけに制限する
• 臨床的鬱病を治療するのに用いられる三環系抗鬱薬（商品名ImavateとTofranil）
• 『sunshine-roof』は『sunroof（サンルーフ）』に対する英国の用語である