17.（常）微分方程模型

一、概述

含导数或微分的方程就是微分方程

导数的最高阶就是几阶方程

特解：不含任意常数的

通解：含相互独立任意常数个的

初值条件：确定通解中任意常数的条件

二、matlab求微分方程解析解

dsolve求解析解函数

实际上大部分应用问题都只能求出数值解

solver（求解器）求解数值解

共7个函数（ode45和ode15s最常用）

90%都可以用ode45 龙格-库塔方法

什么叫刚性和非刚性？ 刚性的意思就是变化的非常剧烈的函数

一阶微分方程组比较好求，但是高阶微分方程组必须转化为一阶微分方程才能求数值解

具体方法就是新引入变量，类似于下图这样

例：求解导弹追踪问题

！！！！快速修改所有变量的方法：shift+enter

三、应用

1.人口模型预测

马尔萨斯模型

就是这样的模型

阻滞增长模型

就是把r看成x 的线性函数，r（x）是减函数

fplot是绘制表达式函数

其中xx是关于x 的函数

2.捕食者猎物模型

高中生物（）

种群相互竞争模型

σ注意一下，是两个种群之间消耗食物的倍数

种群相互依存模型

有三种形式

σ是提供养分相对于另一方消耗事物的倍数

出现这个黄字说明时间取的不合适

可以把时间区间改小一点(一个种群能独立生存一个不能独立生存）

食饵与捕食者模型

这些模型不要局限于生物学中

3.传染病模型

SI模型

一个matlab小技巧：一个变量如何在主函数和函数里使用：用全局变量

global +变量名先在主函数里定义 然后在函数里再次声明

只考虑易感者和感染者之间的转换

拓展

1.某种措施使β降低

2.考虑出生率和自然死亡率和疾病死亡率

SIS模型

感染者被治好了且恢复后还会被感染

SIR模型

有康复者，且不会感染

拓展：康复率变化

考虑疾病死亡率

SIRS模型

拓展：康复者一部分会成为易感者，另一部分免疫

SEIR模型

加入潜伏者E，潜伏者转换为感染者有一个速率

潜伏者有/无传染性