2.3重要：MATLAB常微分方程(组)的数值解法介绍.pptVIP

2.3 常微分方程（组）的数值解法;知识要点;微分方程在化工模型中的应用;Matlab常微分方程求解问题分类;Matlab求解常微分方程初值问题方法;Matlab求解初值问题函数;odefile;odefile的编写---常微分方程;odefile的编写---常微分方程组;高阶微分方程odefile的编写;解算指令的使用方法;常微分方程初值问题解算指令比较;ode解算指令的选择(1);ode45和ode23的比较-1;ode45和ode23的比较-2;ode解算指令的选择(2);ode解算指令的选择(2);function demo1 figure ode23s(@fun,[0,100],[0;1]) hold on, ode45(@fun,[0,100],[0;1]) %-------------------------------------------------------------------------- function f=fun(x,y) dy1dx = 0.04*(1-y(1))-(1-y(2)).*y(1)+0.0001*(1-y(2)).^2; dy2dx = -1e4*dy1dx + 3000*(1-y(2)).^2; f = [dy1dx; dy2dx];;间歇反应器中串联－平行复杂反应系统（例4-1）;三个串联的CSTR等温反应器（例4-3）;Matlab求解边值问题方法：bvp4c函数;Matlab求解边值问题的基本指令;初始解生成函数：bvpinit();sol＝bvp4c（odefun，bcfun，solinit，options，p1，p2，…） 输入参数： odefun是计算导数的M函数文件。该函数的基本形式为：dydx＝odefun(x,y,parameters,p1,p2，…)，在此，自变量x是标量，y，dydx是列向量。其基本形式与初值问题的odefile形式相同 bcfun是计算边界条件下残数的M函数文件。其基本形式为：res＝bcfun(ya,yb,parameters,p1,p2,……)，文件输入参量ya,yb是边界条件列向量,分别代表y在a和b处的值。res是边界条件满足时的残数列向量。注意：例如odefun函数的输入参量中包含若干“未知”和“已知”参数，那么不管在边界条件计算中是否用到，它们都应作为bcfun的输入参量。;sol＝bvp4c（odefun，bcfun，solinit，options，p1，p2，…） 输入参数： 输入参量options是用来改变bvp4c算法的控制参数的。在最基本用法中，它可以缺省，此时一般可以获得比较满意的边值问题解。如需更改可采用bvpset函数，使用方法同odeset函数。 输入参量p1，p2等表示希望向被解微分方程传递的已知参数。如果无须向微分方程传递参数，它们可以缺省。;sol＝bvp4c（odefun，bcfun，solinit，options，p1，p2，…） 输出参数： 输出变量sol是一个结构体 sol.x是指令bvp4c所采用的网格节点； sol.y是y(x)在sol.x网点上的近似解值； sol.yp是y(x)在sol.x网点上的近似解值； sol.parameters是微分方程所包含的未知参数的近似解值。;边值问题的求解-1;边值问题的求解-2;边值问题的求解-3;化工数学模型中的微分方程