《化工数值计算》课件例题123.pdf

《化工数值计算》课件例题 《化工数值计算》课件例题 第一章 复习思考题 计算方法的主要研究对象是什么？ 非线性方程求根、插值法、曲线拟合、数值积分、常微分方程（组）、初值问题求解、线性非 及线性方程组求解 算法中常用的基本方法 1）离散化方法 用差商代替导数、差分代替微分等，将连续性的数学问题转化为离散问题 2 ）逼近方法 用简单函数的值近似代替求解困难或形式未知的复杂函数的值 3 ）迭代法 用一个固定公式反复计算，对较为粗糙的根的近似值进行加工直到满足精度要求的 方法。 误差为什么不可避免，用什么标准衡量近似值的准确程度？ 误差——近似值与真解（准确解）之差可用于衡量数值解对准确解的“近似”程度 误差来源 （数值计算方法只涉及截断误差和舍入误差） 1）有限过程代替无限过程的误差——截断误差； 2 ）由于计算机位数有限引起的误差——舍入误差 3 ）对实际问题进行简化造成的数学模型与实际问题的误差：模型误差 4 ）用实验方法观测模型参数的误差：观测误差 （都是客观存在的） 什么叫绝对误差、相对误差和有效数字？ 相对误差：误差与真值的比值，与绝对误差相比，相对误差更能准确地反映近似值的准确程度 有效数字：定义：若近似值 x 的绝对误差限是某位数字的半个单位，从该位到 x 的第一位非 零数字共有 n 位，就说 x 有 n 位有效数字。 根据上课PPT 整理，仅供学习交流 - 1 - 《化工数值计算》课件例题 第二章 非线性方程求根 两类方法： 1）根的搜索：二分法、逐步搜索法等 2 ）迭代法：直接迭代法、牛顿法(切线法) 、弦截法（割线法）。 3 -5 f (x ) x  x 1 0 例 题2-1 求方程 在区间（1，1.5）内的实根，要求精度为10 。 解：在给定区间内有一个实根，用二分法求解 (1) 输入求根区间[a,b]即[1,1.5]，计算f a 和f b ； (2) 取中点 c= (a+b) /2 （二分区间），计算f c ； (3) 判断是否有 f| c|  ，是则以 c 作为近似根，输出c，中止计算，否则进行下一步； (4) 判断是否有 f cf a0 ，是则令a=c，否则令b=c，回到步骤 (2) 。 计算结果为（MATLAB 程序及注释见教材 p.7 ） x=1.3247 (方程的根) m=16 (二分次数) 例 题 2-2 在298K 下，化学反应 2OF =O +2F 的平衡常数为0.410 atm，如在298K 下将OF 2 2 2 2 通入容器，当t=0 时为1 atm，问最后总压是多少？取计算精度为10-3 。 解：写出数学方程式 假设是理想气体，可知 2OF = O +2F 2 2 2 设氧的分压为 p ，平衡时有 1-2p p 2p