欢迎光临散文网 会员登陆 & 注册

数学物理方法笔记2

2023-03-24 11:08 作者:皆是概率  | 我要投稿

传染病问题(Logistic模型)

一个区域有M只老鼠,其中N只患上了传染病。他们可以通过接触将病传染给健康的老鼠。问任意时刻患上传染病的老鼠有多少?

解:设任意t时刻病鼠和健康老鼠的数目分别为u和v,则

u%20%2B%20v%20%3D%20M%20%5Cquad%20(1.1.3)

一只病鼠能接触多少老鼠?

当区域内老鼠越多,能接触到的老鼠也就越多。假定这个比例数是%5Cbeta,那么能接触到%5Cbeta%20M只老鼠。

那么这%5Cbeta%20M只老鼠中,健康老鼠的比例为%5Cfrac%7Bv%7D%7BM%7D,那么一只病鼠能平均接触到%5Cbeta%20v只健康老鼠。

当前一共有u只病鼠,所以

病鼠数目的变化率正比于乘积uv,即

%5Cfrac%7Bdu%7D%7Bdt%7D%20%3D%20%5Cbeta%20u%20v%20%5Cquad%20(1.1.4)

这里非线性项uv刻画老鼠的接触性传染,比例系数。利用式(1.1.3),v = M - u,可以将方程(1.1.4)修改为

%5Cfrac%7Bdu%7D%7Bdt%7D%20%3D%20%5Cbeta%20u(M%20-%20u)

%5Cfrac%7Bdu%7D%7Bdt%7D%20%3D%20%5Cbeta%20Mu%20-%20%5Cbeta%20u%5E2

%5Cfrac%7Bdu%7D%7Bdt%7D%20%3D%20%5Calpha%20u%20-%20%5Cbeta%20u%5E2%20%20%5Cquad%20(1.1.5)

其中,。现在求解方程(1.1.5),首先将它写成

%5Cfrac%7Bdu%7D%7Bdt%7D%20%3D%20%5Calpha%20u%20-%20%5Cfrac%7B%5Calpha%7D%7BM%7Du%5E2

%5Cfrac%7Bdu%7D%7Bdt%7D%20%3D%20(u%20-%20%5Cfrac%7Bu%5E2%7D%7BM%7D)%5Calpha

%5Cfrac%7Bdu%7D%7Bdt%7D%20%3D%20(%5Cfrac%7BuM%7D%7BM%7D%20-%20%5Cfrac%7Bu%5E2%7D%7BM%7D)%5Calpha

%5Cfrac%7Bdu%7D%7Bdt%7D%20%3D%20(%5Cfrac%7BuM%20-%20u%5E2%7D%7BM%7D)%5Calpha

%5Cfrac%7Bdu%7D%7Bdt%7D%20%3D%20(%5Cfrac%7Bu(M%20-%20u)%7D%7BM%7D)%5Calpha

为了以后凑积分方便,交换分子上的方向。

%5Cfrac%7Bdu%7D%7Bdt%7D%20%3D%20-(%5Cfrac%7Bu(u-M)%7D%7BM%7D)%5Calpha

du%20%3D%20-(%5Cfrac%7Bu(u-M)%7D%7BM%7D)%5Calpha%20dt

Mdu%20%3D%20-(u(u-M))%5Calpha%20dt

%5Cfrac%7BMdu%7D%7Bu(u%20-%20M)%7D%20%3D%20-%20%5Calpha%20dt%20%5Cquad%20(1.1.6)

的形式,然后两边取定积分

M%20%5Cint_N%5E%7Bu(t)%7D%20%20%5Cfrac%7Bdu%7D%7Bu%20(u-M)%20%7D%20%3D%20-%5Calpha%5Cint_0%5Etdt%20%5Cquad%20(1.1.7)

利用积分公式

%5Cint%5Cfrac%7Bdu%7D%7B(u%2Ba)(u%2Bb)%7D%20%3D%20%5Cfrac%7B1%7D%7Bb-a%7Dln%5Cfrac%7Ba%2Bu%7D%7Bb%2Bu%7D%20(a%20%5Cnot%3D%20b)%20%5Cquad%20(1.1.8)

先看左边的部分:

M%5Cint%5Cfrac%7Bdu%7D%7B(u%2B0)(u%2B(-M))%7D%20%3DM%20%5Cfrac%7B1%7D%7B(-M)-0%7Dln%5Cfrac%7B0%2Bu%7D%7B(-M)%2Bu%7D

%3D%20-ln%5Cfrac%7Bu%7D%7B(u-M)%7D%0A%0A%3Dln%5Cfrac%7Bu-M%7D%7Bu%7D

M%20%5Cint_N%5E%7Bu(t)%7D%20%20%5Cfrac%7Bdu%7D%7Bu%20(u-M)%20%7D%20%3D%20ln%5Cfrac%7Bu-M%7D%7Bu%7D%20-%20ln%5Cfrac%7BN-M%7D%7BN%7D

%3D%20ln%5Cfrac%7Bu-M%7D%7Bu%7D%20%5Cfrac%7BN%7D%7BN-M%7D%20%3D%20ln%5Cfrac%7B(u-M)N%7D%7Bu(N-M)%7D

右边部分

-%5Calpha%5Cint_0%5Etdt%20%3D%20-%5Calpha%20t

左右两边结合

ln%5Cfrac%7B(u-M)N%7D%7Bu(N-M)%7D%20%3D%20-%5Calpha%20t

%5Cfrac%7B(u-M)N%7D%7Bu(N-M)%7D%20%3D%20exp(-%5Calpha%20t)

(u-M)N%20%3D%20u(N-M)%20exp(-%5Calpha%20t)

uN%20-%20MN%20%3D%20uNexp(-%5Calpha%20t)%20-%20uM%20exp(-%5Calpha%20t)

uN%20-%20uNexp(-%5Calpha%20t)%20%2B%20uM%20exp(-%5Calpha%20t)%20%3D%20%20MN

u%5BN%20-%20Nexp(-%5Calpha%20t)%20%2B%20M%20exp(-%5Calpha%20t)%5D%20%3D%20%20MN

u%20%3D%20%20%5Cfrac%7BMN%7D%7BN%20-%20Nexp(-%5Calpha%20t)%20%2B%20M%20exp(-%5Calpha%20t)%7D

u%20%3D%20%20%5Cfrac%7BMN%7D%7BN%5B1%20-%20exp(-%5Calpha%20t)%20%2B%20%5Cfrac%7BM%7D%7BN%7D%20exp(-%5Calpha%20t)%5D%7D

u%20%3D%20%20%5Cfrac%7BM%7D%7B1%20-%20exp(-%5Calpha%20t)%20%2B%20%5Cfrac%7BM%7D%7BN%7D%20exp(-%5Calpha%20t)%7D

u%20%3D%20%20%5Cfrac%7BM%7D%7B1%20%2B%20%5Cfrac%7BM%7D%7BN%7D%20exp(-%5Calpha%20t)%20-%20exp(-%5Calpha%20t)%20%7D

式(1.1.7)给出

u(t)%20%3D%20%5Cfrac%7BM%7D%7B1%2B(%5Cfrac%7BM%7D%7BN%7D%20-%201)%20exp(-%5Calpha%20t)%7D%20%5Cquad%20(1.1.9)

式(1.1.9)作为非线性方程(1.1.5)的解给出任意时刻病鼠的数目。我们看出,它正是众所周知的生物学中的“生长曲线”,如图1.1所示。


数学物理方法笔记2的评论 (共 条)

分享到微博请遵守国家法律