一类常系数线性差分方程的特解探究

第21卷第z期

2009年6月甘肃科学学报JoumalofGansusciencesV01.21No.2Jun.2009

一类常系数线性差分方程的特解探究

李自珍,龚东山

(兰州大学数学与统计学院,甘肃兰州730000)

摘要:针对一类常系数线性差分方程,运用特征函数法和比较系数法,得到了方程特解的显式表达.当方程非齐次项∥P。(五)中多项式P。(惫).一A(A为非零常数)时,可采用特征函数法得到方

程的一个公式化特解;当P。(五)一d。五“+d—l愚一1+…+d。(d。≠O)时,可采用比较系数法来得到

方程的一个特解.该方法简单易行,特解形式直观,避免了以前方法计算量过大的不足.

关键词:差分方程;非齐次;特解;特征函数法;比较系数法

0175.7中图分类号:文献标志码:A文章编号:1004—0366(2009)02一0001一03

ResearchontheSpecialSoIutionofaKindof

LinearDifferenceEquationwithCOnstantCoefficient

LIZi—zhen,GONGDong-shan

(Sc^oozD,M&咖聊m£ic54行dS缸£i5tics,L口竹柚o“U斑抛仃缸y,L口n加D“730000,ClIli,阮)

Abstr北t:

pressionhtllisarticle,weu辩methodsofspecialSolutionofachmcteriSticf吼ctiona11d唧aring∞eff.cientconstaIlttodeducetheex-o“hekindoflineardiff盯enceequation而thcoeffide此WhenpolynoIIlial

funmionis

solution0ff乙(忌)=A(Aisn0忏zeroconstant)iIlthenon.homogeneousitem∥P。(忌),method0fch牡解te五sticused;W1lenP0@)=厶五”+akl忌一1+…+磊(a乙≠O),、代candeducet11efomulaoft11especial

theequationf撒nnlethodofcomparing∞ef矗ci蚰t.ThetworIlethodsa、roidingtIledefectofso瑚屺01dII】Iethods

arewKchcalculateexcessively

Keywords:simple柚ddoable,andthespecialsolutionllasanintuidomsticexpressioll.differenceequation;non—homogeneous;specialsolution;methodofcharacteristicfunction;methodofcomparingcoefficient

咒阶常系数线性差分方程

Py(尼)+口。矿1y(忌)+…+口,卜l毋(忌)+口。y(忌)一,(七),

E[E卜1y(愚)]一E[y(惫+i一1)]一y(忌+i)(z=o,1,…,竹).

若,(忌)不恒为o,则称式(1)是竹阶常系数非齐次线性差分方程,并称方程(1)其中五∈(口,口+1,…,6—1,6)cZ,行∈N,口。,n2,…,%为常数,且口。≠0;E为位移算子,满足Fy(忌)=

E玎y(忌)+口lE”1y(愚)+…+口,广l毋(忌)+口。3,(志)=O

为式(1)对应的齐次方程.(2)

依据线性差分方程的通解结构理论知,式(1)的通解可以表示为式(1)的一个特解与齐次式(2)的通解之和[1’2].关于式(2)通解的计算,一般可采用待定参数法,先将差分方程化为代数方程,再求解这个代数方程即可[3].因此,要得到式(1)的通解,只需找到它的一个特解即可.

计算式(1)特解的常用方法有常数变异法‰5。、z_变换法(也称离散Laplace变换)[6。、算子法‘71及升阶法[8]等.其中采用常数变异法与z-变换法时,往往需要借助复杂的和分运算与逆变换z_1,而采用算子法和收稿日期:2008—09—09基金项目s国家自然科学基金项目(30470298).

2甘肃科学学报2009年第2期升阶法时,其计算量较大,且升阶法仅适用于带差分算子△的方程,从而限制了这些方法的使用与推广.

我们试图利用特征函数法‘引,针对非齐次项厂(忌)=∥P。(忌)(p为非零常数,P。(志)为班次多项式)的差分方程

py(忌)+口l矿1y(惫)+…+口,l毋(愚)+口。了(最)=∥P。(志)

得到它的一个公式化特解.

1(3)定义与引理

引理l若A‘是方程(2)的解,则有

.:l”+口1.:I”1+…+口,rlA十口。=O.(4)

定义l称式(4)为方程(2)的特征方程.并称

F(.;【)=A“+口1.;I”1+…+口,1A+口。(5)

为方程(2)的特征函数.

引理2特征函数(5)的导数满足

£;卑一c飘一‘+co。口1.:【一卜・+…+co。口,f_1.:I+a口,;,z—o,1,…,,1.

定义2(6)利用特征函数(5)以及特征函数的导数(6)求解方程的方法称为特征函数法.

2.1主要结论非齐次项^(七)=A(A为非零常数)的情形

此时,式(3)变为

Py(忌)+口1矿1y(忌)+…+4,r1勋(愚)+口。y(惫)一却‘.

定理l(特征函数法)常系数线性差分方程(7)存在特解(7)夕=杰如H,2而忽≯一’y(8)L芍,其中Z(o≤Z≤行)表示特征根卢的重数.

证明①若卢不是特征式(4)的根,则有Z一0且F(∞(芦)=F(卢)≠0.

为了得到式(7)的一个特解夕,可令夕=卑‘(B为待定常数).将夕(志)、毋(五)、…、F夕(志)代人式(7),整理后有卑‘[,+口l∥。1+…+口,r-产+口。]=舡‘.由于∥为非零常数,故产‘≠o,化简上式后有B一≯了iii]兰}干而一F各・

此时式(7)存在特解夕=卑‘一F各∥一声害茜户‘.

②若口是特征方程(4)的z重特征根(1≤z≤竹),则有F(卢)+,q)=…一FH’(户)一。且pq)≠o.由线性齐次差分方程的通解结构理论知,∥,如‘,…,忌H∥是式(2)的z个线性无关解,且是么‘不是该齐次方程的解.令式(7)有特解夕一擞么‘(B为待定常数).将夕(忌)、毋(志)、…、P夕(曼)代人式(7),消去因子∥并整理后有B一尹F会丽・此时式(7)存在特解:

推论1

存在特解方程夕一耳耸万正么H,1≤z≤兀(9)E,l了(足)+口lE”1y(矗)+…+口,l毋(惫)+口。y(惫)=A,

夕。声焉正‘,忙o,1’..‘・”

其中Z(O≤£≤竹)表示特征根卢一1的重数.

证明式(9)是式(7)卢一l的特殊情况.将p=1代入式(8)可得式(10)成立.

2.2d(10)非齐次项n(七)=以足。+d一.七一‘+…+巩(以。d一。。…。磊为常数)的情形此时,式(3)变为

第21卷李自珍等:一类常系敷线性差分方程的特解探究

Fy(五)+口。矿1y(矗)+…+口。y(五)一[d。志4+d—l五一1+…+d。]矿.

定理2(比较系数法)式(11)存在特解形如(11)

夕=[6m五“+址1七一+…+60]五2弘H,z—o,1,…,卯

其中z(o≤Z≤竹)表示特征根卢的重数,k,6一。,…,60为常数.

证明①若卢不是特征式(4)的根,则有Z=O且F(严)≠O.(12)

为了得到式(11)的一个特解夕,可令夕=[6。愚”+6一。忌一1+…+6。]∥(6。,靠。,…,60为待定系数).将夕(五)、印(忌)、…、_E,I多(忌)代入式(11),消除非零因子矿后有

6私”[(五+挖)“+口。(五+n一1)“+…+口,rl(五+1)“+口。惫“]+

址1卢一[(志+,1)一1+口l(惫+以一1)一1+…十口,1(五+1)一1+口。忌一1]+…+

61止(愚+n)+口1(志十刀一1)+…+口,l(忌+1)+口。忌]+60[1+口1+…+口,1+口。]=

d。+dl奄+…十d。惫4.

利用特征函数导数式(6),整理上式可得

6。c::IF(产)惫”+[6,lc:::rlF(Ifl)+6mc扣F7(p)]七一1+

[%靠zFq)+‰艮。∥7Q)+kc:(∥P(产)+∥7q))]五卅十…一磊+西忌+…+姒“.

通过比较系数可得到一个关于6。,6一’..・,60线性方程组.解此方程组可依次得出6。,6一,,…,6。.从而求得式(11)的一个特解.

②若口是特征式(4)的Z重特征根(1≤Z≤行),则有

F(岸)=F7(口)一…=F‘H’(岸)=O,且F‘o(岸)≠0.

令式(11)的一个特解夕一[60“+6一.五一1+…+60]忌么H(k,6一l,..・,60为待定系数),将夕(惫)、E箩(愚)、…、P夕(忌)代人式(11).类似①中计算思路,先消除非零因子∥,再通过比较系数可得到一个关于6m,6,。,..・,60线性方程组,从而可求出式(11)的一个特解形如式(12).

例1求方程E3y(忌)一3E2y(忌)+3毋(愚)一y(忌)一100的一个特解

例2

解卢一1为特征方程的三重根,且F(3’(1)一6.由推论1可得方程的特解为夕=未尝%忌3=警忌3.求方程Egy(五)一6E2y(惫)+12毋(五)一8y(五)=志・2‘的一个特解弘=2是特征方程的三重根,且,(2)一6.由定理2,可令方程的特解为夕一(6。矗+60)忌32H,其中6。,60为待定常数.将夕(惫)、坊(意)、E2夕(足)、E3夕(忌)代人原方程有246。七十366l+660一惫,比较系数有{;:6li三。;。,解这个二元线性方程组,可得6。一去,60一一{.因此所述差分方程存在特解:i366:+66。;o,解这个二元线性方程组,可得6-一麦,60一一专‘因此所述差分方程存在特解:

夕=(去忌一{)矗32H一(南五一矗)忌32‘.

参考文献:

[1]

[2]sNElaydi.AnIntfoductionDifferencetoEquations[M].NewYorklSpringer.Verlag,1996.Press,1991.KeUeywG.Petef∞nAC.DifferenceEquationswithAppli∞tions[M].NewYork:A∞demic

[3]阮炯.差分方程和常微分方程[M].上海:复旦大学出版社.1991.

[4]王联、王慕秋.常差分方程[M].乌鲁木齐t新疆大学出版社.1991.

[5]saberElaydi.AnIntroductionto.DifferenceEquations[M].NewYork:spriIlger-Verlag,2005.

[6]张广.高英.差分方程的振动理论[M].北京:高等教育出版杜,2001.

[7]黄兆良.算子法在常系数线性差分方程中的应用[J].武汉船舶职业技术学院学报,2008。7(1):26・28.

[8]王炳兴.王海敏.常系数非齐次线性差分方程的特解[J].工科数学.1999,15(4):16l-165.

[9]龚东山,刘岳巍,牛富俊.特征函数在高阶常微分方程特解计算中的应用[J].吉林师范大学学报(自然科学版),2008,29(4),8—10.作者简介:

李自珍(1944一)男.河南省许昌人,1969年毕业于兰州大学数力系,理学博士.现任兰州大学数学与统计学院教授,博士研究生导师,研究方向为数学生态学.

一类常系数线性差分方程的特解探究

作者:

作者单位:

刊名:

英文刊名:

年,卷(期):

被引用次数:李自珍, 龚东山, LI Zi-zhen, GONG Dong-shan兰州大学,数学与统计学院,甘肃,兰州,730000甘肃科学学报JOURNAL OF GANSU SCIENCES2009,21(2)0次

参考文献(9条)

1. S N Elaydi An Introduction to Difference Equations 1996

2. Kelley W G. Peterson A C Difference Equations with Applications 1991

3. 阮炯 差分方程和常微分方程 1991

4. 王联. 王慕秋 常差分方程 1991

5. Saber Elaydi An Introduction to Difference Equations 2005

6. 张广. 高英 差分方程的振动理论 2001

7. 黄兆良 算子法在常系数线性差分方程中的应用[期刊论文]-武汉船舶职业技术学院学报 2008(01)

8. 王炳兴. 王海敏 常系数非齐次线性差分方程的特解 1999(04)

9. 龚东山. 刘岳巍. 牛富俊 特征函数在高阶常微分方程特解计算中的应用[期刊论文]-吉林师范大学学报(自然科学版) 2008(04)

相似文献(10条)

1.期刊论文 莫国良 用升阶法求常系数线性非齐次差分方程的特解 -高等数学研究2005,8(4)

利用升阶法可求解常系数线性非齐次差分方程的特解.

2.学位论文 龚东山 几类常差分方程精确解的研究 2009

差分方程研究的主要内容包括两个问题:差分方程的精确解和方程解的定性分析。其中对于解的定性分析研究,一般以差分模型为基础,具体讨论差分方程的稳定性、有界性、振动性、渐近性、周期性和概周期性等问题,相关的研究在过去几十年中取得了许多重要的成果。然而对于差分方程精确解的研究则相对滞后。尽管许多学者在寻找差分方程的精确解方面取得了一定的进展,但仍有许多结论还不太成熟:在线性若分方程中,有些论证还不是很完善,一般函数的和分不易得到;常系数齐次线性方程的多重特征根所对应特解的线性无关性并非显然;变系数线性方程的主要结论还停留在解的结构与形式上,缺乏一个普适性的方法;在非线性差分方程中,即使是一阶的离散黎卡提方程,也很难得到精确解。

本文重点讨论了定义在整数集Z或它的子集D上的几类常差分方程精确求解方法与解的显式表达。所得的主要结论如下:

I)对于线性齐次差分方程

Eny(k)+a1(k)En-1y(k)+…+an-1(k)Ey(k)+an(k)y(k)=0(1)

可以得到方程(1)的解结构,指出其通解为y(k)=c1y1(k)+c2y2(k)+…+cnyn(k)的形式,其中c1,c2,…,cn是独立的任意常数

,y1(k),y2(k),…,yn(k)是方程(1)的一个基本解组。通过运用参数待定法(阮炯,2002),论证了常系数线性齐次差分方程的多重特征根所对应的全部特解存在线性无关性,并推导出该差分方程的通解表达式。

II)对于线性非齐次差分方程

Eny(k)+a1(k)En-1y(k)+…+an-1(k)Ey(k)+an(k)y(k)=f(k)(2)

证明了方程(2)的通解等于它的一个特解与相应的齐次方程(1)的通解之和。如果已知齐次方程的一个基本解组,可运用常数变异法与函数和分计算,得到非齐次方程的一个特解的形式表达。

1.对于常系数线性非齐次差分方程,利用特征函数法(李自珍,龚东山,2009)得到非齐次项.厂(七)为多项式、指数函数、三角函数、多项式与指数函数的乘积、对数函数以及它们的线性组合时的公式化特解。该方法简便易行,克服了常数变异法(Paul Mason Batchelder,1927)、比较系数法(SaberElaydi,2005)及拉普拉斯变换法(张广,张高英,2001)等传统方法计算工作量过人的缺陷,且特解形式非常直观。

2.对于几类可精确求解的变系数线性差分方程,利用构造函数法,将某些变系数方程化为常系数差分方程;通过引入变上限定和分,给出了一阶变系数线性差分方程的通解:运用函数积分法,得到了系数为线性函数时方程的通解;利用观察法找到方程的一个特解,并以此特解为基础,得出二阶变系数线性齐次方程的通解;借助降阶法,当差分方程的系数函数满足一定条件时,将复合差分方程问题化成若干个一阶变系数线性差分方程的求解问题。

III)对于线性差分方程组得到了该方程组有解的一个充要条件,加强了王联(1991)关于一阶线性差分方程组与高阶线性差分方程同解的充分性结论,完善了线性差分方程精确解的理论体系,并利用线性代数与矩阵理论,推导出常系数线性差分方程组精确解的显式表达。

IV)对于两类非线性差分方程

通过引入独立通解(龚东山,2008),得到了(4)与(5)相应齐次方程的通解表示。研究表明:这两类齐次差分方程的通解由若干个独立通解共同构成,且独立通解的个数与差分方程的阶数无关;当非齐次项为某些特殊函数时,给出了方程的若干个特解,且特解的个数与方程的次数有关,并论证了这些特解的线性相关性。

3.期刊论文 陈景年. CHEN Jing-nian 二阶非齐次线性差分方程为极限圆型的判定 -山东大学学报(理学版)2005,40(3)

提出了二阶非齐次线性差分方程为极限圆型的定义,并给出了几个充分性的判定方法.

4.期刊论文 胡劲松 求非齐次差分方程特解的推广 -重庆师范学院学报(自然科学版)2003,20(3)

教材中用"待定系数"法介绍了一、二阶常系数线性非齐次差分方程在f(x)=dxPm(x)时特解的求法.本文将该方法推广,讨论了当

f(x)=dx[Ps(x)coswx+Pn(x)siwx]时常系数线性非齐次差分方程特解的求法.

5.学位论文 陈庆华 无标度网络的建模分析与度分布计算方法 2006

本文研究无标度网络。我们选择无标度网络的建模分析、度分布计算方法以及相关性等作为主要的研究方向。本论文系统深入地研究了这些问题。 首先,研究了无标度网络的模型构造与模型分析,侧重于揭示现实网络的演化机制,构建适合现实网络的演化模型。

其次,研究了度分布的计算问题。根据马尔可夫链理论,我们提出了一种新的度分布数值计算方法.

第三,研究了无标度网络的相关性问题,重点讨论了BA网络的度相关,给出了BA网络的联合度分布.

本文取得了以下几个创新成果:

(1)构建了四个无标度网络的演化模型,即模型3.5.1至模型3.5.4(见第3.5节).在模型3.5.1中,我们考虑了网络的局部相互作用,即内部边和重新连接等。在模型3.5.2中,我们提出了一种反择优删除连线的演化机制。根据平均场方法,计算了这两个网络的度分布P(k),它们都是幂律分布。

在模型3.5.3和模型3.5.4中,我们首次提出了网络的对数增长,这是一种新的演化机制。特别指出,我们无法求出这两个网络度分布P(k)的解析表达式。我们利用自己提出的马氏链方法给出了度分布P(k)的数值计算(见第四章),数值结果表明这两个度分布都具有幂律尾部.

(2)用随机过程的观点研究复杂网络,发现了无标度网络与马氏链之间的内在联系。对于BA模型,任意给定一个结点i,设k(t)表示它在t时刻的度数,我们证明了随时间变化的度数序列{K(t),t=i,i+1,…)是一个非齐次马氏链,给出了具有时间相依的一步转移概率矩阵

P(t+1),i=1,2,…。

(3)根据马氏链理论,由转移概率矩阵P(t+1)可以给出网络在t时刻的度分布P(七,t)的矩阵运算表达式。因为矩阵P(t+1)具有特殊的简单结构,利用矩阵运算性质,我们提出了一种度分布数值计算的新方法,简称为马氏链方法.

应用马氏链方法研究了BA模型和三个加速增长网络模型,其中BA网络度分布的数值计算与原有的解析解和数值模拟进行比较,三种结果十分接近(见图4.1.1)。特别地,对于两个具有对数增长的(有向)网络模型,用原有的解析方法无法得到度分布的表达式,我们进行了度分布的数值计算,数值结果表明这两个度分布都具有幂律尾部(见图4.3.1和图4.4.1),这两个系统都演化成无标度网络。

(4)研究了无标度网络的相关性问题,求出了BA网络的联合度分布。应用率方程方法和二维母函数性质,我们给出了BA网络相邻点对的联合度分布P(k,l)的公式(5.2.1)。应用平均场方法和顺序统计量性质,我们也给出了BA网络任意点对的联合度分布P(k,k)的公式(5.3.1)。这两个联合度分布都证明了BA网络具有结点的度相关特征.

6.期刊论文 邓勇. DENG Yong 常系数非齐次线性递归数列求特解的简易方法 -达县师范高等专科学校学报

2006,16(5)

利用数列的差分将常系数非齐次线性递归数列转化为常系数非齐次线性差分方程,从而得到一种求常系数非齐次线性递归数列特解的简易方法.

7.学位论文 陈景年 二阶差分方程为极限圆型的判定 1999

该文主要研究方程二阶齐次线性差分方程和二阶非齐次线性差分方程为极限圆型的条件,给出了方程二阶齐次线性差分方程为极限圆型或极限点型的充要或充分条件,改进了已有的结果,并给出了非齐次方程二阶非齐次线性差分方程为极限圆型的定义和两个充分性的判定结果。文章内容分为三部分:1。引言,2。齐次方程的情形,3。非齐次方程的情形。

8.期刊论文 杨继明. 蔡炯辉 常系数非齐次线性微分方程组初值问题的求解公式 -宝鸡文理学院学报(自然科学版)2002,22(1)

分别给出了常系数非齐次线性微分方程组和常系数非齐次线性差分方程组在给定的初始条件下的求解公式.

9.学位论文 苑运良 几类偏差分方程的定性研究 2005

本论文分别讨论了几类时滞偏差分方程解的振动性和两类时滞偏差分方程的正解不存在性,同时对非齐次线性偏差分方程的稳定性进行了研究。 对于常系数时滞偏差分方程,讨论了解的振动性质,得到了两个方程所有解都振动的充分条件,所得结果推广了已有的一些结论。对于具有连续变量的两类时滞偏差分方程,研究了解的振动性,利用函数的单调性,建立了方程的振动准则。

另外研究了一类多时滞偏差分方程的正解不存在性,给出了方程不存在正解的若干充分条件,同时研究了一类非线性时滞偏差分方程的正解不存在性,得到了两个新的充分性判据。

最后讨论了一类非齐次变系数线性偏差分方程的稳定性,获得了关于方程稳定的几个充分条件,所得结论推广了相关文献的已有结果。

10.期刊论文 陈旭东 非齐次差分方程解一类行列式 -科技信息(科学·教研)2007,""(26)

本文利用几种常见形式的非齐次差分方程来求解一类具有递归关系的行列式,并给出这类特殊行列式的求解方法和技巧.

本文链接:http://d.wanfangdata.com.cn/Periodical_gskxxb200902001.aspx

授权使用:中共汕尾市委党校(zgsw),授权号:cd481015-a2c7-4e4d-ab93-9dcc01177755

下载时间:2010年8月8日

第21卷第z期

2009年6月甘肃科学学报JoumalofGansusciencesV01.21No.2Jun.2009

一类常系数线性差分方程的特解探究

李自珍,龚东山

(兰州大学数学与统计学院,甘肃兰州730000)

摘要:针对一类常系数线性差分方程,运用特征函数法和比较系数法,得到了方程特解的显式表达.当方程非齐次项∥P。(五)中多项式P。(惫).一A(A为非零常数)时,可采用特征函数法得到方

程的一个公式化特解;当P。(五)一d。五“+d—l愚一1+…+d。(d。≠O)时,可采用比较系数法来得到

方程的一个特解.该方法简单易行,特解形式直观,避免了以前方法计算量过大的不足.

关键词:差分方程;非齐次;特解;特征函数法;比较系数法

0175.7中图分类号:文献标志码:A文章编号:1004—0366(2009)02一0001一03

ResearchontheSpecialSoIutionofaKindof

LinearDifferenceEquationwithCOnstantCoefficient

LIZi—zhen,GONGDong-shan

(Sc^oozD,M&咖聊m£ic54行dS缸£i5tics,L口竹柚o“U斑抛仃缸y,L口n加D“730000,ClIli,阮)

Abstr北t:

pressionhtllisarticle,weu辩methodsofspecialSolutionofachmcteriSticf吼ctiona11d唧aring∞eff.cientconstaIlttodeducetheex-o“hekindoflineardiff盯enceequation而thcoeffide此WhenpolynoIIlial

funmionis

solution0ff乙(忌)=A(Aisn0忏zeroconstant)iIlthenon.homogeneousitem∥P。(忌),method0fch牡解te五sticused;W1lenP0@)=厶五”+akl忌一1+…+磊(a乙≠O),、代candeducet11efomulaoft11especial

theequationf撒nnlethodofcomparing∞ef矗ci蚰t.ThetworIlethodsa、roidingtIledefectofso瑚屺01dII】Iethods

arewKchcalculateexcessively

Keywords:simple柚ddoable,andthespecialsolutionllasanintuidomsticexpressioll.differenceequation;non—homogeneous;specialsolution;methodofcharacteristicfunction;methodofcomparingcoefficient

咒阶常系数线性差分方程

Py(尼)+口。矿1y(忌)+…+口,卜l毋(忌)+口。y(忌)一,(七),

E[E卜1y(愚)]一E[y(惫+i一1)]一y(忌+i)(z=o,1,…,竹).

若,(忌)不恒为o,则称式(1)是竹阶常系数非齐次线性差分方程,并称方程(1)其中五∈(口,口+1,…,6—1,6)cZ,行∈N,口。,n2,…,%为常数,且口。≠0;E为位移算子,满足Fy(忌)=

E玎y(忌)+口lE”1y(愚)+…+口,广l毋(忌)+口。3,(志)=O

为式(1)对应的齐次方程.(2)

依据线性差分方程的通解结构理论知,式(1)的通解可以表示为式(1)的一个特解与齐次式(2)的通解之和[1’2].关于式(2)通解的计算,一般可采用待定参数法,先将差分方程化为代数方程,再求解这个代数方程即可[3].因此,要得到式(1)的通解,只需找到它的一个特解即可.

计算式(1)特解的常用方法有常数变异法‰5。、z_变换法(也称离散Laplace变换)[6。、算子法‘71及升阶法[8]等.其中采用常数变异法与z-变换法时,往往需要借助复杂的和分运算与逆变换z_1,而采用算子法和收稿日期:2008—09—09基金项目s国家自然科学基金项目(30470298).

2甘肃科学学报2009年第2期升阶法时,其计算量较大,且升阶法仅适用于带差分算子△的方程,从而限制了这些方法的使用与推广.

我们试图利用特征函数法‘引,针对非齐次项厂(忌)=∥P。(忌)(p为非零常数,P。(志)为班次多项式)的差分方程

py(忌)+口l矿1y(惫)+…+口,l毋(愚)+口。了(最)=∥P。(志)

得到它的一个公式化特解.

1(3)定义与引理

引理l若A‘是方程(2)的解,则有

.:l”+口1.:I”1+…+口,rlA十口。=O.(4)

定义l称式(4)为方程(2)的特征方程.并称

F(.;【)=A“+口1.;I”1+…+口,1A+口。(5)

为方程(2)的特征函数.

引理2特征函数(5)的导数满足

£;卑一c飘一‘+co。口1.:【一卜・+…+co。口,f_1.:I+a口,;,z—o,1,…,,1.

定义2(6)利用特征函数(5)以及特征函数的导数(6)求解方程的方法称为特征函数法.

2.1主要结论非齐次项^(七)=A(A为非零常数)的情形

此时,式(3)变为

Py(忌)+口1矿1y(忌)+…+4,r1勋(愚)+口。y(惫)一却‘.

定理l(特征函数法)常系数线性差分方程(7)存在特解(7)夕=杰如H,2而忽≯一’y(8)L芍,其中Z(o≤Z≤行)表示特征根卢的重数.

证明①若卢不是特征式(4)的根,则有Z一0且F(∞(芦)=F(卢)≠0.

为了得到式(7)的一个特解夕,可令夕=卑‘(B为待定常数).将夕(志)、毋(五)、…、F夕(志)代人式(7),整理后有卑‘[,+口l∥。1+…+口,r-产+口。]=舡‘.由于∥为非零常数,故产‘≠o,化简上式后有B一≯了iii]兰}干而一F各・

此时式(7)存在特解夕=卑‘一F各∥一声害茜户‘.

②若口是特征方程(4)的z重特征根(1≤z≤竹),则有F(卢)+,q)=…一FH’(户)一。且pq)≠o.由线性齐次差分方程的通解结构理论知,∥,如‘,…,忌H∥是式(2)的z个线性无关解,且是么‘不是该齐次方程的解.令式(7)有特解夕一擞么‘(B为待定常数).将夕(忌)、毋(志)、…、P夕(曼)代人式(7),消去因子∥并整理后有B一尹F会丽・此时式(7)存在特解:

推论1

存在特解方程夕一耳耸万正么H,1≤z≤兀(9)E,l了(足)+口lE”1y(矗)+…+口,l毋(惫)+口。y(惫)=A,

夕。声焉正‘,忙o,1’..‘・”

其中Z(O≤£≤竹)表示特征根卢一1的重数.

证明式(9)是式(7)卢一l的特殊情况.将p=1代入式(8)可得式(10)成立.

2.2d(10)非齐次项n(七)=以足。+d一.七一‘+…+巩(以。d一。。…。磊为常数)的情形此时,式(3)变为

第21卷李自珍等:一类常系敷线性差分方程的特解探究

Fy(五)+口。矿1y(矗)+…+口。y(五)一[d。志4+d—l五一1+…+d。]矿.

定理2(比较系数法)式(11)存在特解形如(11)

夕=[6m五“+址1七一+…+60]五2弘H,z—o,1,…,卯

其中z(o≤Z≤竹)表示特征根卢的重数,k,6一。,…,60为常数.

证明①若卢不是特征式(4)的根,则有Z=O且F(严)≠O.(12)

为了得到式(11)的一个特解夕,可令夕=[6。愚”+6一。忌一1+…+6。]∥(6。,靠。,…,60为待定系数).将夕(五)、印(忌)、…、_E,I多(忌)代入式(11),消除非零因子矿后有

6私”[(五+挖)“+口。(五+n一1)“+…+口,rl(五+1)“+口。惫“]+

址1卢一[(志+,1)一1+口l(惫+以一1)一1+…十口,1(五+1)一1+口。忌一1]+…+

61止(愚+n)+口1(志十刀一1)+…+口,l(忌+1)+口。忌]+60[1+口1+…+口,1+口。]=

d。+dl奄+…十d。惫4.

利用特征函数导数式(6),整理上式可得

6。c::IF(产)惫”+[6,lc:::rlF(Ifl)+6mc扣F7(p)]七一1+

[%靠zFq)+‰艮。∥7Q)+kc:(∥P(产)+∥7q))]五卅十…一磊+西忌+…+姒“.

通过比较系数可得到一个关于6。,6一’..・,60线性方程组.解此方程组可依次得出6。,6一,,…,6。.从而求得式(11)的一个特解.

②若口是特征式(4)的Z重特征根(1≤Z≤行),则有

F(岸)=F7(口)一…=F‘H’(岸)=O,且F‘o(岸)≠0.

令式(11)的一个特解夕一[60“+6一.五一1+…+60]忌么H(k,6一l,..・,60为待定系数),将夕(惫)、E箩(愚)、…、P夕(忌)代人式(11).类似①中计算思路,先消除非零因子∥,再通过比较系数可得到一个关于6m,6,。,..・,60线性方程组,从而可求出式(11)的一个特解形如式(12).

例1求方程E3y(忌)一3E2y(忌)+3毋(愚)一y(忌)一100的一个特解

例2

解卢一1为特征方程的三重根,且F(3’(1)一6.由推论1可得方程的特解为夕=未尝%忌3=警忌3.求方程Egy(五)一6E2y(惫)+12毋(五)一8y(五)=志・2‘的一个特解弘=2是特征方程的三重根,且,(2)一6.由定理2,可令方程的特解为夕一(6。矗+60)忌32H,其中6。,60为待定常数.将夕(惫)、坊(意)、E2夕(足)、E3夕(忌)代人原方程有246。七十366l+660一惫,比较系数有{;:6li三。;。,解这个二元线性方程组,可得6。一去,60一一{.因此所述差分方程存在特解:i366:+66。;o,解这个二元线性方程组,可得6-一麦,60一一专‘因此所述差分方程存在特解:

夕=(去忌一{)矗32H一(南五一矗)忌32‘.

参考文献:

[1]

[2]sNElaydi.AnIntfoductionDifferencetoEquations[M].NewYorklSpringer.Verlag,1996.Press,1991.KeUeywG.Petef∞nAC.DifferenceEquationswithAppli∞tions[M].NewYork:A∞demic

[3]阮炯.差分方程和常微分方程[M].上海:复旦大学出版社.1991.

[4]王联、王慕秋.常差分方程[M].乌鲁木齐t新疆大学出版社.1991.

[5]saberElaydi.AnIntroductionto.DifferenceEquations[M].NewYork:spriIlger-Verlag,2005.

[6]张广.高英.差分方程的振动理论[M].北京:高等教育出版杜,2001.

[7]黄兆良.算子法在常系数线性差分方程中的应用[J].武汉船舶职业技术学院学报,2008。7(1):26・28.

[8]王炳兴.王海敏.常系数非齐次线性差分方程的特解[J].工科数学.1999,15(4):16l-165.

[9]龚东山,刘岳巍,牛富俊.特征函数在高阶常微分方程特解计算中的应用[J].吉林师范大学学报(自然科学版),2008,29(4),8—10.作者简介:

李自珍(1944一)男.河南省许昌人,1969年毕业于兰州大学数力系,理学博士.现任兰州大学数学与统计学院教授,博士研究生导师,研究方向为数学生态学.

一类常系数线性差分方程的特解探究

作者:

作者单位:

刊名:

英文刊名:

年,卷(期):

被引用次数:李自珍, 龚东山, LI Zi-zhen, GONG Dong-shan兰州大学,数学与统计学院,甘肃,兰州,730000甘肃科学学报JOURNAL OF GANSU SCIENCES2009,21(2)0次

参考文献(9条)

1. S N Elaydi An Introduction to Difference Equations 1996

2. Kelley W G. Peterson A C Difference Equations with Applications 1991

3. 阮炯 差分方程和常微分方程 1991

4. 王联. 王慕秋 常差分方程 1991

5. Saber Elaydi An Introduction to Difference Equations 2005

6. 张广. 高英 差分方程的振动理论 2001

7. 黄兆良 算子法在常系数线性差分方程中的应用[期刊论文]-武汉船舶职业技术学院学报 2008(01)

8. 王炳兴. 王海敏 常系数非齐次线性差分方程的特解 1999(04)

9. 龚东山. 刘岳巍. 牛富俊 特征函数在高阶常微分方程特解计算中的应用[期刊论文]-吉林师范大学学报(自然科学版) 2008(04)

相似文献(10条)

1.期刊论文 莫国良 用升阶法求常系数线性非齐次差分方程的特解 -高等数学研究2005,8(4)

利用升阶法可求解常系数线性非齐次差分方程的特解.

2.学位论文 龚东山 几类常差分方程精确解的研究 2009

差分方程研究的主要内容包括两个问题:差分方程的精确解和方程解的定性分析。其中对于解的定性分析研究,一般以差分模型为基础,具体讨论差分方程的稳定性、有界性、振动性、渐近性、周期性和概周期性等问题,相关的研究在过去几十年中取得了许多重要的成果。然而对于差分方程精确解的研究则相对滞后。尽管许多学者在寻找差分方程的精确解方面取得了一定的进展,但仍有许多结论还不太成熟:在线性若分方程中,有些论证还不是很完善,一般函数的和分不易得到;常系数齐次线性方程的多重特征根所对应特解的线性无关性并非显然;变系数线性方程的主要结论还停留在解的结构与形式上,缺乏一个普适性的方法;在非线性差分方程中,即使是一阶的离散黎卡提方程,也很难得到精确解。

本文重点讨论了定义在整数集Z或它的子集D上的几类常差分方程精确求解方法与解的显式表达。所得的主要结论如下:

I)对于线性齐次差分方程

Eny(k)+a1(k)En-1y(k)+…+an-1(k)Ey(k)+an(k)y(k)=0(1)

可以得到方程(1)的解结构,指出其通解为y(k)=c1y1(k)+c2y2(k)+…+cnyn(k)的形式,其中c1,c2,…,cn是独立的任意常数

,y1(k),y2(k),…,yn(k)是方程(1)的一个基本解组。通过运用参数待定法(阮炯,2002),论证了常系数线性齐次差分方程的多重特征根所对应的全部特解存在线性无关性,并推导出该差分方程的通解表达式。

II)对于线性非齐次差分方程

Eny(k)+a1(k)En-1y(k)+…+an-1(k)Ey(k)+an(k)y(k)=f(k)(2)

证明了方程(2)的通解等于它的一个特解与相应的齐次方程(1)的通解之和。如果已知齐次方程的一个基本解组,可运用常数变异法与函数和分计算,得到非齐次方程的一个特解的形式表达。

1.对于常系数线性非齐次差分方程,利用特征函数法(李自珍,龚东山,2009)得到非齐次项.厂(七)为多项式、指数函数、三角函数、多项式与指数函数的乘积、对数函数以及它们的线性组合时的公式化特解。该方法简便易行,克服了常数变异法(Paul Mason Batchelder,1927)、比较系数法(SaberElaydi,2005)及拉普拉斯变换法(张广,张高英,2001)等传统方法计算工作量过人的缺陷,且特解形式非常直观。

2.对于几类可精确求解的变系数线性差分方程,利用构造函数法,将某些变系数方程化为常系数差分方程;通过引入变上限定和分,给出了一阶变系数线性差分方程的通解:运用函数积分法,得到了系数为线性函数时方程的通解;利用观察法找到方程的一个特解,并以此特解为基础,得出二阶变系数线性齐次方程的通解;借助降阶法,当差分方程的系数函数满足一定条件时,将复合差分方程问题化成若干个一阶变系数线性差分方程的求解问题。

III)对于线性差分方程组得到了该方程组有解的一个充要条件,加强了王联(1991)关于一阶线性差分方程组与高阶线性差分方程同解的充分性结论,完善了线性差分方程精确解的理论体系,并利用线性代数与矩阵理论,推导出常系数线性差分方程组精确解的显式表达。

IV)对于两类非线性差分方程

通过引入独立通解(龚东山,2008),得到了(4)与(5)相应齐次方程的通解表示。研究表明:这两类齐次差分方程的通解由若干个独立通解共同构成,且独立通解的个数与差分方程的阶数无关;当非齐次项为某些特殊函数时,给出了方程的若干个特解,且特解的个数与方程的次数有关,并论证了这些特解的线性相关性。

3.期刊论文 陈景年. CHEN Jing-nian 二阶非齐次线性差分方程为极限圆型的判定 -山东大学学报(理学版)2005,40(3)

提出了二阶非齐次线性差分方程为极限圆型的定义,并给出了几个充分性的判定方法.

4.期刊论文 胡劲松 求非齐次差分方程特解的推广 -重庆师范学院学报(自然科学版)2003,20(3)

教材中用"待定系数"法介绍了一、二阶常系数线性非齐次差分方程在f(x)=dxPm(x)时特解的求法.本文将该方法推广,讨论了当

f(x)=dx[Ps(x)coswx+Pn(x)siwx]时常系数线性非齐次差分方程特解的求法.

5.学位论文 陈庆华 无标度网络的建模分析与度分布计算方法 2006

本文研究无标度网络。我们选择无标度网络的建模分析、度分布计算方法以及相关性等作为主要的研究方向。本论文系统深入地研究了这些问题。 首先,研究了无标度网络的模型构造与模型分析,侧重于揭示现实网络的演化机制,构建适合现实网络的演化模型。

其次,研究了度分布的计算问题。根据马尔可夫链理论,我们提出了一种新的度分布数值计算方法.

第三,研究了无标度网络的相关性问题,重点讨论了BA网络的度相关,给出了BA网络的联合度分布.

本文取得了以下几个创新成果:

(1)构建了四个无标度网络的演化模型,即模型3.5.1至模型3.5.4(见第3.5节).在模型3.5.1中,我们考虑了网络的局部相互作用,即内部边和重新连接等。在模型3.5.2中,我们提出了一种反择优删除连线的演化机制。根据平均场方法,计算了这两个网络的度分布P(k),它们都是幂律分布。

在模型3.5.3和模型3.5.4中,我们首次提出了网络的对数增长,这是一种新的演化机制。特别指出,我们无法求出这两个网络度分布P(k)的解析表达式。我们利用自己提出的马氏链方法给出了度分布P(k)的数值计算(见第四章),数值结果表明这两个度分布都具有幂律尾部.

(2)用随机过程的观点研究复杂网络,发现了无标度网络与马氏链之间的内在联系。对于BA模型,任意给定一个结点i,设k(t)表示它在t时刻的度数,我们证明了随时间变化的度数序列{K(t),t=i,i+1,…)是一个非齐次马氏链,给出了具有时间相依的一步转移概率矩阵

P(t+1),i=1,2,…。

(3)根据马氏链理论,由转移概率矩阵P(t+1)可以给出网络在t时刻的度分布P(七,t)的矩阵运算表达式。因为矩阵P(t+1)具有特殊的简单结构,利用矩阵运算性质,我们提出了一种度分布数值计算的新方法,简称为马氏链方法.

应用马氏链方法研究了BA模型和三个加速增长网络模型,其中BA网络度分布的数值计算与原有的解析解和数值模拟进行比较,三种结果十分接近(见图4.1.1)。特别地,对于两个具有对数增长的(有向)网络模型,用原有的解析方法无法得到度分布的表达式,我们进行了度分布的数值计算,数值结果表明这两个度分布都具有幂律尾部(见图4.3.1和图4.4.1),这两个系统都演化成无标度网络。

(4)研究了无标度网络的相关性问题,求出了BA网络的联合度分布。应用率方程方法和二维母函数性质,我们给出了BA网络相邻点对的联合度分布P(k,l)的公式(5.2.1)。应用平均场方法和顺序统计量性质,我们也给出了BA网络任意点对的联合度分布P(k,k)的公式(5.3.1)。这两个联合度分布都证明了BA网络具有结点的度相关特征.

6.期刊论文 邓勇. DENG Yong 常系数非齐次线性递归数列求特解的简易方法 -达县师范高等专科学校学报

2006,16(5)

利用数列的差分将常系数非齐次线性递归数列转化为常系数非齐次线性差分方程,从而得到一种求常系数非齐次线性递归数列特解的简易方法.

7.学位论文 陈景年 二阶差分方程为极限圆型的判定 1999

该文主要研究方程二阶齐次线性差分方程和二阶非齐次线性差分方程为极限圆型的条件,给出了方程二阶齐次线性差分方程为极限圆型或极限点型的充要或充分条件,改进了已有的结果,并给出了非齐次方程二阶非齐次线性差分方程为极限圆型的定义和两个充分性的判定结果。文章内容分为三部分:1。引言,2。齐次方程的情形,3。非齐次方程的情形。

8.期刊论文 杨继明. 蔡炯辉 常系数非齐次线性微分方程组初值问题的求解公式 -宝鸡文理学院学报(自然科学版)2002,22(1)

分别给出了常系数非齐次线性微分方程组和常系数非齐次线性差分方程组在给定的初始条件下的求解公式.

9.学位论文 苑运良 几类偏差分方程的定性研究 2005

本论文分别讨论了几类时滞偏差分方程解的振动性和两类时滞偏差分方程的正解不存在性,同时对非齐次线性偏差分方程的稳定性进行了研究。 对于常系数时滞偏差分方程,讨论了解的振动性质,得到了两个方程所有解都振动的充分条件,所得结果推广了已有的一些结论。对于具有连续变量的两类时滞偏差分方程,研究了解的振动性,利用函数的单调性,建立了方程的振动准则。

另外研究了一类多时滞偏差分方程的正解不存在性,给出了方程不存在正解的若干充分条件,同时研究了一类非线性时滞偏差分方程的正解不存在性,得到了两个新的充分性判据。

最后讨论了一类非齐次变系数线性偏差分方程的稳定性,获得了关于方程稳定的几个充分条件,所得结论推广了相关文献的已有结果。

10.期刊论文 陈旭东 非齐次差分方程解一类行列式 -科技信息(科学·教研)2007,""(26)

本文利用几种常见形式的非齐次差分方程来求解一类具有递归关系的行列式,并给出这类特殊行列式的求解方法和技巧.

本文链接:http://d.wanfangdata.com.cn/Periodical_gskxxb200902001.aspx

授权使用:中共汕尾市委党校(zgsw),授权号:cd481015-a2c7-4e4d-ab93-9dcc01177755

下载时间:2010年8月8日


相关文章

  • 二阶变系数齐次微分方程通解的求法
  • ((高等数学研究1:?,($$" 一类二阶变系数齐次微分方程通解的求法 李永利! 摘要!!桑改莲!(河南质量工程职业学院基础部!河南平顶山!!"#$$$)作为文[%]两种情形的统一推广,给出一类二阶变系数线性微分方程的通 ...查看


  • 常系数线性差分方程的微分解法
  • [***********][***********][***********]77生的整体素质也起到了重要的作用:符合国家教育部提出的; 重在参与:重在普及:扩大受益面 >以竞赛促教改 撰写论文的能力? 查阅资料的能力:分析问题? 解 ...查看


  • 常微分方程差分解法.入门.多解法
  • 毕业论文 题 目 抛物型方程的差分解法 学 院 专 业 信息与计算科学 班 级 学 生 王丹丹 学 号 指导教师 王宣欣 二〇一二 年 五 月二十五日 摘 要 偏微分方程的数值解法在数值分析中占有重要的地位,很多科学技术问题的数值计算包括了 ...查看


  • 河南专升本数学复习的七个知识点你清楚么?
  • 河南专升本数学复习的七个知识点你清楚么? 在河南专升本数学复习过程中,最需要注意的这七个知识点:极限.等价无穷小代换,函数连续性,可导性和可微性,参数估计,级数,微分方程,随机变量的数字特征,一维随机变量函数的分布详细内容下面马上会为大家讲 ...查看


  • 用matlab实现线性常系数差分方程的求解
  • 数字信号处理课程设计 题目: 试实现线性常系数差分方程的求解 学院: 专业: 班级: 学号: 组员: 指导教师: 题目:用Matlab 实现线性常系数差分方程求解 一. 设计要求 1. 2. 3. 掌握线性常系数差分方程的求解 熟练掌握Ma ...查看


  • 10-345二阶线性微分方程
  • 10.3 可降阶的高阶微分方程 一. y ( n ) = f ( x ) 型的微分方程 二. y ′′ = f ( x , y ′) 型的微分方程 三. y ′′ = f ( y , y ′) 型的微分方程 第一种类型 y( n) = f ...查看


  • 计量经济学考点整理
  • 第一章 计量经济学定义:统计学.经济理论和数学三者的结合.正经济学中,我们用数学的函数概念表达对经济变量间的关系的看法. 计量经济学模型建立的步骤:一.理论模型的设计 二.样本数据的收集 三.模型参数的估计 四.模型的检验 计量经济学模型成 ...查看


  • 有限差分法.边界元法和离散元法
  • 有限差分法 已经发展的一些近似数值分析方法中,最初常用的是有限差分法,它可以处理一些相当困难的问题.但对于几何形状复杂的边界条件,其解的精度受到限制,甚至发生困难.作为60年代最重要的科技成就之一的有单元法.在理论和工程应用上都_得到迅速发 ...查看


  • 差分方程及其应用
  • 差分方程及其应用 在经济与管理及其它实际问题中,许多数据都是以等间隔时间周期统计的.例如,银行中的定期存款是按所设定的时间等间隔计息,外贸出口额按月统计,国民收入按年统计,产品的产量按月统计等等.这些量是变量,通常称这类变量为离散型变量.描 ...查看


热门内容