二阶线性常微分方程的幂级数解法.docVIP

二阶线性常微分方程的幂级数解法 从微分方程学中知道，在满足某些条件下，可以用幂级数来表示一个函数。因此，自然想到，能否用幂级数来表示微分方程的解呢？ 例1、求方程的通解 解：设 为方程的解，这里是待定常系数，将它对微分两次，有 将，的表达式代入方程，并比较的同次幂的系数，得到 ， 或一般的可推得 ， ， 其中，是任意的，因而代入设的解中可得： 这个幂级数的收敛半径是无限大的，因而级数的和（其中包括两个任意常数及）便是所要求的通解。 求方程的满足初值条件及的解。 解 设级数为方程的解。首先，利用初值条件，可以得到 ， ， 因而 将，，的表达式带入原方程，合并的各同次幂的项，并令各项系数等于零，得到 因而 最后得 , , 对一切正整数成立。 将的值代回就得到 这就是方程的满足所给初值条件的解。 是否所有方程都能按以上方式求出其幂级数解？或者说究竟方程应该满足什么条件才能保证它的解可用幂级数来表示呢？级数的形式怎样？其收敛区间又如何?这些问题，在微分方程解析理论中有完满的解答，但因讨论时需要涉及解析函数等较专门的知识，在此我们仅叙述有关结果而不加证明，若要了解定理的证明过程，可参考有关书籍。 考虑二阶齐次线性微分方程 及初值条件及的情况。 不失一般性，可设 ，否则，我们引进新变量，经此变换，方程的形状不变，在这时对应于的就是了，因此，今后我们总认为。 定理10 若方程中系数和都能展成的幂级数，且收敛区间为，则方程有形如 的特解，也以为级数的收敛区间。 在上两例中方程显然满足定理的条件，系数，和可看作是在全数轴上收敛的幂级数，故方程的解也在全数轴上收敛。但有些方程，例如阶贝赛尔方程 这里为非负常数，不一定是正整数，（）在此，，显然它不满足定理10 的条件，因而不能肯定有形如的特解。但它满足下述定理11的条件，从而具有别种形状的幂级数解。 定理11 若方程中系数，具有这样的性质，即和均能展成的幂级数，且收敛区间为，若，则方程有形如 即 的特解，是一个特定的常数，级数也以为收敛区间。若，或更一般的，，但，则引入记号，，则 ， 这里，而仍为待定常数。 求解阶贝赛尔方程。 解 将方程改写成 ， 易见，它满足定理11的条件（和均能展成的幂级数，且收敛区间为），且，按展成的幂级数收敛区间为，由定理11，方程有形如 的解，这里，而和是待定常数，将代入：中，得 ， 把同幂次项归在一起，上式变为 令各项的系数等于0，得一系列的代数方程 因为，故从解得的两个值 和 先考虑时方程的一个特解，这时我们总可以从以上方程组中逐个地确定所有的系数。把代入以上方程组，得到 ， 或按下标为奇数或偶数，我们分别有 从而求得 一般地 将各代入得到方程的一个解 既然是求的特解，我们不妨令 其中函数定义如下： 当0时，；当0且非整数时，由递推公式定义。 具有性质 ; 为正整数 而变为 注意到函数的性质，即有 是由贝塞尔方程定义的特殊函数，称为阶贝赛尔函数。 因此，对于阶贝塞尔方程，它总有一个特解。为了求得另一个与线性无关的特解，我们自然想到，求时方程的形如 的解，我们注意到只要不为非负整数，像以上对于时的求解过程一样，我们总可以求得 使之满足中的一系列方程，因而 是的一个特解。此时，若令 则变为 称为阶贝赛尔函数。 利用达朗贝尔判别法不难验证级数和（在中）都是收敛的，因此，当不为非负整数时，和都是方程的解，而且是线性无关的，因为它们可展为由的不同幂次开始的级数，从而它们的比不可能是常数。于是方程的通解可写为 这里，是任意常数。此情形的和称为第一类贝塞尔函数。 例8 求方程的通解。 解 引入新变量，我们有 ， 将上述关系代入院方程，得到 ， 这是，的贝塞尔方程，由例7可知，方程的通解可表为 ， 代回原来变量，就得到原方程的通解 其中是任意常数。 第二宇宙速度计算 作为这一节的应用，我们计算发射人造卫星的最小速度，即所谓第二宇宙速度。在这个速度你下,物体将摆脱地球的引力,向地球一样绕着太阳运行,成为人造卫星. 让我们首先建立物体垂直上抛运动的微分方程.以和分别表示地球和物体的质量.按牛顿万有引力定律,作用于物体的引力(空气阻力忽略不计)为 这里表示地球的中心和物理体重心之间的距离，为万有引力常数。因为，物体运动规律应满足下面的微分方程 或 这里的负号表示物体的加速度是负的。 设地球半径为，物理发射速度为，因此，当物体刚刚离开地球表面时，我们有，即应取初值条件为 方程不显含自变量,应用4.3.1（可降阶的一些方程类型）的方法，把方程降阶成为一阶方程 解得 注意到这时初值条件为 因而