数学建模（5）微分方程的建立与求解

数学建模（5）微分方程的建立与求解

特点： 常与优化问题相结合
可以用微分方程建模的情况： 涉及改变、变化、增加、衰变、边际、追赶的问题；考虑变化率的问题
一般形式分为隐式和显式，区别类似于函数中的函数与隐函数的关系，隐式微分方程中所求的第n项一般会出现与等式两边。
显式第n项的结果由之前的第1步，第2步，…第n-1步逐步迭代推导出，但迭代量不能太大，易发散
隐式第n项无法利用之前的第1步，第2步，…第n-1步逐步迭代推导出，必须一步求解

分类

常微分方程
ODE，是未知函数中只含有一个自变量的微分方程

偏微分方程
PDE，是未知函数是多元函数的微分方程

基本思想

一定要有背景
利用数学物理力学等学科中的定理或实验规律建立微分方程

微元法：
利用已知定理与规律寻找与微元之间的关系式，对象是微元而不是函数与导数

步骤：
1.微元法建立微分方程
2.确定定解条件（初边值条件）0
3.求解：解析解——数值解——图形解
4.讨论分析

解析解有限，一般是利用数值方法来求近似解

注：

打公式的时候dx中的d一定要用印刷体，即d要立着，不可以使用斜体。其中x，t可以使用斜体，但是d必须使用印刷体，要正直

Eular折线法

核心思想：差商代替微商
原理近似于泰勒公式：分割求解区间-差商代替微商

迭代思想

向前欧拉迭代法

向后欧拉迭代法

梯形公式

改进欧拉迭代法

其中后项欧拉算法和梯形公式为隐形算法，向前欧拉公式和改进的欧拉算法为显性算法。欧拉算法易计算，精度低，梯形公式经度高，不易求解。综合两者得到改进欧拉算式，是主要的求解方法。

实现代码可见连接：
https://blog.csdn.net/susanliuliu28/article/details/78804759

Runge_kutta法

用的较多的为四阶Runge_kutta法
算法精度高对应ode工具箱为ode45

MATLAB求解微分方程问题

请进传送门：

matlab求解微分问题 https://blog.csdn.net/qq_43649786/article/details/98636095