matlab之常微分方程(ODE)求解

问题描述：已知理想全混釜的初始进料条件，并知道各物料的动力学，求解足够长时间后全混釜中各物质的浓度

常微分方程：只包含一个自变量的微分方程是常微分方程（Ordinary differential equations,ODE）

其求解分为初值问题和边值问题。

dsolve：以直观的方式求取解析解，便于理解

Y = dsolve(‘eq1,eq2,…’ , ’cond1,cond2,…’ , ’Name’)

采用dsolve表达理想全混釜反应器中的反应过程，其中一种物质的量关系如下：

-dC/dt=f(C,P01,P02,M1,M2) C=C0-∫Rc dt’

具体求解过程参考了：

关于求解微分方程——初学Matlab里的 ODE求解器_HFUTdumpling的博客-CSDN博客_matlab默认求解器

常微分方程

Y = dsolve(‘eq1,eq2,…’ , ’cond1,cond2,…’ , ’Name’)

其中，‘eq1,eq2,…’:表示微分方程或微分方程组;

        ’cond1,cond2,…’:表示初始条件或边界条件;

        ‘Name’:表示变量。没有指定变量时，matlab默认的变量为t；

dsolve 求解微分方程

求解微分方程：

在命令行输入: dsolve('Dy=3*x^2','x')

计算用例：自变量是t

-dC/dt=f(C,P01,P02,M1,M2) C=C0-∫Rc dt

描述第一个式子代码如下，第二个式子实际上包含在边界条件里

dsolve('DC = KC','C(0)=0.1','t')

常