变系数线性常微分方程的求解
张慧敏，数学计算机科学学院
摘 要：众所周知，所有的常系数一阶、二阶微分方程都是可解的，而变系数
二阶线性微分方程却很难解，至今还没有一个普遍方法。幂级数解法是一个非常有效的方法，本文重点讨论二阶变系数线性常微分方程的解法，从幂级数解法、降阶法、特殊函数法等方面探究了二阶微分方程的解法，简单的介绍了几种高阶微分方程的解法，并讨论了悬链线方程等历史名题。
关键词：变系数线性常微分方程；特殊函数；悬链线方程；幂级数解法
Solving linear ordinary differential equations with variable
coefficients
Huimin Zhang , School of Mathematics and Computer Science
Abstract：As we know, all of ordinary differential equations of first, second order
differential equations with constant coefficients are solvable. However, the linear differential equations of second order with variable coefficients are very difficult to solve. So far there is not a universal method. The method of power-series solution is a very efficient method. This article focuses on solving linear ordinary differential equations of second order with variable coefficients, and exploring the solution of in terms of power-series solution, the method of reducing orders, the method of special functions. Also, this paper applies the above methods to solve several linear differential equations of higher order and especially discusses the famous catenary equation.
Key words：Linear ordinary differential equations with variable coefficients;
Special Functions; catenary equation; Power Series Solution.
I
前言
随着科学的发展和社会的进步，常微分方程在越来越多的领域内有着重要的作用，例如化学，生物学，自动控制，电子技术等，都提出了大量的微分方程问题，同样在社会科学的领域也存在着微分方程问题。
此外，微分方程与数学的其他分支的关系也是非常密切的，他们往往互相联系，互相促进，例如几何学就是常微分方程理论的丰富源泉之一和有力工具，对微分方程的发展产生了深刻的影响。反过来，微分方程进一步发展的需要，也推动着其他数学分支的发展。
众所周知，所有的常系数一阶、二阶微分方程都是可解的，而变系数二阶线性微分方程却很难解，除了近似解法外，至今还没有一个普遍方法。因此，变系数二阶线性微分方程的求解在微分方程理论之中有着十分重要的地位，寻求一种简便的计算方法是完全有必要的。
1
第一部分 二阶线性微分方程的解法探究
一、幂级数解法
⑴一般微分方程的幂级数解法
二阶变系数齐次线性微分方程的求解问题可归结为寻求它的一个非零解。由于方程的系数是自变量的函数，我们不能用之前的代数方法去求解。但是，从微积分学中知道，在满足某些条件下，可以用幂级数来表示一个函数。因此，自然想到，能否用幂级数来表示微分方程的解，下面先以两个例子来探讨一下。
【2】
例1.1.1 求方程y''?2xy'?4y?0的满足初值条件y(0)?0及y'(0)?1的解。
解 设
y?a0?a1x?a2x2???anxn?? (1.1)
为方程的解。首先，利用初值条件，可以得到a0?0,a1?1，因而 y?x?a2x?a3x?????anx????
23ny'?1?2a2x?3a3x2?????nanxn?1????
n?2 y''?2a2?3?2a3x?????n(n?1)anx????
将y,y',y''的表达式代入原方程，合并x的各同次幂的项，并令各项系数等于零，得到
a2?0,a3?1,a4?0,???,an?因而
2an?2,??? n?11111 a5?,a6?0,a7??,a8?0,a9?,???
2!63!4!最后得 a2k?1?对一切正整数k成立。
将ai(i?0,1,2,???)的值代回（1.1）就得到
111??,a2k?0, k(k?1)!k!x5x2k?1????????? y?x?x?2!k!342kxx2x2?????????)?xe, ?x(1?x?2!k! 2
这就是方程的满足所给初值条件的解。 在上例中方程显然满足定理的条件，系数?2x和?4可看作是在全数轴上收敛的幂级数，故方程的解也在全数轴上收敛。但有些方程却未必，例如n阶贝塞尔方程
d2ydy22?x?x?ny?0 x2dxdx2??(1.2)
1n2这里n为非负常数,不一定是正整数。在此p?x??,q?x??1?2.
xx⑵n阶贝塞尔方程
【2】
例1.1.2 求解n阶贝塞尔方程（1.2）。 解 将方程改写成
d2y1dyx2?n2?y?0, 2?dxxdxx2易见，xp?x??1,x2q?x??x2?n2,按x展成的幂级数收敛区间为???x???,从而方程有形如
y??akx??k
k?0?(1.3)
的解，这里a0?0，而ak和?是待定常数。将（1.3）代入（1.2）中，得
x2????k????k?1?axkk?1????k?2?x????k?akxk?1???k?1?x?n?22??axk???k?0，
k?0把x同次幂项归在一起，上式变为
????k????k?1?????k??nakx2k?0????k??akx??k?2?0.
k?0?令各项的系数等于零，得一系列的代数方程
?a0?2?n2?0,??2 ?a1???1??n2?0,
?22??a??k?n?ak?2?0,k?2,3,??k??????(1.4)
?因为a0?0,故从（1.4）的第一个方程解得?的两个值??n和???n. 先考虑??n时方程（1.2）的一个特解。这时我们总可以从（1.4）中逐个地确定所有的系数ak.把??n代入（1.4），得到a1?0,ak??或按下标为奇数或偶数，我们分别有
3
ak?2,k?2n?k?k?2,3,???
变系数线性常微分方程的求解
