微分方程模型

浅夏110

微分方程建模是数学建模的重要方法，因为许多实际问题的数学描述将导致求解微 分方程的定解问题。把形形色色的实际问题化成微分方程的定解问题，大体上可以按以 下几步：
4 ?; R* I  o* q3 w2 u/ z
. k5 v; b  _0 H: C- y4 X% E, r( H1. 根据实际要求确定要研究的量(自变量、未知函数、必要的参数等)并确定坐标系。

2. 找出这些量所满足的基本规律(物理的、几何的、化学的或生物学的等等)。

+ z. V# P. I2 H3 p2 g! x; G3. 运用这些规律列出方程和定解条件。
3 {, U. T: _% R; S4 _
列方程常见的方法有：

（i）按规律直接列方程
! C( P; A- C* y
在数学、力学、物理、化学等学科中许多自然现象所满足的规律已为人们所熟悉， 并直接由微分方程所描述。如牛顿第二定律、放射性物质的放射性规律等。我们常利用 这些规律对某些实际问题列出微分方程。

（ii）微元分析法与任意区域上取积分的方法
+ L3 Y0 g5 R) c; \) h
自然界中也有许多现象所满足的规律是通过变量的微元之间的关系式来表达的。对 于这类问题，我们不能直接列出自变量和未知函数及其变化率之间的关系式，而是通过 微元分析法，利用已知的规律建立一些变量（自变量与未知函数）的微元之间的关系式， 然后再通过取极限的方法得到微分方程，或等价地通过任意区域上取积分的方法来建立 微分方程。
/ M. G0 w' a  H* `7 X/ ?9 D
6 W8 `2 C4 c( B4 ~. d2 p（iii）模拟近似法

- |0 f" Q; H  W" N在生物、经济等学科中，许多现象所满足的规律并不很清楚而且相当复杂，因而需 要根据实际资料或大量的实验数据，提出各种假设。在一定的假设下，给出实际现象所 满足的规律，然后利用适当的数学方法列出微分方程。 在实际的微分方程建模过程中，也往往是上述方法的综合应用。不论应用哪种方法， 通常要根据实际情况，作出一定的假设与简化，并要把模型的理论或计算结果与实际情 况进行对照验证，以修改模型使之更准确地描述实际问题并进而达到预测预报的目的。 本章将利用上述方法讨论具体的微分方程的建模问题。
3 F% w. a8 H; A
. _) b  ]% Y" k9 D§1 发射卫星为什么用三级火箭
采用运载火箭把人造卫星发射到高空轨道上运行，为什么不能用一级火箭而必须用 多级火箭系统？ 下面通过建立运载火箭有关的数学模型来回答上述问题。 火箭是一个复杂的系统，为了使问题简单明了，我们只从动力系统和整体结构上分 析，并且假设引擎是足够强大的。
* F5 \8 F3 e( t3 s+ J& t
1.1 为什么不能用一级火箭发射人造卫星
  }+ J; ], L, m% \下面用三个数学模型回答这个问题

1.1.1 卫星进入 600km 高空轨道时，火箭必须的最低速度 首先将问题理想化，假设：

$ r7 u  K5 R1 b（i）卫星轨道是以地球中心为圆心的某个平面上的圆周，卫星在此轨道上以地球 引力作为向心力绕地球作平面匀速圆周运动；
7 h; g# f9 d: T
（ii）地球是固定于空间中的一个均匀球体，其质量集中于球心；

（iii）其它星球对卫星的引力忽略不计。 建模与求解：设地球半径为 R ，质量为 M ；卫星轨道半径为 r ，卫星质量为m 。
1 j! I. p4 V2 n) a+ j7 O
. B: ~4 L. F* l' W根据假设（ii）和（iii），卫星只受到地球的引力，由牛顿万有引力定律可知其引力大小为



1 M. E( N0 Y/ S1 K) ]即要把卫星送入离地面 600km 高的轨道，火箭的末速度最低应为 7.6km/s。

+ Q3 |6 E0 s3 I( q3 I1.1.2 火箭推进力及升空速度
2 ^5 Z2 {0 d: I( F火箭的简单模型是由一台发动机和一个燃料仓组成。燃料燃烧产生大量气体从火箭 末端喷出，给火箭一个向前的推力。火箭飞行要受地球引力、空气阻力、地球自转与公 转等的影响，使火箭升空后作曲线运动。为使问题简化，假设：

$ g0 z8 a+ h* q" @（i）火箭在喷气推动下作直线运动，火箭所受的重力和空气阻力忽略不计。

（ii）在t 时刻火箭质量为m(t) ，速度为v(t) ，且均为时间t 的连续可微函数；

& {2 o$ _4 G9 a- B& J3 i（iii）从火箭末端喷出气体的速度（相对火箭本身）为常数u 。

: |, r0 q- K& h7 C建模与分析：由于火箭在运动过程中不断喷出气体，使其质量不断减少，在 (t,t + Δt) 内的减少量可由台劳展式表示为


* P4 u. z! W/ O' }2 F


$ [4 Y( }, Z' B0 ?+ k1.1.3 一级火箭末速度上限

- {2 N9 A. w/ i: Y" y7 g


因此，用一级火箭发射卫星，在目前技术条件下无法达到相应高度所需的速度。

1.2 理想火箭模型
从前面对问题的假设和分析可以看出：火箭推进力自始至终在加速着整个火箭，然 而随着燃料的不断消耗，所出现的无用结构质量也在随之不断加速，作了无用功，故效 益低，浪费大。

+ X! E- _5 s6 j9 u8 ^所谓理想火箭，就是能够随着燃料的燃烧不断抛弃火箭的无用结构。下面建立它的 数学模型。

假设：在(t,t + Δt) 时段丢弃的结构质量与烧掉的燃料质量以α 与1−α 的比例同 时进行。

建模与分析：由动量守恒定律，有

2 M  _2 L; s$ |0 G3 u4 N& B. o2 D
4 d( T9 @. ^) m* Z4 Q. _& G
1.3 多级火箭卫星系统
理想火箭是设想把无用结构质量连续抛弃以达到最佳的升空速度，虽然这在目前的 技术条件下办不到，但它确为发展火箭技术指明了奋斗目标。目前已商业化的多级火箭 卫星系统便是朝着这种目标迈进的第一步。多级火箭是从末级开始，逐级燃烧，当第i 级燃料烧尽时，第i +1级火箭立即自动点火，并抛弃已经无用的第i 级。我们用 表 示第i 级火箭质量， 表示有效负载。为了简单起见，先作如下假设:
% ], w; V/ a# e) \



: k" o) d" r9 q9 k8 O, c
实际上，由于受技术条件的限制，采用四级或四级以上的火箭，经济效益是不合算 的，因此采用三级火箭是最好的方案。
7 H" `9 [4 t! a9 z0 C* s
3 Q5 E3 R1 R5 h* R1.4 最佳结构设计
' E, `2 N3 |- Q; P8 I6 T/ F% p
+ t/ G* L! b- g! ]& x  n( ]1 F& a


: a1 l  V- x; N5 U% I§2 人口模型
5 F! i2 N$ x2 c" Z* s7 d2.1 问题提出
据考古学家论证，地球上出现生命距今已有 20 亿年，而人类的出现距今却不足 200 万年。纵观人类人口总数的增长情况，我们发现：1000 年前人口总数为 2.75 亿。经过 漫长的过程到 1830 年，人口总数达 10 亿，又经过 100 年，在 1930 年，人口总数达 20 亿；30 年之后，在 1960 年，人口总数为 30 亿；又经过 15 年，1975 年的人口总数是 40 亿，12 年之后即 1987 年，人口已达 50 亿。 我们自然会产生这样一个问题：人类人口增长的规律是什么？如何在数学上描述这 一规律。

2.2 Malthus 模型
1789 年，英国神父 Malthus 在分析了一百多年人口统计资料之后，提出了 Malthus 模型。 模型假设
3 U; [' P- Z5 X) m" G( K
! h6 q) [4 F" {9 c2 H（i）设 x(t) 表示t 时刻的人口数，且 x(t) 连续可微。

（ii）人口的增长率 r 是常数（增长率=出生率—死亡率）。
2 y5 g6 S: k  |6 b' G3 A2 x2 l
（iii）人口数量的变化是封闭的，即人口数量的增加与减少只取决于人口中个体的 生育和死亡，且每一个体都具有同样的生育能力与死亡率.

建模与求解

由假设，t 时刻到t + Δt 时刻人口的增量为



根据 1700～1961 年间世界人口统计数据，我们发现这些数据与（16）式的计算结 果相当符合。因为在这期间地球上人口大约每 35 年增加 1 倍，而（16）式算出每 34.6 年增加 1 倍。 但是，当人们用（15）式对 1790 年以来的美国人口进行检验，发现有很大差异。

利用（16）式对世界人口进行预测，也会得出惊异的结论：当 t = 2670 年时 x(t) = 4.4× ，即 4400 万亿，这相当于地球上每平方米要容纳至少 20 人。

显然，用这一模型进行预测的结果远高于实际人口增长，误差的原因是对增长率 r 的估计过高。由此，可以对 r 是常数的假设提出疑问。

; E7 g2 y* I! i' u' t: _5 F2.3 阻滞增长模型（Logistic 模型）
6 a8 u/ P  i$ b2 R' M6 u) F如何对增长率 r 进行修正呢？我们知道，地球上的资源是有限的，它只能提供一定 数量的生命生存所需的条件。随着人口数量的增加，自然资源、环境条件等对人口再增 长的限制作用将越来越显著。如果在人口较少时，我们可以把增长率 r 看成常数，那么 当人口增加到一定数量之后，就应当视 r 为一个随着人口的增加而减小的量，即将增长 率 r 表示为人口 x(t) 的函数 r(x) ，且 r(x) 为 x 的减函数。

模型假设
) ]. u# j. x5 H
/ b7 ~1 K/ ?: i$ t0 V
- R. b0 a- m+ V1 X- p) ]

2 S; @6 t( [, {% w
实际数据都能较好地吻合，在 1930 年之后，计算与实际偏差较大。原因之一是 60 年代 的实际人口已经突破了假设的极限人口 ，由此可知，本模型的缺点之一就是不易确 定 。

2.4 模型推广
* Y( ]2 @  b$ c3 C1 S' V可以从另一个角度导出阻滞增长模型，在 Malthus 模型上增加一个竞争项   ( b > 0)  ，它的作用是使纯增长率减少。如果一个国家工业化程度较高，食品供应 较充足，能够提供更多的人生存，此时b 较小；反之b 较大，故建立方程


% }9 c- ^. j3 A+ F6 X
参数a 和b 可以通过已知数据利用 Matlab 中的非线性回归命令 nlinfit 求得。
  n: w- K/ E. H2 H8 v+ H
§3 战争模型
早在第一次世界大战期间，F. W. Lanchester 就提出了几个预测战争结局的数学模 型，其中包括作战双方均为正规部队；作战双方均为游击队；作战的一方为正规部队， 另一方为游击队。后来人们对这些模型作了改进和进一步的解释，用以分析历史上一些 著名的战争，如二次世界大战中的美日硫黄岛之战和 1975 年的越南战争。 影响战争胜负的因素有很多，兵力的多少和战斗力的强弱是两个主要的因素。士兵 的数量会随着战争的进行而减少，这种减少可能是因为阵亡、负伤与被俘，也可能是因为疾病与开小差。分别称之为战斗减员与非战斗减员。士兵的数量也可随着增援部队的 到来而增加。从某种意义上来说，当战争结束时，如果一方的士兵人数为零，那么另一 方就取得了胜利。如何定量地描述战争中相关因素之间的关系呢？比如如何描述增加士 兵数量与提高士兵素质之间的关系。

2 A; M% r2 F$ t8 R1 N$ Y& m+ ^3.1 模型一  正规战模型
. [1 A( K+ d; S$ ?1 q4 a, A$ P模型假设
6 J" f3 [" z5 l/ g
" ]! t  a: R6 t8 [


' @6 D8 x3 D; k
" P# k' D: @/ x; S# O

( R; a# ?4 E2 c; i; y以上是研究双方之间兵力的变化关系。下面将讨论每一方的兵力随时间的变化关 系。

( E- n9 c$ l9 q! f  l

) M; w6 M" b- P! {- i' V, @5 c5 ^

) {0 S. d/ ~; A1 w/ ^/ y, X. p3.2 模型二 游击战模型
模型假设

) M, j- g* M* K% l+ u& ~
. o! }& |. z% S

! Z4 l" p0 P& e3 B! P/ l
% m+ d8 J- y9 {
' f1 \0 h, B' y1 @' `& r
7 L% N* `, N9 h6 Y
3 N& J& Q9 E* b) f3.3 模型三 混合战模型
模型假设

1 d0 T8 U: b) S" g; X/ q（i） x 方为游击队， y 方为正规部队。

6 D( Y. J4 k( Y（ii）交战双方均无战斗减员与增援。 模型与求解 借鉴模型一与二的思想，可得


, `- Y5 Q+ A" m8 c
4 R% U7 Z% I5 Z# q) }一般来说，正规部队以火力强而见长，游击队以活动灵活，活动范围大而见长。这 可以通过一些具体数据进行计算。

+ ]  x/ Y7 z: R3 b7 g& c+ Q" p+ @" u
0 W' S4 q) t. e) ?; s
美国人曾用这个模型分析越南战争。根据类似于上面的计算以及四五十年代发生在 马来亚、菲律宾、印尼、老挝等地的混合战争的实际情况估计出，正规部队一方要想取 胜必须至少投入 8 倍于游击部队一方的兵力，而美国至多只能派出 6 倍于越南的兵力。 越南战争的结局是美国不得不接受和谈并撤军，越南人民取得最后的胜利。

3.4 模型四 一个战争实例
# t  B* a6 M7 y6 ]J. H. Engel 用二次大战末期美日硫黄岛战役中的美军战地记录，对正规战争模型进 行了验证，发现模型结果与实际数据吻合得很好。 硫黄岛位于东京以南 660 英里的海面上，是日军的重要空军基地。美军在 1945 年 2 月开始进攻，激烈的战斗持续了一个月，双方伤亡惨重，日方守军 21500 人全部阵亡 或被俘，美方投入兵力 73000 人，伤亡 20265 人，战争进行到 28 天时美军宣布占领该 岛，实际战斗到 36 天才停止。美军的战地记录有按天统计的战斗减员和增援情况。日 军没有后援，战地记录则全部遗失。

用 A(t) 和 J (t) 表示美军和日军第t 天的人数，忽略双方的非战斗减员，则
' D. V# T) v. \
) D# A' G/ v1 l2 T

6 K* i8 w. I  s9 U' o- t
2 ?& H6 k7 c1 G' F
————————————————
版权声明：本文为CSDN博主「wamg潇潇」的原创文章，遵循CC 4.0 BY-SA版权协议，转载请附上原文出处链接及本声明。
2 h6 W% p, O9 d. Q, b" \" [原文链接：https://blog.csdn.net/qq_29831163/article/details/89702947

! N2 E; z3 {5 J/ w" S) }5 i/ v