応用数学 - 定数係数2階線形常微分方程式

概要

線形微分方程式とは、構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle f_1(x), \ f_2(x), \ \cdots, \ f_n(x), \ g(x)} をxのみの関数とする。
この時、以下の形の微分方程式をN階線形常微分方程式という。

${\frac {d^{n}y}{dx^{n}}}+f_{1}(x){\frac {d^{(n-1)}y}{dx^{(n-1)}}}+\cdots +f_{(n-1)}(x){\frac {dy}{dx}}+f_{n}(x)y=g(x)\quad \cdots \quad (1)$

上式において、$g(x)=0$ の時を同次方程式、$g(x)\neq 0$ の時を非同次方程式と呼ぶ。

線形常微分方程式には、連立1次方程式の解の構造との類似性があり、線形代数とのアナロジーがある。
定数係数2階線形常微分方程式は、電磁気学、動力学、量子力学、振動現象等の記述に現れる。

解の存在と一意性

定理（存在定理） :
xのみの関数 $f_{1}(x),\ f_{2}(x),\ \cdots \ f_{n}(x),\ g(x)$ が、区間Iで連続とする。
この時、I内の点 $x=a$ における以下の初期条件のもとで、以下のN階線形常微分方程式の解は区間Iでただ1つ存在する。

初期条件 ：
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle y(a) = b_0, \quad \frac{dy(a)}{dx} = b_1, \quad \cdots , \quad \frac{d^{(n - 1)} y(a)}{dx^{(n - 1)}} = b_{(n - 1)}}


N階線形常微分方程式 ：
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \frac{d^n y}{dx^n} + f_1(x) \frac{d^{(n - 1)} y}{dx^{(n - 1)}} + \cdots + f_{(n - 1)}(x) \frac{dy}{dx} + f_n(x)y = g(x) \quad \cdots \quad (1)}

存在定理により、N階線形常微分方程式の解の存在が保障され、当該初期条件を満たすものがただ1つに決まる。
1階常微分方程式の初期条件は1つであり、N階常微分方程式では、N個の初期条件がないと解は一意に定まらない。

存在定理により、当該問題を解くための数値解析プログラムを作成することに対する正当性が保障される。

定数係数2階線形常微分方程式(同次方程式)の解全体の集合

以下の定数係数2階線形常微分方程式(同次方程式)について考える。
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \frac{d^2y}{dx^2} + f_1(x)\frac{dy}{dx} + f_2(x)y = 0 \quad \cdots \quad (*)}

上式の解関数全体の集合をVとする。
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle V = \left \{ y \quad \Big \vert \quad \frac{d^2y}{dx^2} + f_1(x)\frac{dy}{dx} + f_2(x)y = 0 \right \}}

次に、集合Vは、実数上のベクトル空間（2次元線形空間）となること（すなわち、Vは平面ベクトル全体と同じ構造を持つこと）を記述する。

定理1 ：
解関数 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle y1, \, y2 \, \in V}
 を考える。
この時、次の関数u, vも集合Vの要素である。
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle (1) \ u = y_1 + y_2}

構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle (2) \ v = ky_1 \quad (k : \mbox{ 実 数 } )}

定理2 （重ねあわせの原理） :
同次方程式 : 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \frac{d^2y}{dx^2} + f_1(x)\frac{dy}{dx} + f_2(x)y = 0}
 の解である2つの関数y1, y2を考える。
この時、次の形も同次方程式の解である。
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle y = C_1 y_1 + C_2 y_2 \qquad (C1, \ C2 : \mbox{ 実 数 } )}

定理1の（1）は、ベクトルにおける和に関する性質（ベクトルの和はベクトルである）に対応する。
定理1の（2）は、ベクトルのスカラー倍の性質に対応する。

定理2は、定理1から導かれる。
定理2は、線形常微分方程式が「線形」と呼ばれる由来である。

定理1および定理2により、同次方程式の解はベクトル（線形代数）とのアナロジーで扱うことができるということが分かる。

定理3 ：
定数係数2階線形常微分方程式(同次方程式) 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \frac{d^2 y}{dx^2} + f_1(x) \frac{dy}{dx} + f_2(x) y = 0}
 の解全体の集合Vは、実数上の線形空間を定める。

定理3は、同次方程式の解全体の集合は、線形空間（ベクトル空間）を定める。
すなわち、同次方程式の解関数は、ベクトルと同様に扱うことができる。（定理1を参照）

同次方程式の解は、線形代数とのアナロジーで扱うことができる。

解関数の線形独立性

定義 ：
定数係数2階線形常微分方程式 (同次方程式) 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \frac{d^2 y}{dx^2} + f_1(x) \frac{dy}{dx} + f_2(x)y = 0}
 の解となる関数y1, y2が、
ある区間Iにおいて、少なくとも1つは0でない実数k1, k2を用いて以下のように記述できる時、関数y1とy2はIで線形従属（1次従属）であるという。

また、2つの関数y1, y2が線形従属でない時、線形独立（1次独立）であるという。

構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle k_1 y_1 + k_2 y_2 = 0}

上の定義は、ベクトルの線形従属・線形独立の定義とのアナロジーである。

線形同次微分方程式の解関数の線形従属性・線形独立性を判定するには、次の定理を使用する。

定理4 ：
定数係数2階線形常微分方程式 (同次方程式) 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \frac{d^2 y}{dx^2} + f_1(x) \frac{dy}{dx} + f_2(x)y = 0}
 の解関数y1, y2について、
次のロンスキー行列式 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle W [y_1, y_2]}
 を考える。

この時、考察中の区間Iにおいて以下が成立する。

構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle y_1, \, y_2 \mbox{ が 線 形 独 立 } \iff W[y_1, \, y_2] \ne 0}


構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle W[y_1, \, y_2] = \begin{vmatrix} y_1 & y_2 \\ \frac{dy_1}{dx} & \frac{dy_2}{dx} \end{vmatrix}}

定理4により、ロンスキー行列式を使用すると、簡単な演算のみで線形独立性が判定できる。
(1) 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle y_1, \, y_2 \mbox{ が 線 形 独 立 } \iff W[y_1, \, y_2] \ne 0}
(2) 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle y_1, \, y_2 \mbox{ が 線 形 従 属 } \iff W[y_1, \, y_2] = 0}

同次方程式の解

定理5 ：
集合Vを定数係数2階線形常微分方程式 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \frac{d^2 y}{dx^2} + f_1(x) \frac{dy}{dx} + f_2(x) y = 0}
 の解全体の集合とする。
この時、線形独立な2つの関数 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle y_1,\, y_2 \, \in \, V}
 が存在し、任意の関数 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle y \, \in \, V}
 は、y1とy2の線形結合 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle y = C_1 y_1 + C_2 y_2 \quad (C_1, \, C_2 : \mbox{ 実 数 } )}
で、ただ1通りに表される。

上記のような線形独立な解の組 (y1, y2) を基本解という。
基本解は、線形代数でいう基底ベクトルに対応する。

定理の意味 ：
基本解 y1, y2が見つかれば、他の全ての解（一般解）yもそれらを使って、以下のように表すことができる。
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle y = C_1 y_1 + C_2 y_2 (C_1, C_2 : \mbox{ 実 数 } )}

基本解 (y1, y2) で全ての解関数yを表現できる。

非同次方程式の解

定理6 ：
定数係数2階線形常微分方程式 (非同次方程式) 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \frac{d^2y}{dx^2} + f_1(x) \frac{dy}{dx} + f_2(x) y = g(x) \quad \cdots \quad (1)}
 の1つの特殊解をvとする。
この時、上式(1)の解yは、定数係数2階線形常微分方程式 (同次方程式) 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \frac{d^2y}{dx^2} + f_1(x) \frac{dy}{dx} + f_2(x) y = 0 \quad \cdots \quad (2)}
 の基本解 y1, y2を使用して、
次式のように記述できる。
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle y = (C_1 y_1 + C_2 y_2) + v \quad (C_1, C_2 : \mbox{ 実 数 } )}