数学模型差分方程【精选】
- 格式:ppt
- 大小:1.23 MB
- 文档页数:52
差分方程模型一. 引言数学模型按照离散的方法和连续的方法,可以分为离散模型和连续模型。
1. 确定性连续模型1) 微分法建模(静态优化模型),如森林救火模型、血管分支模型、最优价格模型。
2) 微分方程建模(动态模型),如传染病模型、人口控制与预测模型、经济增长模型。
3) 稳定性方法建模(平衡与稳定状态模型),如军备竞赛模型、种群的互相竞争模型、种群的互相依存模型、种群弱肉强食模型。
4) 变分法建模(动态优化模型),如生产计划的制定模型、国民收入的增长模型、渔业资源的开发模型。
2. 确定性离散模型1) 逻辑方法建模,如效益的合理分配模型、价格的指数模型。
2) 层次分析法建模,如旅游景点的选择模型、科研成果的综合评价模型。
3)图的方法建模,如循环比赛的名次模型、红绿灯的调节模型、化学制品的存放模型。
4)差分方程建模,如市场经济中的蛛网模型、交通网络控制模型、借贷模型、养老基金设置模型、人口的预测与控制模型、生物种群的数量模型。
随着科学技术的发展,人们将愈来愈多的遇到离散动态系统的问题,差分方程就是建立离散动态系统数学模型的有效方法。
在一般情况下,动态连续模型用微分方程方法建立,与此相适应,当时间变量离散化以后,可以用差分方程建立动态离散模型。
有些实际问题既可以建立连续模型,又可建立离散模型,究竟采用那种模型应视建模的目的而定。
例如,人口模型既可建立连续模型(其中有马尔萨斯模型Malthus、洛杰斯蒂克Logistic模型),又可建立人口差分方程模型。
这里讲讲差分方程在建立离散动态系统数学模型的的具体应用。
二. 差分方程简介在实际中,许多问题所研究的变量都是离散的形式,所建立的数学模型也是离散的,譬如,像政治、经济和社会等领域中的实际问题。
有些时候,即使所建立的数学模型是连续形式,例如像常见的微分方程模型、积分方程模型等。
但是,往往都需要用计算机求数值解。
这就需要将连续变量在一定的条件下进行离散化,从而将连续型模型转化为离散型模型。
差分方程对连续型变量而言,我们常常回导致到微分方程的问题. 对离散型变量将导致一类的问题.一、差分的定义定义 设)(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 .二、差分方程定义 设是含有未知函数差分的等式,称为差分方程。
差分方程建模举例差分方程建模方法的思想与与一般数学建模的思想是一致的,也需要经历背景分析、确定目标、预想结果、引入必要的数值表示(变量、常量、函数、积分、导数、差分、取最等)概念和记号、几何形式(事物形状、过程轨迹、坐标系统等),也就是说要把事物的性态、结构、过程、成分等用数学概念、原理、方法来表现、分析、求解。
当然,由于差分方程的特殊性,首先应当把系统或过程进行特别分解,形成表现整个系统的各个部分的离散取值形式,或形成变化运动过程的时间或距离的分化而得到离散变量。
然后通过内在的机理分析,找出变量所能满足的平衡关系、增量或减量关系及规律,从而得到差分方程。
另外,有时有可能通过多个离散变量的关系得到我们关心的变量的关系,这实际上建立的是离散向量方程,它有着非常重要的意义。
有时还需要找出决定变量的初始条件。
有时还需要将问题适当分成几个子部分,分别求解。
模型1 种群生态学中的虫口模型:在种群生态学中,考虑像蚕、蝉这种类型的昆虫数目的变化 ,他的变化规律是:每年夏季这种昆虫成虫产卵后全部死亡,第二年春天每个虫卵孵化成一个虫子。
建立数学模型来表现虫子数目的变化规律。
模型建立:假设第n 年的虫口数目为n P ,每年一个成虫平均产卵c 个(这个假设有点粗糙,应当考虑更具体的产卵分布状况),则有:n n cP P =+1,这是一种简单模型;如果进一步分析,由于成虫之间会有争斗以及传染病、天敌等的威胁,第n+1年的成虫数会减少,如果考虑减少的主要原因是虫子之间的两两争斗,由于虫子配对数为)1(21-n n p p 221n p ≈,故减少数应当与它成正比,从而有: 21n n n bP cP P -=+这个模型可化成:)1(1n n n x x x -=+λ,这是一阶非线性差分方程。
这个模型的解的稳定性可以用相应一阶差分方程的判断方法来获得。
如果还考虑其它的影响成虫孵卵及成活的因素的定量关系,这个模型在此基础上仍可进一步改进,更加符合实际情形。
数学建模课程设计实验报告题目:差分阻滞增长模型问题研究姓名:班级:日期:摘要该文以生物数量增长的预测为例,建立了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 阶常系数齐次线性差分方程。