数模学习笔记-微分方程模型
模型介绍
当我们描述实际对象的某些特性随时间(或空间)而演变的过程、分析它的变化规律、预测它的未来性态、研究它的控制手段时,通常要建立对象的动态微分方程模型。微分方程大多是物理或几何方面的典型问题,假设条件已给出,只需用数学符号将已知规律表示出来,即可列出方程。
微分方程在各个领域中都有着广泛应用,比较典型的有传染病的预测模型、经济增长预测模型、兰彻斯特(Lanchester)战争预测模型、人口的预测模型、药物在体内的分布与排除预测模型、烟雾的扩散与消失预测模型。
在高数中我们对微分方程有一定认识,但是把学习中心放在了复杂微分方程的求解上,但是忽略了如何建立微分方程的问题。其基本规律随着时间的增长趋势呈指数形式,根据变量的个数建立微分方程。微分方程模型的建立基于相关原理的因果预测法。
它们的不同之处在于建立微分方程的过程一般会涉及到变化率的概念(往往和时间相关)。例如,运动学中的速度是位移相对时间的变化率,加速度是速度相对时间的变化率,人口学中人口增长率是人口相对时间的变化率。
该方法的优点是短、中、长期的预测都合适,既能反映内部规律以及事物的内在关系,也能分析两个因素的相关关系,精度相应较高,另外对模型的改进也比较容易理解和实现。该方法的缺点是虽然反映的是内部规律,但由于方程的建立是以局部规律的独立性假定为基础,故做中长期预测时,偏差有点大,而且微分方程的解比较难以遇到。
建立微分方程
下面介绍一些常见的问题说明微分方程建模的基本方法:
eg1. 放射性物质的衰变:
物理学家发现,放射性物质衰变的速度和物质的量成正比。
eg2. 人口增长问题(生物数量增长问题):
观察统计发现,人口增长率和现有人口数成正比。
eg3. 战争问题:
交战双方的伤亡影像因素虽然比较复杂,但是大致与武器杀伤力、防御力、双方兵力对比因素相关。
eg4. 万有引力定理:
开普勒发现地球沿着以太阳为焦点的椭圆轨迹运动,运动速度满足单位时间内扫过面积相等。
eg5. 传染病问题:
医学经验表明,传染病有两个关键因素,易感人群、带病人群。而传染往往通过两种人群的接触导致,不仅与易感人群和带病人群有关,还与他们的接触有关。
微分方程的解法
解析方法:
对于一些简单的微分方程,可以借助一些数学技巧解出,例如我们在高数中学习到的:可分离变量的方程、齐次方程、一阶线性微分方程、二阶常系数微分方程等。
数值方法:
能通过解析方法求解的微分方程毕竟只是少数,对于从实际问题中提取的大量微分方程,数值方法是必不可少的手段。
微分方程的解包含两部分,一是满足方程的一个函数族(或者一条曲线族),我们称其为通解;二是满足特定条件的解,我们称其为特解。
解析方法
数值方法
对于大量的微分方程,只能得到其数值解。一般而言,得到的解是方程的一个特解的近似。求解微分方程数值解的方法很多,例如:欧拉法、龙格-库塔法等。其基本思想是通过已知点得到函数值,并用函数值代替一个小区间上函数的导数,得到在该区间上的一条直线,并用该直线作为方程特解的近似。
传染病问题的建模和分析
模型1:已感染人数i(t)
模型2:SI模型
模型3:SIS模型
模型4:SIR模型
经济增长模型
道格拉斯(Douglas)生产函数
资金与劳动力的最佳分配(静态模型)
经济(生产率)增长的条件(动态模型)
Bernoulli方程
人口增长预测与控制模型
人口发展方程
连续型人口发展方程
模型求解
讨论
兰彻斯特作战模型
参考文献
以上笔记内容基于学习:
[1] 数学建模算法与应用/司守奎, 孙兆亮主编. —2版. —北京:国防工业出版社,2015.4
[2] https://blog.csdn.net/weixin_44949135/article/details/107741248
仅供参考学习