主題
Search

尤拉微分方程

一般的非齊次微分方程由下式給出

x^2(d^2y)/(dx^2)+alphax(dy)/(dx)+betay=S(x),
(1)

齊次方程為

x^2y^('')+alphaxy^'+betay=0
(2)
y^('')+alpha/xy^'+beta/(x^2)y=0.
(3)

現在嘗試將方程從

y^('')+p(x)y^'+q(x)y=0
(4)

轉換為常係數方程

(d^2y)/(dz^2)+A(dy)/(dz)+By=0
(5)

透過使用線性二階常微分方程的標準變換。比較 (3) 和 (5)，函式 p(x) 和 q(x) 是

p(x)=alpha/x=alphax^(-1)
(6)
q(x)=beta/(x^2)=betax^(-2).
(7)

令 B=beta 並定義

z=B^(-1/2)intsqrt(q(x))dx
(8)
=beta^(-1/2)intsqrt(betax^(-2))dx
(9)
=intx^(-1)dx
(10)
=lnx.
(11)

那麼 A 由下式給出

A=(q^'(x)+2p(x)q(x))/(2[q(x)]^(3/2))B^(1/2)
(12)
=(-2betax^(-3)+2(alphax^(-1))(betax^(-2)))/(2(betax^(-2))^(3/2))beta^(1/2)
(13)
=alpha-1,
(14)

這是一個常數。因此，該方程變為具有常係數的二階常微分方程

(d^2y)/(dz^2)+(alpha-1)(dy)/(dz)+betay=0.
(15)

定義

r_1=1/2(-A+sqrt(A^2-4B))
(16)
=1/2[1-alpha+sqrt((alpha-1)^2-4beta)]
(17)
r_2=1/2(-A-sqrt(A^2-4B))
(18)
=1/2[1-alpha-sqrt((alpha-1)^2-4beta)]
(19)

和

a=1/2(1-alpha)
(20)
b=1/2sqrt(4beta-(alpha-1)^2).
(21)

解是

y={c_1e^(r_1z)+c_2e^(r_2z) (alpha-1)^2>4beta; (c_1+c_2z)e^(az) (alpha-1)^2=4beta; e^(az)[c_1cos(bz)+c_2sin(bz)] (alpha-1)^2<4beta.
(22)

用原始變數 x 表示，

y={c_1|x|^(r_1)+c_2|x|^(r_2) (alpha-1)^2>4beta; (c_1+c_2ln|x|)|x|^a (alpha-1)^2=4beta; |x|^a[c_1cos(bln|x|)+c_2sin(bln|x|)] (alpha-1)^2<4beta.
(23)

Zwillinger (1997, p. 120) 給出了另外兩種型別的方程，也稱為尤拉微分方程，

y^'=+/-sqrt((ay^4+by^3+cy^2+dy+e)/(ax^4+bx^3+cx^2+dx+e))
(24)

(Valiron 1950, p. 201) 和

y^'+y^2=alphax^m
(25)

(Valiron 1950, p. 212)，後者可以用貝塞爾函式求解。

另請參閱

尤拉無粘性運動方程

使用 探索

參考文獻

Valiron, G. The Geometric Theory of Ordinary Differential Equations and Algebraic Functions. Brookline, MA: Math. Sci. Press, 1950.Zwillinger, D. Handbook of Differential Equations, 3rd ed. Boston, MA: Academic Press, p. 120, 1997.

在 中引用

尤拉微分方程

請這樣引用

Weisstein, Eric W. "Euler Differential Equation." 來自 Web 資源. https://mathworld.tw/EulerDifferentialEquation.html

主題分類