差分方程模型
- 格式:pdf
- 大小:127.50 KB
- 文档页数:8
数学建模中的差分方程模型数学建模是一种将实际问题转化为数学模型并寻求与之相连的数学方法的学科,不仅仅在理论研究上有很大的应用,也在实际生活中有着广泛的应用。
在各种数学模型中,差分方程模型也是一种很重要的模型。
本文将结合实例,介绍差分方程模型的定义、建立、求解以及应用。
差分方程模型定义差分方程模型是一种通过离散化的方法,将连续时间问题转化为离散时间问题,来描述变量随时间的变化规律的数学模型。
这种数学模型以时间为自变量,以某个状态量为因变量,由一定的关系式组成。
例如: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$$此时,我们便成功地建立了描述简谐振动的差分方程模型。
差分方程模型求解对差分方程模型求解通常有两种方法:一种是使用递推公式进行求解,另一个方法是使用其它数学方法,如拉普拉斯变换或离散傅立叶变换等。
第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 是一条上升的曲线。
差分方程模型与生长率模型下载提示:该文档是本店铺精心编制而成的,希望大家下载后,能够帮助大家解决实际问题。
文档下载后可定制修改,请根据实际需要进行调整和使用,谢谢!并且,本店铺为大家提供各种类型的实用资料,如教育随笔、日记赏析、句子摘抄、古诗大全、经典美文、话题作文、工作总结、词语解析、文案摘录、其他资料等等,想了解不同资料格式和写法,敬请关注!Download tips: This document is carefully compiled by this editor. I hope that after you download it, it can help you solve practical problems. The document can be customized and modified after downloading, please adjust and use it according to actual needs, thank you!In addition, this shop provides you with various types of practical materials, such as educational essays, diary appreciation, sentence excerpts, ancient poems, classic articles, topic composition, work summary, word parsing, copy excerpts,other materials and so on, want to know different data formats and writing methods, please pay attention!差分方程模型与生长率模型在数学和生物学领域中都扮演着重要的角色。
差分方程模型matlab差分方程模型在数学和工程领域中具有重要的应用。
它是描述动态系统行为的一种数学模型,通常由一系列离散时刻的状态变量和状态转移方程组成。
MATLAB作为一种功能强大的数值计算软件,为差分方程模型的建模和求解提供了便捷的工具和环境。
本文将介绍差分方程模型在MATLAB中的使用方法和应用场景。
首先,我们将探讨差分方程模型的基本原理和概念,然后详细介绍在MATLAB中的建模步骤和求解技巧。
最后,我们会给出一些在实际问题中使用差分方程模型的案例,并展示其在系统分析、控制和优化等方面的优势。
差分方程模型是描述离散系统行为的数学模型,常用于描述在给定时间步长下变量之间的关系。
它与连续时间的微分方程模型相对应,但在很多情况下,离散系统更符合实际情况。
差分方程模型可以描述许多系统,例如电路、金融市场、人口增长等。
在MATLAB中建立差分方程模型需要以下步骤:1. 定义变量:首先需要确定模型涉及的状态变量,然后在MATLAB 中声明这些变量。
可以使用向量或矩阵表示多个变量。
2. 构建状态转移方程:差分方程模型通过状态转移方程描述系统变量在不同时间步长之间的变化规律。
在MATLAB中,可以使用循环或矩阵运算构建状态转移方程。
3. 设定初值条件:差分方程模型通常需要给定初始条件,即在 t=0 时刻各个变量的值。
在MATLAB中,可以使用向量或矩阵存储初始条件。
4. 求解差分方程:在MATLAB中可以使用函数或求解器来求解差分方程模型。
常用的函数包括 `solve`、`ode45`、`ode15s`等,它们可以根据模型的具体特点选择合适的求解方法。
在实际应用中,差分方程模型在系统分析、控制和优化等方面具有广泛的应用。
例如,在系统分析中,可以通过建立差分方程模型来预测系统的行为和变化趋势。
在控制问题中,差分方程模型可以描述系统动态行为,从而设计和优化控制策略。
在优化问题中,差分方程模型可以作为约束条件或目标函数进行求解。
第八讲 差分方程模型一、差分方程介绍规定t 只取非负整数。
记为变量在t 点的取值,则称t y y t t t y y y −=Δ+1为的一阶向前差分,简称差分,称Δ为的二阶差分。
类似地,可以定义的阶差分。
t y t t t t t y t t y y y y y y +−=Δ−Δ=ΔΔ=+++12122)(t y t y n t ny Δ由及的差分给出的方程称为的差分方程,其中含的最高阶差分的阶数称为该差分方程的阶。
差分方程也可以写成不显含差分的形式。
例如,二阶差分方程也可改写成t y t 、t y t y t y 02=+Δ+Δt t t y y y 012=+−++t t t y y y 。
满足一差分方程的序列称为差分方程的解。
类似于微分方程情况,若解中含有的独立常数的个数等于差分方程的阶数时,称此解为该差分方程的通解。
若解中不含任意常数,则称此解为满足某些初值条件的特解。
t y 称如下形式的差分方程)(110t b y a y a y a t n t n t n =+++−++L (1) 为阶常系数线性差分方程,其中是常数,n n a a a ,,,10L 00≠a 。
其对应的齐次方程为0110=+++−++t n t n t n y a y a y a L (2)容易证明,若序列与均为(2)的解,则也是方程(2)的解,其中为任意常数。
若是方程(2)的解,是方程(1)的解,则也是方程(1)的解。
)1(t y )2(t y )2(2)1(1t tt y c y c y +=21,c c )1(t y )2(t y )2()1(t t t y y y +=方程(1)可用如下的代数方法求其通解: (I )先求解对应的特征方程(3)00110=+++−a a a n nL λλ(II )根据特征根的不同情况,求齐次方程(2)的通解。
(i )若特征方程(3)有n 个互不相同的实根n λλ,,1L ,则齐次方程(2)的通解为t n n t c c λλ++L 11 (为任意常数)n c c ,,1L (ii )若λ是特征方程(3)的重根,通解中对应于k λ的项为t k k tc c λ)(11−++L ,),,1(k i c i L =为任意常数。
数学建模方法之差分方程模型差分方程模型是数学建模中常用的一种方法,它基于差分方程来描述问题,并用差分方程来求解问题。
所谓差分方程,是指用差分代替微分的方程,它是一种离散的模型。
在实际问题中,很多情况下,并不能直接通过微分方程来描述问题,而差分方程模型则可以通过离散化的方法来近似地描述问题。
差分方程模型的优点之一是可以适用于离散化的数据,对于实际问题的离散化模型建立是非常有帮助的。
差分方程模型的另一个优点是可以通过数值方法来求解,不需要进行繁琐的解析推导,因此适用于复杂问题的求解。
差分方程模型的基本形式为:yn+1 = fn(yn, yn-1, ..., yn-k)其中,yn表示第n个时刻的解,fn是一个给定的函数,表示通过前k个时刻的解来计算第n+1个时刻的解。
这个方程是离散的,通过已知的初始条件来逐步递推获得结果。
差分方程模型的适用范围非常广泛,可以用于描述和预测各种动态过程。
例如,差分方程模型可以用来描述人口增长模型、生态系统模型、传染病模型等等。
在这些例子中,差分方程模型可以通过已知的数据和初始条件来预测未来的发展趋势。
差分方程模型的建立步骤主要包括以下几个方面:1.确定问题的描述和目标:明确问题的背景和目标,确定需要建立差分方程模型的原因和用途。
2.确定模型的变量和参数:根据实际问题,确定需要用到的变量和参数。
3.确定差分方程的形式和函数:根据问题的特点和要求,选择合适的差分方程形式和函数。
这部分需要结合实际问题和数学方法来确定。
4.确定初始条件和边界条件:确定差分方程模型的初始条件和边界条件。
这部分是求解差分方程的前提条件。
5.差分方程的求解和分析:通过数值方法求解差分方程,得到数值解,并对结果进行分析和解释。
差分方程模型数学建模讲座一、关于差分方程模型简单的例子1. 血流中地高辛的衰减地高辛用于心脏病。
考虑地高辛在血流中的衰减问题以开出能使地高辛保持在可接受(安全而有效)的水平上的剂量处方。
假定开了每日0.1毫克的剂量处方,且知道在每个剂量周期(每日)末还剩留一半地高辛,则可建立模型如下:设某病人第n 天后血流中地高辛剩余量为n a , 则1.05.01+=+n n a a (一阶非齐次线性差分方程) n n n n a a a a 5.01−=−=∆+2. 养老金问题对现有存款付给利息且允许每月有固定数额的提款, 直到提尽为止。
月利息为1℅,月提款额为1000元,则可建模型如下:设第n 月的存款额为n a ,则100001.11−=+n n a a (一阶非齐次线性差分方程)3. 兔子问题(Fibonacci 数)设第一月初有雌雄各一的一对小兔,假定两月后长成成兔,同时(即第三个月)开始,每月初产雌雄各一的一对小兔, 新增小兔也按此规律繁殖,设第n 月末共有n F 对兔子,则建模如下:==+=−−12121F F F F F n n n (二阶线性差分方程初值问题) 34232143212211F F F F F F F F F F ≠+=+ 注意上月新生的小兔不产兔(因第n 月末的兔子包括两部分, 一部分上月留下的为1−n F , 另一部分为当月新生的,而新生的小兔数=前月末的兔数)4.车出租问题A ,B 两地均为旅游城市,游客可在一个城市租车而在另一个城市还车。
A , B 两汽车公司需考虑置放足够的车辆满足用车需要,以便估算成本。
分析历史记录数据得出:n x : 第n 天营业结束时A 公司的车辆数 n y :第n 天营业结束时B 公司的车辆数则 +=+=++n n n nn n y x y y x x 7.04.03.06.011 (一阶线性差分方程组)(问题模型可进一步推广)二 差分方程的解法类比: 差分方程是数列间关系; 微分方程是函数间关系 定义1. 形如02211=++++−−−k n k n n n a b a b a b a L 的差分方程,称为}{n a 的k 阶常系数线性齐次差分方程,其中i b 为常数,0≠k b 且k n ≥. 02211=++++−−k k k k b x b x b x L 称为差分方程的特征方程,其根称为特征根。
定理1(单根情形)若特征方程恰有k 个相异的特征根k x x x ,,,21L ,则差分方程的通解为 nk k n n n x c x c x c a +++=L 2211.例1 求解兔子问题 ==+=−−12121F F F F F n n n解:差分方程的特征方程为 012=−−x x特征根 251,25121−=+=x x 通解为 nnn c c F−+ +=25125121, 由初始条件 121==F F 得: = −++=−++)2(1251251)1(1251251222121c c c c 解得 −==515121c c 故−− +=nn n F 25121551定理 2 (重根情形) 若特征方程的相异特征根为12,,,t x x x L , 重数依次为t m m m ,,,21L , 其中k m m m t =++L 21,则差分方程的通解为:121211111211212222()()m m n n n m m a c c n c n x c c n c n x −−=+++++L L112()t t m n t t tm t c c n c n x −++++L L(定理1包含在定理2之中)定理 3 若差分方程的特征方程的特征根出现一对共轭虚根,−=+=ivu x ivu x 21 和相异的2−k 个根k x x ,,3L , 则差分方程的通解为:nk k nnnn x c x c n c n c a ++++=L 3321sin cos θρθρ.定义2 形如)(2211n f a b a b a b a k n k n n n =++++−−−L(k b b b ,,,21L 为常数,0≠k b ,0)(≠n f ,k n ≥) 的差分方程为k 阶常系数线性非齐次差分方程。
称011=+++−−k n k n n a b a b a L为其对应的齐次方程。
定理4 非齐次差分方程的通解等于对应齐次差分方程的通解加上非齐次方程的特解,即n nn a a a +=*,其中*n a 为通解, n a 为特解.例2(地高辛问题)解 1.05.01+=+n n a a齐次特征方程 05.0=−λ, 齐次方程通解 n nc a )5.0(*=. 设特解为D a n =, 代入1.05.0+=D D 得0.2D =,于是所求通解为2.0)5.0(+=nnc a 例3 (养老金问题)解法1 100001.11−=+n n a a齐次特征方程001.1=−λ, 齐次方程通解nnc a )01.1(*=. 设特解为 D a n =, 代入代入原方程得 100000=D . 通解为 100000)01.1(+=n n c a .解法2(化齐) 100001.11−=−+n n a a ,100001.11−=−−n n a a .相减得001.101.211=+−−+n n n a a a , 特征方程001.101.22=+−λλ.0)1)(01.1(=−−λλ,通解为1(1.01)n n a c c =+,代入原方程得1000)01.1(01.101.11111−=−−⋅+++c c c c n n , 1100000c =.故 100000)01.1(+=nnc a . 例4 求非齐次差分方程nn n n a a a 24421=+−−−的通解. 解法1 齐次特征方程0442=+−λλ, 二重根2=λ,对应齐次方程的通解为n n n n c c a 2221*⋅+=. 因nn f 2)(=中, 2 是2 重根,故设特解为nn n A a 22⋅⋅=代入得12A =, 故通解为1221222−⋅+⋅+=n n n n n n c c a 方法2(化齐) nn n n a a a 24421=+−−−,132122)44(2−−−−⋅=+−n n n n a a a相减得12361280n n n n a a a a −−−−+−=,特征方程0812623=−+−λλλ 特征根2=λ为三重根,通解为nn n n n c n c c a 2222321⋅+⋅+=. 代入原方程得312c =,故1221222−⋅+⋅+=n n n n n n c c a三 差分方程的平衡点及稳定性1.一阶线性方程b aa a n n =+−1的平衡点及稳定性 平衡点可由b x a x =+解得abx +=10(相当于0,x a n n =∀的那种点).当初始条件00x a =, 则 0,x a n n =∀.若对任何初始条件,都有n →∞时, 0x a n →,则称平衡点0x 是稳定的,否则称为不稳定的。
一阶方程的通解aba c a nn ++−=1)(,因此1||<a 时才是稳定的. 地高辛问题1.05.01+=−n n a a ,通解2.0)5.0(+=nn c a , 0.2是平衡点,且是稳定的。
就是说,不管初始值如何,若干天以后,血中地高辛剩留量接近0.2.养老金问题 100001.11−=−n n a a ,通解100000)01.1(+=n n c a , 510是平衡点,不稳定.若01000000=⇒=c a , 则100000,=∀n a n若5010,0,a c >> 则n a →+∞若5010,0,a c << 则n a →−∞2. 二阶方程的平衡点及稳定性只须讨论齐次方程021=+−−−n n n ba aa a ; 对非齐次方程d ba aa a n n n =+−−−21(d 为常数)可作线性变换e a b n n −=化成齐次方程,稳定性相同.齐次方程通解nn n x c x c a 2211+=,平衡点为0, 21,x x 是互异特征根(或重根),当∞→n 仅当1||,1||21<<x x 才是稳定的。
3. n 阶齐次次方程组平衡点O 的稳定性.)(n a 为n 维列向量, A 为n n ×阵。
齐次线性差分方程组 0)1()(=−+n a n a ,平衡点O 稳定的条件是A 的所有特征根||1i λ<。
4.求解n 阶齐次线性差分方程组方法: 仿照线性微分方程组解的法,注意二者的区别 1t nn n n d xx x ce d ta a a c λλλλ− ==== .汽车出租问题 全类完似的问题:选民下一次选举的投票趋势110.60.30.40.7n n nn n n x x y y x y ++=+ =+n n X X=+7.04.03.06.01 0.60.30.40.7A=汽车出租问题解法120.60.3||(0.6)(0.7)0.12 1.30.300.40.7A E λλλλλλλ−−==−−−=−+=−3.0,121==λλ 特征根互异对11=λ, 设= =′b a b a X nn 1 , 由0)(1= −b a E A λ得 03.04.0=+−b a取 = 4.03.0b a ,=′4.03.0n X 对3.01=λ, 设=′′b a X nn 3.0 由0)(2=−b a E A λ得 03.03.0=+b a取 −= 11b a ,−=′′113.0n nX 于是−=213.04.03.03.0c c X n n n 通解为 −=+=n n n n c c y c c x 3.04.03.03.02121汽车出租问题解法2(理论解法)+=+=++n n n n n n y x y y x x 7.04.03.06.011 =7.04.03.06.0A ,B AT T = =−3.00011,1−=TBTA , 1110n n n X TBT X TB T X −−−== 令TBAT =. 即= 3.00017.04.03.06.021212121ββααββαα 即 = ++++2121221122113.03.07.04.07.04.03.06.03.06.0ββααβαβαβαβα,取−=14.013.0T 得通解 −= −=21213.04.03.03.03.000114.013.0c c c c X n n n n 平衡点是稳定的, 平衡点由+=+=y x y yx x 7.04.03.06.0 即y x3.04.0=A 、B 两地投放车辆比为3 : 4。