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

第五讲常微分方程数值解 化工学院软件应用教科组 2006-10 本章知识要点 数值计算 常微分方程初值问题 常微分方程边值问题 MATLAB 微分方程求解常微分方程的相关函数 ode45 ode23 bvp4c 微分方程在化工模型中的应用 间歇反应器的计算 活塞流反应器的计算 全混流反应器的动态模拟 定态一维热传导问题 逆流壁冷式固定床反应器一维模型 固定床反应器的分散模型 Matlab常微分方程求解问题分类 初值问题： 定解附加条件在自变量的一端 一般形式为： 初值问题的数值解法一般采用步进法，如Runge-Kutta法 边值问题： 在自变量两端均给定附加条件 一般形式： 边值问题可能有解、也可能无解，可能有唯一解、也可能有无数解 边值问题有3种基本解法 迭加法 打靶法 松弛法 Matlab求解常微分方程初值问题方法 将待求解转化为标准形式，并“翻译”成Matlab可以理解的语言，即编写odefile文件 选择合适的解算指令求解问题 根据求解问题的要求，设置解算指令的调用格式 Matlab求解初值问题函数 指令 含义 指令 含义 解 算 ode23 普通2-3阶法解ODE odefile ODE文件格式 ode45 普通4-5阶法解ODE 选 项 odeset 创建、更改ODE选项的设置 ode113 普通变阶法解ODE odeget 读取ODE选项的设置 ode23t 解