WJF8-9差分方程简介
差分方程简介

以上三个定理揭示了n阶齐次及非齐次线性差 分方程的通解结构, 它们是求解线性差分方程非常 重要的基础知识.
我们只探讨一阶常系数线性差分方程的解法.
二、 一阶常系数线性差分方程 一阶常系数线性差分方程的一般形式为
yt 1 ayt f ( t ),
其中 a 0 为常数,f ( t ) 为已知函数. 当 f (t ) 0 时,称方程
Pt Pt 1 ( St 1 Dt 1 )
( 为常数),
即
Pt [1 (b d )]Pt 1 (a c ).
定义2 含有未知函数差分或未知函数几个时期值的方 程就称为差分方程. 例如
F ( x , yt , yt 1 , , yt n ) 0,
它对应的齐次方程
yt n a1 yt n1 an1 yt 1 an yt 0
的通解与它自己本身的一个特解之和,即通解等于
Y C1 y1 ( t ) C 2 y2 ( t ) C n yn ( t ) y* ( t ),
* 其中 y ( t ) 是它自己本身的一个特解.
2
(C ) 0;
(Cyt ) C ( yt );
3
4
(ayt bzt ) a ( yt ) b( zt );
( yt zt ) zt 1yt yt zt yt 1zt zt yt ;
yt z t yt yt z t z t 1 yt yt 1 z t . 5 zt zt 1 zt zt 1 zt
2 yt 2 yt 3t .
定义3 如果一个函数代入差分方程后,方程两边 恒等,则称此函数为差分方程的解.
差分方程介绍

例如,如认为第一季度的销售量大体按线性增长,可设销售量
(1) yk = ak + b
(1) (1) yk = 1.3k + 9.5, y6 = 17.3
得到
缺点:数据少,用回归分析不好。改用差分方程
yk = a1 yk −1 + a2 yk = a1 yk −1 + a2 yk − 2 + a3 或者 用二阶差分, yk = a1 yk −1 + a2 yk − 2 + a3
和最小二乘法,使 最小,求出
∑[ y
3
5
k
− (a1 yk −1 + a2 yk − 2 + a3 a3 = −8, y6 = 21, y7 = 19
上述为预测各年第一季度销售量而建立的二阶差分方程,虽然其系数与前 上述为预测各年第一季度销售量而建立的二阶差分方程,虽然其系数与前5 年第一季度的统计数据完全吻合,但用于预测时预测值与事实不符。 年第一季度的统计数据完全吻合,但用于预测时预测值与事实不符。凭直 第六年估计值明显偏高,第七年销售量预测值甚至小于第六年。 觉,第六年估计值明显偏高,第七年销售量预测值甚至小于第六年。稍作 分析,不难看出,如分别对每一季度建立一差分方程, 分析,不难看出,如分别对每一季度建立一差分方程,则根据统计数据拟 合出的系数可能会相差甚大,但对同一种商品,这种差异应当是微小的, 合出的系数可能会相差甚大,但对同一种商品,这种差异应当是微小的, 故应根据统计数据建立一个共用于各个季度的差分方程。 故应根据统计数据建立一个共用于各个季度的差分方程。 为此, 为此,将季度编号为
∗
只取一次项近似为: 只取一次项近似为: (5)是(4)的近似线性方程,x ∗ 也是 ( 5 ) ) )的近似线性方程, 的平衡点, 的平衡点,关于线性方程平衡点稳定的条 件上面已给出。 件上面已给出。
差分方程

第七节 差分方程对连续型变量而言,我们常常回导致到微分方程的问题. 对离散型变量将导致一类的问题.一、差分的定义定义 设)(x y y =是一个函数, 自变量从x 变化到x +1, 这时函数的增量记为)()1(x y x y y x -+=∆, 我们称这个量为)(x y 在点x 步长为1的一阶差分,简称为)(x y 的一阶差分. 为了方便我们也记)(),1(1x y y x y y x x =+=+,即x x x y y y -=∆+1.称x x x x x x x x y y y y y y y y +-=---=∆∆+++++121122)()()(为)(x y 二阶差分,简记为x y 2∆.同样记)(2x y ∆∆为x y 3∆,并称为三阶差分.一般记)(1x n x n y y -∆∆=∆,称为n 阶差分.且有i n x i ni i n x ny C y -+=-=∆∑)1(0. 性质: 当a,b,C 是常数, y x 和z x 是函数时,(1) Δ(C )=0;(2) Δ(Cy x )= C Δ(y x );(3) Δ(ay x + b z x )= a Δy x + b Δ z x ;(4) Δ(y x z x )= z x+1Δy x +y x Δ z x = y x+1Δz x +z x Δy x ;(5) 1111++++∆-∆=∆-∆=⎪⎪⎭⎫ ⎝⎛∆x x x x x x x x x x x x x x z z z y y z z z z y y z z y .例 已知),0(≠=x x y x α求Δ(y x ).解 Δ(y x )= ααx x -+)1(.特别, 当n 为正整数时, Δ(y x )= i n n i i n x C-=∑1, 阶数降了一阶.推论 若m, ,n 为正整数时, m,> n P(x)为n 次多项式,则0)(=∆x P m .例 已知),10(≠<=a a y x x 求Δ(y x ).解 Δ(y x )= )1(1-=-+a a a a x x x .二、差分方程定义 设是含有未知函数差分的等式,称为差分方程。
差分方程介绍

yt a(1 b) yt1 abyt2 G (4.23)
(4.23)式是一个二阶常系数差分方程,其特征方程为
2 a(1 b) ab 0 ,相应特征根为
1 a2(1 b)2 ab 1 4
(4.24)
成立时才是稳定的。 (4.24)式可用于预报经济发展趋势。
现用待定系数法求方程 (4.23)的一个特解
代入(4.23)式,得
C G 1a
y。t 令 yt C
故当(4.24)式成立时,差分方程 (4.23)的通解为
yt
t (C1 cost
C2 sint )
G 1a
其中ρ为 1,2 的模,ω为其幅角。
例如,若取
a
1 ,4
b
1 2
反之若ab商品紧缺易引起顾客抢购该商品供售市场易造成混乱如果生产者对市场经济的蛛网模型有所了解为了减少因价格波动而造成的经济损失他应当提高自己的经营水平不应当仅根据上一周期的价格来决定现阶段的生产量
§4.4 差分方程建模
一、差分方程简介 以t 表示时间,规 定t只取非负整数。t=0表示第一周期初, t=1表示第二周期初等。 记yt 为变量y在时刻t 时的取值,则
(步三) 求的非通齐解次,则方非程齐(4次.1方5)的程一(4个.15特)的解通y解t.若为yt为yt 方 程yt(4.16)
求非齐次方程(4.15)的特解一 般要用到 常数变易法,计算较繁。 对特殊形式 的b(t)也可使用 待定 系数法。
例4.13 求解两阶差分方程 yt2 yt t
在市场经济中,对每一商品事实上存在着两个不同的 函数: (1)供应函数x=f(P),它是价格P的单增函数,其曲 线称为供应曲线。 (2)需求函数x=g(P),它是价格P的单降函数,其 曲线称为需求曲线,供应曲线与需求曲线的 形状如图所示。
3.4.差分方程简介

故原方程的通解为
(2)方程对应的特征方程为 λ 1= 0 ,其特征根为 λ =1 ,于是齐次方
* 程的通解为 yi = C 。设方程的特解为: yi = Acos π i + Bsin π i 2 2
将其代入原方程可解得 A = B = 1 2 故原方程的通解为
yi = C 1 (cos π i + sin π i) 2 2 2
λn + Pλn1 ++ P 1λ + P = 0 1 n n
(3.4.3)
方程(3.4.3)称为(3.4.2)的特征方程。若 λ1, λ2 ,, λn 是(3.4.3)的 n 个不同的根,则 Y (i) = λ1 ,Y2 (i) = λ2 ,,Yn (i) = λn 就是(3.4.2)的 n 个 1
r 1+ P ++ P 1 n
是稳定的条件与对应的齐次方程(3.4.2)完全相同。
此外,对于 n 维向量 yi 和 n× m 常数矩阵 A 构成的一阶线性差分方程组
yi +1 + Ayi = 0
其 平 衡点 稳 定的 条 件是 矩阵 A 的特 征 根
λi ( i =1,2,, n) 均有 λi <1 。 即均在复平面上的
(3)
若 Y1 (i) ,…,Yn (i) 是方程(3.4.2) n 个线性
无关的解,则它们的线性组合 C1Y (i) ++ CnYn (i) 1 就是 (3.4.2)的通解。 Y1 (i) ,…, Yn (i) 称为(3.4.2) 的一组基本解。 (4) 若 C1Y (i) ++ CnYn (i) 是 (3.4.2) 的通解, y *(i) 是 1 非齐次方程(3.4.1)的一个特解,则
微分方程和差分方程简介

返 回
(二)建立数值解法的一些途径
设 xi 1 xi h, i 0,1,2, n 1, 可用以下离散化方法求解微分方程: y' f(x,y) y(x0 ) y0
1、用差商代替导数 若步长h较小,则有
y ' ( x) y ( x h) y ( x ) h
解 首先分离变量 ,得
g ( y )dy
f ( x ) dx C
2 例1 求微分方程 y 3x y的通解。
1 2 dy 3 x dx y 两端积分,得 即 ln y x 3 C1 y e
x 3 C1
或y e e
C1
x3
因 e C1 仍是任意常数,令其为C,则所求得通解为 y Ce
二、常见的微分方程的类型及其解法:
1.一阶微分方程
y f ( x, y )
常用的解法:分离变量法
形如
dy f ( x) g ( y ) dx P ( x) P2 ( y ) dx Q1 ( x)Q2 ( x) 0 1
的方程均为可分离变量 的微分方程。
对(2)式两端分别积分,便可得到微分方程的通解 其中C为任意常数。
例1 求
解
d2y
2
dx du 1 u 2 的通解. dt
0 应表达为:D2y=0.
输入命令:dsolve('Du=1+u^2','t')
结
果:u = tg(t-c)
例 2 求微分方程的特解.
d 2 y dy 2 4 29 y 0 dx dx y (0) 0, y ' (0) 15
对马尔萨斯人口模型的解作进一步分析, 当 t 时,x(t ) ,表明人口将无限增长。马 尔萨斯人口论的核心内容是:人口按几何级数 增长,而生活资料则按算术级数增长,两者的 矛盾必会给人类社会进步造成障碍。马尔萨斯 并不认为: 解决人口过剩和生活资料匮乏两 者之间的矛盾,只有通过失业、饥饿、犯罪甚 至战争等方式来自发调节。使用消极手段来遏 制人口增长,这是人们对马尔萨斯人口论的一 种误解。
差分方程讲解

解 特征方程为
2 4 + 16 = 0.
方程的根为
1,2 2 2 3i , r 4, .
3
原方程的通解为
y x C1 cos x C 2 sin 3 3
x x4 .
代入初始条件 y0=0, y1=1得
C1 cos 0 C 2 sin 0 40 0, 1 C1 cos C 2 sin 4 1, 3 3
其中 B 为待定系数.
例11 求差分方程 yx+2 3yx+1 + 2yx = 2x的一个特解.
解 对应的齐次方程的特征方程为 方程的根为
2 3 + 2 = 0. 1 = 1, 2 = 2,
因为 q = 2 =2, 设特解为 y Bx 2 x ,
x
代入原方程, 得 B(x+2)2x+23B(x+1)2x+1+2Bx2x = 2x, 1 B , 2 1 x x 所求特解为 yx x 2 x 2 . 2
设特解的待定式为 m y x B0 B1 x Bm x (a 1)
或
(6)
(7)
y x ( B0 B1 x Bm x m ) x (a 1)
其中B0 , B1 , , Bm为待定系数.
例5 求差分方程 yx+1 2yx = 3x2 的一个特解.
为二阶差分, 记为2 yx, 即
2 yx = (yx) = yx+2 2 yx+1 + yx
同样可定义三阶差分3yx, 四阶差分4yx, 即
3yx = (2yx), 4yx = (3yx) .
差分方程简介

它的通解是 y x Cx A ( A 是任何实常数). ( 3) y x Pn1 ( x ) ( n 1次多项式) 通解 y x Pn ( x ) ( n次多项式 )
4 n y x 0
通解 y x是n 1次多项式.
二、一阶常系数线性差分方程
形如: y x 1 ay x f ( x ) 齐次方程: y x 1 ay x 0
y x ( x 2 ) ( x 1)2 x 2 2 x 1
2 y x 2 ( x 2 ) (2 x 1) 2( x 1) 1 ( 2 x 1) 2
3 y x ( 2 y x ) ( 2) 2 2 0
x n x( x 1)( x 2)( x n 1) , x 0 1
例2 设 求 x n
解: x n ( x 1)n x n
( x 1) x( x 1)( x 1 n 1) x( x 1)( x n 1) [( x 1) ( x n 1)]x( x 1)( x n 2) nx n1
2. 差分方程 有某种商品 t 时期的供给量St与需求 一个例子: 量Dt都是这一时期价格Pt 的线性函数:
St a bPt (a , b 0) , Dt c dPt (c, d 0)
设 t 时期的价格Pt由 t –1时期的价格 Pt 1与供给量 及需求量之差 St 1 Dt 1 按如下关系确定.
Pt Pt 1 ( St 1 Dt 1 )
( 为常数),
即
Pt [1 (b d )]Pt 1 (a c )
这样的方程就是差分方程.
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
例2 设 求 ∆x (n)
解: ∆x (n) = ( x + 1)(n) − x(n)
= ( x + 1) x( x − 1)L( x + 1 − n + 1) − x( x − 1)L( x − n + 1) = [( x + 1) − ( x − n + 1)]x( x − 1)L( x − n + 2) = nx(n−1)
简单差分方程的解
(1) ∆yx = 0
yx+1 − yx = 0 ⇒ yx+1 = yx = yx−1 = L = y0 = C
(2) ∆yx = C
是任何实常数) 它的通解是 yx = A ( A 是任何实常数).
yx+1 − yx = C ⇒ yx+1 = yx + C= yx−1 + 2C = L = y0 + xC
= y x+2 − 2 y x+1 + y x = ∆2 y x 同样定义 ∆3 y x = ∆(∆2 y x ) , ∆4 y x = ∆(∆3 y x ) , LL
∆
二阶及二阶以上的差分统称为高阶差分 二阶及二阶以上的差分统称为高阶差分. 高阶差分
例1 求 ∆( x2 ) ,
∆2 ( x 2 ) ,
解∶齐次方程的通解为 Pt = A[1 − λ (b + d )]t
a ≠ 1 设 ~t = k 则有 k − [1 − λ (b + d )]k = λ (a + c ) p a+c 即 k= b+d ~ = a+c 特解∶ 方程 特解∶ pt b+d
例3 求 ∆λx ;
∆lnx
解:∆λx = λx+1 − λx = (λ − 1)λx 补充差分的性质: 补充差分的性质: 1. ∆cux = c∆ux (c为常数)
1 = ln( x + 1) − lnx = ln 1 + x
2. ∆ ( u x + v x ) = ∆ u x + ∆ v x
(ii) f ( x) = cbx (其中 c, b ≠ 1 均为常数 其中 均为常数)
y b ≠ a 设 ~x = kbx 代入方程 kb x+1 − akbx = cb x C k= 即 b−a ~ = c bx 的特解为∶ 故方程 的特解为∶ y x b − a ~ = kxbx y b=a
x 为二阶差分方程, 如 y x+2 − 2 y x+1 − y x = 3 为二阶差分方程,
它等价于 ∆2 y x − 2 y x = 3 x 定义3 如果一个函数代入差分方程后, 定义 如果一个函数代入差分方程后,方程两边恒 则称此函数为该差分方程的解. 等,则称此函数为该差分方程的解 对差分方程附加的一定条件称为初始条件 初始条件. 对差分方程附加的一定条件称为初始条件 满足初 特解. 始条件的解称为特解 始条件的解称为特解 如果差分方程的解中含有相 互独立的任意常数的个数恰好等于方程的阶数, 互独立的任意常数的个数恰好等于方程的阶数,则 称它为差分方程的通解 通解. 称它为差分方程的通解
是任何实常数) 它的通解是 yx = Cx + A ( A 是任何实常数). ( 3) ∆y x = Pn−1 ( x ) ( n − 1次多项式 ) 通解 y x = Pn ( x ) ( n次多项式 )
(4) ∆ x =0 y
n
次多项式 通解 y x是 n − 1次多项式 .
二、一阶常系数线性差分方程
∆3 ( x 2 )
解:设 y x = x 2 ,那么
∆y x = ∆( x 2 ) = ( x + 1)2 − x 2 = 2 x + 1 ∆2 y x = ∆2 ( x 2 ) = ∆(2 x + 1) = [2( x + 1) + 1] − (2 x + 1) = 2
∆3 y x = ∆(∆2 y x ) = ∆(2) = 2 − 2 = 0
8.9
差分方程简介
差分、 差分、差分方程基本概念 一阶常系数线性差分方程 二阶常系数线性差分方程 小结
一、差分方程的基本概念
dy 连续变化的时间内, 来刻画; 连续变化的时间内,变量 y 的变化速度用 来刻画; dt 但有时,变量要按一定的离散时间取值. 但有时,变量要按一定的离散时间取值. ∆y 来刻画变化速度. 来刻画变化速度. 常取在规定的时间区间上的差商 ∆t
整理得 (− B0 + B1 + B2 ) + (− B1 + 2B2 ) x − B2 x = 3 x B0 = −9 , − B0 + B1 + B2 = 0 比较系数 − B1 + 2B2 = 0 解得 B1 = −6 , B = −3 − B2 = 3 2 y 即特解 ~x = −9 − 6 x − 3 x 2
Pt = Pt −1 − λ ( St −1 − Dt −1 )
为常数), ( λ 为常数),
即
Pt − [1 − λ (b + d )]Pt −1 = λ (a + c)
这样的方程就是差分方程. 这样的方程就是差分方程
定义2 含有自变量、 定义 含有自变量、未知函数以及未知函数差分的 方程称为差分方程. 方程称为差分方程 方程中含有未知函数差分的最高阶数称为差分 方程的阶. 方程的阶 n 阶差分方程的一般形式为: 阶差分方程的一般形式为: H( x , yx , ∆yx , ∆2 yx , L, ∆n yx ) = 0 将 ∆y x = y x+1 − y x
, 定义1 定义 设函数 y = f ( x) 记为 yx
取遍非负整数时,函数值可以排成一个数列: 当 x 取遍非负整数时,函数值可以排成一个数列:
y0 , y1,L , yx ,L L L
则差 yx+1 − yx 称为 yx 的差分,也称为一阶差分, 的差分,也称为一阶差分, 记为 ∆yx ,即 ∆yx = yx+1 − yx . 二阶差分 一阶差分的差分 ∆(∆yx ) = ∆( yx+1 − yx ) = ( yx+2 − y x+1 ) − ( yx+1 − yx )
x
kx + k − kx = C 即 k = C 故方程的特解为: y 故方程的特解为:~x = Cx
例4
yx+1 − 2 yx = 3 x 2的通解. 的通解. 求差分方程
x 解 yx+1 − 2 yx = 0 的通解为Yx = A2 yx+1 − 2 yx = 3 x 2 的特解为 Qa = 2 ≠ 1, 设 ~ = B + B x + B x 2 代入方程,则有 yx 代入方程, 0 1 2 B0 + B1 ( x + 1) + B2 ( x + 1)2 − 2B0 − 2B1 x − 2B2 x 2 = 3 x 2
1 5 所以方程的特解为 y x = 2 2
x
x
1 5 1 差分方程的通解为 y x = + A 2 2 2
x
x
例8 本节引例的差分方程为
Pt +1 − [1 − λ (b + d )]Pt = λ (a + c )
a, b, c, d > 0,λ ≠ 0
设
x
代入方程 即
k( x + 1)b x+1 − akxbx = cb x
y x = cxbx−1 方程 特解∶ 特解∶
c k ( x + 1)b − bkx = c ⇒ kb = c ⇒ k = b
1 5 的通解. 例7 求差分方程 y x+1 − 2 y x = 2 的通解 x 1 1 解∶ y x+1 − y x = 0 的通解为 Y = A x 2 2 x 1 5 5 Q ≠ 设 ~x = k 代入方程 y 2 2 2 x+1 x x 1 1 5 5 5 1 5 k − k = ⇒ k − k = 1 ⇒ k = 2 2 2 2 2 2 2
yx+1 − yx = 3 x 2 + x + 4 的特解. 的特解. 例6 求差分方程 解 设特解为 ~ = x( B + B x + B x 2 ) 代入原方程得 y
x 0 1 2
3B2 x 2 + (2B1 + 3B2 ) x + ( B0 + B1 + B2 ) = 3 x 2 + x + 4 3B2 = 3 B0 = 4 比较系数得 2B1 + 3B2 = 1 ⇒ B1 = −1 B + B + B = 4 B = 1 0 1 2 2 ~ = x(4 − x + x 2 ) 差分方程的特解为 yx
形如: 形如: yx+1 − ayx = f ( x) 齐次方程: 齐次方程: yx+1 − ayx = 0
(a ≠ 0,常数) 常数) (a ≠,常数) 0 常数)
(1)
的方程称为一阶常系数线性差分方程 的方程称为一阶常系数线性差分方程. 一阶常系数线性差分方程 (2) 1.齐次方程通解 yx+1 − ayx = 0 1.齐次方程通解 yx+1 = a, 通解为 y x = Aa x( A 是任何实常数). a≠0 是任何实常数) yx 例3 yx+1 − 3 yx = 0 是任何实常数) y x = A3 x ( A 是任何实常数). 通解: 通解: