数学建模之差分方程

  • 格式:pdf
  • 大小:323.08 KB
  • 文档页数:8

差分方程

对连续型变量而言,我们常常回导致到微分方程的问题. 对离散型变量将导

致一类的问题.

一、差分的定义

定义 设)(xyy是一个函数, 自变量从x变化到x+1, 这时函数的增量记为

)()1(xyxyy

x, 我们趁这个量为)(xy在点x步长为1的一阶差分,

简称为)(xy的一阶差分. 为了方便我们也记)(),1(

1xyyxyy

xx

,即

xxxyyy

1.

称xxxxxxxxyyyyyyyy

121122)()()(为)(xy二阶差

分,简记为xy2.

同样记)(2

xy为xy3,并称为三阶差分.

一般记)(1

xn

xnyy,称为n阶差分.且有inxin

ii

nxnyCy



)1(

0.

性质: 当a,b,C是常数, yx和zx 是函数时,

(1) Δ(C)=0;

(2) Δ(Cyx)= CΔ(yx);

(3) Δ(ayx+ b zx)= aΔyx+ bΔ zx ;

(4) Δ(yx zx)= zx+1Δyx+yx Δ zx = yx+1Δzx+zx Δyx;

(5)

111

1













xxxxxx

xxxxxx

xx

zzzyyz

zzzyyz

zy

.

例 已知),0(xxy

x求Δ(yx).

解 Δ(yx)= xx)1(.

特别, 当n为正整数时, Δ(yx)= inn

ii

nxC



1, 阶数降了一阶. 推论 若m, ,n为正整数时, m,> n P(x)为n次多项式,则0)(xPm.

例 已知),10(aayx

x求Δ(yx).

解 Δ(yx)= )1(1aaaaxxx.

二、差分方程

定义 设是含有未知函数差分的等式,称为差分方程。

它的一般形式为0),,,,(

1

nxxxyyyxF或0),,,,(

xn

xxyyyxG,

其中F, G是表达式,x是自变量. 使等式成立自变量的取值范围称为该方程的定

义域. 的0),,,,(

1

nxxxyyyxF的方程,也称为n阶差分方程. n为方程的阶.

形如

)()()()(

110xfyxayxayxa

xnnxnx

 (14-7-1)

称为n阶线性差分方程.. 0)(xf时为齐次的. 0)(xf为非齐次的.

差分方程是含有未知函数及其导数的方程, 满足该方程的函数称为差分方

程的解.对于一阶差分方程来说,它的含有一个任意常数的解,称为此微分方程

的通解.一般来说,对于n阶差分方程,其含有n个互相独立的任意常数的解称

为差分方程的通解.不含有任意常数的解称为差分方程的特解.同微分方程一样

椰油初值问题. 初值条件也有如下情形: 一阶的如: 00yy

xxx

.二阶的如:

00yy

xxx

,00yy

xxx

等等.

对于线性差分方程的解的结构有如下结论.

定理 如果)(

1xyy和)(

2xyy都是方程(14-7-1)的解,则对任意常数C1, C2,

)()(

2211xyCxyC也是方程(14-7-1)的解.

定理 设0)(

0xa ,)()2()1(,.....,,n

xxxyyy是

0)()()(

110

xnnxnxyxayxayxa 的n个线性无关的特解,则)()2(

2)1(

1.....n

xnxxxyCyCyCy是它的通解.

定理 设0)(

0xa ,)()2()1(,.....,,n

xxxyyy是齐次方程

0)()()(

110

xnnxnxyxayxayxa

的n个线性无关的特解,*

xy是非齐次方程

)()()()(

110xfyxayxayxa

xnnxnx



的一个特解,则*)()2(

2)1(

1.....

xn

xnxxxyyCyCyCy是非齐次方程的通解.

定理 设,)1(

xy是方程 )()()()(

1110xfyxayxayxa

xnnxnx

的解,

)2(

xy是方程 )()()()(

2110xfyxayxayxa

xnnxnx

的解,则

)2()1(

xxxyyy是方程

)()()()()(

21110xfxfyxayxayxa

xnnxnx

的解.

本书着重研究一阶和二阶常系数的差分方程.

三、一阶常系数的差分方程

一阶常系数的差分方程是)(

1xfpyy

xx

 (常数p≠0).

(a)当0)(xf,设x

xry是其齐次方程的解, 即 01xxprr,

所以 r=p . 那么01xxprr有通解x

xCpy(C为任意常数)

例 求差分方程023

1

xxyy的通解.

解 事实上原方程是0

32

1

xxyy所以其通解为x

xCy





32

(C为任意

常数)..

(b)当0)(xf,用待定系数法求其特解. (i) 如果)()(xPxf

n(n次多项式),则非齐次方程为 )(

1xPpyy

nxx

.

若 p=1, 即 )(

1xPyy

nxx

, 那么xy可以是n+1次多项式.,相减时常数项

和最高次数相被消去, 所以可以设][2

210n

nxxbxbxbbxy, 代入方

程后,比较系数确定nbbbb,,,,

210便得到一个特解.

若 p≠1, 最高次数相不可能被消去, 所以可以设有特解

n

nxxbxbxbby2

210, 同样代入方程后,比较系数确定

nbbbb,,,,

210便得到一个特解..

(ii) 如果)()(xPxf

nx()(xP

n是n次多项式,λ是常数),则非齐次方程为

)(

1xPpyy

nx

xx

.

为了求之一个特解,分两步: 第一步, 令 xx

xzy,代入方程得

)(

11xPzpz

nx

xx

xx

,

它等价于)(

1xPpzz

nxx

.

第二步, 用(i)的方法.

总之,对这种情况,可以直接设其特解为)(2

210n

nsx

xxbxbxbbxy,

其中当p≠λ时, s=0 , 当p=λ时, s=1 .

例 求差分方程 x

xxyy273

1

 的通解.

解 显然其齐次方程的通解为x

xCy3(C为任意常数).

设其特解为x

xby2, 所以有xxxbb272321, 从而得b=-7.

因此,原方程的通解为xx

xCy273.

四、二阶常系数的差分方程 这里讨论的是这样的方程: )(

12xfqypyy

xxx

 (p ,q是常数). 先给

结论 .

定理 x

xry是方程 0

12

xxxqypyy (16-7-2)

的解的充分必要条件r为方

程 02qprr (16-7-3)

的根 (读者自己证明).

(16-7-3))称为原方程的特征方程. 下面分步讨论.

(a)当0)(xf,

如果 042qp, 即其特征方程有两个不同实根,记为21,rr. 注意到

xxrr21,是线性无关的, 所以(16-7-2)有通解xx

xrCrCy2211, (21,CC是任意常

数).

如果042qp, 即其特征方程有两个相同实根,记为

221p

rr.,可以

验证xxp

xp









2,

2是(16-7-2)的线性无关的特解. 所以

x

xp

xCCy





2)(

21(21,CC是任意常数)是(16-7-2)的通解.

如果 042qp,因 p, q是实数, 即其特征方程有两互为共轭的复根, 记为24

22pqip

,记为0),sin(cosii . 可以验证

)sin(),cos(xxxx是(16-7-2)的线性无关的特解. 所以

))sin()cos((

21xCxCyx

x(21,CC是任意常数)是(16-7-2)的通解 .

例 求034

12

xxxyyy的通解.

解 其特征方程0342rr, 有根 -1, -3 . 原方程有通解