化工计算方法化工计算方法方程求根.pptVIP

2 非线性方程求根 2.1 二分法(一种简单而直观的搜索求根方法 ) 设方程f?x?=0在区间 [a,b]内有单实根， 反复对分，可得到一系列有根区间 多个实根的情况怎么办？ 2.1.2 应用举例例 题2-1 求方程 在区间（1，1.5）内的实根，要求精度为10-5 。 解：本题在给定区间内有一个实根，用二分法求解步骤如下 例题 2-2 在298K下，化学反应 2OF2=O2+2F2 的平衡常数为0.410 atm，如在298K下将OF2通入容器，当t=0时为1atm，问最后总压是多少？取计算精度为10－3。 解：首先写出求解问题的数学方程式 假设是理想气体，由反应的化学计量式可知 2OF2= O2+2F2 例题 2-2 在298K下，化学反应 2OF2=O2+2F2 的平衡常数为0.410 atm，如在298K下将OF2通入容器，当t=0时为1atm，问最后总压是多少？取计算精度为10－3。 例题 2-3 求非线性方程 sin(x)=0 在区间（-4，7）内的全部实根，要求精度为10-5。 （取 h = 0.2） 例题 2-4 在在容积VR=10m3的全混反应器中进行绝热一级不可逆反应。已知数据如下 反应放热速率物料温度变化速率 2.2 直接迭代法（简单迭代） 2.2 直接迭代法（简单迭代） 迭代过程收敛情况 迭代：用固定公式逐次逼近方程的根 直接迭代法(2) 例题2-5 用迭代法求以下方程在1.5 附近的 根，要求精度10-5。 例题 2-6 在9.33 atm，300.2 K时，容器中充以2mol 氨气，试求该容器的容积，要求精度为10－2。已知氨气的范德华常数 。 解: 范德华方程 2.3 牛顿法——恒收敛且收敛速度快的迭代方法 如果 x0 为 f（x）=0 的一个近似根，函数 f（x）在点 x0 附近可展开成泰勒级数 牛顿迭代过程的几何解释 例题 2-7 用牛顿法求方程 f(x) = ex - 3x= 0，在 0.5 附近的根，要求精度为10-5。 例题 2-8 简单蒸馏时，某时刻釜底残液量与低沸点组分组成x之间有如下关系式 例题 2-9 对多元体系中的组分i，其饱和蒸气压pi ，可用安托因(Antoine)公式表示（式中t为泡点） 据热力学知识，多元体系中，与液相摩尔分率xi平衡的气相分率 yi 在泡点时，有 2.4 弦 截 法 牛顿法缺点：对复杂函数，无法计算导数时不便使用。 避开导数计算的方法：差商代替导数 弦截法的几何意义 例题 2-10 用弦截法求方程 * * 例：化工热力学中压缩因子Z的计算式 (关于Z的3次多项式) 超越方程 形式为 ，解x? 称为方程的根，或称函数f?x?的零点 如何找到方程的根 x? ？ 两类方法 根的搜索：二分法、逐步搜索法等 迭代法：直接迭代法、牛顿法(切线法)、弦截法（割线法）。 计算机只能找到近似根，本章介绍用以上主要算法找近似根的原理及应用。 # 0 y x a b f(a) f(b) x* 必有f ?a? f ?b? 0 二分过程进行到什么时候为止？ x0 步骤： (1)从x0=1/2(a+b)，对分区间 (2) 检查 f?x0?与f?a? 的符号是否相同 相同时令 a1= x0, b1= b； 不相同时令 a1= a, b1= x0 ∴区间为(a1, b1) f(x0) a1 b1 x1 （3）重复⑴步，直到结束 (a，b)，(a1，b1)，(a2，b2)，?，(ak，bk)，长度依次减半 搜索法求根的基本思想：通过逐步缩小有根区间找到近似根 k→∞时，二分过程的极限必为方程的根 x? 。 有根区间长度为（a k+1，b k+1 ）时，a k与b k的中点是 x k=1/2（a k + b k） 以x k为近似根，则（b k+1 - a k+1 ）是 x 的绝对误差限 有根区间(ak, bk) 的长度 计算时预先给出允许误差 ?，区间长度（b－a） ?时， 近似根 x 满足f (x)=0，此时x与 x?的误差不超过允许误差? |x －x?| ? 计算精度 划分小区间，对各个小区间实施二分求根 y x a b x1 x2 x3 x4 x5 将[a,b]分割成小区间，保证每一区间只有一个实根 从a开始搜索根。先检查小区间内是否有根，无根转向下一区间；有根则用二分法求出 依次对每一区间进行搜索和检查，直至求出全部实根。 h 步