微分方程模型

合集 - 数学建模算法(27)

1.整数规划2024-07-072.非线性规划2024-10-113.非线性规划之飞行管理问题2024-10-114.遗传算法与直接搜索2024-10-115.图论基础：最小路径问题2024-07-076.网络最大流量问题2024-07-077.最小费用最大流问题2024-07-058.旅行商（TSP）问题2024-07-05

9.微分方程模型2024-07-09

10.数据预处理2024-10-1111.插值方法2024-07-0612.拟合算法2024-07-0813.数据初始化2024-10-1114.层次分析法2024-07-1315.模糊综合评价2024-07-1416.CRITIC指标客观赋权方法2024-10-1117.Topsis评价法2024-10-1118.熵权法2024-10-1119.灰色关联度分析2024-10-1120.聚类分析2024-10-1121.主成分分析2024-10-1122.因子分析2024-10-1123.判别分析2024-10-1124.典型相关分析2024-10-1125.对应分析2024-10-1126.多维标度法2024-10-1127.亮度变化与空间滤波2024-10-11

收起

matlab求解微分方程

%dsolve('equations','initial conditions')
%使用
syms y(x); %定义y是关于x的函数
eq = y-diff(y,x) == 2*x %将方程付给eq
dsolve(eq)%求通解
dsolve(eq,y(0) == 3)%求特解

对多元微分方程

$$\{\begin{array}{l}\frac{dx}{dt}=y\\ \frac{dy}{dt}=-x\end{array}$$

syms x(t),y(t);
eq1 = diff(x,t)==y;
eq2 = diff(y,t)==-x;
eqs = [eq1 eq2];
conds = [x(0)==0 y(0)==0];
[xsol,ysol] = dsolve(eqs, conds);

对于大多数微分方程，很难求解解析解，一般求数值解。

数值解一般是通过将函数离散化,例如$\frac{df}{dx}=\frac{f(n+1)-f(n)}{\mathrm{\Delta }x}$,然后将$\mathrm{\Delta }x\mathrm{理}\mathrm{解}\mathrm{为}\mathrm{步}\mathrm{长}\mathrm{通}\mathrm{过}\mathrm{迭}\mathrm{代}\mathrm{求}\mathrm{解}.$

ODE系类函数（Ordinary Differential Equation）

求数值解常用ode系列函数（常微分方程）

ode函数分为两类，求解刚性和非刚性问题

两类问题的分类大体上是根据变化的陡峭程度来分的。

非刚性ode求解命令
求解器solver	功能	说明
ode45	一步算法：4、5阶龙格库塔方程：累计截断误差(Δx)^5	大部分尝试的首选算法
ode23	一步算法：2、3阶龙格库塔方程：累计截断误差(Δx)^3	适用于精度较低的情形
ode113	多步算法：Adams	计算时间比ode45短

刚性ode求解命令
求解器solver	功能	说明
ode23t	梯形算法	适度刚性情形
ode15s	多步法：Gear’s反向数值微分：精度中等	若ode45失效时，可以尝试使用
ode23s	一步法：2阶Rosebrock算法：精度低	当精度较低时，计算时间比ode15s短
ode23tb	梯形算法：精度低	当精度较低时，计算时间比ode15s短

函数使用见常微分方程的函数：常微分方程 - MATLAB & Simulink - MathWorks 中国