数学建模之差分方程
- 格式: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(1aaaaxxx.
二、差分方程
定义 设是含有未知函数差分的等式,称为差分方程。
它的一般形式为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)(
0xa ,)()2()1(,.....,,n
xxxyyy是
0)()()(
110
xnnxnxyxayxayxa 的n个线性无关的特解,则)()2(
2)1(
1.....n
xnxxxyCyCyCy是它的通解.
定理 设0)(
0xa ,)()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是其齐次方程的解, 即 01xxprr,
所以 r=p . 那么01xxprr有通解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
nxxbxbxbby2
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为方
程 02qprr (16-7-3)
的根 (读者自己证明).
(16-7-3))称为原方程的特征方程. 下面分步讨论.
(a)当0)(xf,
如果 042qp, 即其特征方程有两个不同实根,记为21,rr. 注意到
xxrr21,是线性无关的, 所以(16-7-2)有通解xx
xrCrCy2211, (21,CC是任意常
数).
如果042qp, 即其特征方程有两个相同实根,记为
221p
rr.,可以
验证xxp
xp
2,
2是(16-7-2)的线性无关的特解. 所以
x
xp
xCCy
2)(
21(21,CC是任意常数)是(16-7-2)的通解.
如果 042qp,因 p, q是实数, 即其特征方程有两互为共轭的复根, 记为24
22pqip
,记为0),sin(cosii . 可以验证
)sin(),cos(xxxx是(16-7-2)的线性无关的特解. 所以
))sin()cos((
21xCxCyx
x(21,CC是任意常数)是(16-7-2)的通解 .
例 求034
12
xxxyyy的通解.
解 其特征方程0342rr, 有根 -1, -3 . 原方程有通解