数学建模第六章代数方程与差分方程模型模板
差分方程模型(讲义)
差分方程模型一.引言数学模型按照离散的方法和连续的方法，可以分为离散模型和连续模型。
1. 确定性连续模型1）微分法建模（静态优化模型），如森林救火模型、血管分支模型、最优价格模型。
2）微分方程建模（动态模型），如传染病模型、人口控制与预测模型、经济增长模型。
3）稳定性方法建模（平衡与稳定状态模型），如军备竞赛模型、种群的互相竞争模型、种群的互相依存模型、种群弱肉强食模型。
4）变分法建模（动态优化模型），如生产计划的制定模型、国民收入的增长模型、渔业资源的开发模型。
2. 确定性离散模型1）逻辑方法建模，如效益的合理分配模型、价格的指数模型。
2）层次分析法建模，如旅游景点的选择模型、科研成果的综合评价模型。
3）图的方法建模，如循环比赛的名次模型、红绿灯的调节模型、化学制品的存放模型。
4）差分方程建模，如市场经济中的蛛网模型、交通网络控制模型、借贷模型、养老基金设置模型、人口的预测与控制模型、生物种群的数量模型。
随着科学技术的发展，人们将愈来愈多的遇到离散动态系统的问题，差分方程就是建立离散动态系统数学模型的有效方法。
在一般情况下，动态连续模型用微分方程方法建立，与此相适应，当时间变量离散化以后，可以用差分方程建立动态离散模型。
有些实际问题既可以建立连续模型，又可建立离散模型，究竟采用那种模型应视建模的目的而定。
例如，人口模型既可建立连续模型(其中有马尔萨斯模型Malthus、洛杰斯蒂克Logistic 模型)，又可建立人口差分方程模型。
这里讲讲差分方程在建立离散动态系统数学模型的的具体应用。
差分方程简介在实际中，许多问题所研究的变量都是离散的形式，所建立的数学模型也是离散的，譬如，像政治、经济和社会等领域中的实际问题。
有些时候，即使所建立的数学模型是连续形式，例如像常见的微分方程模型、积分方程模型等。
但是, 往往都需要用计算机求数值解。
这就需要将连续变量在一定的条件下进行离散化，从而将连续型模型转化为离散型模型。
数学建模差分方程PPT课件
的方程都是差分方程。 方程中所含未知函数角标的最大值与最小值的差数称为差分
方程的阶。 若一个函数代入差分方程后，方程两端恒等，则称此函数为
差分方程的解。如果解中所含相互独立的任意常数的个数等于方 程的阶数，则称该解为差分方程的通解。满足初始条件的解称为 特解。
• 第一阶段: w(k)每周减1千克, c(k)减至下限10000千卡
w (k)w (k1)1 w ( k 1 ) w ( k ) c ( k 1 ) w ( k )
c(k1) 1[w(k)1] w (k)w (0)k
c(k1) w (0) 1(1k)
1 8000
0.025
120 200 k 00Cm 10000 k 10
2 x k 2 x k 1 x k 2 ( 1 ) x 0 , k 1 , 2 ,
二阶线性常系数差分方程
x0为平衡点 研究平衡点稳定，即k, xkx0的条件
模型的推广 2 x k 2 x k 1 x k 2 ( 1 ) x 0
方程通解
xk
c1
k 1
c2
k 2
(c1, c2由初始条件确定)
相当于70千克的人每天消耗2000千卡 ~ 3200千卡；
3）运动引起的体重减少正比于体重，且与运动 形式有关；
4）为了安全与健康，每周体重减少不宜超过1.5 千克，每周吸收热量不要小于10000千卡。
减肥计划
某甲体重100千克，目前每周吸收20000千卡热量， 体重维持不变。现欲减肥至75千克。
1 2 k 是（3）的 k 重根，则只要将 Y1 (i),Y2 (i),,Yk (i) 换为
差分方程模型ppt课件
P1( x1, y1 ), P2 ( x2 , y1), P3 ( x2 , y2 ), P4 ( x3, y2 ),
图上的箭头表示求出 Pk 的次序，由图知
lim
k
Pk
(
x,
y)
P0
(
x0
,
y0
)
即市场经济趋于稳定。
14
并不是所有的需求g 函数和供 应函数都趋于稳定，若给定
其中含 的最yt高阶差分的阶数称为该差分方程的阶数。
差分方程也t 可以写成不显含差分的形式，例如二阶差分
方程
2 yt yt yt可以0 写成
yt2 yt1 yt 0
2
满足一阶差分方程的序列 yt 称为差分方程的解，若 解中含有独立的常数的个数等于差分方程的阶数时，称 此解为该差分方程的通解。
(3) 下一时段的商品数量由上一时段的商品价格决定，
xk 1 g( yk )
称为供应函数，由于价格越高可导致产量越大，所以可以假 设供应函数是一个单调递增的函数。
12
3、模型求解
在同一坐标系中同时做出 供应函数和需求函数的图形 ，设两条曲线相交于 P0 (x0, y0 ) 则 P0为平衡点。因为此时
t3
最小。根据这一方程可以迭代求解以后各年第一 季度销售
23
量的预测值 y6 21, y7 19,。第7年销售量预测值居然小于第 6年的，稍作分析，不难看出，如分别对第一季度建立差分 方程，则根据统计数据拟合出的系数可能会相差甚大，但对 同一种商品，这种差异应当是微小的，故应根据统计数据建 立一个共用于各个季度的差分方程，为此，将季度编号为 t 1,2, 20，令 yt a1 yt4 a2 yt8 a3，利用全体数据来拟合 求拟合得到最好的系数。即求 a1, a2 , a3使得
差分方程模型
差分方程模型
周家全
对连续型变化的问题而言, 常常可建立微分方程模型. 而对离散状态转移的问题, 则可建立差分方程模型. 差分方 程与常微分方程有很多类似的性质和结论．首先引入差分的 概念．
1 差分定义及其性质
定义 设函数 y = y(x) 在等距节点 xi = x0 + ih ( i = 0,1, , n)
对于一般的差分方程 xn+2 + axn+1 + bxn = f 来讲, 其平衡 点的稳定性问题可以同样给出. 二阶方程的上述结果可以推
广到 n 阶线性差分方程, 即稳定平衡点的条件是特征根： n
次代数方程的根 λi (i = 1, 2, , n) 均有| λi |< 1.
4 经济学中的蛛网模型
1. 提出问题 在自由竞争的社会中, 很多领域会出现循环波动的现象. 在经济领域中, 可以从自由集市上某种商品的价格变化看到 如下现象：在某一时期, 商品的上市量大于需求, 引起价格 下跌, 生产者觉得该商品无利可图, 转而经营其它商品；一
解
Δf (0) = f (0.5) − f (0) = 0.75 ,
-2-
洛阳理工学院数学建模竞赛培训教案
Δf (0.5) = f (1) − f (0.5) = 1.25
周家全
Δ2 f (0)= Δ(Δf (0)) = Δf (0.5) − Δf (0) = 1.25 − 0.75 = 0.5
计算较多点的差分可按差分表进行, 容易看出表中每一 个需要计算的差分值分别等于其左侧的数减去左上侧的 数．每个点 xi 处的各阶差分位于与主对角线平行的斜线上．
(I) 先求解对应的特征方程
a0λn + a1λn−1 + + a0 = 0
7.数学建模-差分方程法
pt 发生动态等幅振荡;
ab t ) p* (5) 当 0 < ab < 2 , pt ( A1 sin kt A2 cos kt)( 2 ab ab t 1 ( ) 为衰减因子 2 2
pt → p*
( t → + ∞ ) , pt 动态发展趋于稳定 .
5．差分形式的生物数量 ic（阻滞增长）模型及其稳定性研究 描述生物生长受到环境约束的微分方程模型是 Logistic（阻滞增 长）模型 。其形式是 ： y
0
这时还贷公司需要还清银行的债务的时限变为:
b ln b ry0 x 503.5 ( 半月) 21年 . ln(1 r )
这表明还贷公司只用 21 年就可还清银行的债务, 由此 , 还贷公司赚 了购房人 一年的钱: 24 × 316 = 7584 ( 元 ) . 故问题 (2) 的解答是 : 此方案对还贷公司而言是有利可图的 。
模型II . 模型假设: (1) t 时刻的商品价格 pt 是商品数量 xt 的直线下降函数: pt = pM - a xt ; (2) 这一时期的商品数量 xt 是前两个时期的商品价格 pt-1 与 pt-2 的 算术平均值的直线上升函数（企业对市场的分析、判断应更成 b( pt 1 pt 2 ) 熟一些）: 模型建立:
p ( 0 ) = p0 ，p(1) = p1 ( 初始价格 ) . （二阶线性常系数差分方程）
r1, 2
ab ab(ab 8) 4
p M axm p* 1 ab
(2) 当 ab = 8 时,
ab t pt ( A1 A2 t )( ) p * ( A1 A2 t )(2) t p * 4 ab t ) p* (3) 当 ab < 8 时, pt ( A1 sin kt A2 cos kt)(
第六章微分方程与数学建模1
第一节 微分方程 第二节 微分方程在数学建模中的应用
第一节 微分方程
一、微分方程的基本概念 二、一阶微分方程 三、一阶线性微分方程及可降阶的高阶 微分方程 四、二阶常系数线性微分方程
一、微分方程的基本概念
1. 引例
例 一曲线通过点(1,2),且在该曲线上任一点
M ( x, y)处的切线的斜率为2 x,求这曲线的方程.
dp dy
原方程化为
p
dp dy
f
(
y,
p)
设方程
p dp dy
f (y, p)
的通解为
y
p
(x,C1)
则利用
分离变量法得原方程的通解为
(
dy y,C1)
x
C2
四、二阶常系数线性微分方程
二阶常系数线性微分方程的一般形式
二阶常系数齐次线性微分方程的形式
如果 f (x) 不恒等于 0，则方程
y1 e( i ) x , y2 e( i )x ,
重新组合
1
y1
( 2
y1
y2 )
ex cos x,
y2
1 2i
( y1
y2 )
ex sin x,
得齐次方程的通解为
y ex (C1 cos x C2 sinx).
2.二阶常系数非齐次线性微分方程
解 这是齐次方程，可以化成如下形式
dx x dy
x2 y2 x
y
y
x y
2
1
令
差分方程模型
设特解为 an D 代入 D 0.5D 0.1 得 D 0.2 ， 于是所求通解 an c(0.5) n 0.2 例3 (养老金) 解: 齐次特征方程 设特解 an D
an1 1.01an 1000
1.01 0,
* an c(1.01) n.
代入原方程得 D 100000
例 4 求非齐次差分方程
* 对应齐次方程的通解为 an c1 2n c2 n 2n
的通解
f (n) 2 中, 2 是2 重根, 设特解为
n
an A n 2 2 n
n 2 n1
代入
得 A 1 2 方法2 (化齐) :
故通解为 an c1 2 c2 n 2 n 2
Fn Fn 1 Fn 2 F1 F2 1
解：差分方程的特征方程为 x 2 x 1 0 特征根
x1
n
1 5 1 5 , x2 2 2
n
1 5 1 5 Fn c1 c2 2 2
n
2(an1 4an2 4an3 ) 2 2n1 相减得 an 6an1 12an2 8an3 0 特征方程 3 62 12 8 0 特征根 2 为三重根, 通解为：
an 4an1 4an2 2n
an c1 2n c2 n 2n c3n 2 2n
x k b1 x k 1 b2 x k 2 bk 0
称为差分方程的特征方程，其根称为特征根。 定理1（单根）若特征方程恰有k个相异的特 x1 , x2 ,, x 征根 ， k 则差分方程的通解为
an c x c x ck x
第六章 代数与建模
第六章 代数与建模§1 特征值与特征向量在层次分析法中的应用模型一、层次分析模型的基本原理层次分析法的本质就是多因素的权重确定方法。
设有n 个因素，确定它们相对某个事情、标准之下各自的重要性大小、各自所占的地位的大小的量化。
先分析一种特殊情况之下的权重的存在与确定方法：（1）单位重量大小的石头分解模型：将一个单位重量的石头分成n 块小石头C 1,C 2,…，C n精确称出它们的重量w 1,…w n,则向量就是每一个小石头在整个石头中的权重。
这是非常形象化的概念与意义。
Tn w w w w ),,(21L =为了从总体上显示这些权重以及其重要的性质，权重之间的比较应当是非常有价值的一种数学模型，因此将它们进行两两比较，形成矩阵：⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=n n n n n n w w w w w w w w w w w w w w w w w w A /////////212221212111L M M M M L L 如果记：，则有：Tn w w w w ),,(21L =nw Aw =，这表明权重向量是矩阵A 的对应于特征值n 的特征向量。
w 矩阵A 满足，这表明：位置处的元素为所在的行和所在的列上，在对称位置即列、行数相同处的两个元素的乘积，称这种矩阵A 为一致矩阵。
jkij ik a a a .=).(k i 一致矩阵A 的基本性质:1、A 有唯一非零特征值n 2、A 的任何一列向量都是对应于特征根n 的特征向量。
这表明：真正精确的权重，进行两两对比所形成的矩阵，必是一致矩阵。
其基本特点是：有唯一非零特征值n ；A 的任何一列向量都是对应于特征根n 的特征向量，并且在所有特征向量中，归一化的即分量之和为1的特征向量就是相应的权向量。
（2）一般情况下，如果对于各因素的权重不易具体确定时，可以先形成相互的对比量化结果，即权重的对比，这种对比显示了相对分析的模式，也是单一的描述到直接进行关系分析与显示思想的进展。