清风数学建模——模型学习微分方程模型

微分方程及微分方程模型的基础概念

Matlab求微分方程的解析解

Matlab求一阶微分方程的数值解

Matlab求高阶微分方程的数值解

例题2：刚性问题

例题3：导弹追击问题

微分方程模型之：人口预测模型

微分方程模型之：捕食者-猎物模型

微分方程模型之：人口预测模型微分方程模型之：种群相互竞争模型、相互依存模型

微分方程模型之：传染病SI模型

非原创，看清风数学建模视频课做的笔记，建议去B站观看清风数学建模的课程或者购买正课~

基础概念：
常微分方程：y’+4y-2=0；导数的方程，一个变量
偏微分方程：偏导数的方程，多个变量

Matlab求微分方程的解析解：dsolve(‘方程1’,‘方程2’，…，‘方程n’,‘初始条件’,‘自变量’)；一阶微分方程；含未知系数的微分方程；高阶微分方程及初始条件；微分方程组

1、Matlab在表示微分方程时，用字母D表示微分，D2,D3分别表示二阶、三阶微分，后面跟的是要求解的因变量；
2、必须要指定自变量，要不然Matlab默认自变量是t，若实际自变量不是t，则计算结果是不正确的；
3、初始值可以不给，不给求出来的就是通解
4、乘号* 不能省略
%% 求解y-y’=2x
clear;clc
dsolve(‘y-Dy=2*x’,‘x’)
% 新版本的Matlab可能会淘汰上面的这种写法，但目前就用这种写法吧

%% 求解含有未知系数的微分方程 y-y’=ax
clear;clc
dsolve(‘y-Dy=a*x’,‘x’)

%% 求解二阶微分方程，已知初始条件 y’‘+4y’+29y=0 且 y(0)=0 y’(0)=15
clear;clc
dsolve(‘D2y+4Dy+29y=0’,‘y(0)=0,Dy(0)=15’,‘x’)

%% 求解微分方程组

clear;clc
[x,y,z]=dsolve(‘Dx=2x-3y+3z+t’,'Dy=4x-5y+3z+t’,‘Dz=4x-4y+2*z+t’,‘t’)

Matlab求一阶微分方程的数值解：
微分方程的数值解函数如下： [x,y]=solver(‘f’,ts,x0,options)
1、x代表自变量； y代表函数值； solver代表求解函数，总共有七种，如下图所示。ts=[t0,tf]代表自变量的初值和终值；x0表示函数的初始值； f是待求解的微分方程需要用函数文件编写，且要写成标准形式，如下图所示
并用@进行调用，或者用单引号；options = odeset(‘reltol’,1e-4,‘abstol’,1e-8) 用来设置相对误差和绝对误差，相对误差是reltol，绝对误差是abstol
2、使用Matlab求解数值解时，高阶微分方程组必须转换为一阶微分方程组；
3、优先使用ode45;解不