微分方程模型

微分方程模型

重点：

车间空气清洁问题、减肥问题、单种群增长问题与多物种相互作用问题等数学模型的建立过程与所使用的方法

要求： 1．进一步理解建模基本方法与基本建模过程，掌握平衡原理与微元法在建模中的用法．

所谓平衡原理是指自然界的任何物质在其变化的过程中一定受到某种平衡关系的支配.注意发掘实际问题中的平衡原理是从物质运动机理的角度组建数学模型的一个关键问题.就象中学的数学应用题中等量关系的发现是建立方程的关键一样.

微元法是指在组建对象随着时间或空间连续变化的动态模型时，经常考虑它在时间或空间的微小单元变化情况，这是因为在这些微元上的平衡关系比较简单，而且容易使用微分学的手段进行处理．这类模型基本上是以微分方程的形式给出的．

例1 设警方对司机饮酒后驾车时血液中酒精含量的规定为不超过80%(mg/ml). 现有一起交通事故,在事故发生3个小时后,测得司机血液中酒精含量是56%(mg/ml), 又过两个小时后, 测得其酒精含量降为40%(mg/ml),试判断: 事故发生时,司机是否违反了酒精含量的规定? 解：模型建立

设)(t x 为时刻t 的血液中酒精的浓度, 则依平衡原理时间间隔],[t t t ∆+内, 酒精浓度的改变量t t x x ∆⋅∝∆)(, 即

t t kx t x t t x ∆-=-∆+)()()(

其中k >0为比例常数, 式前负号表示浓度随时间的推移是递减的, 遍除以t ∆, 并令0→∆t , 则得到

,d d kx t x

-= 且满足40)5(,56)3(==x x 以及0)0(x x =.

模型求解

容易求得通解为kt c t x -=e )(, 代入0)0(x x =,得到

kt x t x -=e )(0

则)0(0x x =为所求. 又由,40)5(,56)3(==x x 代入0)0(x x =可得

17.040

56e 40

e 56

e 25030=⇒=⇒⎩⎨⎧==--k x x k k

k 将17.0=k 代入得 25.93e 5656e 17.03017.030≈⋅=⇒=⨯⨯-x x >80 故事故发生时,司机血液中的酒精浓度已超出规定.

2．理解种群的相互关系模型的建立原理与结论． ∙ 马尔萨斯模型 模型假设

（1）初始种群规模已知（设为N 0），种群数量非常大，世代互相重叠，因此种群的数量可以看作是连续变化的；

（2）种群在空间分布均匀，没有迁入和迁出（或迁入和迁出平衡）；

（3）种群的出生率和死亡率为常数，即不区分种群个体的大小、年龄、性别等. （4）环境资源是无限的. 确定变量和参数 :t 自变量，t t N :)(时刻的种群密度， :b 出生率，:d 死亡率.

模型的建立与求解

由上述假设，单种群增长模型与马尔萨斯人口模型极为类似，于是使用完全相同的建模过程易得

)(:)()(d )

(d t rN t N d b t t N =-= （3.1） 满足初始条件0)0(N N =的解为

.e e )(0)(0rt t d b N N t N ==-

于是有

,)(lim ,,0+∞=>>+∞

→t N d b r t 则有即

,)(lim ,,00N t N d b r t ===+∞

→则有即

,0)(lim ,,0=<<+∞

→t N d b r t 则有即

在种群生长的初期，种群规模较小，有足够的生存空间、足够的食物，彼此间没有利益冲突.但随着种群规模的逐渐扩大，对有限的空间、食物和其他生存必须条件的种内竞争越来越激烈，这必然影响种群的出生率和死亡率，从而降低实际增长率，因而在上述模型中假设出生率、死亡率为常数，资源无限不尽合理.

∙ 罗捷斯蒂克模型

完全类似于人口模型的分析知道，种群的增长模型为

⎪⎩

⎪

⎨

⎧=-=,)0(),1(0N N K N rN dt dN

（3.2） 其中r 是种群的固有（N =0时）增长率，K 是环境的最大容纳量.

方程（3.2）既是变量可分离方程，又是贝努利型方程.容易求得其解为

00

)()(N e N K KN t N rt +-=

- （3.3）

3．会建立较为简单的相关实际问题的数学模型．

例2 在凌晨1时警察发现一具尸体, 测得尸体温度是29︒C, 当时环境温度是21︒C .

一小时后尸体温度下降到27︒C , 若人的正常体温是37︒C , 估计死者的死亡时间.

解 运用牛顿冷却定律T ')(T T out -=-α, 得到它的通解为 )(0out out T T T T -+=t α-e , 这里0T 是当0=t 时尸体的温度, 也就是所求的死亡时间时尸体的温度, 将题目提供的参数代入:

⎩⎨⎧=-+=-++--27e

)2137(2129

e )2137(21)

1(t t αα 解得: 168e =-t α 和 166e )1(=+-t α 则34e =α 求得:

)(409.2)

12(,2877.0h Ln t ≈-

=≈α

α 这时求得的t 是死者从死亡起到尸体被发现所经历的时间, 因此反推回去可推测死

者的死亡时间大约是前一天的夜晚10:35.

例3 设某种动物头数的变化服从Logistic 规律．在正常情况下净相对增长率为a 1，环境容许的极限头数为N 1．假设当头数增加到Q （Q < N 1）时瘟疫流行，使净相对增长率为a 2，极限头数降为N 2（N 2< Q ），于是头数下降．当降至q （q >N 2）时，瘟疫停止，恢复正常．试建立这种情况下动物头数的模型，并讨论在瘟疫影响下动物头数的周期性变化，周期与哪些因素有关．

解 由题中条件知，动物头数x (t )应满足：

⎪⎪⎩

⎪⎪⎨⎧-=-=瘟疫流行时

正常时

)1(~d ~

d )1(d d 22

1

1N x x a t x N x x a t x

解得

⎪

⎪⎪⎩

⎪

⎪⎪⎨

⎧-+=-+=----瘟疫流行时

正常时

)(2

2

)

(111201e )1(1)(~e )1(1)(t t a t t a Q N N t x q N N t x

其中10,t t 分别为开始转入正常的时刻和开始转入瘟疫流行的时刻，由

Q q

N

N t x t t a =-+=

--)

(11

01e )1(1)(

解得 )

()(ln 11110Q N q q N Q a t t --=

-