求常系数线性常微分方程特解的有限递推法.docVIP

求一类常系数线性常微分方程特解的 有限递推法* 方有康** 北京航空航天大学北海学院 Email: fykfyk2004@ *收入教育数学深圳会议论文集；被写入教材《高等数学新讲》(2008.10)；发表于《数学的认识与实践》(2009.09） **作者简介：教育数学学会个人常务理事。北京航空航天大学北海学院教授。 德国Albert-Ludwigs-University(Freiburg)大学应用数学博士。 摘 要 : 对于非齐次项为多项式，指数函数，正(余)弦函数，或它们的乘积形式的常系数线性常微分方程，本文提出了求其特解的有限递推法.它方法统一，计算简洁，便于编程，能解决高阶问题,能在有限步内得出方程的解析特解，因而优于目前广泛采用的待定系数法. 关键词：常系数线性常微分方程；递推法；待定系数法 MR(2000) 主题分类 11D41 中图分类 O175.11 问题的提出 本文求如下的n阶常系数线性常微分方程的特解: SKIPIF 1 0 (1.1) SKIPIF 1 0 SKIPIF 1 0 SKIPIF 1 0 (1.2) SKIPIF 1 0 式中 SKIPIF 1 0 ; SKIPIF 1 0 =1; SKIPIF 1 0 ∈R; SKIPIF 1 0 ∈c; SKIPIF 1 0 . 在科学技术上,这是一类很重要的微分方程. 目前国内外的高等数学,工程数学或常微分方程的教科书中都采用待定系数法来求这类微分方程的特解.这需要首先求出所对应的齐次方程的特征根（因为要确定特征方程各特征根的重数），再按 SKIPIF 1 0 是特征根的重数 SKIPIF 1 0 （当(1.1)或(1.2)的右端不含三角函数时）设特解为 SKIPIF 1 0 ； 当(1.1)或(1.2)的右端含三角函数时，按 SKIPIF 1 0 是特征根的重数 SKIPIF 1 0 设特解为 SKIPIF 1 0 . 然后计算的各阶导数，再把它们代回原方程中, 比较方程两边同类项的系数，得出一个线性方程组,然后再解这个线性方程组才能求出的各项待定系数.整个求解过程繁琐，计算量大.特别是方程的右端是高次多项式和三角函数及指数函数的乘积时，求导过程中的各阶导数表达式的项数将以几何级数的方式急剧增加，得出的线性方程组也是一个大型的线性方程组.当方程阶数较大时，待定系数法更是显得无能为力.本文提出的有限递推法无需先求出所对应的齐次方程的特征根，无需设定特解的形式且具有计算简洁，方法统一，便于编程，能解决高阶问题和能在有限步内得出方程的解析特解的特点,很好地解决了待定系数法所遇到的困难. 2．主要成果 让我们首先来解决方程右边仅为多项式的情况， 设方程 SKIPIF 1 0 (2.1) 式中为 SKIPIF 1 0 的次多项式.不失一般性,我们规定上式中 SKIPIF 1 0 ≠0.因为若不然,我们可以令方程左边非零最低阶项为新的 SKIPIF 1 0 SKIPIF 1 0 ,其余各项(可以为零)按 SKIPIF 1 0 的导数的阶数由低到高分别为 SKIPIF 1 0 ,且方程右边不变，对这样得出的新方程我们称之为原方程的降阶方程.求出降阶方程的特解再积分,就得出原方程的特解.例如要求 SKIPIF 1 0 的特解,我们先求出 SKIPIF 1 0 的特解,再对其积分就得出原方程的特解了. 定理1 令方程(2.1)左端导数阶数最小(非零)项为 SKIPIF 1 0 记 SKIPIF 1 0 ,则当 SKIPIF 1 0 时(当 SKIPIF 1 0 时,方程的特解显然为 SKIPIF 1 0 ), 方程的特解 SKIPIF 1 0 可由如下的递推公式最多在 SKIPIF 1 0 步内推出: SKIPIF 1 0 , (2.2) 对 SKI