主题
Search

二阶常微分方程


形如以下的常微分方程

 y^('')+P(x)y^'+Q(x)y=0.
(1)

当 x=x_0 为有限值时,这类方程在以下条件下具有奇点: (a) 如果当 P(x)Q(x) 发散,但当 x->x_0 时,(x-x_0)P(x)(x-x_0)^2Q(x) 保持有限,则 x_0 称为正则奇点或非本质奇点。(b) 如果 P(x) 的发散速度比 (x-x_0)^(-1) 快,以至于当 x->x_0(x-x_0)P(x)->infty,或者 Q(x) 的发散速度比 (x-x_0)^(-2) 快,以至于当 x->x_0(x-x_0)^2Q(x)->infty,则 x_0 称为非正则奇点或本质奇点。

为了研究方程 (1) 在无穷远处的奇点,我们进行替换 x=z^(-1),因此 dx=-z^(-2)dz,得到

 (dy)/(dx)=-z^2(dy)/(dz)
(2)
(d^2y)/(dx^2)=-z^2d/(dz)(-z^2(dy)/(dz))
(3)
=-z^2(-2z(dy)/(dz)-z^2(d^2y)/(dz^2))
(4)
=2z^3(dy)/(dz)+z^4(d^2y)/(dz^2).
(5)

那么 (3) 变为

 z^4(d^2y)/(dz^2)+[2z^3-z^2P(z^(-1))](dy)/(dz)+Q(z^(-1))y=0.
(6)

情况 (a):如果

alpha(z)=(2z-P(z^(-1)))/(z^2)
(7)
beta(z)=(Q(z^(-1)))/(z^4)
(8)

x=+/-infty (z=0) 处保持有限,则该点是常点。情况 (b):如果 alpha(z) 的发散速度不超过 1/z,或者 beta(z) 的发散速度不超过 1/z^2,则该点是正则奇点。情况 (c):否则,该点是非正则奇点。

Morse 和 Feshbach (1953, pp. 667-674) 给出了按奇点类型分类的二阶常微分方程的规范形式和解法。

对于特殊类型的线性二阶常微分方程,可变系数可以转换为常系数。给定一个具有可变系数的二阶线性 ODE

 (d^2y)/(dx^2)+p(x)(dy)/(dx)+q(x)y=0.
(9)

定义函数 z=y(x),

(dy)/(dx)=(dz)/(dx)(dy)/(dz)
(10)
(d^2y)/(dx^2)=((dz)/(dx))^2(d^2y)/(dz^2)+(d^2z)/(dx^2)(dy)/(dz)
(11)
((dz)/(dx))^2(d^2y)/(dz^2)+[(d^2z)/(dx^2)+p(x)(dz)/(dx)](dy)/(dz)+q(x)y=0
(12)
(d^2y)/(dz^2)+[((d^2z)/(dx^2)+p(x)(dz)/(dx))/(((dz)/(dx))^2)](dy)/(dz)+[(q(x))/(((dz)/(dx))^2)]y
(13)
=(d^2y)/(dz^2)+A(dy)/(dz)+By=0.
(14)

如果 AB 不是 x 的函数,这将具有常系数。但是,对于 q(x)>=0,我们可以通过将 z 定义为以下形式,将 B 设置为任意常数

 z=B^(-1/2)int[q(x)]^(1/2)dx.
(15)

那么

(dz)/(dx)=B^(-1/2)[q(x)]^(1/2)
(16)
(d^2z)/(dx^2)=1/2B^(-1/2)[q(x)]^(-1/2)q^'(x),
(17)

A=(1/2B^(-1/2)[q(x)]^(-1/2)q^'(x)+B^(-1/2)p(x)[q(x)]^(1/2))/(B^(-1)q(x))
(18)
=(q^'(x)+2p(x)q(x))/(2[q(x)]^(3/2))B^(1/2).
(19)

因此,方程 (◇) 变为

 (d^2y)/(dz^2)+(q^'(x)+2p(x)q(x))/(2[q(x)]^(3/2))B^(1/2)(dy)/(dz)+By=0,
(20)

只要满足以下条件,它就具有常系数

 A=(q^'(x)+2p(x)q(x))/(2[q(x)]^(3/2))B^(1/2)=[constant].
(21)

消除常数,得到

 A^'=(q^'(x)+2p(x)q(x))/([q(x)]^(3/2))=[constant].
(22)

因此,对于 A^' 为常数的常微分方程,解可以通过求解具有常系数的二阶线性 ODE 得到

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

对于 z,其中 z 的定义如上。

一般形式的线性二阶齐次微分方程

 y^('')+P(x)y^'+Q(x)y=0
(24)

可以转换为标准形式

 z^('')+q(x)z=0
(25)

通过使用以下替换消除一阶项

 lny=lnz-1/2intP(x)dx.
(26)

那么

(y^')/y=(z^')/z-1/2P(x)
(27)
(yy^('')-y^('2))/(y^2)=(zz^('')-z^('2))/(z^2)-1/2P^'(x)
(28)
(y^(''))/y-((y^')/y)^2=(z^(''))/z-(z^('2))/(z^2)-1/2P^'(x)
(29)
(y^(''))/y=[(z^')/z-1/2P(x)]^2+(z^(''))/z-(z^('2))/(z^2)-1/2P^'(x)
(30)
=(z^('2))/(z^2)-(z^')/zP(x)+1/4P^2(x)+(z^(''))/z-(z^('2))/(z^2)-1/2P^'(x),
(31)

因此

(y^(''))/y+P(x)(y^')/y+Q(x)=-(z^')/zP(x)+1/4P^2(x)+(z^(''))/z-1/2P^'(x)+P(x)[(z^')/z-1/2P(x)]+Q(x)
(32)
=(z^(''))/z-1/2P^'(x)-1/4P^2(x)+Q(x)=0.
(33)

因此,

 z^('')+[Q(x)-1/2P^'(x)-1/4P^2(x)]z=z^('')(x)+q(x)z=0,
(34)

其中

 q(x)=Q(x)-1/2P^'(x)-1/4P^2(x).
(35)

如果 Q(x)=0,则微分方程变为

 y^('')+P(x)y^'=0,
(36)

可以通过乘以以下项来求解

 exp[int^xP(x^')dx^']
(37)

得到

 0=d/(dx){exp[int^xP(x^')dx^'](dy)/(dx)}
(38)
 c_1=exp[int^xP(x^')dx^'](dy)/(dx)
(39)
 y=c_1int^x(dx)/(exp[int^xP(x^')dx^'])+c_2.
(40)

对于非齐次二阶常微分方程,其中 x 项不出现在函数 f(x,y,y^') 中,

 (d^2y)/(dx^2)=f(y,y^'),
(41)

v=y^',则

 (dv)/(dx)=f(v,y)=(dv)/(dy)(dy)/(dx)=v(dv)/(dy).
(42)

因此,一阶 ODE

 v(dv)/(dy)=f(y,v),
(43)

如果是线性的,可以作为线性一阶 ODE 求解 v。一旦解已知,

 (dy)/(dx)=v(y)
(44)
 int(dy)/(v(y))=intdx.
(45)

另一方面,如果 y 从 f(x,y,y^') 中缺失,

 (d^2y)/(dx^2)=f(x,y^'),
(46)

v=y^',则 v^'=y^(''),方程简化为

 v^'=f(x,v),
(47)

如果是线性的,可以作为线性一阶 ODE 求解 v。一旦解已知,

 y=intv(x)dx.
(48)

如果已知齐次版本的通解,则可以求解非齐次常微分方程,在这种情况下,可以使用参数变分法找到特解。特别地,非齐次二阶常微分方程的特解 y^*(x)

 y^('')+p(x)y^'+q(x)y=g(x)
(49)

可以使用参数变分法找到,由以下方程给出

 y^*(x)=-y_1(x)int(y_2(x)g(x))/(W(x))dx+y_2(x)int(y_1(x)g(x))/(W(x))dx,
(50)

其中 y_1(x)y_2(x) 是无源方程的齐次解

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

W(x) 是这两个函数的 Wronskian 行列式


另请参阅

阿贝尔微分方程恒等式, 伴随, 一阶常微分方程, 常微分方程, 二阶常微分方程的第二解, 待定系数法, 参数变分法 在 MathWorld 课堂中探索此主题

使用 Wolfram|Alpha 探索

参考文献

Arfken, G. "A Second Solution." §8.6 in Mathematical Methods for Physicists, 3rd ed. Orlando, FL: Academic Press, pp. 467-480, 1985.Boyce, W. E. and DiPrima, R. C. Elementary Differential Equations and Boundary Value Problems, 4th ed. New York: Wiley, 1986.Morse, P. M. and Feshbach, H. Methods of Theoretical Physics, Part I. New York: McGraw-Hill, pp. 667-674, 1953.

在 Wolfram|Alpha 中被引用

二阶常微分方程

请引用为

Weisstein, Eric W. "Second-Order Ordinary Differential Equation." 来自 MathWorld——Wolfram Web 资源。 https://mathworld.net.cn/Second-OrderOrdinaryDifferentialEquation.html

学科分类