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.差分方程简介

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 )
这样的方程就是差分方程.

差分方程初步介绍

差分方程初步介绍
用Δ表示时差分的最高阶数)称为差分方程的阶。
例如 yt 5 yt 3 yt 2 0是 3 阶差分方程。
2 yt yt 0是 1 阶差分方程。
经济学中常用到的是以下标表示的差分方程。
第七章 方程类经济数学模型
返 回
上一张
下一张
退 出
2012年10月21日星期日
7
定义 若将函数yt=f(t)代入差分方程使其成为恒等式,则 称函数yt=f(t)为此差分方程的解。 对n阶差分方程,若它的解中含有n个(独立的)任意常数,
( yt zt ) z t y t y t z t z t 1 z t
第七章 方程类经济数学模型
返 回
上一张
下一张
退 出
2012年10月21日星期日
6
二、差分方程的基本概念
定义 含有自变量、未知函数及其差分的方程称为(常)差分
方程,出现在差分方程中的yi的下标的最大值与最小值之差(或
返 回 上一张 下一张 退 出
2012年10月21日星期日
2
定义
设函数yt=f(t)的定义域为Z+,其自变量t(通常表示时
间)取正整数数得到数列:
y0 f (0) y1 f (1) yt f (t )
称yt+1-yt为函数yt的一阶差分,记为 y t .
第七章 方程类经济数学模型
若a 1 ,则设特解为 t t (b0 t y
b1t
m 1
例 求方程yt 1 2 yt t的通解。 答案 yt C 2 t 1 t
t 0, 1, 2,
例 求方程yt 1 yt 2t的通解。 答案
yt C t t

差分方程

差分方程

第九节 差分方程迄今为止,我们所研究的变量基本上是属于连续变化的类型. 但在经济管理或其它实际问题中,大多数变量是以定义在整数集上的数列形式变化的,银行中的定期存款按所设定的时间等间隔计息,国家财政预算按年制定等等. 通常称这类变量为离散型变量. 对这类变量,我们可以得到在不同取值点上的各离散变量之间的关系,如递推关系等. 描述各离散变量之间关系的数学模型称为离散型模型. 求解这类模型就可以得到各离散型变量的运行规律. 本节将介绍在经济学和管理科学中最常见的一种离散型数学模型—差分方程.内容分布图示★ 引言 ★ 差分的概念 ★ 例1-5 ★ 差分方程的概念 ★ 例6★ 例7★ 一阶常系数线性齐次差分方程 ★ 一阶常系数线性非齐次差分方程★ 例9-14 ★ 例15★ 例16★ 二阶常系数线性差分方程★ 二阶常系数线性齐次差分方程的通解★ 例17 ★ 例18★ 例19 ★ 二阶常系数线性非齐次差分方程的特解★ 例20-23 差分方程在经济学中的应用★ 模型1 ★ 模型2★模型3★ 内容小结 ★ 课堂练习 ★ 习题8-9 ★ 返回内容要点:一、 差分的概念与性质一般地,在连续变化的时间范围内,变量y 关于时间t 的变化率是用dtdy来刻画的;对离散型的变量y ,我们常取在规定的时间区间上的差商ty∆∆来刻画变量y 的变化率. 如果选择1=∆t ,则)()1(t y t y y -+=∆ 可以近似表示变量y 的变化率. 由此我们给出差分的定义.定义1 设函数).(t y y t = 称改变量t t y y -+1为函数t y 的差分, 也称为函数t y 的一阶差分, 记为t y ∆, 即t t t y y y -=∆+1 或 )()1()(t y t y t y -+=∆.一阶差分的差分称为二阶差分t y 2∆, 即t t t t y y y y ∆-∆=∆∆=∆+12)(.2)()(12112t t t t t t t y y y y y y y +-=---=+++++类似可定义三阶差分, 四阶差分,……),(),(3423t t t t y y y y ∆∆=∆∆∆=∆一般地,函数t y 的1-n 阶差分的差分称为n 阶差分,记为t n y ∆,即t n t n t ny y y 111-+-∆-∆=∆i n t inni i y C -+=∑-=0)1( 二阶及二阶以上的差分统称为高阶差分.差分的性质:(1) t t y C Cy ∆=∆)( );(为常数C (2) ;)(t t t t z y z y ∆±∆=±∆ (3);)(1t t t t t t z y y z z y ∆+∆=⋅∆+ (4)t t t t t t t t z z z y y z z y ⋅∆-∆=⎪⎪⎭⎫ ⎝⎛∆+1 ).0(≠t z二、差分方程的概念定义2 含有未知函数t y 的差分的方程为差分方程. 差分方程的一般形式:0),,,,,(2=∆∆∆t n t t t y y y y t F或 .0),,,,,(21=+++n t t t t y y y y t G差分方程中所含未知函数差分的最高阶数称为该差分方程的阶. 差分方程的不同形式可以互相转化.定义3 满足差分方程的函数称为该差分方程的解.如果差分方程的解中含有相互独立的任意常数的个数恰好等于方程的阶数, 则称这个解为该差分方程的通解.我们往往要根据系统在初始时刻所处的状态对差分方程附加一定的条件,这种附加条件称为初始条件, 满足初始条件的解称为特解.定义4 若差分方程中所含未知函数及未知函数的各阶差分均为一次的, 则称该差分方程为线性差分方程.线性差分方程的一般形式是)()()()(1111t f y t a y t a y t a y t n t n n t n t =+++++--++其特点是t n t n t y y y ,,,1 +++都是一次的. 三、一阶常系数线性差分方程一阶常系数线性差分方程的一般形式为)(1t f Py y t t =-+ (9.1) 其中, P 为非零常数, )(t f 为已知函数. 如果,0)(=t f 则方程变为01=-+t t Py y (9.2)方程(9.2)称为一阶常系数线性齐次差分方程, 相应地,方程(9.1)称为一阶常系数线性非齐次差分方程.一阶常系数线性齐次差分方程的通解 一阶常系数线性非齐次差分方程定理1 设t y 为方程(9.2)的通解,*t y 为方程(9.1)的一个特解, 则*t t t y y y +=为方程(9.1)的通解.(1)C t f =)( (C 为非零常数)(2)t Cb t f =)( (C , b 为非零常数且1≠b ) 四、二阶常系数线性差分方程二阶常系数线性差分方程的一般形式:)(12t f by ay y t t t =++++ (9.9)其中b a ,均为常数, 且,0≠b )(x f 是已知函数. 当0)(=x f 时, 方程(9.9)变为012=++++t t t by ay y (9.10)方程(9.10)称为二阶常系数线性齐次差分方程,相应地,方程(9.9)称为二阶常系数线性非齐次差分方程.定理2 设t y 为方程((9.10)的通解, *t y 为方程(9.9)的一个特解, 则*t t t y y y +=为方程(9.9)的通解.二阶常系数线性齐次差分方程的通解特征方程 02=++b a λλ (9.11) 二阶常系数线性非齐次差分方程的特解和通解仅考虑方程(9.9)中的)(x f 取某些特殊形式的函数时的情形.(1))()(t P x f m =(其中)(t P m 是t 的m 次多项式), 方程(9.9)具有形如)(*t R t y m k t =的特解, 其中)(t R m 为t 的m 次待定多项式.五、 差分方程在经济学中的应用采用与微分方程完全类似方法,我们可以建立在经济学中的差分方程模型,下面举例说明其应用.1.“筹措教育经费”模型某家庭从现在着手, 从每月工资中拿出一部分资金存入银行, 用于投资子女的教育, 并计算20年后开始从投资账户中每月支取1 000元, 直到10年后子女大学毕业并用完全部资金. 要实现这个投资目标, 20年内要总共筹措多少资金? 每月要在银行存入多少钱? 假设投资的月利率为0.5%, 为此, 设第t 个月, 投资账户资金为,t a 每月存资金为b 元, 于是20年后, 关于,t a 的差分方程模型为1000)005.1(1-=+t t a a (9.11)且.,00120x a a ==二、价格与库存模型本模型考虑库存与价格之间的关系设)(t P 为第t 个时段某类产品的价格, )(t L 为第t 个时段的库存量. L 为该产品的合理库存量. 一般情况下, 如果库存量超过合理库存, 则该产品的售价要下跌, 如果库存量低于合理库存, 则该产品售价要上涨, 于是有方程)(1t t t L L k P P -=-+ (9.13)其中k 为比例常数.三、国民收入的稳定分析模型本模型主要讨论国民收入与消费和积累之间的关系问题.设第t 期内的国民收入t y 主要用于该期内的消费t G , 再生产投资t I 和政府用于公共设施的开支G (定为常数), 即有G I C y t t t ++= (9.17)又设第t 期的消费水平与前一期的国民收入水平有关, 即)10(1<<=-A Ay C t t (9.18)第t 期的生产投资应取决于消费水平的变化, 即有)(1--=t t t C C B I (9.19)由方程(9.17), (9.18), (9.19)合并整理得G BAy y B A y t t t =++---21)1( (9.20)于是, 对应A , B , G 以及,,0y y 可求解方程, 并讨论国民收入的变化趋势和稳定性.例题选讲:差分的概念与性质例1(讲义例1)设,2t y t =求 ).(),(),(32t t t y y y ∆∆∆例2(讲义例2)设.1),1()2)(1()0()(=+---=t n t t t t t n 求)(n t ∆.例3(讲义例3)求t t t y 32⋅=的差分. 例4 设,22t t y += 求.,,32t t t y y y ∆∆∆ 例5 试改变差分方程023=∆+∆t t y y 的形式. 差分方程的概念例6(讲义例4)试确定下列差分方程的阶..735)2(;0)1(15423=+=+-++--+t t t t t y y y y y例7(讲义例5)指出下列等式哪一个是差分方程, 若是, 进一步指出是否为线性方程..432)2(;33)1(12=+-+=∆-++t t t t t t y y y a y y一阶常系数线性差分方程例8(讲义例6)求差分方程031=-+t t y y 的通解. 例9(讲义例7)求差分方程231-=-+t t y y 的通解.例10(讲义例8)求差分方程tt t y y ⎪⎭⎫⎝⎛=-+233211在初始条件50=y 时的特解.例11(讲义例9)求差分方程2134t y y t t =-+的通解. 例12 求差分方程t y y t t πsin 341=++的通解.例13 求差分方程 t y y t t 231+=-+满足初始条件50=y 的特解. 例14(讲义例10)求差分方程t t t t y y 4221+=++的通解.例15 设某产品在时期t 的价格, 供给量与需求量分别为t t S P ,与),2,1,0( =t D t . 1当121+=t t P S , t t t t D S P D =+-=- 3,5421时, 求证(1) 由 3,2,1推出差分方程.221=++t t P P (2) 已知0P , 求上述差分方程的解.例16(讲义例11)在农业生产中, 种植先于产出及产品出售一个适当的时期, t 时期该产品的价格t P 决定着生产者在下一时期愿意提供市场的产量t t P S ,1+还决定着本期该产品的需求量,t Q 因此有1,-+-=-=t t t t dP c S bP a Q (a , b , c , d 均为正的常数)求价格随时间变动的规律. 二阶常系数线性差分方程例17(讲义例12)求差分方程04312=--++t t t y y y 的通解. 例18(讲义例13)求差分方程04412=++++t t t y y y 的通解. 例19(讲义例14)求差分方程04212=+-++t t t y y y 的通解.例20 求差分方程 12212=-+++t t t y y y 的通解及0,010==y y 的特解. 例21(讲义例15)求差分方程t y y y t t t =-+++4312的通解. 例22(讲义例16)求差分方程t t t t y y y 23212⋅=++++的通解. 例23 求差分方程tt t t y y y ⎪⎭⎫⎝⎛-=++++214112的通解. 差分方程在经济学中的应用课堂练习1.求差分方程21t y y t t =-+的通解.2.求差分方程t y y y t t t =-+++4312的通解.3.求差分方程t t t t y y y 57612⨯=-+++的通解.。

差分方程基础知识

差分方程基础知识

yt yt y
* t
C APt , 1 P A Ct ,
其中, A为任意常数,且当
P 1 时,

P 1 时,
C A y0 A1 , 1 P
A y0 A 1 .
例5 求差分方程
解 由于
yt 1 3 yt 2 的通解.
,故原方程的通解为
2.一阶常系数线性非齐次差分方程的通解
定理 设
yt
为齐次方程的通解, yt 为非齐次方程的一个
*
特解,则
yt yt yt* 为非齐次方程的通解.
y t 1 P y t 0
* * 证明 由题设,有 yt 1 Pyt f (t ) ,及
将这两式相加得 ( y t 1 yt*1 ) P ( y t yt* ) f (t ) ,即
其特点是
yt n , yt n 1,, yt
都是一阶的.
三 、一阶常系数线性差分方程 一阶常系数差分方程的一般方程形式为
yt 1 Pyt f (t )
其中 则方程变为
f (t ) 为已知函数.如果 f (t ) 0 P 为非零常数,
yt 1 Pyt 0
f (t ) 0 时方程
yt t ( n) t (t 1)(t 2) (t n 1) ,则
( n)
yt (t 1)
.
t
( n)
(t 1)t (t 1) (t 1 n 1)
t (t 1) (t n 2)(t n 1)
( n 1)
y3 Py2 P y0
t
, yt Pyt 1 P y0 .

数学建模差分方程问题

数学建模差分方程问题

数学建模差分方程问题数学建模是运用数学方法解决现实问题的一种方法。

而差分方程是数学建模中常用的一种数学工具,用于描述离散时间的动态系统。

本文将介绍差分方程的基本概念和应用,并以一个实际问题为例进行论述。

一、差分方程概述差分方程是一种用差分代替导数的方程,适用于离散时间的动态系统建模。

差分方程常用于描述离散时间下的变量变化规律,包括时序数据和动态优化等问题。

差分方程可以通过迭代求解来获得系统的演化过程。

二、差分方程的类型差分方程可分为线性差分方程和非线性差分方程两种类型。

线性差分方程的形式为:y(n+1) = a*y(n) + b*y(n-1)其中,y(n)表示第n个时间点的变量值,a和b为常数。

非线性差分方程的形式更加复杂,可以包含更多的项和参数,例如:y(n+1) = a*y(n)^2 + b*y(n-1) + c*n其中,y(n)^2表示y(n)的平方,c*n表示变量与时间的乘积。

三、差分方程的应用差分方程广泛应用于各个领域的实际问题,在科学研究、工程设计和金融市场等方面都有重要的应用价值。

下面以生态系统模型为例,来介绍差分方程的具体应用。

生态系统模型是生态学领域中的重要问题之一。

考虑一个简化的生态系统,由捕食者和被捕食者两个物种组成。

假设捕食者的数量为x,被捕食者的数量为y。

捕食者的增长速率与被捕食者的数量成正比,而被捕食者的减少速率与捕食者的数量成正比。

则可以建立如下差分方程模型:x(n+1) = x(n) + a*x(n)*y(n)y(n+1) = y(n) - b*x(n)*y(n)其中,a和b为模型的参数,表示捕食者与被捕食者之间的相互作用强度。

通过迭代求解这个差分方程模型,可以得到生态系统中捕食者和被捕食者数量的变化趋势。

四、差分方程的求解方法差分方程的求解可以通过数值方法进行。

常见的有欧拉法和龙格-库塔法等。

这些方法可以将差分方程转化为计算机程序进行求解,得到系统的近似解。

五、差分方程与其他数学工具的关系差分方程与微分方程是数学建模中常用的两种数学工具。

什么叫差分方程

什么叫差分方程

什么叫差分方程?给我举几个例子呗§1 基本理论1. 差分2. 任意数列{xn },定义差分算子Δ如下:Δxn=xn+1-xn对新数列再应用差分算子,有Δ2xn=Δ(Δkxn).性质性质1 Δk(xn+yn)=Δkxn+Δkyn性质2 Δk(cxn)=cΔkxn性质3 Δkxn=∑(-1)jCjkXn+k-j性质4 数列的通项为n的无限次可导函数,对任意k>=1,存在η,有Δkxn=f(k)(η) 差分方程定义8。

1 方程关于数列的k阶差分方程:xn-a1xn-1-a2xn-2-……aBxn-k=b(n=k,k+1,……)其中a1,a2,------ak 为常数,ak≠0. 若b=0,则该方程是齐次方程关于λ 的代数方程λk-a1λk-1-------ak-1λ-ak=0为对应的特征方程,根为特征值。

1.实验内容与练习2.1 插分例1 Xn={n3},求各阶差分数列:xn △xn △2xn △3xn △4xn1 7 12 6 08 19 18 6 027 37 24 6 064 61 30 6125 91 36216 127343可见,{n3},三阶差分数列为常数数列,四阶为0。

练习1 对{1},{n},{n2},{n4},{n5}, 分别求各阶差分数列。

练习2 {C0n-1}{C1n-1}{C2n-1},{C4n-1},分别求各阶差分数列.{Xn}的通项为n的三次函数,Xn=a3n3+a2n2+a1n+a0证明它为常数数列。

证明由Xn=a3n3+a2n2+a1n+a0可直接计算。

定理8。

1 若数列的通项是关于n 的k次多项式,则k 阶差分数列为非零数列,k+1阶差分数列为0。

练习3 证明定理8。

1 。

定理8。

2 若{Xn}的k 阶插分为非零常数列,则{Xn}是n的k次多项式,练习4 根据插分的性质证明定理8。

2例2。

求∑i3例3例4解设Sn=∑i3 表Sn △Sn △2Sn △3Sn △4Sn △5Sn1 8 19 18 6 09 27 37 24 6 036 64 61 30 6 0100 125 91 36 6 0225 216 127 42441 343 169784 5121296设Sn=a4n4+a3n3+a2n2+a1n+a0, s1=1,s2=9,s3=36,s4=100,s5=225,得a0=0, a1=0, a2=1/4, a3=1/2, a4=1/4.所以,Sn=(1/4)n4+(1/2)n3+(1/4)n2.练习{Xn}的通项Xn为n的k次多项式,证明∑xi为n的k+1次多项式;求∑i4.由练习2 {Crn-1}可得。

第l六节差分与差分方程的概念

第l六节差分与差分方程的概念
第六节 差分与差分方程的概念
机动
目录
上页
下页
返回
结束
当变量被认为是离散或间断地变化而不是连续或瞬时 地变化时,差分方程就适合表示这些变化之间的关 系,而微分方程就不适合。
在企业管理和经济分析中,差分方程常常是有用的。
下面介绍差分与差分方程的一些概念并介绍简单的差
分方程的解法。这对我们研究和解决一些实际问题是 颇有益处的。
yx 4( x2 ) (2) 4(2x 1) 0 8x 4
2 yx yx2 2 yx1 yx
2 4( x 2)2 2 2 4 x 1 2 4 x 2 2 8
或者
2 yx yx 8x 4 8
2 yx 1 ;当 n 2, yx x 则 yx 2x 1.
一般地,若 y xn ,则
k n k yx ( x 1) n x n Cn x k 1 n
机动
目录
上页
下页
返回
结束
例3 解
已知 y 4x 2 , 求 yx 和 yx .
2 2 3
(yx ) yx 1 yx ( yx 2 yx 1 ) ( yx 1 yx ) yx 2 2 yx 1 yx
称为函数 y y ( x)的二阶差分,记作 2 yx ,即
2 yx yx2 2 yx1 yx
实际上二阶差分是一阶差分的差分。
yx5 4 yx3 3 yx2 2 0
就是 ( x 5) ( x 2) 3 阶的,而不是5阶。事实上,作
x x 2 变换,便可得未知函数 yx 的差分方程
yx3 4 yx1 3 yx 2 0

差分方程初步

差分方程初步
返回 上页 下页
返回 上页 下页
定理3(齐次线性差分方程通解结构定理) 如 果 y1(t),y2(t),…,yn(t) 是 齐 次 线 性 差 分 方 程
yt+n+a1yt+n1 +a2yt+n2 +…+an1yt+1+anyt=0的n个线性无关 的特解,则方程 的通解为:
yA(t)=A1y1(t)+A2y2(t)+…+Anyn(t), 其中A1,A2,…,An为n个任意(独立)常数.
=
A( 1 )t 2
+
1 3
2t +1
A = A 2 为任意常数. 3
返回 上页 下页
2.待定系数法求特解
情形Ⅰ f(t)为常数.
方程变为yt+1+ayt=b, a,b均为非零常数.
试以 yt = (为待定常数)形式的特解代入方程得 +a =(1+a) =b.
当a≠-1时,可求得特解
由差分方程的通解来确定它的特解,需要给出确定 特解的定解条件.n阶差分方程F(t,yt,yt+1,…,yt+n)=0常 见的定解条件为初始条件.
y0=a0, y1=a1,…,yn1=an1, 这里a0,a1,a2,…,an1均为已知常数.
返回 上页 下页
只要保持差分方程中的时间滞后结构不变,无论对 t提前或推后一个相同的等间隔值,所得新方程与原方程 是等价的,即二者有相同的解.例如,方程
三、
返回 上页 下页
返回 上页 下页
一般地,k阶差分(k为正整数)定义为
Dk yt = D(Dk1 yt )
= Dk1 yt +1 Dk1 yt

差分方程性质

差分方程性质

Lax等价定理 讨论了差分问题的相容性、收敛性和稳定性。 相容性是收敛性的必要条件 ; 还发现 , 稳定性与 收敛性有一定的联系。Lax等价定理就是阐述相 容性、收敛性和稳定性三者之间的关系的。 Lax等价定理:对一个适定的线性微分问题及一 个与其相容的差分格式,如果该格式稳定则必收 敛,不稳定必不收敛。
i ui U ( xi , t j )
j j
称为离散化误差。如果当Δ x->0、Δ t->0 时 , 离散化误差的某种范数‖e‖趋近于零 , 即:
lim
x , t
0
则说明此差分格式是收敛的,即此差分格式的 解收敛于相应微分问题的解,否则不收敛。与相 容性类似,收敛又分为有条件收敛和无条件收 敛。
因此除了必需要求差分格式能逼近微分方程和定解条件外还进一步要求差分格式的解精确解与微分方程定解问题的解精确解是一致即当步长趋于零时要求差分格式的解趋于微分方程定解问题的解
截Байду номын сангаас误差和相容性
以FTBS格式为例
u t
j i
u C x
j i
0
j+1
u ij u ij1 u ij 1 u ij c 0 t x
i-1 i
j
u t
j i
1 2u 2 t 2
j i
2 u C u j j 2 t (t 2 ) C x ( x )0 i 2 i x 2 x
u t
j i
u j 1 2 u C i x 2 t 2
j i
C 2u j 2 2 t x ( t , x )0 2 i 2 x
离散误差与收敛性 所谓相容性,是指当自变量的步长趋于零时,差 分格式与微分问题的截误差的范数是否趋于零, 从而可看出是否能用此差分格式来逼近微分问 题。然而方程是物理问题的数学表达形式,其 目的是为了借助数学的手段来求问题的解。因 此,除了必需要求差分格式能逼近微分方程和 定解条件外,还进一步要求差分格式的解(精确 解)与微分方程定解问题的解(精确解)是一致 的。

数理经济学03-微分方程与差分方程

数理经济学03-微分方程与差分方程

数理经济学03-微分方程与差分方程微分方程与差分方程微分方程与差分方程简介本章简单地介绍微分方程、差分方程的一些基本概念和稳定性概念。

§ 2.1 微分方程的基本概念微分方程的定义及其阶在许多实际和理论问题中,需要寻找变量之间的函数关系。

一般来说,变量之间的函数关系很难直接求出,然而,根据以知条件,往往可以得到一个自变量、未知函数与它的导数之间的关系式。

因此,希望利用以知的函数与它的导数之间的关系式,去求出这个函数本身。

为此,给出下列描述性的定义:定义含有未知函数和未知函数各阶导数的等式称为微分方程。

在该等式中,若未知函数及其导数是一元函数,就称该微分方程是常微分方程。

若未知函数是多元函数,且该等式中所含的导数是偏导数,则称该微分方程是偏微分方程。

本章仅介绍常微分方程。

在下面,“微分方程” 一词,均是指常微分方程。

微分方程的一般形式是F(x,y,y , ,y(n)) 0其中,X是自变量,y是X的函数,y, ,y(n)是y对x的各阶导数。

微分方程的解、通解、特解和初始条件若函数(可以是显函数,也可以是隐函数) y y(x)满足该微分方程,即将y y(x),y y (x),,y(n)y(n)(x )代入到微分方程F(x,y,y, , y(n)) 0,能使等式成为恒等式,则称这个函数y y(x)是这个微分方程的解。

例假设曲线在点x处的切线斜率是2x。

求满足这一条件的所有曲线。

解:根据导数的几何意义,有y 2x 这是一个一阶微分方程。

两边同时积分,有y dx 2xdx x2 c 所以,该微分方程的解是y x2 c 由于一个函数对应平面上的一条曲线,故也常常称微分方程F(x, y, y , , y(n)) 0的解y y(x)是该微分方程的积分曲线。

上例的积分曲线如图2.1所示。

从图中可以看到,该微分方程有无穷多条积分曲线,并且,所有的积分曲线都可以通过其中的某一条积分曲线沿纵轴平行移动而得到。

一般来说,若一个微分方程有解,则它有无穷多个解,且这些解的图象互相平行。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
x( n) = x( x − 1)( x − 2)L( x − n + 1) , x( 0) = 1
例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

次多项式 通解 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 是任何实常数). 通解: 通解:
相关文档
最新文档