Matlab与化工数值计算-第5讲常微分方程数值解22409.pptVIP

第五讲常微分方程数值解 化工学院软件应用教科组 2006-10 本章知识要点 Matlab常微分方程求解问题分类 Matlab求解初值问题函数 odefile 所谓的odefile实际上是一个Matlab函数文件，一般作为整个求解程序的一个子函数，表示ode求解问题 Matlab提供了odefile的模板，采用type odefile命令显示其详细内容，然后将其复制到脚本编辑窗口，在合适的位置填入所需内容 一般而言，对于程序通用性要求不高的场合，只需将原有模型写成标准形式，然后“翻译”成Matlab语言即可 odefile的编写规定 常微分方程组odefile的编写 高阶微分方程odefile的编写 解算指令的输出控制 解算指令的options选项 边值问题求解指令：bvp4c () 边值问题求解指令：bvp4c () 常微分方程边值问题的应用 已知在球形氧化铝催化剂进行环己烷脱氢反应，催化剂直径为5mm，操作温度为700K，在此温度下，k＝4s-1，D=0.05cm2/s。计算催化剂颗粒内环己烷的浓度分布及催化剂的有效因子。 边值问题的求解 function Cha5demo8 lmb=15; solinit = bvpinit(linspace(0,pi,10),[1;1],lmb); opts=bvpset(Stats,on); sol = bvp4c(@ODEfun,@