数模(差分方程模型)
数学建模中的差分方程模型
数学建模中的差分方程模型数学建模是一种将实际问题转化为数学模型并寻求与之相连的数学方法的学科，不仅仅在理论研究上有很大的应用，也在实际生活中有着广泛的应用。
在各种数学模型中，差分方程模型也是一种很重要的模型。
本文将结合实例，介绍差分方程模型的定义、建立、求解以及应用。
差分方程模型定义差分方程模型是一种通过离散化的方法，将连续时间问题转化为离散时间问题，来描述变量随时间的变化规律的数学模型。
这种数学模型以时间为自变量，以某个状态量为因变量，由一定的关系式组成。
例如：y(n+1)=ay(n)+b，式子中y(n)代表第n时刻系统状态，y(n+1)代表第n+1时刻系统状态，a和b为常数。
差分方程模型建立建立差分方程模型的关键是将实际问题中的连续变化离散化。
一般情况下，对于所建立的模型，首先要确定它的思路和范围，然后根据实际情况，确定差分方程的形式。
此外，还需要进行参数的估计和参数变化的分析，以及对模型精确性的验证。
以物理学中的简谐振动为例，建立一个差分方程模型描述其运动，即一个质点在回复力作用下以简谐运动形式振动。
设t为时间，y为质点的位移，v为质点的速度，a为质点的加速度，则有：$$y=n\Delta y \\v=\dfrac{y(n+1)-y(n-1)}{2\Delta t} \\a=\dfrac{y(n+1)-2y(n)+y(n-1)}{(\Delta t)^2}$$其中n为时间步长，$\Delta t$为时间间隔。
我们利用受力平衡的原理，即简谐振动中的$F=-ky$得到：$$\dfrac{y(n+1)-2y(n)+y(n-1)}{(\Delta t)^2} = -\dfrac{k}{m}y(n)$$将$\alpha=\dfrac{k}{m}$带入上式得到：$$y(n+1)-2(1+\alpha)y(n)+y(n-1) = 0$$此时，我们便成功地建立了描述简谐振动的差分方程模型。
差分方程模型求解对差分方程模型求解通常有两种方法：一种是使用递推公式进行求解，另一个方法是使用其它数学方法，如拉普拉斯变换或离散傅立叶变换等。
差分方程模型的基本概念
预测经济趋势
通过建立差分方程模型，可以对 未来的经济趋势进行预测，帮助 决策者制定相应的经济政策。
评估经济政策
差分方程模型可以用来评估不同 经济政策的实施效果，为政策制 定者提供参考依据。
在物理学中的应用
描述振动现象
差分方程模型可以用来描述物体的振动规律，如弹簧振荡、单摆 等。
预Байду номын сангаас波动传播
在声学和波动理论中，差分方程模型可以用来描述波动传播的规 律，如声波、电磁波等。
可以采用动态模型来反映数据的变化趋势，减少时间滞后的影 响。
可以利用大数据技术来处理大规模的数据集，提高模型的预测 精度和稳定性。
可以尝试优化参数估计方法，例如采用全局优化算法或贝叶斯 推断等方法，以提高参数估计的准确性和稳定性。
THANKS FOR WATCHING
感谢您的观看
确定差分关系
根据时间序列数据的特性，确定合适的差分关系，以描述数据的变化规律。差分关系通常表示为变量在不同时间 点的变化量或变化率。
建立差分方程模型
根据变量和参数建立模型
根据确定的变量和参数，建立差分方程模型，以描述变量的变化规律。
验证模型的适用性
建立差分方程模型后，需要验证模型的适用性，确保模型能够准确描述实际问题的变化规律。
Python
使用Python的数值计算库，如NumPy和 SciPy，求解差分方程。
Mathematica
使用Mathematica的符号计算和数值计算功 能求解差分方程。
04 差分方程模型的应用
在经济学中的应用
描述经济周期
差分方程模型可以用来描述经济 活动的周期性变化，如经济增长、 通货膨胀、就业率等的时间序列 数据。
第三章差分方程模型 ppt课件
输入必要信息 轻击鼠标即得
单利和复利 两种计算利息的基本方式
单利 ~１万元存5年定期, 年利率4.75%, 到期后本 息(本金加利息)：10000(1+0.04755)=12375元.
复利 ~１万元存1年定期, 年利率为3%, 到期不取则 自动转存, 5年后本息：10000 (1+0.03)5=11593元.
3. 差分方程模型
• 差分方程的基本类型及求解 3.1 贷款购房 3.2 管住嘴迈开腿 3.3 物价的波动 3.4 动物的繁殖与收获 3.5 中国人口增长预测——全国大学生
数学建模竞赛2007年A 题
差分方程的基本类型及求解
xk~未知变量x在时段k的数值(k=0,1,2, …)
1. 一阶线性常系数差分方程 xk 1 axk b, x0已知，k 0,1,2,
• 由x0, x1按照方程递推地计算x2, x3,…
•
求解公式
xk
c11k
c2k2
b 1 a1 a2
,
k 0,1,2,
1, 2～特征根 2 a1 a2 0 ~ 特征方程
c1, c2 ～常数, பைடு நூலகம்始值x0, x1代入求解公式确定.
1, 2<1
k→∞,
xk
x
1
b a1 a2
~稳定平衡点
3. 线性常系数差分方程组
x1(k), x2(k),, xn(k) ~n个未知变量在时段k的数值
x1(k 1) a11x1(k) a12x2 (k) a1n xn (k) b1 x2 (k 1) a21x1(k) a22x2 (k) a2n xn (k) b2 xn (k 1) an1x1(k) an2x2 (k) ann xn (k) bn
经济数学建模
经济数学模型
问题分析
商品信息传播一般有两个途径：
消费者外部信息:广告、亲眼看到商品等。
消费者内部信息:部分人使用并有所评价,使周围
人了解到有关产品信息。 由于是耐用消费品，所以一般不会重复购买， 故产品累计销售量可以认为是购买者人数。
经济数学模型
建模与求解
设K为潜在的消费者总数。
n(t)为t 时刻购买该产品的人数，在 [ t , t+Δt ]
经济数学模型
设厂家获得的利润为u, 每台电视机的生产成本 为c,销售价格为p,销售量为x, 则利润函数为 u = (p - c) x (3) 问题变化为在条件(1)(2)下求解利润函数的最大值。 构造拉格朗日函数
L( x, p, c, , ) ( p c) x
( x Me
说明： 在最优销售价格p*的表达式中含有待定的规模 参数k、价格系数a。为了确定电视机的最优销售价 格，必须预先给出这些参数。
S2 C3 S1 C4
经济数学模型
不妨假设市场容量函数为
M (t ) (1 e
, 0为常数
t
)
M(t)
t
0
t
将上式代入（3），再代入（1），得
dS1 t AS1 Be C dt
经济数学模型
解得
S1 (t ) K1e
同理
At
K 2e
经济数学模型
第二章 微分方程与差分方程模型
经济数学模型
模型一 利率模型
经济数学模型
设年利率为r,初始资金量为S0,n年后资金量为Sn
一、单利模型
n年后的本利和为
Sn S （ 0 1 nr）
差分方程模型
第7章 差分方程模型在第三章中我们看到，动态连续模型用微分方程方法建立，与此相应，当时间变化离散后，可以用差分方程方法建立。
有些实际问题既可建立连续模型，又可建立离散模型。
本章将主要讨论差分方程模型。
7.1市场经济中的蛛网模型在自由贸易市场上你注意过这样的现象吗：一个时期以来某种消费品如肉的上市量远大于需求，由于销售不畅导致价格下跌，生产者发现养猪赔钱，于是转而经营其他农副业。
这一段时间猪肉上市量就会大减，供不应求将导致价格上涨。
生产者看到有利可图，有重操旧业。
这样在下一个时期会重现供大于求、价格下降的局面、在这种没有外界干预的情况下，这种现象将如此循环下去。
在完全自由竞争的市场经济中上述现象通常是不可避免的，因为商品的价格是由消费者的需求关系决定的，商品数量越多价格越底，而下一时期商品的数量有生产者的需求关系决定的，商品的价格越低生产的数量越少，这样的需求和供应关系决定了市场经济中商品的价格和数量必然是振荡的。
在现实世界里这样的振荡会出现不同的形式，有的振荡渐小趋向平稳，有的则振幅越来越大，没有外界如政府的干预，将导致经济崩溃。
本节先用图形方法建立所谓“蛛网模型”，对上述现象进行分析，给出市场经济区域稳定的条件，再利用差分方程建模，对结果进行解释，并讨论当市场经济不稳定时政府可以采取什么样的干预措施，最后对上述模型做适当推广。
蛛网模型 记第k 时段商品的数量为k x ，价格为k y ，1,2,3,k = 这里我们把时间离散化为时段，1个时段相当于商品的1个生产周期，如蔬菜、水果是一个种植周期，肉类是牲畜的饲养周期。
同一时段商品的价格k y 取决于数量k x ，设()k k y f x = (1)它反映消费者对这种商品的需求关系，称需求函数，因为商品的数量越多价格越低，所以在图1中用一条下降曲线f 表示它，f 称需求曲线。
下一时段商品的数量1k x +由上一时段价格k y 决定，设1()k k x h y +=,或1()k k y g x += (2)这里g 是h 的反函数，反映生产者的供应关系，称供应函数，因为价格越高生产量越大，所以在图中供应曲线g 是一条上升的曲线。
差分方程模型
差分方程模型
周家全
对连续型变化的问题而言, 常常可建立微分方程模型. 而对离散状态转移的问题, 则可建立差分方程模型. 差分方 程与常微分方程有很多类似的性质和结论．首先引入差分的 概念．
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
数学建模之差分方程
差分方程对连续型变量而言，我们常常回导致到微分方程的问题. 对离散型变量将导致一类的问题.一、差分的定义定义 设)(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 .二、差分方程定义 设是含有未知函数差分的等式，称为差分方程。
差分方程
第八章 差分方程模型在经济与管理及其它实际问题中，许多数据都是以等间隔时间周期统计的。
例如，银行中的定期存款是按所设定的时间等间隔计息，外贸出口额按月统计，国民收入按年统计，产品的产量按月统计等等。
这些量是变量，通常称这类变量为离散型变量。
描述离散型变量之间的关系的数学模型成为离散型模型。
对取值是离散化的经济变量，差分方程是研究他们之间变化规律的有效方法。
下面介绍差分方程的基本概念、解的基本定理及其解法，与微分方程的基本概念、解的基本定理及其解法非常类似，可对照微分方程的知识学习本章内容。
函数的差分对离散型变量，差分是一个重要概念。
下面给出差分的定义。
设自变量t 取离散的等间隔整数值：，，，， 210±±=t t y 是t 的函数，记作)(t f y t =。
显然，t y 的取值是一个序列。
当自变量由t 改变到1+t 时，相应的函值之差称为函数)(t f y t =在t 的一阶差分，记作t y ∆，即)()1(1t f t f y y y t t t -+=-=+∆。
由于函数)(t f y t =的函数值是一个序列，按一阶差分的定义，差分就是序列的相邻值之差。
当函数)(t f y t =的一阶差分为正值时，表明序列是增加的，而且其值越大，表明序列增加得越快；当一阶差分为负值时，表明序列是减少的。
按一阶差分的定义方式，我们可以定义函数的高阶差分。
函数)(t f y t =在t 的一阶差分的差分为函数在t 的二阶差分，记作t y 2∆，即)()()(11212t t t t t t t t y y y y y y y y ---=-==++++∆∆∆∆∆t t t y y y +-=++122。
依次定义函数)(t f y t =在t 的三阶差分为t t t t t t t y y y y y y y ∆∆∆∆∆∆∆∆+-=-==+++12212232)(t t t t y y y y -+-=+++12333。
数学建模差分方程模型
yk
x k 1 bk(1 x x k) (2 )
记br1 一阶(非线性)差分方程
(1)的平衡点y*=N
(2)的平衡点 x* r 11 r1 b
讨论 x* 的稳定性
补充知识
一阶非线性差分方程 xk1f(xk)(1)的平衡点及稳定性 (1)的平衡点 x*——代数方程 x=f(x)的根 (1)的近似线性方程 x k 1 f(x * ) f(x * )x k ( x * )( 2 ) 稳定性判断 x*也是(2)的平衡点
需求函数不变 y k y 0 (x k x 0 ) 2 x x x 2 ( 1 ) x , k 1 , 2 ,
k 2 k 1 k
0
二阶线性常系数差分方程
x0为平衡点 研究平衡点稳定，即k, xkx0的条件
模型的推广 2 x k 2 x k 1 x k 2 ( 1 ) x 0
• 运动(内容同前) C 80 0 0 .00 2 78 5 16(千 80 )
3 差分形式的阻滞增长模型
连续形式的阻滞增长模型 (Logistic模型)
x(t) ~某种群 t 时刻的数量(人口)
x (t)rx(1 x) N
t, xN, x=N是稳定平衡点(与r大小无关)
离散
yk ~某种群第k代的数量(人口)
y
g
需求曲线变为水平 y0 以行政手段控制价格不变
0
2. 使 尽量小，如 =0 y
供应曲线变为竖直
靠经济实力控制数量不变
0
f
x g
f
x0
x
模型的推广 生产者管理水平提高 xk1h(yk)
• 生产者根据当前时段和前一时 段的价格决定下一时段的产量。
xk1
h
第二讲 简介_差分方程模型与案例
第二讲简介_差分方程模型与案例1.1 差分方程的基本定义差分方程是在离散时段上描述现实世界中变化过程的数学模型。
现实中的问题通常是连续变化的，但我们常常只能在离散的时间点上对其进行观测和描述。
为了表述这一类的数学模型，我们引入了差分方程的方法。
1.2 一阶线性常系数差分方程一阶线性常系数差分方程的一般形式差分方程的平衡点差分方程的解平衡点稳定的条件1.3高阶线性常系数差分方程n 阶线性常系数差分方程的一般形式称方程011...0k n k n n k a x a x a x ++-+++= (*)为对应的齐次方程设方程(*)有形如kk x λ=的解,则有 特征方程特征根(假定有k 个不同的实根,其它情形参见资料《差分方法建模理论与案例》)平衡点差分方程的解平衡点稳定的条件所有特征值的模均小于11.4非线性差分方程建模案例：题目1 濒危物种的自然演变和人工孵化问题：Florida 沙丘鹤属于濒危物种，生态学家估计它在自然环境下，年平均增长率为 -3.24%。
如果在某自然保护区内开始有100只鹤，每年人工孵化5只鹤放入该保护区，建立描述其数量变化规律的模型，并作数值计算。
模型及其求解记第k 年沙丘鹤的数量为x k ，自然环境下的年平均增长率为r ，记a =1+r ，每年孵化的数量为b ，则第k +1年鹤的数量为,2,1,0,1,1=+=+=+k r a b ax x k k模型分析讨论时间充分长以后沙丘鹤数量的变化趋势, 即k →∞时x k 的极限状态。
,2,1,0,11)1(010=--+=++++=-k a ab x a aa b x a x kkk k k当a <1即r <0时x k →x =b /(1-a )。
MATLAB演示计算计算并作图，程序如下：子程序：function x=exf11(x0,n,r,b)% 建立名为exf11的函数M文件，x0,n,r,b可调节a=1+r;x=x0; % 赋初值for k=1:nx(k+1)=a*x(k)+b; % 按照（3）迭代计算end主程序:clc; clf; clear all% x0:初始值；r=-0.0324; b:人工孵化数x0=100;n=20;b=5; r=-0.0324;% 给定x0,n,r，b,调用exf11计算k=(0:n)';y=exf0201(x0,n,r,b);plot(k,y,'r*-'),title('Florida 沙丘鹤数量变化趋势');% 在图上做标记（运行结果显示）题目2 汽车租赁公司的运营问题:一家汽车租赁公司在3个相邻的城市运营，为方便顾客起见公司承诺，在一个城市租赁的汽车可以在任意一个城市归还。
数学建模中的差分法
西北大学数学系
三 常微分方程向差分方程转化（数值解）
1 Euler 方法 求初值问题的近似解。
先把自变量所在的区间 n 等分；
dx
dt
f (t, x)
x(t0 ) x0
t1 t0 t t2 t1 t x f (t, x)t
tn tn1 t
1 1, 2 1
时，方程（4）的平衡点是稳定的。
非齐次线性方程（5）的稳定性可转化为齐次方 程（4）来研究。
xk2 a1xk1 a2 xk b,
(5)
对于n阶线性方程平衡点稳定的条件是特征根
i 1,(i 1,2,n,)
西北大学数学系
3 一阶非线性差分方程
xk1 f (xk )
由于
xk 1
b 1 a
a( xk
b) 1 a
0,
(1) k 0,1,2,
方程（1）平衡点的稳定性问题可转化为下面 方程零点的稳定性。
xk1 axk 0, k 0,1,2, (2)
方程（2）的解可表示为
xk (a)k x0 , k 1,2,
可得到下面的稳定性结论。
(6)
平衡点 x 通过求解方程 x f (x)
而得到。
研究稳定性的方法之一是研究其对应的线性部 分的稳定性。
将方程（6）的右端在 x 点作泰勒展开只取 一次项， （6）近似为
xk1 f (x )( xk x ) f (x )
(7)
x 也是（7）平衡点。
西北大学数学系
xk1 axk b, k 0,1,2,,
(1)
满足方程 x ax b 的解，称为上方程的平衡点。
-数学建模-差分方程--混沌现象的起源
r 1 1 (2)的平衡点 x r 1 b
*
讨论 x* 的稳定性
补充知识
一阶非线性差分方程 xk 1 f ( xk ) (1) 的平衡点及稳定性 (1)的平衡点 x*——代数方程 x=f(x)的根
* * * x f ( x ) f ( x )( x x ) (2) (1)的近似线性方程 k 1 k
差分方程模型
市场经济中的蛛网模型
混沌-----差分形式的阻滞增长 模 型
市场经济中的蛛网模型
供大于求 价格下降
数量与价格在振荡 增加产量 价格上涨 供不应求
减少产量
现 象
问 题
描述商品数量与价格的变化规律 商品数量与价格的振荡在什么条件下趋向稳定
当不稳定时政府能采取什么干预手段使之稳定
蛛网模型
xk~第k时段商品数量；yk~第k时段商品价格
t, xN, x=N是稳定平衡点(与r大小无关)
离散 形式
yk ~某种群第k代的数量(人口)
yk yk 1 yk ryk (1 ), k 1,2, N y*=N 是平衡点 若yk=N, 则yk+1,yk+2,…=N
讨论平衡点的稳定性，即k, ykN ？
离散形式阻滞增长模型的平衡点及其稳定性
k x x ( ) ( x1 x0 ) xk 1 x0 ( xk x0 ) k 1 0
1 ( 1 / )
xk x0 xk
P0稳定 K f K g P0不稳定 K f K g
1 ( 1 / )
方程模型与蛛网模型的一致
二阶线性常系数差分方程 x0为平衡点 研究平衡点稳定，即k, xkx0的条件
数学建模差分方程
数学建模课程设计实验报告题目：差分阻滞增长模型问题研究姓名：班级：日期：摘要该文以生物数量增长的预测为例，建立了logistic 阻滞增长的微分方程模型，并把它离散化而得相应的差分方程模型。
将logistic 阻滞增长模型的差分形式进行简化并对简化差分形式进行迭代求解。
做出随固有增长率的变化， 按logistic 阻滞增长模型的差分形式增长的序列{k y }收敛、 2倍周期、 4倍周期......直至一片混乱的图形。
以参数b 为横坐标、 序列{k y }的收敛点为纵坐标，用数学软件模拟展示了这一简单差分方程从收敛、分叉、倍周期收敛进入混沌现象的过程。
为部分工程领域的混沌现象的研究提供了模拟方法。
关键词：logistic 模型 分叉 倍周期收敛 混沌现象一：问题分析生物数量在增长过程中，由于环境因素与自然资源的作用，受到阻滞。
此时，其增长率呈现递减趋势。
基于此的logistic 模型可以对类似问题进行分析。
但是，现实对象的活动一般都是具有周期性的，所以采用离散化的时间比采用连续的时间更为方便，于是采用差分形式的离散模型。
对于平衡点的稳点问题，我们知道，logistic 模型中x*=N 是稳定平衡点，x*=0不是稳定平衡点，那么对于差分形式的离散模型)1(1k k k x bx x -=+，k=0,1,2，... 是否还具有同样的性质?以下，我们将对模型从平衡点和稳定性的角度进行分析并借助计算机对倍周期收敛、分岔和混沌的现象进行分析；二：模型假设（1）自然资源，环境条件等对生物的增长起着阻滞作用，并随着数量的增加阻滞作用越来越大。
（2）自然资源与黄精条件所容纳的最大生物数量，现有生物数量和固有增长率已知。
（3）阻滞作用体现在对增长率的影响上，使得增长率随着生物数量表的增加而下降。
（4）所研究该对象每年有固定的周期性活动。
三：模型建立与图解模型建立设当前（即 t=0时） 生物数量为0x ,固有增长率为r 生物数量为x 。
差分方程数学建模分析
规 律 和运 算关 系等式 ，建 立起 差分 方 程 。 １２ 对事 物 系统进 行 划分 ，划分 成若干 子 系统 ，在 每个 子 系统 ． 中 引入恰 当的 变量 或 向量 ，然 后分析 建 立起 子过 程 间的这 种 量 的关 系等式 ，从而 建立 起 差分 方程 ， 着重 强调 的是在 此过 程 中 ，时段 或子 系统 的 划分方 式是 非 常非 常 重要 的 ，应 当结 合 已有 的信 息和 分 析条 件 ，从 多种可 选方 式 中挑 选 易于 分 析 、针对 性 强的 划分 ，同时 ，对 划 分后 的时 段或 子过 程 ，引入 哪些 变量 或 向量 都是 至关 重要 的 ，要 仔 细分 析 、选择 ，尽 量 扩大 对 过程 或 系统 的数 量感 知范 围 ，包 括对 已有的 、 已知的 若干 量 进行 结 合运 算 、取 最运 算等 处理 方式 ， 目的 是 建立 起 简洁 、深 刻 、易于 求解 分析 的 差分方 程 。 ２ 模 型 举 例
Ｙ ＝ｆ ｈｙ ） ［（ 】
这就 是 两个 差分 方程 ， 属一 阶非 线性 差分 方程 。
于是 ２．１ ｘ＋＋
即 ２． ＋ ｘ＋ ２
一＋ａ ， ＝（＋ａ ）ｏ ｌ ｐｘ １ ￣ Ｘ
＋＋ ｌ ：（＋ ）。 １
３几 何模 型 分析
为了表现出两个变量 Ｘ 和 的变化过程， 我们可以借助已有的函 数 ｆ 和 ｇ ， 过对应 关系 的几何表现 把点列 （ ， ，和 ｎ 在 通 Ｙ） ＋Ｙ） ｌ 坐标 系 中描 绘 出来 ，进而 分析 它们的变化规 律 、 势、找稳 定点． 中 趋 其 （ｎＹ ） Ｘ， （ ） Ｙ） Ｘ ｇ ＋ ） Ｘ ， ＝（ｎｆ ｘ） （ ， ：（ｎ （ ） 将 点 列 Ｐ（ Ｉ 。 Ｐ （２Ｙ）Ｐ （３Ｙ ）Ｐ （４Ｙ ） … 一 接起来 ， １Ｘ， ） ２Ｘ １ ３Ｘ ３ ｘ ，，… Ｙ， ， ， 连 就会形成 象蛛 一样的折线，这个图形被称作为蛛网模型。
第三章_差分方程模型
第三章 差分方程模型§1、 差分方程设有未知序列{}k y ，称0),,,;(1=++n k k k y y y k F （1）为n 阶差分方程。
若有)(k y y k =，满足0))(,),1(),(;(=++n k y k y k y k F则称)(k y y k =是差分方程（1）的解，包含n 个任意常数的解称为（1）的通解， 当110,,,-n y y y 为已知时，称其为（1）的初始条件，通解中的任意常数都由初始条件确定后的解称为（1）的特解。
[例1] 设第一月初有雌雄各一的一对小兔，假定两月后长成成兔，同时即第三月开始每月初产雌雄各一的一对小兔，新增小兔也按此规律繁殖。
设第k 月末共有k y 对兔子，试建立关于k y 的差分方程。
[解] 因为第2+k 月末的兔子包括两部分，一部分为上月留下的，另一部分为当月新生的，而由题设当月生的小兔数等于前月末的兔数，所以有⎩⎨⎧==+=++1,01012y y y y y k k k 这是著名的裴波那契数列。
[例2] 汉诺塔问题将k 个大小不同的圆盘依其半径大小依次套在桩A 上，大的在下，小的在上。
现将此k 个盘移到空桩B 或C 上，但要求一次只能移动一个盘且移动过程中，始终保持大盘在下，小盘在上，移动过程中桩A 也可利用。
设移动k 个盘的次数为k y ，试建立k y 的差分方程。
[解] 先将桩A 上的k 个大小不同的圆盘按题设要求移到C 上，这需要移动k y 次，再将A 上的最大盘移到B 上，这需要移动一次，最后将C 上的k 个盘按要求移到B 上，这又需要移动k y 次。
所以，差分方程为⎩⎨⎧=+=+01201y y y k k§2、 差分方程的解法一.常系数线性齐次差分方程形如 0110=+++-++k n n k n k y a y a y a ——（1）其中n a a a ,,,10 为常数，且0,00≠≠n a a ，称为n 阶常系数齐次线性差分方程。
数学建模方法之差分方程模型
数学建模方法之差分方程模型差分方程模型是数学建模中常用的一种方法，它基于差分方程来描述问题，并用差分方程来求解问题。
所谓差分方程，是指用差分代替微分的方程，它是一种离散的模型。
在实际问题中，很多情况下，并不能直接通过微分方程来描述问题，而差分方程模型则可以通过离散化的方法来近似地描述问题。
差分方程模型的优点之一是可以适用于离散化的数据，对于实际问题的离散化模型建立是非常有帮助的。
差分方程模型的另一个优点是可以通过数值方法来求解，不需要进行繁琐的解析推导，因此适用于复杂问题的求解。
差分方程模型的基本形式为：yn+1 = fn(yn, yn-1, ..., yn-k)其中，yn表示第n个时刻的解，fn是一个给定的函数，表示通过前k个时刻的解来计算第n+1个时刻的解。
这个方程是离散的，通过已知的初始条件来逐步递推获得结果。
差分方程模型的适用范围非常广泛，可以用于描述和预测各种动态过程。
例如，差分方程模型可以用来描述人口增长模型、生态系统模型、传染病模型等等。
在这些例子中，差分方程模型可以通过已知的数据和初始条件来预测未来的发展趋势。
差分方程模型的建立步骤主要包括以下几个方面：1.确定问题的描述和目标：明确问题的背景和目标，确定需要建立差分方程模型的原因和用途。
2.确定模型的变量和参数：根据实际问题，确定需要用到的变量和参数。
3.确定差分方程的形式和函数：根据问题的特点和要求，选择合适的差分方程形式和函数。
这部分需要结合实际问题和数学方法来确定。
4.确定初始条件和边界条件：确定差分方程模型的初始条件和边界条件。
这部分是求解差分方程的前提条件。
5.差分方程的求解和分析：通过数值方法求解差分方程，得到数值解，并对结果进行分析和解释。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
Fibonacci 数列
数学建模
问题
13世纪意大利著名数学家Fibonacci在他的著作《算盘书》 中记载着这样一个有趣的问题：
一对刚出生的幼兔经过一ຫໍສະໝຸດ 月可长成成兔，成兔再经过一个月后可以繁殖出一对幼兔. 若不计兔子的死亡数，问一年之 后共有多少对兔子？
月份 0 1 2 3 4 5 6 7 …
幼兔 1 0 1 1 2 3 5 8 … 成兔 0 1 1 2 3 5 8 13 … 总数 1 1 2 3 5 8 13 21 …
数学建模 将兔群总数记为 fn, n=0,1,2,…，经过观察可以发现，数列{fn} 满足下列递推关系：
f0 = f1 =1, fn+2 = fn+1 + fn , n=0,1,2,…
这个数列称为Fibonacci数列. Fibonacci数列是一个十分有趣 的数列，在自然科学和数学领域中都有着广泛的应用.
Fibonacci数列的一些实例. 1. 蜜蜂的家谱 2. 钢琴音阶的排列 3. 树的分枝 4. 杨辉三角形
数学建模
日常的经济问题中的差分方程模型
1. 银行存款与利率
假如你在银行开设了一个1000元的存款账户，银行的年利 率为7%. 用an表示n年后你账户上的存款额，那么下面的数列 就是你每年的存款额：
（7.1）
a0 xnt a1xnt1 ... an xt 0
（7.2）
容易证明，若序列
x (1) t
与
x(2) t
均为方程（7.2）的解，则
xt c1xt(1) c2 xt(2)
也是方程（7.2）的解，其 中c1、c2为任意常数，这说明， 齐次方程的解构成一个 线性空间（解空间）。
此规律对于（7.1）也成立。
...
Cnnn
(C1,…,Cn为任意常数)
情况2 若λ 是特征方程（7.3）的k重根，通解中对应
于λ的项为 (C1 L Ck nk1 ) n
C i为任意常数，i=1,…,k。
情况3 若特征方程（7.3）有单重复根 a i
通解中对应它们的项为 C1 n cosn C2 n sinn
2an= an+1- an=(an+2-an+1)-(an+1-an)=an+2-2an+1+an 一般地，可以定义n阶差分.
差分方程 an+1= an+b的解
数学建模
定理1 一阶线性差分方程 an+1= an+b 的通解是：
nb c,
1,
an
cn
b
1
,
1.
定理2 对一阶线性差分方程 an+1= an+b, 若 | |<1, 则 an无限趋近于平衡解 b/(1- ) (收敛型不动点); 若 | |>1, 则 an逐渐远离平衡解 b/(1- ) (发散型不动点).
a0=c, an+1=an+b, n=0,1,2,3,… 称此类递推关系为一阶线性差分方程. 当b=0时称为齐次差分方
程，否则称为非齐次差分方程.
定义1 对任意数列A={a1,a2,…,an,…}，其差分算子定义如下： a1=a2-a1, a2=a3-a2,… an=an+1-an, …
定义2 对数列A={a1,a2,…,an,…}，其一阶差分的差分称为二 阶差分, 记为2A=(A). 即：
数学建模
第七章 差分方程模型
7.1 差分方程基本知识 7.2 市场经济中的蛛网模型 7.3 减肥计划——节食与运动 7.4 差分形式的阻滞增长模型 7.5 按年龄分组的种群增长
数学建模
7.1 差分方程基本知识
• 1、差分方程： 差分方程反映的是关于离散变量的取值 与变化规律。通过建立一个或几个离散变量取值所满 足的平衡关系，从而建立差分方程。
设贷款额为a0，每月还贷额为x，月利率为r，第n个月后的欠 款额为an，则
a0=200000, a1=(1+r)a0-x, a2=(1+r)a1-x, ……
an=(1+r)an-1-x, n=1,2,3,…
一阶线性差分方程
数学建模
在上述模型中，给出了an+1与an之间的递推公式. 将它们写成
统一的形式：
数学建模
方程（7.1）可用如下的代数方法求其通解：
（步一）先求解对应的特征方程
a0n a1n1 ... an 0
（7.3）
（步二）根据特征根的不同情况，求齐次方 程(7.2）的通解 情况1 若特征方程（7.3）有n个互不相同的实根
,…1 , ，n则齐次方程（7.2）的通解为
C ， n 11
设n年后教育基金总额为an，每年向银行存入x元，依据复利 率计算公式，得到家庭教育基金的数学模型为：
a0=x, an+1=(1+r)an+x, n=0,1,2,3,…
3 . 抵押贷款
数学建模
小李夫妇要购买二居室住房一套，共需30万元. 他们已经筹 集10万元，另外20万元申请抵押贷款. 若贷款月利率为0.6%， 还贷期限为20年，问小李夫妇每月要还多少钱？
• 差分方程就是针对要解决的目标，引入系统或过程中 的离散变量，根据实际背景的规律、性质、平衡关系， 建立离散变量所满足的平衡关系等式，从而建立差分 方程。通过求出和分析方程的解，或者分析得到方程 解的 特别性质（平衡性、稳定性、渐近性、振动性、 周期性等），从而把握这个离散变量的变化过程的规 律，进一步再结合其他分析，得到原问题的解。
数学建模
a0 (t)xtn a1(t)xtn1 .... an (t)xt 0
则被称为方程对应的 齐次线性差分方程 。
若所有的 ai(t)均为与t无关的常数，则称其为 常系数差分 方程，即n阶常系数线性差分方程可分成
a0 xnt a1xnt1 ... an xt b(t)
的形式，其对应的齐次方程为
a0, a1, a2, a3, …, an,… 设r为年利率，由于an+1=an+r an, 因此存款问题的数学模型 是：
a0=1000, an+1=(1+r)an, n=1,2,3,…
2. 家庭教育基金
数学建模
从1994年开始，我国逐步实行了大学收费制度. 为了保障子女 将来的教育费用，小张夫妇从他们的儿子出生时开始，每年向 银行存入x元作为家庭教育基金. 若银行的年利率为r，试写出第 n年后教育基金总额的表达式. 预计当子女18岁入大学时所需的 费用为100000元，按年利率3%计算，小张夫妇每年应向银行存 入多少元?