Matlab与化工数值计算第5讲常微分方程数值解(.pptxVIP

Matlab与化工数值计算第5讲常微分方程数值解(2)

Contents本讲内容概览从经典方法到现代求解器，系统掌握常微分方程数值解法的核心脉络。01回顾与衔接：欧拉法的局限02高阶单步法：龙格-库塔方法族03线性多步法与预测-校正格式04刚性方程问题与隐式方法05MATLABODE求解器体系详解06化工过程综合案例实战

CHAPTER01回顾与衔接：欧拉法的局限从一阶精度出发，理解高阶方法的必要性

NUMERICALMETHODS欧拉法核心原理回顾欧拉法以一阶泰勒展开为基础，用差商替代导数实现离散递推。其公式简洁但精度仅为O(h)，在化工复杂非线性系统中累积误差显著，需更小步长补偿精度不足，计算效率低下。欧拉法切线逼近曲线几何原理示意01基本公式：y(n+1)=y(n)+h·f(xn,yn)，基于前向差商近似导数，几何意义为沿切线方向线性外推一个步长02误差阶数：局部截断误差O(h2)、整体截断误差O(h)，属于一阶方法，步长减半时整体误差仅减半03显式与隐式：显式欧拉计算简单但稳定性差，隐式欧拉需迭代求解但稳定性更好04化工应用局限：反应动力学方程含指数项和强非线性，一阶精度在合理步长下难以满足工程要求

NUMERICALMETHODS·ERRORANALYSIS化工实例：一级反应的欧拉法误差分析以一级不可逆反应dCA/dt=-kCA为典型场景，欧拉法在合理