创建微分方程模型思路与方法

一、 寻找变化率具线性关系的等式

基于数学方法进行分类的模型中,微分方程模型是非常重要的一类(其他的还有几何模型、图论模型、概率模型、最优控制模型、规划论模型、排队论模型和马氏链模型等)。但是怎样建立,却并不那么自明。回忆一下初中解应用题的时候,我们建立方程都是找变量间的一个等式。根据这个思路,微分方程的等式可以这样寻找。

因变量的变化率 = 自变量的线性关系

用形式化语言可以抽象表达为:


图片.png

方程左边使用的是增量符号Δ,这是便于建模时好理解(一般情况下都是如此,对这个方法熟悉后可以直接使用“率”来理解),毕竟有


图片.png

X 可以是一个常量,也可以是一个变量,还可以是一个向量,花体字 ℒ 表示一个线性关系。

这就是我们要建立微分方程的一个基本思路和出发点。以下通过两个🌰来说明怎样寻找等式关系。

二、 速度和人口的两个栗子

(一) 匀速直线运动

假设位移为 s ,时间为 t ,匀速直线运动就是相同的时间内位移相同。也就是
图片.png

取极限后得到

图片.png

这是一个非常简单的可分离变量的微分方程,其通解为 s = kt + C 。C是任意常数,由于 t = 0 时,s = 0 ,这是初始条件,带入上式得到特解 s = kt ,如果把k换成v,就是我们熟悉的 s = vt

你可能会说,这么简单的公式我早就会了,还需要你来说?且慢,一个可以从另外一个角度去看看我们习以为常的一些事情有别的理解。另外,也是为理解以下更复杂模型打一个基础。

对于没有阻力的自由落体运动,该怎样求其s和t之间的关系呢?

提示:根据实验,相同时间内,速度的变化量相同。于是,方程可以建立为:

图片.png

解出上述公式后,将 s = vt 带入,就能得到s 和 t 的平方之间的线性关系模型。

(二) 人口模型

1. 人口增长率为常数

基本假设 : 人口自然增长率(出生率-死亡率) r 是常数, x(t) ~时刻t的人口,于是可以得到下列等式。
(t 时刻 + Δt时间)后人口变化数量 = t 时刻人口数量人口增长率Δt
即:x(t+Δt)-x(t) = x(t)rΔt, 移项, [x(t+Δt)-x(t)]/ Δt = x(t)r ,

图片.png

X(0) = x_0 ,方程求解后得到
图片.png

②就是马尔萨斯在1798年提出的著名的指数增长模型,其实,只要是可以表示为等比数列的关系,都可以通过这个方式建立微分方程。由于资源有限,环境承载能力有上限,故以上模型在短期、人口压力不大的时候可以使用,但是长期一定有一个增长的极限,资源、环境等外部因素对人口增长有阻滞作用人口,反应到增长率是一个随时间变化的量。

2. 人口增长率为变量

当到了某个临界点后,人口数量越多,各种资源越来越难以满足人的需求,出生率一定会下降,死亡率假定不变情况下,人口自然增长率r 是x 的减函数,为了计算简单,我们还是假定 r 和 x 存在线性关系的减函数。
假设: r(x) = r – sx …③ (r,s >0) r~固有增长率(x很小时)
x_M ~人口容量(资源、环境能容纳的最大人口数量),此时,人口不再增加,即 r(x_M) = 0, 于是 s = r/ x_M. , 带入③式得到


图片.png

将④带入①得


图片.png

求解⑤得到
图片.png

图形如下
图片.png

图 1 阻滞模型
⑥式为著名的logistic 模型,广泛应用于生物、科技、社会等各个领域,这个模型在动力学系统中也叫寻的模型。

三、 一种流行病模型简单理解

按照单一群体和复合群体两种基于群体进行的测算,使用的都是如下几类基本模型,具备S(易感者)、E(潜伏者)、I(感染者)和R(抵抗者)的全部或者部分对象人群,以这四类人群为基础,可以选择全部或者部分构成不同的模型类型,如适合普通流感的SI模型、适合急性传染病的SIR模型等。

(一) SEIR模型的基本构成

单一群体模型采用宏观视角建模,关注整个人群状态的变化,一般采用微分方程描述。最流行的单一群体模型是仓室模型,所有处于相同状态的人构成一个仓室,随着状态的变化,人员在仓室之间移动,仓室模型的基本假设是每个人都在相同人群中均匀混合,接触是瞬时的,并与历史无关。每个仓室的人口数量足够大,在传染病流行过程中,感染率和恢复率是常数。具体的抽象模型见下图:


图片.png

图 2 传染病一般发病形式对应的专业路线
表 1 符号说明


图片.png

(一) SEIR模型理解

SEIR模型由⑦、⑧、⑨、⑩四个方程构成的一个方程组和⑪这样一个约束条件组成。方程组中,左边表示不同人群的变化率。∂是求偏导数符号,表示变化量的符号∆趋于0时,就使用∂来表示。如:抵抗者增加的数量/天可以表示为 ΔR/Δt 。当 ΔR→0 与 Δt→0 时,就得到了 ∂R/∂t, 也就是抵抗者增加率。下面分别对四个方程进行说明。


图片.png

⑦式:易感者S增加的数量/天= -感染概率β(易感者接触感染者后转变为潜伏者的概率)×当天易感者S数量×感染者占比(感染者I/总人数N)。SI/N 表示易感者S接触的人中恰好是感染者I(假设接触的人中只有感染者I具有传染性,其他三类S、E、R都不具备传染性)概率。负号表示易感者每天减少数量。


图片.png

⑧式:潜伏者E增加的数量/天=当天由易感者转化为潜伏者的数量-转化率σ×当天潜伏者E数量。


图片.png

⑨式:感染者I增加的数量/天=转化率σ×当天潜伏者数量-治愈率γ×当天感染者数量


图片.png

⑩式:抵抗者R增加的数量/天=治愈率γ×当天感染者数量I
以上四个方程结合图2传染病一般发病形式对应的专业路线就更清晰了。
s.t N=S+E+I+R……⑪
只有极少部分(偏)微分方程有解析解,SEIR模型也不能得到解析解,通过数值算法求解。具体为用差商代替导数,将微分方程离散化后成为差分方程,利用4阶龙格-库塔(R-K)方法插值求解差分方程,这个过程就是 RK 算法。在R中编程计算还是很方便的,绘图效果也不错。

四、 微分方程模型常见应用场景

(一) 静态优化模型

存储模型,价格模型,消费者均衡模型、万有引力定律的发现

(二) 动态模型

人口模型、传染病模型、经济增长模型、药物在体内的分布与排除、烟雾的扩散与消失

 描述对象特征随时间(空间)的演变过程
 分析对象特征的变化规律
 预报对象特征的未来性质
 研究控制对象特征的手段

(三) 稳态模型

种群的竞争模型,军备竞赛模型、正规战与游击战、香烟过滤嘴的作用

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 213,928评论 6 493
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 91,192评论 3 387
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 159,468评论 0 349
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 57,186评论 1 286
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 66,295评论 6 386
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 50,374评论 1 292
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 39,403评论 3 412
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 38,186评论 0 269
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 44,610评论 1 306
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 36,906评论 2 328
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 39,075评论 1 341
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 34,755评论 4 337
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 40,393评论 3 320
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 31,079评论 0 21
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,313评论 1 267
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 46,934评论 2 365
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 43,963评论 2 351

推荐阅读更多精彩内容