微分方程模型建模实例

合集下载

微分方程(组)模型

③
(2) 方程③是一阶线性微分方程，通解为②当n>0时，有特解y=0.
求微分方程（组）的解析解命令: dsolve(‘方程1’, ‘方程2’,…‘方程n’, ‘初始条件’, ‘自变量’) 符号说明：在表达微分方程时，用字母D表示求微分， D2、D3等表示求2阶、3阶等微分。任何D后所跟的字母为因变量，自变量可以指定或由系统规则选定为确省。 d2y
方法：
• 规律分析法：根据相关学科的定理或定律、规律（这些涉及到某些函数变化率）建立微分方程模型，如曲线的切线性质. • 微元分析法:应用一些已知规律和定律寻求微元之间的关系式. • 近似模拟法：在社会科学、生物学、医学、经济学等学科的实际问题中，许多现象的规律性不清楚，常常用近似模拟的方法建立微分方程模型.
4.符号说明
• • • • • • • a---某人每天在食物中摄取的热量 b---某人每天用于新陈代谢（及自动消耗）的热量 α ---某人每天从事工作、生活每千克体重必需消耗的热量 β---某人每天从事体育锻炼每千克体重消耗的热量 w---体重(单位：千克) w0---体重的初始值 t---时间（单位：天）
若Q(x)≡0,则称为一阶线性齐次方程，一阶线性微分方程通解为 P ( x ) dx P ( x ) dx ② y ( x) e ( Q( x)e dx C )
从而可得
dz (1 n) P ( x) z (1 n)Q ( x) dx
dz dy (1 n) y n dx dx
一、微分方程模型二、微分方程的数学形式三、微分方程（组）的MATLAB解法四、减肥的数学模型五、人口增长数学模型六、兰彻斯特（Lanchester）作战模型七、硫磺岛战役案例

微分方程模型的建立与求解

微分方程模型的建立与求解微分方程是自然界中许多现象的数学描述，通过建立微分方程模型可以更好地理解和预测各种现象。

本文将介绍微分方程模型的建立与求解方法。

一、微分方程模型的建立微分方程通常用来描述系统内部的变化规律，要建立微分方程模型，首先需要根据具体问题分析系统的特点，确定影响系统变化的因素，并建立相关的数学表达式。

以一个简单的弹簧振子系统为例，假设弹簧的位移为x(t)，弹簧的弹性系数为k，质量为m，外力为f(t)，则可以建立微分方程模型：$$ m\\frac{{d^2x}}{{dt^2}} + kx = f(t) $$二、微分方程模型的求解1. 解析解法对于一些简单的微分方程，可以通过解析的方法求解。

例如，对于一阶线性微分方程：$$ \\frac{{dy}}{{dx}} + P(x)y = Q(x) $$可以通过积分因子的方法求解。

2. 数值解法对于复杂的微分方程或无法求得解析解的情况，可以借助数值方法进行求解。

常用的数值解法包括欧拉方法、龙格-库塔法等，通过逐步迭代逼近真实解。

3. 计算机模拟借助计算机编程，可以通过数值方法对微分方程进行求解，这在实际工程和科学研究中非常常见。

利用计算机程序，可以模拟出系统的运行状态，观察系统的响应特性。

三、实例分析以简单的振动系统为例，通过建立微分方程模型并利用数值方法进行求解，可以分析系统的振动特性。

通过调节参数值，可以观察到系统振动的变化规律，为系统设计和控制提供重要参考。

结论微分方程模型的建立与求解是数学建模中的重要一环，通过适当的模型建立和求解方法，可以更好地了解和预测系统的行为。

在实际应用中，需要综合运用解析方法、数值方法和计算机模拟，以全面分析和解决问题。

以上是关于微分方程模型的建立与求解的介绍，希望对读者有所帮助。

微分方程方法建模概述及举例

微分方程方法建模概述及举例微分方程是数学中的一个重要分支，广泛应用于各个领域，特别是自然科学和工程学科中的建模问题。

本文将概述微分方程方法建模的基本思路，并通过举例说明其在实际问题中的应用。

1.问题抽象化：首先需要将实际问题抽象成一个或一组微分方程。

通过观察问题的物理过程和规律，了解问题中的变量、因果关系以及其演化过程。

将这些信息用数学语言表示出来，通常是通过建立数学模型来描述问题。

2.建立微分方程：基于问题的抽象化模型，我们可以建立相应的微分方程。

根据物理规律和描述问题演化的数学关系，确定方程中的变量、常数和系数。

对于复杂问题，可能需要引入附加的假设和近似，以简化问题求解。

3.求解微分方程：通过求解微分方程，可以得到问题的数学解。

求解方法包括解析解和数值解两种。

解析解通常是通过变量分离、常数变易、积分变换等方法，求得方程的具体解析形式。

数值解则是通过数值计算方法，如欧拉法、龙格-库塔法等，近似计算出微分方程的解。

4.模型验证和分析：将求得的数学解与实际问题进行比较和分析，验证模型的有效性和准确性。

通过对模型进行敏感性分析和参数优化，对模型进行改进和完善。

现在我们来通过两个实际问题的建模例子，进一步说明微分方程方法的应用。

1.指数增长模型问题：假设一个生物种群遵循指数增长规律，种群数量在一段时间内以固定比率增加。

已知在初始时刻，种群数量为100只，经过3个小时后，种群数量增加到了1000只。

求解该问题。

解答：我们可以建立如下的微分方程模型:dy/dt = k * y其中，y表示种群数量，t表示时间，k为增长率。

根据已知条件，当t=0时，y=100；当t=3时，y=1000。

将这些条件代入微分方程，就可以求解得到k的值。

然后再根据k的值，求解出种群数量y随时间t的变化。

2.弹簧振动模型问题：一个弹簧系统在无外力作用下，其振动满足以下微分方程：m* d^2y/dt^2 = -k * y，其中m为弹簧的质量，k为弹簧的劲度系数。

微分方程模型-伪造名画案

知识回顾 Knowledge Review
祝您成功！
设 t 时刻的原子数为N (t) ，则有
dN N
dt
为物质的衰变常数。
初始条件
N t t0
N0
N (t)
N e (tt0 ) 0
t
t0
1
ln
N0 N
t
t0
1
ln
N0 N
半衰期 T 1 ln 2
碳-14 T 5568 年
镭-226
T 1600 年
铀-238 T 45亿年铅-210 T 22年
微分方程建模案例—— 范. 梅格伦伪造名画案
第二次世界大战比利时解放后，荷兰保安机关开始搜捕纳粹分子的合作者，发现一名三流画家H.A.Vanmeegren 曾将17世纪荷兰著名画家Jan.Vermeer的一批名贵油画盗卖给德寇，于1945年5月29日通敌罪逮捕了此人。
Vanmeegren被捕后宣称他从未出卖过荷兰的利益，所有的油画都是自己伪造的，为了证实这一切，在狱中开始伪造Vermeer的画《耶稣在学者中间》。当他的工作快完成时，又获悉他可能以伪造罪被判刑，于是拒绝将画老化，以免留下罪证。
但是，许多人还是不相信其余的名画是伪造的，因为， Vanmeegren在狱中作的画实在是质量太差，所找理由都不能使怀疑者满意。直到20年后，1967年，卡内基梅隆大学的科学家们用微分方程模型解决了这一问题。
原理著名物理学家卢瑟夫（Rutherford）指出：
物质的放射性正比于现存物质的原子数。
（2）钋的半衰期为138天容易测定，铅210的半衰期为22年，对要鉴别的300多年的颜料来说，每克白铅中每分钟钋的衰变数与铅210的衰变数可视为相等。

微分方程建模案例

微分方程建模案例微分方程是数学中的一种重要工具，它被广泛应用于各个领域的建模和问题求解中。

下面将以一个具体的案例来介绍微分方程建模的过程，并通过求解微分方程来解决实际问题。

案例：生物种群的增长模型在生态学中，研究生物种群的增长是一个重要的课题。

种群的增长速度与种群中的个体数量有关。

如果种群中个体数量增加的速度与当前个体数量成正比，可以建立如下的微分方程模型：$$\frac{dN}{dt} = rN$$其中，$N$表示种群的个体数量，$t$表示时间，$r$表示增长的速率。

这个微分方程描述了种群个体数量随时间变化的规律。

解：首先，我们需要求解上述微分方程，得到种群个体数量随时间的函数关系。

这是一个一阶线性常微分方程，我们可以使用分离变量的方法求解。

将微分方程变形为：$$\frac{dN}{N} = rdt$$将方程两边同时积分，得到：$$\int \frac{dN}{N} = \int rdt$$经过积分运算，得到：$$\ln N = rt + C$$其中，$C$为积分常数。

进一步求解，得到：$$N = e^{rt + C}$$根据初始条件，当$t=0$时，$N=N_0$，其中$N_0$为初始种群个体数量。

代入初始条件，解得$C=\ln N_0$，将其代入上述方程，得到最终的解：$$N = N_0e^{rt}$$这个解描述了种群个体数量随时间的增长情况。

接下来，我们来解决一个具体的问题，一个兔子种群的增长情况。

假设初始时刻兔子种群中有100只兔子，增长速率$r=0.02$，那么该种群在未来的10个月内，兔子的数量会如何变化？根据上面的微分方程解，代入初始条件$N_0=100$，$r=0.02$，$t=10$，得到：$$N=100e^{0.02t}$$将$t=10$代入上述方程，可以得到10个月后兔子种群的个体数量：所以，10个月后的兔子种群中大约有122只兔子。

通过这个模型，我们可以预测种群在未来的增长情况，并在实践中应用于生态学、环境保护等领域，为实际问题的决策提供参考。

微分方程预测模型实例

微分方程预测模型实例引言微分方程是数学中的重要概念，用于描述自然界中的各种变化和现象。

它在物理学、工程学、经济学等领域都有广泛应用。

在本文中，我们将介绍微分方程预测模型的概念和实例，以帮助读者更好地理解和应用这一方法。

什么是微分方程预测模型？微分方程预测模型是一种利用已知条件和规律，通过建立微分方程来预测未来变化的方法。

它基于数学原理和统计学方法，通过对已有数据进行拟合和分析，得出一个能够描述系统行为的微分方程，并利用该方程进行未来的预测。

微分方程预测模型的应用微分方程预测模型广泛应用于各个领域，下面我们以经典案例为例介绍其中两个：1. 成长模型成长模型是一类常见的微分方程预测模型。

它通常用于描述人口、生物群体等在时间上的增长情况。

以人口增长为例，我们可以假设人口增长率与当前人口数量成正比，即：dPdt=kP其中，P表示人口数量，k为比例常数。

这是一个一阶线性常微分方程，可以通过求解得到人口数量随时间的变化情况。

通过拟合已有的人口数据，我们可以得到合适的k值，并利用该方程进行未来人口数量的预测。

2. 热传导模型热传导模型是另一个常见的微分方程预测模型。

它通常用于描述物体内部温度随时间和空间的变化情况。

以一维热传导为例，我们可以假设物体内部温度变化率与温度梯度成正比，即：∂T ∂t =α∂2T∂x2其中，T表示温度，α为热扩散系数。

这是一个二阶偏微分方程，可以通过求解得到物体内部温度随时间和空间的变化情况。

通过拟合已有的温度数据和边界条件，我们可以得到合适的α值，并利用该方程进行未来温度分布的预测。

微分方程预测模型实例下面我们以一维热传导模型为例，介绍微分方程预测模型的具体实现步骤。

步骤一：收集数据首先，我们需要收集已有的温度数据。

假设我们有一个金属棒，长度为L，初始时刻t=0时，金属棒上各点的温度分布已知。

步骤二：建立微分方程根据热传导模型的假设，我们可以建立如下的一维热传导方程：∂T ∂t =α∂2T∂x2其中，T(x,t)表示金属棒上某点处的温度，α为热扩散系数。

3.1微分方程模型-微分方程的几个简单实例

微分方程模型浙江大学数学建模实践基地§3.1 微分方程的几个简单实例在许多实际问题中，当直接导出变量之间的函数关系较为困难，但导出包含未知函数的导数或微分的关系式较为容易时，可用建立微分方程模型的方法来研究该问题，本节将通过一些最简单的实例来说明微分方程建模的一般方法。

在连续变量问题的研究中，微分方程是十分常用的数学工具之一。

例1（理想单摆运动）建立理想单摆运动满足的微分方程，并得出理想单摆运动的周期公式。

从图3-1中不难看出，小球所受的合力为mgsin θ，根据牛顿第二定律可得：sin ml mg θθ=- 从而得出两阶微分方程：0sin 0(0)0,(0)g l θθθθθ+==⎪=⎧⎪⎨⎩ （3.1）这是理想单摆应满足的运动方程（3.1）是一个两阶非线性方程，不易求解。

当θ很小时，sin θ≈θ，此时，可考察（3.1）的近似线性方程：00(0)0,(0)g l θθθθθ+==⎧=⎪⎨⎪⎩ （3.2）由此即可得出2g T l π=（3.2）的解为: θ(t )=θ0cosωtg l ω=其中当时,θ(t )=04T t =42g T l π=故有M Q P mgθl 图3-1（3.1）的近似方程例2我方巡逻艇发现敌方潜水艇。

与此同时敌方潜水艇也发现了我方巡逻艇，并迅速下潜逃逸。

设两艇间距离为60哩，潜水艇最大航速为30节而巡逻艇最大航速为60节，问巡逻艇应如何追赶潜水艇。

这一问题属于对策问题，较为复杂。

讨论以下简单情形：敌潜艇发现自己目标已暴露后，立即下潜，并沿着直线方向全速逃逸，逃逸方向我方不知。

设巡逻艇在A 处发现位于B 处的潜水艇，取极坐标，以B 为极点，BA 为极轴，设巡逻艇追赶路径在此极坐标下的方程为r =r (θ)，见图3-2。

B AA1dr ds dθθ图3-2由题意，，故ds =2dr 2ds dr dt dt =图3-2可看出，222()()()ds dr rd θ=+故有:2223()()dr r d θ=即:3rdr d θ=（3.3）解为：3r Ae θ=（3.4）先使自己到极点的距离等于潜艇到极点的距离然后按（3.4）对数螺线航行，即可追上潜艇。

数学建模-微分方程模型-饮酒驾车问题

和 x0 ，将体重 70kg 的某人在快速喝下 2 瓶啤酒之后一段时间内他血液中酒精含量的
测量值进行处理后，得到附录 1 所示的 y0 0 时的一组数据，并采用非线性最小二乘法拟合算法对系数进行求解，得出参数如下。 x0 5193
=2.00796
=0.1855
同时可以看到，每瓶啤酒含酒精量为 2596.5mg。所以，得出的血液中酒精含量关于时间的函数如下。
0.1855 t e 2.00756t ) 2860.78604(e y (t ) 0.1855( t 6) 2860.8028e 2.00756(t 6) 3800.7595e
0t 6 6 t 12
利用 matlab 对以上模型进行求解。图 3 大李血液中酒精含量随时间变化图像
y (t ) ( y0 +5721.57208)e 0.1855t 5721.57208e 2.00796t
拟合效果如图。图 1 函数的拟合效果
图 2 残差分析图
残差分析图
600 500 400 300 200 100 0 10 11 12 13 14 15 0.5 1.5 2.5 3.5 0.25 ‐100 ‐200 ‐300 ‐400 残差 0.75 4.5 16 1 2 3 4 5 6 7 8 9
时刻为 t 时胃肠道中的酒精含量。
y (t ) 时刻为 t 时血液中的酒精含量。

胃肠道中的酒精进入血液的转移率与胃肠道中酒精量的比值。血液中的酒精的排除率与血液中酒精量的比值。
五、模型的建立与求解
5.1 问题一根据题目叙述，大李的实际情况符合快速饮酒的模型。为了确定函数中的系数，

微分方程模型案例分析

微分方程模型案例分析-------传染病传播的数学模型张清华由于人体的疾病难以控制和变化莫测，因此医学中的数学模型较为复杂。

医学中的数学模型分为两大类：传染病传播的数学模型和疾病数学模型。

以下仅讨论传染病的传播问题。

人们将传染病的统计数据进行处理和分析，发现在某一民族或地区，某种传染病传播时，每次所涉及的人数大体上是一常数。

这一现象如何解释呢？关于这个问题，医学工作者试图从医学的不同角度进行解释都得不到令人满意的解释。

最后由于数学工作者的参与，在理论上对上述结论进行了严格的证明。

同时又由于传染病数学模型的建立，分析所得结果与实际过程比较吻合，这个现象才得到了比较满意的解释。

传染病传播所涉及的因素很多，如传染病人的多少，易受传染者的多少，传染率的大小，排除率的大小，人口的出生和死亡等。

如果还要考虑人员的迁入与迁出，潜伏期的长短以及预防疾病的传播等因素的影响，那么传染病的传播就变得非常复杂。

如果一开始就把所有的因素考虑在内，那么将陷入多如乱麻的头绪中不能自拔，倒不如舍去众多的次要因素，抓住主要因素，把问题简化，建立相应的数学模型。

将所得结果与实际比较，找出问题，修改原有假设，再建立一个与实际比较吻合的模型。

下面由简单到复杂将建模的思考过程作一个示范，读者可以从中得到很好的启发。

1 模型一假设(1)，每个病人在单位时间内传染的人数是常数K 0；假设(2)，一人得病后，经久不愈，并在传染期内不会死亡。

记i t ()表示t 时刻病人数，K 0表示每个病人单位时间内传染的人数，i i ()00=，即最初有i 0个传染病人。

则在∆t 时间内增加的病人数为i t t i t K i t t ()()()+-=∆∆0于是得微分方程⎪⎩⎪⎨⎧==00)0()()(i i t i K dt t di (1)，其解为 i t i e k t ()=00结果表明：传染病的传播是按指数函数增加的。

这个结果与传染病传播初期比较吻合，传染病传播初期，传播快，被传染人数按指数函数增长。

微分方程建模个例

A1
C
C1
分析：1.追击开始后，大家将进入正方 A 形里面，距离将变小，由于追击的规则及四个人速度和方向的假定，四人还是在某个正方形的顶点上。 2.会不会出现四个人绕一个圆循环追？不会！距离会不断缩小最后到一点，就是正方形的中心。追击曲线是四条指向 D1 中心的螺旋线（可能绕中心几周） 3.坐标架怎么建？ D O点在中心，直角坐标架。
2H g
2.二氧化碳的吸收
空气通过盛有CO2的吸收剂的圆柱形器皿，已知它吸收CO2的量与 CO2的浓度及吸收层的厚度成正比，今有含CO28%的空气通过厚度为10cm的吸收层后浓度为2%，求：（1）若吸收层变为30cm厚，出口浓度是多少？（2）要使出口浓度为1%，应该设多厚的吸收层？解：记吸收层厚度为d，等分n份，每小层d/n厘米。入口浓度为8%，在每小层看吸收量，第一层后被吸收量为： kd k8%d/n，含量变为： 8%(1)
v0t y x(0) 0 y , 就是曲线的切向量， 1 x y (0) 0
Q(1,v0t) 模型里y(t)，x(t)都是t的函数，但是三个变量不好处理，注意我们要求的是y(x)。 P(x,y) O 1 x
(1 x) y y v0t实现了变量t的分离
再建立一个y(t)，x(t)，t的关系：t时间里导弹已飞行的距离是可求的。 x 1 y2 dx 5v0t (1 x) y y v0t , x0 0, y0 0
v r (0) 2 2 , (2r cos dx cos dr r sin d dx r sin cos d , , y r sin dy sin dr r cos d dy r cos sin dr d 1 sin cos dx dr r r cos r sin dy

4-微分方程建模实例——Malthus模型与Logistic模型-课件PPT

23
于是，
N0 N (t)e (tt0 ) r[e (tt0 ) 1].
若此画是真品，t - t0 ≈ 300 (年) . 从而可求出 λN0 的近似值. 对油画《在埃牟斯的门徒》具体计算如下：
N0 N (t)e300 r[e300 1]
由于半衰期: T ln 2 ,
于是， ln 2 .
4.1. 人口增长模型 4.2. 赝品的鉴定 4.3. 耐用新产品的销售速度问题 4.4. 传染病模型
1
4.1 人口增长模型
世界人口增长概况
年
1625 1830 1930 1960 1974 1987 1999
人口(亿) 5 10 20 30 40 50 60
中国人口增长概况
年 1908 1933 1953 1964 1982 1990 1995 2000 人口(亿) 3.0 4.7 6.0 7.2 10.3 11.3 12.0 13.0
马尔萨斯(1766～1834) Malthus，Thomas Robert
4
模型假设： • 人口增长率 r 是常数. • 人口的数量本应取离散
值，但由于人口数量一般较大，为建立微分方程模型，可以将人口数量看作连续变量，甚至允许它为可微变量，由此引起的误差将是十分微小的.
5
模型构成：
设 x(t) 表示 t 时刻的人口，有
16
• 六十年后，美国记者、专栏作家乔纳森·洛佩兹(Jonathan Lopez)出版了《制造维米尔的人》(The man who made Vermeers) 一书. 在书中，洛佩兹表达了对那个时代荷兰人民的体谅：“荷兰人对米格伦的态度并非不可理解. 在二战中，这个国家遭遇了残酷的羞辱，光复也是在盟国的帮助下完成. 米格伦给了未能主宰自身命运的荷兰人内心深处想要得到的东西. 而对于 ‘欺骗’这种事情，他又是太熟谙了.”

微分方程型建模实例题

一个数学问题都可以用不同的方法来求解的，不同的方法做出来效果不同，效率也不同。

下面就微分方程模型建模展开建模。

下面给出些微分方程建立模型的实例，供大家参考。

1．一个半球状雪堆，其体积融化的速率与半球面面积S成正比，比例系数k > 0。

设融化中雪堆始终保持半球状，初始半径为R且3小时中融化了总体积的7/8，问雪堆全部融化还需要多长时间？2．从致冰厂购买了一块立方体的冰块，在运输途中发现，第一小时大约融化了1/4 （1）求冰块全部融化要多长时间（设气温不变）（2）如运输时间需要2.5小时，问：运输途中冰块大约会融化掉多少？3．一展开角为α的圆锥形漏斗内盛着高度为H的水，设漏斗底部的孔足够大（表面张力不计），试求漏斗中的水流光需要多少时间？4．容器甲的温度为60度，将其内的温度计移入容器乙内，设十分钟后温度计读数为70度，又过十分钟后温度计读数为76度，试求容器乙内的温度。

5．一块加过热的金属块初始时比室温高70度，20分钟测得它比室温高60度，问：（1）2小时后金属块比室温高多少？（2）多少时间后，金属块比室温高10度？6．设初始时容器里盛放着含净盐10千克的盐水100升，现对其以每分钟3升的速率注入清水，容器内装有搅拌器能将溶液迅时搅拌均匀，并同时以每分钟2升的速率放出盐水，求1小时后容器里的盐水中还含有多少净盐？7．某伞降兵跳伞时的总质量为100公斤（含武器装备），降落伞张开前的空气阻力为0.5v,该伞降兵的初始下落速度为0，经8秒钟后降落伞打开，降落伞打开后的空气阻力约为0.6 试球给伞降兵下落的速度v(t),并求其下落的极限速度。

8．1988年8月5日英国人Mike McCarthy创建了一项最低开伞的跳伞纪录，它从比萨斜塔上跳下，到离地179英尺时才打开降落伞，试求他落地时的速度。

9．证明对数螺线r=A 上任一处的切线与极径的夹角的正切为一常数，（）10．实验证明，当速度远低于音速时，空气阻力正比与速度，阻力系数大约为0.005。

常微分方程数学建模案例分析

常微分方程数学建模案例分析常微分方程是运用微积分中的概念与理论研究变化率的方程。

它是数学建模中常用的方法之一，可用于描述各种实际问题，如经济增长、生物扩散、化学反应等。

本文将通过一个关于人群传染病的数学建模案例，分析常微分方程在实际问题中的应用。

假设地有一种传染病，病毒的传播速度与感染者的接触频率有关。

现在我们要研究传染病的传播速度以及控制措施对传染病传播的影响。

为此，我们可以建立如下的数学模型：设N(t)表示时间t时刻的总人口数，而I(t)表示感染者的人口数，S(t)表示易感者的人口数。

根据该模型，易感者的人数随时间的变化率可表示为：dS/dt = -βSI其中，β表示感染率，即感染者每接触到一个易感者，会使其发病的概率。

感染者的人数随时间的变化率可表示为：dI/dt = βSI - γI其中，γ表示恢复率，即感染者每天被治愈的人数。

总人口数随时间的变化率可以通过易感者和感染者的变化率求和得到：dN/dt = dS/dt + dI/dt通过对该方程进行求解，我们可以得到感染者和易感者的人数随时间变化的解析解。

进一步，我们可以通过调节β和γ来研究不同的传播速度和控制措施对传染病传播的影响。

例如，如果β较大，表示感染率较高，此时传染速度会加快，可能导致传染病扩散的速度加快。

反之，如果β较小，表示感染率较低，传染病传播的速度会减慢。

另外，如果γ较大，表示恢复率较高，此时感染者的人数会快速减少，传染病传播的速度会减慢。

相反，如果γ较小，传染病传播的速度会加快。

通过对这些参数的调节，我们可以研究不同的控制措施对传染病传播的影响。

例如，我们可以通过降低感染率β或增加恢复率γ来减缓传染病传播的速度，从而控制疫情的爆发。

在实际应用中，常微分方程数学建模方法可以用于预测传染病的传播趋势，评估各种干预措施的效果。

此外，还可以通过引入更多的变量和参数，建立更复杂的模型，以更好地解释实际问题。

总之，常微分方程是数学建模中常用的方法之一，可以用于描述各种实际问题，如传染病的传播、经济增长等。

数学建模,第三章-微分方程模型

8小时20分－2小时57分＝5小时23分
即死亡时间大约在下午5:23，因此张某不能被排除在嫌疑犯之外。
理学院
3.2 目标跟踪模型
例1 饿狼追兔问题黑龙现有一直兔子，一只狼，兔子位于狼的正西100米处，假江科设兔子与狼同时发现对方并一起起跑，兔子往正北60米处的技巢穴跑，而狼在追兔子，已知兔子、狼是匀速跑且狼的速度学是兔子的2倍。兔子能否安全回到巢穴? 整理得到下述模型：院解：设狼的行走轨迹为y=f(x),则有：
理பைடு நூலகம்院
本章将通过一些最简单的实例来说明微分方程建模的一般方法。在连续变量问题的研究中，微分方程是十分常用的数学工具之一。
在许多实际问题中，当直接导出变量之间的函数关系较为困难，但导出包含未知函数的导数或微分的关系式较为容易时，可用建立微分方程模型的方法来研究该问题，
黑龙江科技学院数学建模
数学建模
B
60
2 2xf' ' x 1 f' x y' x 0 , y 0 100 x 100 解得狼的行走轨迹为： 100 0 100 (0,h) 0, f' f 假设在某一时刻，兔子跑到处，而狼在 (x,y)处，则有：
理学院
y y0 g e
g
车间空气中CO2浓度y 与时间t的数学模型
黑龙江科技学院数学建模
3.4 学习模型
一般认为，对一项技术工作，开始学得较快，但随着学得越来越多时，内容也越来越复杂，学员学得就会越来越慢。
员学习的速度，则随y的增长而下降。
dy 设y%表示已经掌握了这项工作的百分数， dt

常微分方程在数学建模中的应用(免费版)

常微分方程在数学建模中的应用这里介绍几个典型的用微分方程建立数学模型的例子. 一、人口预测模型由于资源的有限性,当今世界各国都注意有计划地控制人口的增长,为了得到人口预测模型,必须首先搞清影响人口增长的因素,而影响人口增长的因素很多,如人口的自然出生率、人口的自然死亡率、人口的迁移、自然灾害、战争等诸多因素,如果一开始就把所有因素都考虑进去，则无从下手.因此,先把问题简化,建立比较粗糙的模型,再逐步修改,得到较完善的模型.例1（马尔萨斯（Malthus ）模型）英国人口统计学家马尔萨斯（1766—1834）在担任牧师期间,查看了教堂100多年人口出生统计资料,发现人口出生率是一个常数,于1789年在《人口原理》一书中提出了闻名于世的马尔萨斯人口模型,他的基本假设是：在人口自然增长过程中,净相对增长（出生率与死亡率之差）是常数,即单位时间内人口的增长量与人口成正比,比例系数设为r ,在此假设下,推导并求解人口随时间变化的数学模型.解设时刻t 的人口为)(t N ,把)(t N 当作连续、可微函数处理（因人口总数很大,可近似地这样处理,此乃离散变量连续化处理）,据马尔萨斯的假设,在t 到t t ∆+时间段内,人口的增长量为t t rN t N t t N ∆=-∆+)()()(,并设0t t =时刻的人口为0N ,于是|⎪⎩⎪⎨⎧==．，00)(d d N t N rN t N这就是马尔萨斯人口模型,用分离变量法易求出其解为)(00e )(t t r N t N -=,此式表明人口以指数规律随时间无限增长.模型检验：据估计1961年地球上的人口总数为91006.3⨯,而在以后7年中,人口总数以每年2%的速度增长,这样19610=t ,901006.3⨯=N ,02.0=r ,于是)1961(02.09e1006.3)(-⨯=t t N .这个公式非常准确地反映了在1700—1961年间世界人口总数.因为,这期间地球上的人口大约每35年翻一番,而上式断定年增加一倍（请读者证明这一点）．但是,后来人们以美国人口为例,用马尔萨斯模型计算结果与人口资料比较,却发现有很大的差异,尤其是在用此模型预测较遥远的未来地球人口总数时,发现更令人不可思议的问题,如按此模型计算,到2670年,地球上将有36 000亿人口.如果地球表面全是陆地（事实上,地球表面还有80%被水覆盖）,我们也只得互相踩着肩膀站成两层了,这是非常荒谬的,因此,这一模型应该修改.;例2（逻辑Logistic 模型）马尔萨斯模型为什么不能预测未来的人口呢这主要是地球上的各种资源只能供一定数量的人生活,随着人口的增加,自然资源环境条件等因素对人口增长的限制作用越来越显著,如果当人口较少时,人口的自然增长率可以看作常数的话,那么当人口增加到一定数量以后,这个增长率就要随人口的增加而减小.因此,应对马尔萨斯模型中关于净增长率为常数的假设进行修改.1838年,荷兰生物数学家韦尔侯斯特(Verhulst)引入常数m N ,用来表示自然环境条件所能容许的最大人口数（一般说来,一个国家工业化程度越高,它的生活空间就越大,食物就越多,从而m N 就越大）,并假设将增长率等于⎪⎪⎭⎫⎝⎛-m N t N r )(1,即净增长率随着)(t N 的增加而减小,当m N t N →)(时,净增长率趋于零,按此假定建立人口预测模型.解由韦尔侯斯特假定,马尔萨斯模型应改为⎪⎩⎪⎨⎧=⎪⎪⎭⎫⎝⎛-=，，000)(1d d N t N N N N r t N 上式就是逻辑模型,该方程可分离变量,其解为,)(00e 11)(t t r m mN N N t N --⎪⎪⎭⎫ ⎝⎛-+=.下面,我们对模型作一简要分析.（1）当∞→t ,m N t N →)(,即无论人口的初值如何,人口总数趋向于极限值m N ；@（2）当m N N <<0时,01d d >⎪⎪⎭⎫ ⎝⎛-=N N N r t N m ,这说明)(t N 是时间t 的单调递增函数；（3）由于N N N N N r t N m m ⎪⎪⎭⎫ ⎝⎛-⎪⎪⎭⎫ ⎝⎛-=211d d 222,所以当2m N N <时,0d d 22>t N ,t N d d 单增；当2m N N >时,0d d 22<tN ,t N d d 单减,即人口增长率t Nd d 由增变减,在2m N 处最大,也就是说在人口总数达到极限值一半以前是加速生长期,过这一点后,生长的速率逐渐变小,并且迟早会达到零,这是减速生长期；（4）用该模型检验美国从1790年到1950年的人口,发现模型计算的结果与实际人口在1930年以前都非常吻合,自从1930年以后,误差愈来愈大,一个明显的原因是在20世纪60年代美国的实际人口数已经突破了20世纪初所设的极限人口.由此可见该模型的缺点之一是m N 不易确定,事实上,随着一个国家经济的腾飞,它所拥有的食物就越丰富, m N 的值也就越大；(5）用逻辑模型来预测世界未来人口总数.某生物学家估计,029.0=r ,又当人口总数为91006.3⨯时,人口每年以2%的速率增长,由逻辑模型得⎪⎪⎭⎫⎝⎛-=m N N r t N N 1d d 1, 即 ⎪⎪⎭⎫ ⎝⎛⨯-=m N 91006.31029.002.0, 从而得 91086.9⨯=m N ,即世界人口总数极限值近100亿. ）值得说明的是：人也是一种生物,因此,上面关于人口模型的讨论,原则上也可以用于在自然环境下单一物种生存着的其他生物,如森林中的树木、池塘中的鱼等,逻辑模型有着广泛的应用.二、市场价格模型对于纯粹的市场经济来说,商品市场价格取决于市场供需之间的关系,市场价格能促使商品的供给与需求相等(这样的价格称为(静态)均衡价格).也就是说,如果不考虑商品价格形成的动态过程,那么商品的市场价格应能保证市场的供需平衡,但是,实际的市场价格不会恰好等于均衡价格,而且价格也不会是静态的,应是随时间不断变化的动态过程.例3 试建立描述市场价格形成的动态过程的数学模型解假设在某一时刻t ,商品的价格为)(t p ,它与该商品的均衡价格间有差别,此时,存在供需差,此供需差促使价格变动.对新的价格,又有新的供需差,如此不断调节,就构成市场价格形成的动态过程,假设价格)(t p 的变化率tpd d 与需求和供给之差成正比,并记),(r p f 为需求函数,)(p g 为供给函数（r 为参数）,于是()()[]⎪⎩⎪⎨⎧=-=，，0)0(,d d p p p g r p f tpα 其中0p 为商品在0=t 时刻的价格,α为正常数.若设b ap r p f +-=),(,d cp p g +=)(,则上式变为—⎪⎩⎪⎨⎧=-++-=，，0)0()()(d d p p d b p c a t pαα ① 其中d c b a ,,,均为正常数,其解为ca db c a d b p t p t c a +-+⎪⎭⎫ ⎝⎛+--=+-)(0e)(α. 下面对所得结果进行讨论:（1）设p 为静态均衡价格 ,则其应满足0)(),(=-p g r p f ,即d p c b p a +=+-,于是得ca db p +-=,从而价格函数)(t p 可写为。

数学建模实例

我缉私舰雷达发现，在其正西方距c海里处
有一艘走私船正以匀速度a沿直线向北行驶，缉
私舰立即以最大的速度b追赶，若用雷达进行跟
踪，保持船的瞬时速度方向始终指向走私船，试求缉私舰追逐路线和追上的时间。
图2 走私船与缉私舰的位置关系
走私船
R(0,at)
缉私艇 D(x,y)
O
(c,0)
x
几何关系
dy y at tg dx x dy 即 x y at dx
Z取不同值时的浓度C(30)和时间T
Z / m3 C (30) / m3
T / min
5 0.00239 552
10 0.00478 738
15 0.00717 918
20 0.00956 1014
案例5 药物在体内的变化（房室模型）
何为房室系统？
在用微分方程研究实际问题时，人们常常采用一种叫‚房室系统‛的观点来考察问题。根据研究对象的特征或研究的不同精度要求，我们把研究对象看成一个整体（单房室系统）或将其剖分成若干个相互存在着某种联系的部分（多房室系统）。
二、微分方程模型
微分方程的建模步骤
1、翻译或转化：在实际问题中许多表示导数的常用词，如 ‚速率‛、‘增长‛(在生物学以及人口问题研究中)， ‚衰变‛(在放射性问题中)，以及‚边际的‛(在经济学中)等． 2、建立瞬时表达式：根据自变量有微小改变△t时，因变量的增量△W，建立起在时段△t上的增量表达式，令 dW △t →0，即得到的表达式．
微分方程模型
一、微分方程建模简介
二、微分方程模型
三、微分方程案例分析四、微分方程的MATLAB求解五、微分方程综合案例分析
一、微分方程模型简介

3.微分方程模型(缉私艇追击走私船问题)解析

c 2
1 1 r
x c
1r
r
1 ，1
c r1
x
cr r2 1
dy
1
x
r
c
r
dx 2 c x
y(c) 0
当 x 0 时， y 缉私艇不可能追赶上走私船。
3)
r 1 ，
y
1 2
x2 c2 2c
c ln
x c
当 x 0 时， y ，缉私艇不可能追赶上走私船。
微分方程模型实例2——缉私艇追击走私船问题
六. 实验任务
1. 有一只猎狗在B点位置发现了一只兔子在正东北方距离它200米的地方O处，此时兔子开始以8米/秒的速度向正西北方距离为120米的洞口A全速跑去，假设猎狗在追赶兔子的时候始终朝着兔子的方向全速奔跑，用计算机仿真法等多种方法完成下面的实验： (1) 问猎狗能追上兔子的最小速度是多少？ (2) 在猎狗能追上兔子的情况下，猎狗跑过的路程是多少？ (3) 画出猎狗追赶兔子奔跑的曲线图。 (4) 假设在追赶过程中，当猎狗与兔子之间的距离为30米时，兔子由于
t=t+dt; jstx=jstx-b*dt*jstx/sqrt(jstx^2+(a*t-jsty)^2); jstxb=[jstxb,jstx]; jsty=jsty+b*dt*(a*t-jsty)/sqrt(jstx^2+(a*t-jsty)^2); jstyb=[jstyb,jsty]; zscy=a*t; zscyb=[zscyb,zscy]; end zscxb=zeros(length(zscyb)); plot(jstxb,jstyb,zscxb,zscyb,'*')
5.3.2 缉私艇追击走私船问题