主題
Search

引數變易法


對於一個二階常微分方程

 y^('')+p(x)y^'+q(x)y=g(x).
(1)

假設線性獨立解 y_1(x)y_2(x) 已知齊次方程

 y^('')+p(x)y^'+q(x)y=0,
(2)

並尋求 v_1(x)v_2(x) 使得

y^*=v_1y_1+v_2y_2
(3)
y^('*)=(v_1^'y_1+v_2^'y_2)+(v_1y_1^'+v_2y_2^').
(4)

現在,施加附加條件

 v_1^'y_1+v_2^'y_2=0
(5)

因此

y^('*)(x)=v_1y_1^'+v_2y_2^'
(6)
y^(''*)(x)=v_1^'y_1^'+v_2^'y_2^'+v_1y_1^('')+v_2y_2^('').
(7)

y^*y^*^'y^*^('') 代回原方程以獲得

 v_1(y_1^('')+py_1^'+qy_1)+v_2(y_2^('')+py_2^'+qy_2)+v_1^'y_1^'+v_2^'y_2^'=g(x),
(8)

簡化為

 v_1^'y_1^'+v_2^'y_2^'=g(x).
(9)

結合方程 (◇) 和 (9) 並同時求解 v_1^'v_2^',然後得到

v_1^'=-(y_2g(x))/(W(x))
(10)
v_2^'=(y_1g(x))/(W(x)),
(11)

其中

 W(y_1,y_2)=W(x)=y_1y_2^'-y_2y_1^'
(12)

朗斯基行列式,它是僅 x 的函式,因此可以直接積分得到

v_1=-int(y_2g(x))/(W(x))dx
(13)
v_2=int(y_1g(x))/(W(x))dx,
(14)

可以代入以給出特解

 y^*=v_1y_1+v_2y_2.
(15)

推廣到 n 階 ODE,設 y_1, ..., y_n 為齊次 ODE 的解,並選擇 v_1^'(x), ..., v_n^'(x) 使得

 {y_1v_1^'+y_2v_2^'+...+y_nv_n^'=0; y_1^'v_1^'+y_2^'v_2^'+...+y_n^'v_n^'=0; |; y_1^((n-1))v_1^'+y_2^((n-1))v_2^'+...+y_n^((n-1))v_n^'=g(x).
(16)

則特解為

 y^*(x)=v_1(x)y_1(x)+...+v_n(x)y_n(x),
(17)

另請參閱

常微分方程, 二階常微分方程

使用 探索

請引用為

Weisstein, Eric W. "引數變易法。" 來自 Web 資源。 https://mathworld.tw/VariationofParameters.html

學科分類