化工计算方法化工计算方法绪论.pptVIP

化工计算方法 1 绪 论 1.1 什么是计算方法 在工程设计和科学研究中，经常会遇到一些复杂的数学计算问题，对这些计算，仅靠人工计算难以完成，必须借助于计算机求解 数学计算举例 数学计算举例2 例5： 化工原理中对流传热系数的计算公式 解析法求解数学方程的难易程度 计算方法 针对不同问题，将加减乘除等基本运算按一定的顺序和步骤构成适用于计算机的、计算时间较短、占用资源较少的完整的求解过程，称为计算方法，简称算法。 本课程包括的算法有： 非线性方程求根 插值法 曲线拟合 数值积分 常微分方程（组）初值问题求解 线性及非线性方程组求解 这些方法也称为数值计算方法，所得的解称为数值解。 算法中常用的基本方法 1.2 误差和有效数字 误差——近似值与真解（准确解）之差 误差可用于衡量数值解对准确解的“近似”程度 误差是怎样产生的？ 1.2.1 误差来源 (1) 有限过程代替无限过程的误差——截断误差 1.2 误差和有效数字2 (2) 由于计算机位数有限引起的误差——舍入误差 （与计算机和所用数学表达式有关） 例：圆周率 p =3.14159， 等 运算次数较少时舍入误差可以忽略不计，但大量的运算会造成舍入误差的传播和积累，有可能对求解结果造成严重影响 1.2.2 误差和误差限 相对误差与相对误差限 有效数字 有效数字 1.2.4 防止误差危害 有效数字（2） 复习思考题 * * 科学计算中涉及到的复杂计算 用数学解析方法无法求解或者求解困难的计算——如复杂函数积分和微分、微分方程求解、代数方程组求解等 需要反复多次进行或采用人工计算太复杂繁琐的计算——如迭代计算、插值计算等 大量的数据处理——数据筛选和处理、曲线拟合等 例1: (化工热力学）给定温度 T 和压力 P 下实际纯气体的压缩因子Z 的表达式为 Z=PV/RT 由R-K 方程式求 Z 式中 （常数取决于气体性质） 以上两例均为非线性方程求根问题 例2：求解方程 例3:（化学反应工程）某气相反应的管式反应器，管长与转化率的关系为： 已知管径为d=0.1m，求转化率为0.4时的管长是多少？ 例4：函数 f (x) 的形式未知，由实验知 x，f (x)数据 2.83 1.7 3.11 1.8 2.15 1.4 2.35 1.5 2.58 1.6 1.66 1.1 1.54 1.0 1.97 1.81 f (x) 1.3 1.2 x 以上两例积分，难以用解析法求解，需要用计算机求数值解，数值积分问题 。 求积分 由因次分析得 强制对流 自然对流 已知函数形式 问题： 如何从实验数据得到合适的公式？ 如何求参数的值？ 如何对所求公式进行检验？ 这类问题是数据处理及曲线拟合问题。 参数a,b,c,d由实验确定 此外，包括化工在内的工程计算还有函数插值、常微分方程、代数方程组等需要求数值解的问题。 易 几个 极易 一个 基本不可能 代数方程 多个 线 性 基本不可能 难 难 易 不可能 偏微分方程 非常难 常微分方程 非常难 几个 非常难 一个 不可能 代数方程 多个 非 线 性 不可能 不可能 不可能 非常难 不可能 偏微分方程 不可能 常微分方程 解析法求解线性问题，大部分不可能 解析法求解非线性问题，几乎完全不可能 ＃ 迭代法 离散化方法 逼近方法 用简单函数的值近似代替求解困难或形式未知的复杂函数的值，称为逼近函数的方法。 用差商代替导数、差分代替微分等，将连续性的数学问题转化为离散问题处理。 用一个固定公式反复计算，对较为粗糙的根的近似值进行加工直到满足精度要求的方法。 例：指数函数ex 展开为无穷级数，实际计算时，通常只截取前面有限项来计算 截断误差是评价一个计算方法优劣的重要指标 泰勒展开只取前面两项 在建立数学模型时，由于对实际问题进行简化造成的数学模型与实际问题的误差称为模型误差 用实验方法观测模型参数时的误差称为观测误差 模型误差和观测误差都是客观存在的 数值计算方法只涉及截断误差和舍入误差 绝对误差与绝对误差限 用 x* 和x分别表示准确值和近似值，绝对误差为 由于x*通常未知，无法求出e，只能估计其取值范围，即： e 称为近似值 x 的绝对误差限，反映了近似值 x 可信赖的程度，可写为 误差如何表示呢？ 绝对误差限的含义是，近似值x与准确值x*之间的误差不超过e 。 # 相对误差限：满足下式的 e r 是近似值 x 的相对误差限 近似值x的相对误差 (2) 相对误差与相对误差限 相对误差是误差与真值的比值，与绝对误差相比，相对误差更能准确地反映近似值的准确程度。 由于真值实际上无法得到，相对误差通常表示为 1