「応用数学 - 定数係数2階線形常微分方程式」の版間の差分
| 105行目: | 105行目: | ||
<br> | <br> | ||
基本解 (y1, y2) で全ての解関数yを表現できる。<br> | 基本解 (y1, y2) で全ての解関数yを表現できる。<br> | ||
<br> | |||
[[ファイル:Solution of Homogeneous Equations 1.png|フレームなし|中央]] | |||
<br><br> | <br><br> | ||
2023年7月22日 (土) 00:08時点における版
概要
線形微分方程式とは、構文解析に失敗 (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階線形常微分方程式という。
構文解析に失敗 (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)}
上式において、構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle g(x) = 0}
の時を同次方程式、構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle g(x) \ne 0}
の時を非同次方程式と呼ぶ。
線形常微分方程式には、連立1次方程式の解の構造との類似性があり、線形代数とのアナロジーがある。
定数係数2階線形常微分方程式は、電磁気学、動力学、量子力学、振動現象等の記述に現れる。
解の存在と一意性
定理(存在定理) :
xのみの関数 構文解析に失敗 (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)}
が、区間Iで連続とする。
この時、I内の点 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle 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階線形常微分方程式(同次方程式)について考える。
上式の解関数全体の集合をVとする。
次に、集合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{ 実 数 } )}
同次方程式の例題
例題1 :
微分方程式 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \frac{d^2y}{dx^2} - 3 \frac{dy}{dx} + 2y = 0 \quad \cdots \quad (1)}
と関数 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle y_1 = e^x, \quad y_2 = e^{2x}}
について以下の問いに答えよ。
(1) 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle y_1, \, y_2}
は、上式(1)の解であることを示せ。
(2) y1とy2は、任意の区間で線形独立であることを示せ。
(3) y1とy2の線形結合の関数 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle y = 2e^x - e^{2x}}
も上式(1)の解であることを示せ。
例題1 (1)の解答
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle y_1 = e^x, y_2 = e^{2x}}
がともに方程式 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \frac{d^2y}{dx^2} - 3 \frac{dy}{dx} + 2y = 0 \quad \cdots \quad (1)}
を満たすことを示せばよい。
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle y_1 = e^x}
について、構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \frac{dy_1}{dx} = e^x, \quad \frac{d^2 y_1}{dx^2} = e^x}
なので、上式(1)の左辺に代入すると、
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \frac{d^2 y_1}{dx^2} - 3 \frac{dy_1}{dx} + 2 y_1 = e^x - 3 e^x + 2 e^x = 0}
ゆえに、構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle y_1 = e^x}
は、上式(1)の解である。
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle y_2 = e^{2x}}
について、 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \frac{dy_2}{dx} = 2 e^{2x}, \quad \frac{d^2 y_2}{dx^2} = 4 e^{2x}}
なので、上式(1)の左辺に代入して、
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \frac{d^2 y_2}{dx^2} - 3 \frac{dy_2}{dx} + 2 y_2 = 4 e^{2x} - 3 \times 2 e^{2x} + 2 e^{2x} = 0}
ゆえに、構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle y_2 = e^{2x}}
は、上式(1)の解である。
例題1 (2)の解答
ロンスキー行列式 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle W [y_1, \, y_2]}
を計算して、線形独立性を判定する。
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \begin{align} W [ y_1, y_2 ] &= \begin{vmatrix} y_1 & y_2 \\ \frac{d y_1}{dx} & \frac{d y_2}{dx} \end{vmatrix} = \begin{vmatrix} e^x & e^{2x} \\ e^x & 2 e^{2x} \end{vmatrix} \\ &= 2 e^{3x} - e^{3x} \\ &= e^{3x} \end{align} }
ゆえに、構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle W [y1, \, y2]}
は、任意の区間で零関数ではないため、y1とy2とは線形独立である。
例題1 (3)の解答
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle y = 2 e^x - e^{2x}}
が方程式 を満たすことを示せばよい。
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \begin{align} \frac{dy}{dx} &= \frac{d}{dx} (2 e^x - e^{2x}) \\ &= 2 e^x - 2 e^{2x} \end{align} }
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \begin{align} \frac{d^2 y}{dx^2} &= \frac{d2}{dx} (2 e^x - 2 e^{2x}) \\ &= 2 e^x - 4 e^{2x} \end{align} }
これらを、上式(1)の左辺に代入すると、
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \begin{align} (2 e^x - 4 e^{2x}) - 3(2 e^x - 2 e^{2x}) + 2(2 e^x - 2 e^{2x}) &= 2 e^x - 4 e^{2x} - 6 e^x + 6 e^{2x} + 4 e^x - 2 e^{2x} \\ &= (2 e^x - 6 e^x + 4 e^x) + (- 4 e^{2x} + 6 e^{2x}- 2 e^{2x}) \\ &= 0 \end{align} }
ゆえに、 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle y = 2 e^x - e^{2x}}
は、上式(1)の解である。
特性方程式
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \frac{d^2 y}{dx^2} + a \frac{dy}{dx} + by = 0 \quad \cdots \quad (1)}
に対して、以下をこの方程式の特性方程式(characteristic equation)という。
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \lambda^2 + a \lambda + b = 0 \quad \cdots \quad (2)}
上式(2)の特性方程式の解λ1, λ2が分かれば、上式(1)の基本解y1, y2を求めることができる。
したがって、定数係数2階線形常微分方程式(同次方程式)の一般解を求めるには、特性方程式を解けばよいことになる。
特性方程式 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \lambda^2 + a \lambda + b = 0}
の解は、判別式 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle D = a^2 - 4b}
の値により以下の3種類が考えられる。
- 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle D > 0}
のとき
- 相異なる2つの実数解をもつ。
- 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \lambda = \alpha, \, \beta}
- 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle D = 0}
のとき
- 1つの実数解(重解)をもつ。
- 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \lambda = \alpha}
- 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle D < 0}
のとき
- 共役な2つの複素数解をもつ。
- 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \lambda = p \pm q i}
定理 :
定数係数2階線形常微分方程式 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \frac{d^2 y}{dx^2} + a \frac{dy}{dx} + by = 0}
の基本解と一般解は、特性方程式の解の種類により以下の3通り(1)(2)(3)で求められる。
(1) 実数解の場合 :
特性方程式: 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \lambda^2 + a \lambda + b = 0}
特性方程式の相異なる2つの実数解: 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \lambda = \alpha, \, \beta}
基本解: 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle (e^{\alpha x}, \, e^{\beta x})}
一般解: 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle y = C_1 e^{\alpha x} + C_2 e^{\beta x} \quad (C_1, \, C_2 : \mbox{ 任 意 定 数 } )}
(2) 実数重解の場合 :
特性方程式: 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \lambda^2 + a \lambda + b = 0}
特性方程式の1つの実数解(重解): 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \lambda = \alpha}
基本解: 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle (e^{\alpha x}, \, x e^{\alpha x})}
一般解: 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \begin{align}y &= C_1 e^{\alpha x} + C_2 x e^{\alpha x} \\ &= (C_1 + C_2 x) e^{\alpha x} \quad (C_1, \, C_2 : \mbox{ 任 意 定 数 }) \end{align}}
(3) 共役複素数解の場合 :
特性方程式: 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \lambda^2 + a \lambda + b = 0}
特性方程式の共役な2つの複素数解: 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \lambda = p \pm qi}
基本解: 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle (e^{px} \cos{qx}, \, e^{px} \sin{qx})}
一般解: 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \begin{align} y &= C_1 e^{px} \cos{qx} + C_2 e^{px} \sin{qx} \\ &= (C_1 \cos{qx} + C_2 \sin{qx}) e^{px} \quad (C_1, \, C_2 : \mbox{ 任 意 定 数 }) \end{align}}
定数係数2階線形常微分方程式 (同次方程式)の例題
例題 1 :
以下の微分方程式の一般解を求めよ。
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \frac{d^2 y}{dx^2} - 5 \frac{dy}{dx} + 6y = 0}
解答 :
上式の特性方程式は、構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \lambda^2 - 5 \lambda + 6 = 0}
より、特性方程式の解は 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \lambda = 2, \, 3}
(2つの実数解)
したがって、基本解の組は 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle (e^{2x}, \, e^{3x})}
である。
よって、一般解は次式となる。
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle y = C_1 e^{2x} + C_2 e^{3x} \quad (C_1, \, C_2 : \mbox{ 任 意 定 数 })}
定数係数2階線形常微分方程式 (非同次方程式)
定理 :
定数係数2階線形微分方程式 (非同次方程式) 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \frac{d^2 y}{dx^2} + a \frac{dy}{dx} + by = g(x)}
の一般解yは、次式で求めることができる。
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle y = \mbox{ 同 次 方 程 式 の 一 般 解 } + \mbox{ 非 同 次 方 程 式 の 特 殊 解 }}
定数係数2階線形常微分方程式(非同次方程式)の特殊解の求め方として、以下の3つがある。
(1) 定数変化法(variation of constants)
(2) 未定係数法(method of undetermined coefficients)
(3) 記号法(symbolic method)
定数変化法
以下の定理の証明を与える方法を定数変化法という。
定理 :
定数係数2階線形常微分微分方程式 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \frac{d^2 y}{dx^2} + a \frac{dy}{dx} + by = 0}
の基本解を(y1, y2)とする時、
以下のv(x)は、定数係数2階線形常微分方程式 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \frac{d^2 y}{dx^2} + a \frac{dy}{dx} + by = g(x)}
の特殊解である。
ここで、構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle W [y_1, \, y_2]}
は、ロンスキー行列式であり、構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \int}
は原始関数を表す。
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle v(x) = - y_1 \int{ \frac{y_2 \, g(x)}{W [ y1, \, y2 ]}} dx + y_2 \int{ \frac{y_1 \, g(x)}{W [ y_1, \, y_2 ]}} dx}
定数変化法では、g(x)が考察中の区間Iで連続であれば、どんな関数に対しても特殊解v(x)を求めることができる。
ただし、g(x)によっては、原始関数をよく知られた関数で表せない場合もある。
定数変化法の例題1 :
以下の微分方程式の特殊解v(x)を定数変化法で求め、一般解を求めよ。
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \frac{d^2 y}{dx^2} + 3 \frac{dy}{dx} + 2y = x}
同次方程式 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \frac{d^2 y}{dx^2} + 3 \frac{dy}{dx} + 2y = 0}
の基本解は、構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle (e^{-2x}, \, e^{-x})}
、
一般解は、構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle y = C_1 e^{-2x} + C_2 e^{-x} \quad \cdots \quad (C_1,\, C_2 : \mbox{ 任 意 定 数 })}
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle y_1 = e^{-2x}, \, y_2 = e^{-x}, \, g(x) = x}
とおいて、定数変化法の公式を用いて特殊解v(x)を求める。
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \begin{align} W [ y_1, \, y_2 ] &= \begin{vmatrix} e^{-2x} & e^{-x} \\ -2 e^{-2x} & -e^{-x} \end{vmatrix} \\ &= -e^{-3x} + 2e^{-3x} \\ &= e^{-3x} \end{align} }
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \begin{align} v(x) &= - y_1 \int{ \frac{y_2 \, g(x)}{W [ y1, \, y2 ]}} dx + y_2 \int{ \frac{y_1 \, g(x)}{W [ y_1, \, y_2 ]}} dx \\ &= -e^{-2x} \int{\frac{e^{-x} \, x}{e^{-3x}}} dx + e^{-x} \int{\frac{e^{-2x} \, x}{e^{-3x}}} dx \\ &= -e^{-2x} \int{x e^{2x}} dx + e^{-x} \int{x e^x} dx \\ &= -e^{-2x} \left ( \frac{1}{2} x e^{2x} - \frac{1}{2} \int{e^{2x}} \right ) dx + e^{-x} \left ( x e^x - \int{e^x} \right ) dx \\ &= - \frac{1}{2} x + \frac{1}{4} + x - 1 \\ &= \frac{1}{2} x - \frac{3}{4} \\ &= \frac{2x - 3}{4} \end{align} }
したがって、与式の一般解は次式となる。
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle y = C_1 e^{-2x} + C_2 e^{-x} + \frac{2x - 3}{4}}
未定係数法
未定係数法とは
定数係数2階線形常微分方程式 (非同次方程式) 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \frac{d^2 y}{dx^2} + a \frac{dy}{dx} + by = g(x)}
において、
g(x)が特定の形(例えば、xの多項式やeaxの形)の場合、特殊解v(x)の形を推定できる場合がある。
この場合、v(x)の形は、以下の定数係数2階線形常微分方程式 (同次方程式)の特性方程式の解とg(x)の形で決まる。
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \frac{d^2 y}{dx^2} + a \frac{dy}{dx} + by = 0}
このように、g(x)が特定の形の場合に有効な方法が未定係数法である。
ただし、未定係数法は万能な方法ではなく、
定数係数2階線形常微分方程式 (非同次方程式)において、g(x)が以下のような形の場合に有効な方法である。
(1) xの多項式の形。
(2) eaxの形。
(3) 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \sin{\alpha x}}
の形。
(4) の形。
(5) 上記(1)〜(4)の組み合わせでできている形。
以下、上記の場合を考慮して、g(x)の形、特性方程式の解、および特殊解v(x)の形の間の関係を記述する。
特殊解の形
多項式 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle p_n, \, P_n, \, Q_n}
を定める。
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \begin{cases} p_n(x) &= k_n x^n + k_{n - 1} x^{n - 1} + \cdots + k_1 x + k_0 \\ P_n(x) &= A_n x^n + A_{n - 1} x^{n - 1} + \cdots + A_1 x + A_0 \\ Q_n(x) &= B_n x^n + B_{n - 1} x^{n - 1} + \cdots + B_1 x + B_0 \\ \end{cases} }
( 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle n = 0}
の時は、構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle p_n(x), \, P_n(x), \, Q_n(x)}
は定数となる)
このとき、定数係数2階線形常微分方程式 (非同次方程式): 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \frac{d^2 y}{dx^2} + a \frac{dy}{dx} + by = g(x)}
におけるg(x)の形に応じて、
4通りに分類して、特殊解v(x)の形を見る。
- 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle g(x) = p_n(x)}
の場合の特殊解v(x)
- 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \lambda \ne 0}
の時
- 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle v(x) = P_n(x)}
- 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \lambda = 0 \mbox{ ( 単 解 ) }}
の時
- 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle v(x) = x P_n(x)}
- 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \lambda = 0 \mbox{ ( 2 重 解 ) }}
の時
- 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle v(x) = x^2 P_n(x)}
- 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \lambda \ne 0}
の時
- 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle g(x) = p_n(x) e^{\alpha x}}
の場合の特殊解v(x)
- 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \lambda \ne \alpha}
の時
- 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle v(x) = P_n(x) e^{\alpha x}}
- 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \lambda = 0 \mbox{ ( 単 解 ) }}
の時
- 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle v(x) = x P_n(x) e^{\alpha x}}
- 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \lambda = 0 \mbox{ ( 2 重 解 ) }}
の時
- 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle v(x) = x^2 P_n(x) e^{\alpha x}}
- 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \lambda \ne \alpha}
の時
- 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle g(x) = p_n(x) \sin{\beta x}}
または 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle p_n(x) \cos{\beta x}}
の場合
- 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \lambda \ne i \beta}
の時
- 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle v(x) = Pn(x) \cos{\beta x} + Q_n(x) \sin{\beta x}}
- 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \lambda = i \beta}
の時
- 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle v(x) = x \{ Pn(x) \cos{\beta x} + Q_n(x) \sin{\beta x} \}}
- 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \lambda \ne i \beta}
の時
- 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle g(x) = p_n(x) e^{\alpha x} \sin{\beta x}}
または 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle g(x) = p_n(x) e^{\alpha x} \cos{\beta x}}
の場合
- 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \lambda \ne \alpha \pm i \beta}
の時
- 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle v(x) = e^{\alpha x} \{ P_n(x) \cos{\beta x} + Q_n(x) \sin{\beta x} \}}
- 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \lambda = \alpha \pm i \beta}
の時
- 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle v(x) = x e^{\alpha x} \{ P_n(x) \cos{\beta x} + Q_n(x) \sin{\beta x} \}}
- 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \lambda \ne \alpha \pm i \beta}
の時
未定係数法の例題
例題1 :
以下の微分方程式の特殊解v(x)を未定係数法で求め、一般解を求めよ。
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \frac{d^2 y}{dx^2} + 3 \frac{dy}{dx} + 2y = x}
※ 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle g(x) = p_n(x) = k_1 x + k_0 = x}
の形なので、構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle v(x) = P_n(x) = A_1 x + A0}
の形になる。
上式の特性方程式は、構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \lambda^2 + 3 \lambda + 2 = 0}
より、構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \lambda = -2, \, -1}
したがって、上式の基本解は 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle ( e^{-2x}, \, e^{-x} )}
となり、一般解は 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle y = C_1 e^{-2x} + C_2 e^{-x} \quad (C_1, \, C_2 : \mbox{ 任 意 定 数 } )}
次に、非同次方程式 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \frac{d^2 y}{dx^2} + 3 \frac{dy}{dx} + 2y = x}
の特殊解v(x)を求める。
v(x)の形は、g(x)がpn(x)の場合の特殊解v(x)はPn(x)になることを考慮すると、
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle v(x) = A_1 x + A_0 \quad (n = 1 \mbox{ の ケ ー ス } )}
の形になることがわかる。
以下、このA1、および、A0を求める。
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle v(x) = A_1 x + A_0, \quad \frac{dv(x)}{dx} = A_1, \quad \frac{d^2 v(x)}{dx} = 0}
これらを与式 構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \frac{d^2 y}{dx^2} + 3 \frac{dy}{dx} + 2y = x}
の左辺へ代入すると、
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \begin{align} \frac{d^2 v(x)}{dx^2} + 3 \frac{d v(x)}{dx} + 2v(x) &= 0 + 3 A_1 + 2(A_1 x + A_0) \\ &= (2 A_1)x + (3 A_1 + 2 A_0) \\ \end{align} }
これが右辺の多項式xに等しくなるためには、構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle 2 A_1 = 1, \quad 3 A_1 + 2 A_0 = 0}
となる必要がある。
よって、構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle A_0 = - \frac{3}{4}, \quad A_1 = \frac{1}{2}}
したがって、特殊解は以下になる。
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \begin{align} v(x) &= \frac{1}{2} x - \frac{3}{4} \\ &= \frac{2x - 3}{4} \end{align} }
以上より、与式の一般解は以下になる。
構文解析に失敗 (SVG (ブラウザーのプラグインで MathML を有効にできます): サーバー「https://wikimedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle y = C_1 e^{-2x} + C_2 e^{-x} + \frac{2x - 3}{4}}