第五章微分方程模型
微分方程模型的建立与求解
微分方程模型的建立与求解微分方程是自然界中许多现象的数学描述，通过建立微分方程模型可以更好地理解和预测各种现象。
本文将介绍微分方程模型的建立与求解方法。
一、微分方程模型的建立微分方程通常用来描述系统内部的变化规律，要建立微分方程模型，首先需要根据具体问题分析系统的特点，确定影响系统变化的因素，并建立相关的数学表达式。
以一个简单的弹簧振子系统为例，假设弹簧的位移为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）定性理论方法。
建立微分方程模型的方法：1）利用数学、力学、物理、化学等学科中的定理或经过实验检验的规律等来建立微分方程模型。
2）微元分析法利用已知的定理与规律寻找微元之间的关系式，与第一种方法不同的是对微元而不是直接对函数及其导数应用规律3）模拟近似法在生物、经济等学科的实际问题中，许多现象的规律性不很清楚，即使有所了解也是极其复杂的，建模时在不同的假设下去模拟实际的现象，建立能近似反映问题的微分方程，然后从数学上求解或分析所建方程及其解的性质，再去同实际情况对比，检验此模型能否刻画、模拟某些实际现象。
下面我们以生态学模型为例介绍微分方程模型的建立过程： 一. 单种群模型1. 马尔萨斯(Malthus)模型假定只有一个种群，()N t 表示t 时刻生物总数，r 表示出生率，0t 表示初始时刻，则生物总数增长的数学模型为()()()00d ,d (1)t t N t rN t t N t N =⎧=⎪⎨⎪=⎩不难得到其解为()0()0r t t N t N e-=.2. 密度制约模型由马尔萨斯模型知，种群总数将以几何级数增长，显然与实际不符，因为种群密度增大时，由于食物有限，生物将产生竞争，或因为传染病不再按照增长率r 增长，因而有必要修改，在(1)式右端增加一项竞争项。
()()()d (1)(2)d N t N t rN t tK=-其中K 为最大容纳量，可以看出当()N t K =时，种群的规模不再增大。
这个模型就是著名的Logistic 模型，可以给出如下解释：由于资源最多仅能维持K 个个体，故每个个体平均需要的资源为总资源的1K，在t 时刻个体共消耗了总资源的()N t K此时资源剩余()1N t K-，因此Logistic 模型表明：种群规模的相对增长率与当时所剩余的资源份量成正比，这种种群密度对种群规模增长的抑制作用。
培训资料--微分方程模型人口模型等
x0
x0 0
t
人口发展方程
• 年龄分布对于人口预测的重要性 • 只考虑自然出生与死亡，不计迁移
F(r,t) ~ 人口分布函数 (年龄 r的人口） p(r,t) ~ 人口密度函数 N(t) ~ 人口总数 r ( ) ~ 最高年龄
m
F(0,t) 0, F(r ,t) N(t) m p(r, t) F r
人口发展方程
f
(t
)
(t
)r2 r1
h(r,
t
)k
(r,
t
)
p(r,
t
)dr
(t) ~总和生育率——控制生育的多少
h(r, t) ~生育模式——控制生育的早晚和疏密
p(r,t)
p 0
(r
r
t)e (s)ds r t r
,
0
t
r
f (t r)e0(s)ds , t r
p0 (r)
• 正反馈系统
(r,t) p(r,t)dt, dt dr1
p p (r,t) p(r,t) 一阶偏微分方程
r t
p
r
p t
(r, t )
p(r, t )
人口发展方程
p(r,0) p0 (r), r 0 ~已知函数（人口调查）
p(0,
t)
f
(t ),
t0
~生育率（控制人口手段）
p0 (0) f (0) --------相容性条件
b(r,
t)k
(r,
t)
p(r,
t)dr
b(r,t) (t)h(r,t)
0
r2 r1
h(r , t )dr
1
h~生育模式
(t)
常微分方程第五章 微分方程建模案例
第五章微分方程建模案例微分方程作为数学科学的中心学科，已经有三百多年的发展历史，其解法和理论已日臻完善，可以为分析和求得方程的解（或数值解）提供足够的方法，使得微分方程模型具有极大的普遍性、有效性和非常丰富的数学内涵。
微分方程建模包括常微分方程建模、偏微分方程建模、差分方程建模及其各种类型的方程组建模。
微分方程建模对于许多实际问题的解决是一种极有效的数学手段，对于现实世界的变化，人们关注的往往是其变化速度、加速度以及所处位置随时间的发展规律，其规律一般可以用微分方程或方程组表示，微分方程建模适用的领域比较广，涉及到生活中的诸多行业，其中的连续模型适用于常微分方程和偏微分方程及其方程组建模，离散模型适用于差分方程及其方程组建模。
本章主要介绍几个简单的用微分方程建立的模型，让读者一窥方程的应用。
下面简要介绍利用方程知识建立数学模型的几种方法：1．利用题目本身给出的或隐含的等量关系建立微分方程模型这就需要我们仔细分析题目，明确题意，找出其中的等量关系，建立数学模型。
例如在光学里面，旋转抛物面能将放在焦点处的光源经镜面反射后成为平行光线，为了证明具有这一性质的曲线只有抛物线，我们就是利用了题目中隐含的条件——入射角等于反射角来建立微分方程模型的。
2．从一些已知的基本定律或基本公式出发建立微分方程模型我们要熟悉一些常用的基本定律、基本公式。
例如从几何观点看，曲线y=上某点的切线斜率即函数)yy=在该点的导数；力学中的牛顿第二运(x)(xy动定律：maF=，其中加速度a就是位移对时间的二阶导数，也是速度对时间的一阶导数等等。
从这些知识出发我们可以建立相应的微分方程模型。
例如在动力学中，如何保证高空跳伞者的安全问题。
对于高空下落的物体，我们可以利用牛顿第二运动定律建立其微分方程模型，设物体质量为m，空气阻209210力系数为k ，在速度不太大的情况下，空气阻力近似与速度的平方成正比；设时刻t 时物体的下落速度为v ，初始条件：0)0(=v . 由牛顿第二运动定律建立其微分方程模型：2kv mg dtdv m -= 求解模型可得：)1]2(exp[)1]2(exp[+-=mkg t k m kg tmg v 由上式可知，当+∞→t 时，物体具有极限速度：kmg v v t ==∞→lim 1， 其中，阻力系数s k αρ=，α为与物体形状有关的常数，ρ为介质密度，s 为物体在地面上的投影面积。
第五章 微分方程模型讲1
i0
1-1/σ σ
di 1 = −λi[i − (1 − )] σ =λ/ µ dt σ
σ >1
i
σ ≤1
di/dt < 0
i0
0
1-1/σ σ
1 i
i0
0
1 , σ > 1 1 − i(∞ ) = σ 0, σ ≤ 1
t
0
t
接触数σ =1 ~ 阈值
σ >1
σ ≤ 1 ⇒ i (t ) ↓
s i ( s ) = ( s 0 + i0 ) − s + ln σ s0
i
1
1D = {( s ,源自i ) s ≥ 0 , i ≥ 0 , s + i ≤ 1}
D 0
s
1
模型4 模型
相轨线 i ( s ) 及其分析
i
1 D
SIR模型 模型
s i(s) = (s0 + i0 ) − s + ln σ s0
dP dP = kP(10000− P) 把 P t=0 =10, = 100代入微分方程 dt dt t=0
1 得 k= 999 鸟的数量和时间的函数关系为 P =
10000 1+ 999 e
− 10000 t 999
Logistic函数 函数
5.1 传染病模型
问题
• 描述传染病的传播过程 • 分析受感染人数的变化规律 • 预报传染病高潮到来的时刻 • 预防传染病蔓延的手段 • 按照传播过程的一般规律， 按照传播过程的一般规律， 用机理分析方法建立模型 已感染者(the infective) 易感染者 易感染者(the susceptible) 已感染者 移出者(the removed) 移出者
方程模型
例2. 解初值问题
x yd x ( x 2 1 ) d y 0
y( 0 ) 1
dy x 解: 分离变量得 dx 2 y 1 x
两边积分得
即
y x2 1 C
( C 为任意常数 )
由初始条件得 C = 1, 故所求特解为
y x 1 1
数学建模－ 微分方程模型
关晓飞 同济大学数学科学学院
一、什么是微分方程？
最最简单的例子
引例
一曲线通过点（1，2），且在该曲线任一点 若设曲线方程为 y f ( x) ， (1)
M( x ,y )处的切线的斜率为2x，求该曲线的方程。
解
根据导数的几何意义可知未知函数满足关系式：
dy 2x dx
规律。
解
dM 铀的衰变速度就是 M (t ) 对时间t的导数 dt
，
由于衰变速度与其含量成正比，可知未知函数满足 关系式： dM M (1) ( 0) 是衰变系数
dt
且初始条件 M t 0 M0 dM dt 分离变量得 M 对上式两端积分得：ln M t ln c 因此， M (t ) Cet 代入初始条件得
建立模型：设物体在冷却过程中的温度为 T(t)，t≥0，
“T的变化速率正比于T与周围介质的温度差” 翻译为
dT 与 T m 成正比 dt
dT k (T m ), dt T ( 0) 60.
建立微分方程
数学语言
其中参数k >0，m=18. 求得一般解为
ln(T－m)=－k t+c,
或
T m ce
kt
, t 0,
微分方程模型-药物中毒急救模型(附matlab编程)
100μg/ml浓度会出现严重中毒, 200μg/ml浓度可致命.
医生需要判断：
（1）孩子的血药浓度会不会达到100~200 μg/ml？
（2）如果会达到，应采取怎样的紧急施救方案？
调查与分析
口服药物
胃肠道 药量x(t)
转移速率
和x成正比
血液系统
体外
了解药物的半衰期，可以帮助制定合理的用药 间隔，并且可以预知药物在体内的变化趋势。
调查与分析
通常，血液总量约为人体体重的7 % ~8%，体重 50~60 kg的成年人有4000ml左右的血液. 目测这个孩子的体重约为成年人的一半，可认为 其血液总量约为2000ml.
调查与分析
临床施救的办法
• 口服活性炭来吸附药物，可使药物的排除率 增加到原来（人体自身）的2倍.
在t=0瞬间进入胃肠道.
dx x, x(0) 1100
dt
y(t)由吸收而增长的速度是λx，由排除而减少的速度与y(t)
成正比(比例系数μ) , t=0时血液中无药物.
dy x y, y(0) 0
dt
模型求解
dx x, x(0) 1100
dt
药物吸收的半衰期为5 h
x(t) 1100et
2. 血液系统中药物的排除速率与y(t) 成正比，比例系数 μ(>0)，t=0时血液中无药物.
3. 氨茶碱被吸收的半衰期为5 h，排除的半衰期为6 h.
4. 孩子的血液总量为2000 ml.
模型建立
口服药物
胃肠道 药量x(t)
转移率和
x成正比
血液系统
体外
药量y(t)
排除速率 和y成正比
微分方程模型的基本原理
微分方程模型的基本原理微分方程是数学中描述变化的一种重要工具，它能够描述系统中随时间、空间或者其他变量而发生的变化规律。
微分方程模型是一种基于微分方程的数学模型，用于描述各种实际问题的变化过程。
1.变量与变化率的关系：微分方程模型描述了系统中变量随时间的变化率，即变量的导数。
它指出了变量如何随时间而变化，从而提供了数量化的描述。
2.初始条件和边界条件：微分方程模型需要给定初始条件和边界条件，以确定具体的解。
初始条件是在系统起始时给定的变量值，边界条件是在系统边界上给定的限制条件。
这些条件可以是实际问题中必须满足的条件。
3.多变量之间的关系：微分方程模型可以涉及多个变量之间的相互作用。
这些变量可以表示不同的物理量或者变化过程，它们之间的关系可以是线性的、非线性的、常系数的或者变系数的。
这些关系可以通过微分方程进行描述。
4.具体问题的建模过程：微分方程模型的建立需要针对具体问题进行分析和建模过程。
这个过程中需要确定问题中涉及的变量、关系以及边界条件，并将其转化为合适的微分方程模型。
这个过程可以涉及到数学推理、物理实验、统计分析等多个方面。
微分方程模型的应用非常广泛，几乎涉及到各个学科领域。
例如，在物理学中，微分方程模型可以用于描述粒子的运动、电磁场的分布、热传导等问题；在经济学中，微分方程模型可以用于描述市场供需关系、经济增长等问题；在生物学中，微分方程模型可以用于描述生物种群的演化、药物动力学等问题。
微分方程模型的求解方法也非常丰富多样，可以通过数值方法、解析方法、近似方法等进行求解。
数值方法通过将微分方程转化为差分方程，然后采用逼近的方式进行求解。
解析方法通过数学推导和变量分离的方式求得方程的解析解。
近似方法通过针对特定问题的特殊性质，利用适当的近似方法得到问题的近似解。
总之，微分方程模型是一种重要的数学工具，广泛用于各个学科领域中的问题描述和解决。
它通过描述变量与变化率的关系，建立初始条件和边界条件，描述多变量之间的关系等方面，为实际问题提供了准确的数学描述和求解方法。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
解得:
e??t
?
8 16
和
e ?? (t ?1)
?
6 16
则
e?
?
4 3
求得:
? ? 0.2877, t ? ? Ln(21 ) ? 2.409(h) ?
模型建立与求解
求得:
?
?
0.2877, t
?
?
Ln(21 )
?
?
2.409(h)
这时求得的t是死者从死亡起到尸体被发现所 经历的时间, 因此反推回去可推测死者的死亡时 间大约是前一天的夜晚10:35。
分析：考虑圆桶的极限速度
limv(t) ? W ? B ? 527.436? 470.327
t? ?
C
0.08
≈713.86（英尺/秒）＞＞40（英尺/秒）
实际极限速度与圆桶的承受速度相差巨大！
解决思路：避开求t0的难点
令 v(t)=v(y(t)), 其中 y=y(t) 是圆桶下沉深度
将
dv
?
dv dy .
v(300).(试一试)
* 用数值方法求出v(300)的近似值为
v(300)≈45.41（英尺/秒）＞40（英尺/秒）
* 分析 v=v(y) 是一个单调上升函数，而v 增
大,y 也增大,可求出函数y=y(v)
y?
?W g
(W C
?
W ? B W ? B ? Cv C 2 ln W ? B ),
令 v=40(英尺/秒)，g=32.2（英尺/秒），
代入（1）得
dt dy dt
m dv . dy ? W ? B ? Cv , dy dt
或
?
? ?
W
?
v B ? Cv
dv dy
?
g W
,
?? v ( 0 ) ? 0 , y ( 0 ) ? 0 .
两边积分得函数方程：
?
v C
?
W?B C2
ln
W
?B W?
? Cv B
?
gy W
,
若能求出函数v=v(y),就可求出碰撞速度
例3：放射性废物的处理问题 美国原子能委员会（现为核管理委员会）
处理浓缩放射性废物，是将废物放入密封性能 很好的圆桶中，然后扔到水深300英尺的海里. 他们这种做法安全吗？ 联想：安全 、危险 分析：可从各个角度去分析造成危险的因素， 这里仅考虑圆桶泄露的可能.
问题的关键
* 圆桶至多能承受多大的冲撞速度？(40英尺/秒)
算出 y= 238.4 (英尺)＜300（英尺）
问题的实际解答：美国原子能委员会处理 放射性废物的做法是极其危险的，必须改变。
例4：崖高的估算
假如你站在崖顶且身上带着一只具有跑表功 能的计算器，你也许会出于好奇心想用扔下 一块石头听回声的方法来估计山崖的高度， 假定你能准确地测定时间，你又怎样来推算 山崖的高度呢，请你分析一下这一问题。
F阻=0.02m
下一页
铁链下滑示意图
t= 0
t时刻
x(t)
L/5
2x+L/5
x(t)
上一页
例2：刑事侦察中死亡时间的鉴定
问题描述
在凌晨 1时警察发现一具尸体 ,测得尸体温度是 29?C,当时环境温度是 21?C。一小时后尸体温度下 降到 27?C,若人的正常体温是 37?C,估计死者的死 亡时间。
例1：铁链下滑问题
问题： 一条长为 L米质量为 M的链条悬挂在一个钉子 上，初始时，一边长 3/5L，另一边长 2/5L，由静止 启动。分别根据以下情况求出链条下滑的时间： 1、不计摩擦力和空气阻力； 2、阻力为 1/10L的链条重； 3、阻力与速度 v成正比; 4、摩擦力与对钉子的压力成正比，在 v=1时。
求解微分方程有三种方法： 1）求精确解； 2）求数值解（近似解）； 3）定性理论方法。
建立微分方程模型的方法 （1）根据规律列方程
利用数学、力学、物理、化学等学科中的定理或 经过实验检验的规律等来建立微分方程模型。
（2）微元分析法
利用已知的定理与规律寻找微元之间的关系 式，与第一种方法不同的是对微元而不是直 接对函数及其导数应用规律。
这里T0是当t=0时尸体的温度 ,也就是所求的死亡时 间时尸体的温度。
模型建立与求解
T ? Tout ? (T0 ? Tout )e ?? t
将题目提供的参数代入，得
? 21 ? (37 ? 21)e??t ? 29
? ?21
?
(37 ?
21)e ?? (t?1)
?
27
T0=37oC； Tout=21oC； T(t)=29oC；
其中
m
d2y dt2
?
W
?
B
?
D
(1)
m ? W , D ? Cv , dy ? v
g
dt
或
?? dv ? cg v ? g (W ? B), ? dt W W
(2)
??V(0) ? 0.
方程的解为
v(t)
?
W
?
B
(1 ?
? Cg t
e W ),
t? 0
C
计算碰撞速度，需确定圆桶和海底的碰撞时间t0
方法原理
牛顿冷却定律指出：物体在空气中冷却的速 度与物体温度和空气温度之差成正比，现将牛 顿冷却定律应用于刑事侦察中死亡时间的鉴定。
模型建立与求解
设尸体的温度为 T(t)(t从谋杀后计 )，运用牛顿
冷却定律得
速度
dT dt
??
(Tout
? T)
得到它的通解为
T ? Tout ? (T0 ? Tout )e ?? t
* 圆桶和海底碰撞时的速度有多大？ 问题：求这一种桶沉入300英尺的海底时的末 速度.（原问题是什么?） 可利用的数据条件：
圆桶的总重量 W=527.327（磅） 圆桶受到的浮力 B=470.327（磅） 圆桶下沉时受到的海水阻力 D=Cv，C=0.08 可利用牛顿第二定律，建立圆桶下沉位移 满足的微分方程：
（3）模拟近似法
在生物、经济等学科的实际问题中，许多现象 的规律性不很清楚，即使有所了解也是极其复 杂的，建模时在不同的假设下去模拟实际的现 象，建立能近似反映问题的微分方程，然后从 数学上求解或分析所建方程及其解的性质，再 去同实际情况对比，检验此模型能否刻画、模 拟某些实际现象。
5.2 一些简单的微分方程案例
第五章 微分方程模型
5.1 微分方程建模概述 5.2 简单微分方程模型案例 5.3 综合性微分方程模型
传染病模型 古尸断代பைடு நூலகம்
正规战与游击战
5.1 微分方程建模概述
动态 模型
? 描述对象特征随时间 (空间)的演变过程 ? 分析对象特征的变化规律 ? 预报对象特征的未来性态 ? 研究控制对象特征的手段
微分 方程 建模
? 根据函数及其变化率之间的关系确定函数 ? 根据建模目的和问题分析作出简化假设 ? 按照内在规律或用类比法建立微分方程
在研究实际问题时，常常会联系到某些变量的变 化率或导数， 这样所得到变量之间的关系式就是微分 方程模型。微分方程模型反映的是变量之间的间接关 系，因此，要得到直接关系，就得求微分方程。