第32卷第8期 2011年8月 湖南科技学院学报 Journal of Hunan University of Science and Engineering Vl01.32 NO.8 Aug.201 1 ._-.1 Cj 同 阶常系数线性差分方程的解的研究 劳智 (广州大学,广东广州510000;广东海洋大学寸金学院,广东湛江524033) 摘要:在文献【5]中,论文作者将常系数齐次线性差分方程改写为矩阵与向量乘积形式的递推关系,并运用相似矩阵 的理论给出了常系数齐次线性差分方程通解的解析形式。在论文中,则通过引进算子把常系数齐次线性差分方程化为一些式 子之积,再利用算子相关的引理,简便地得到k阶常系数齐次线性差分方程k个线性无关的解,从而得到通解。 关键词:差分方程;常系数差分方程;算子 中图分类号:0155 文献标识码:A 文章编号:1673-2219(2011)08-0019-05 引进算千: =Y “,E Y =E(Ey )=E(y +1)=Y +2,…,E Y =Y¨ 其中:用I表示恒等算子,即ty =Y 。 在解常差分方程时常用到以下算子的有关引理: 引理1 (E一 ,)(E一 1)y =(E一 ,)(E一 ,)y 证明: ( 一 )( 一 1)y = E2一九IEI一九jEl+ jI2、yn =【E 一( 十 )E十 ,]Y =(E一2jI)(E一 ,)y 即算子满足交换律。 引理2(E-M)(Pa(n)2n)= 一1(n) 其中 (,1)为n的 次多项式。 证明:不妨设 Pa(n)=ao,z +口l,毫 +…+口 ln+a , 则(E一 )( (,1) ) 收稿日期:2011一o4—15 作者简介:劳智(1973一),女,广东化州人,广东海洋大学寸金学院讲师,广州大学在职硕士研究生,主要从事常微 分方程理论及其应用研究。 19 =( -M)(ao +以l,l +…+口 _1n+Ⅱ ) :【ao(,l+1)口+Ⅱ1(,l十1) 一十…+d _】(,l+1)+Ⅱ ], + ・(a0 +口1,z 十…+以 _】,l+口 ) “ =(口of(,l+1) 一n 】+口l[(n十1) 一nO,-1】+…+n _l[(,2+1)一 】), = ( ) 由引理2可知:算子(E一 )每作用一次,多项式 (n)降低一次。故: 引理3 (E一 ) (Pa(n)2f1)=0 定理1 如果Y(¨,Y( )是七阶常系数齐次线性差分方程 Y肿 +P lY + 一1+…+PoY =0 (1) 的两个解,那么 =GY +C2y‘ ,( ,C2是任意实数) 也是方程(1)的解。 定理2如果.),(n,),( ,…,y( )是方程(1)的七个线性无关的特解,那么 Yn=Cly( )+C2Y( )+…+C Y( ,(C1,C2,…,C 是任意实数) 就是方程(1)的通解。 定理3方程(1)的特征方程为 +p 一I +...+pI +Po=0 (2)。 根据特征根的情况,可以写出对应解的各项见表3-1: 特征方程(2)的根 方程(1)通解中对应的项 (1)单实根 给出一项C (2) :重实根 给出 项(cl+c2n+...+cjn 一 ) (3)一对单共轭复粮 .2= ±f 给出两项(C1 COSnO+C2 sinnO)r (4)-- ̄crj重共轭复根 1.2 +ilf 给出2j项 li((cc川+c +_1+c2n十・・‘+c.. +c2,I )』月 )c0s,l sinnO}+ l, . 则(1)的通解为Y =cly‘ +C2Y +…+CkY‘ ,其中C1,C2,…,C 是任意实数,且通解Y,l中的每一项都由特征 方程的一个根所对应,其对应情况如表3.I。 证明:由定理3可知,要求方程(1)的通解,就要找出它k个线性无关的解。 20 先找 Y = (3) 形式的解。 若(3)为方程(1)的解,则 +Pk_l ̄,n 一 +…+po2"=0 ( +Pk_l ̄k +…+PO)=0 即 +Pk_l:l,k一 +…+p0=0 (2) 将上式称为方程(1)的特征方程,记为(2)。 设方程(2)可以分解为: ( 一 )嘶( 一 ) …( 一 ) =O (4) 其中口l+ 2+…+ ,,l=k。 引进算子 , Ey =Y肘l,E2y =E(Ey )=E(), +1)=), +2….,Eky Y + Ekyn Yn+ , 且ty =Y 。 则方程(1)可以写成 (E +Pk-1E 一 +・・・+poE)Yn=0 (5) 由方程(4)可知,方程(5)可以写成 (E一 ,)啦(E一221) …(E一2mI) Y =O (6) 由引理1可知,箅子满足交换律。故不妨考虑 ( -,;111)嘶Y =0 (7) 若al=1,即 为单根。此时, (E一 ,) : 州一 ¨=0 即 Y : 是方程(7)的解,也是方程(6)的解。 若al=2,此时由引理1,得 (E一 J) ( ):(E-21I)(E一 ,)(,z ) : 一 ,)【 +1) 州一A・ 】 :(E一 ,) 21 : 一 =O 即 Y = 是方程(7)的解,也是方程(6)的解。 由引理3 (E一 )叶 ( (,z) )=0 可知 ( 一 ,)嘶(n ~ )=0 即 Y =,z 一 是方程(7)的解,也是方程(6)的解。 故,一般的 ,有 个解,分别是 4n, , 。 ,…,na'l-l ̄,nl 因此,有k个线性无关的解 2 一 2 , , ,…, , , , ,…,rta2-1 ̄,…, ,n ,n2 ,…,norm-1 。 故方程),n+ +Pk一1Yn 一I+…+PoYn=0(1)的通解为: cO4.+c1.4"+C?n +…+C ̄r1-1 o ̄1-1 +c +c + c ,z +…+c 2-1,zO ̄2-1 +…+c 0/Ln +c 1,z n+c三,z +…+c m一 ntZm-1 其中c ,c ,cl2,…,c ~,c ,c ,c ,…,c ~,c ,c ,c三,…,c ~是任意实数。 若 为方程(6) 重实根,则 (cl+C2n+..・+c naj-1) 是方程(1)通解的对应项。 若 为方程(6) 重复根,设 = (其中rj, 分别是 的模及辐角) 由于复根应该共轭成对出现,故 -i : , r,e 也是方程(6)的解。 因此, 对应的解为: r;einO)。nr;einO) 2r;e n8i。…,n。i一 r e 8i 对应的解为: r e— 8inr e一 ei 2 r e一 ojai-1 nje-inOj .n,…,n共2 个复解。 由欧拉公式 i8e:cosO+isinO.P一 =cos0一isin0 及定理1可知, (n , 肿 +n,r7e-in01)=n 0n c。s n J 去(n e们L n 4e-,.e):n sin n 』 帆是方程(6)的解。 因此,得到2 』个实数解 c。s 0i,nr ̄cos aj,n 2r]c。s 0 ,…,,z cos 0 , sin 01,nr,sin Oj,n2r;sin 0J,…,,z sin 0J。 故 【(cl十c2,l+.._+cJ )c。s,l +(c川+cj+2n--t--- ̄-I-c2J )sin,l 是方程(1)通解的对应项。 得证。 例2.1 求差分方程 Y +7—9y +6+29y +5—45y +4+55y +3—63y +2+27yn+l一27y =0的通解。 解:该方程的特征方程为 凳一92 29凳一45 55凳一63 272,一27=0 即 ( 一3)。( +1) =0 其根为: ,2,3=3, 45=i,267=-i ..根据定理2_3,得该方程的通解为 (C1+c2n+c3n2)3 +(c4+ )c。s詈+(c6+c7)Sin 其中C1,C2,C3,C4,C5,C6,C7为任意实数。 参考文献: 【l健 士银.一阶常系数差分方程的特解公式【J】.沈阳工程学院学报,2008,(1):91—93. [2】赵士银.一类二阶常系数差分方程特解的简单求法[J]_河南教育学院学报,2007,(9):14.15. 【3]郑桂梅.高等数学【M].长沙:国防科技大学出版社,2008. 【4]祝浩锋.k阶常系数线性差分方程的几种解法【J】.工科数学,1994,(1):16—21. 【5】胡劲松,郑克龙.用矩阵的方法求解常系数齐次线性差分方程[J】.大学数学,2007,(3):133—134 (责任编校:京华) 23