特征函数在高阶常微分方程特解计算中的应用.docVIP

2008年 11月 吉林师范大学学报 (自然科学版 . 4第 4期 Journal of Jilin Normal University (NaturalScience Edition Nov. 2008 收稿日期 :2008 08 04 基金项目 :国家自然科学基金对外交流与合作项目 (40640420072 第一作者简介 :龚东山 (1969 , 男 , 湖北监利人 , 现为兰州大学数学与统计学院讲师 , 博士 . 研究方向 :应用数学科研与教学工作 . 特征函数在高阶常微分方程特解计算中的应用 龚东山 1, 2, 刘岳巍 1, 牛富俊 2 (1. 兰州大学 数学与统计学院 , 兰州 730000; 2. 中国科学院寒区旱区环境与工程研究所 冻土工程国家重点实验室 , 兰州 730000 摘 要 :通过借助特征函数的导数 , 得到了非齐次项为特殊函数情形的一类高阶常微分方程的一个 特解的一种 新的计算方法 . 运用该方法 , 还得到了非齐次项为常见情形时方程的一个特解 . 关键词 :微分方程 ; 特征函数 ; 非齐次 ; 特解 ; 应用 中图分类号 :O 175. 1 文献标识码 :A 文章编号 :1674 3873 (2008 04 0008 031 引言 常系数线性微分方程的一般形式是 y (n +a 1y (n -1 + +a n -1y ! +a n y =g (x (1 其中 a i (i =1, 2, , n 为常数 , g (x 为连续函数 . 若 g (x ? 0, 方程 (1 变成 y (n +a 1y (n -1 + +a n -1y ! +a n y =0(2 称方程 (2 是方程 (1 的齐次方程 . 否则 , 称方程 (1 为非齐次方程 . 依据常系数线性微分方程的通解结构理论 [1] , 知方程 (1 的通解可表示成方程 (2 的通解与 方程 (1 的一个特解之和 . 对于方程 (2 的通解 , 一般可采用 Euler 待定指 数函数法 [1]得到 , 故方程 (1 的通解计算的难点在于 寻找方程 (1 的一个特解 . 常用寻找特解的方法有 :待定系数法、 常数变易法 [2]、 积分法 [3]、 拉普拉斯变 换法 [4]等 . 特别地 , 当 g (x 为特 殊形式 e x [P l (x cos wx +P n (x sin wx ]或 e x P m (x (其中 , P i (x 为 i 次多项式 (i =l, n, m , 为常数 时 , 待定系数法更 常用 . 考虑到线性微分方程的解具有可加性 , 上述两 种结构可归结为统一形式 :f (x =e x P m (x ( 为常 数 . 本文针对微分方程 y (n +a 1y (n -1 + +a n -1y ! +a n y =A (x e x (3 (其中 , a i (i =1, 2, , n 、 为常数 , A (x 为多 项 式 , 试图利用特征函数法 , 得出该类方程的一个特 解 . 2 y (n +a 1y (n -1 + +a n -1y ! +a n y = A(x e x 的一个特解 方程 (2 的特征方程为 n +a 1 n -1 + +a n -1 +a n =0 (4 定义 若记 f ( = n +a 1 n -1 + +a n -1 +a n (5 则称 f ( 为方程 (2 的特征函数 . 对特征函数 f ( 求导 , 有 :f (k ( =P k n # n -k +P k n -1#a 1 n -k -1 + +P k k +1#a n -k -1 +P k k #a n - k (6 (k =0, 1, , n. 特别地 , f (0 ( =f ( . 2. 1 A(x 为常数的情形 A(x 为常数 A 时 , 方程 (3 变为 : y (n +a 1y (n -1 + +a n -1y ! +a n y =A e x (7 其中 , 系数 a i (i =1, 2, , n 、 指数 均为 常 数 . 定理 若 是方程 (4 的 k (k =0, 1, , n 重 特征根 , 则方程 (7 存在一个特解 : y ~ = f (k ( x k e x , (k =0, 1, , n (8 其中 , f (k ( =f (k ( | = , 证明 分 是或不是方程 (4 的特征根两种情 况来讨论 : 1 若 不是方程 (4 的特征根 (即 k =0 , 此时 满足 f ( ? 0. 设方程 (7 的特解为 y ~ =B e x (其中 , B 为待定 常数 . 将 y ~以及 y ~得各阶导数分别代入方程 (7 , 考 虑到 e x ? 0, 有 B = n + 1 n -1 ++ n -1 +