化工计算机数据与图形处理PPT.pptVIP

2.4单变量求解法第29页，共56页，星期日，2025年，2月5日2.4单变量求解法3）最大误差设置从上例可以看出，所得解x＝0.0222并不精确，需要重新设定最大误差：打开“工具”－“选项”－“重新计算”。在该对话框的“最大误差”里填入想要的设置，如1.00E-12。第30页，共56页，星期日，2025年，2月5日2.4单变量求解法3）最大误差设置第31页，共56页，星期日，2025年，2月5日2.4单变量求解法3）最大误差设置分别将最大误差设置为1.00E-5,1.00E-9,1.00E-12,1.00E-15,1.00E-99，用上述单变量求解方法求得相应的x和y值。第32页，共56页，星期日，2025年，2月5日2.5求微分方程数值解1.原理化学反应动力学方程一般用一阶微分方程表示，其通式为：[a,b]是自变量的定义域，f(x,y)为已知函数。第33页，共56页，星期日，2025年，2月5日2.5求微分方程数值解1.原理若一函数y=F(x)代入微分方程，使得式(1)成立，即:则该函数就是微分方程的解。第34页，共56页，星期日，2025年，2月5日2.5求微分方程数值解1.原理解微分方程用到积分，得到的解析式是含一个任意常数C的通解。若有初始条件（初值）：y0＝F(x0)则通解的常数C可以确定，得到特解。第35页，共56页，星期日，2025年，2月5日2.5求微分方程数值解1.原理例如对于一级反应：A→B，其速率方程为：[A]t是反应物A在时间t时的浓度，其通解为：第36页，共56页，星期日，2025年，2月5日2.5求微分方程数值解1.原理若已知t＝0时A的浓度为[A]0，则在初始条件下的特解为：第37页，共56页，星期日，2025年，2月5日2.5求微分方程数值解1.原理并非所有微分方程都有解析解。事实上除了一些简单的基元反应，大多数反应动力学难以得到解析解或解析式很复杂，甚至不存在解析解。此时必须求助于数值解。另一方面反应动力学关心的问题是在t时刻反应体系内各物质的浓度[A]t,[B]t,……，有足够精度的近似解即可。第38页，共56页，星期日，2025年，2月5日2.5求微分方程数值解1.原理数值解在化工“三传一反”的模型化中都广泛应用。数值解有普适性，可用于复杂的微分方程体系及任何初始条件。第39页，共56页，星期日，2025年，2月5日2.5求微分方程数值解2.Euler法常微分方程数值解采用离散方法，即找出一种有效的数值计算方法，计算自变量的离散点：x0，x1，x2，…，xn以及对应的y近似值y0，y1，y2，…，yn。最简单的是Euler法，通常取等间距的x值：x1-x0=x2-x1=…=xn-xn-1=hh称为步长。第40页，共56页，星期日，2025年，2月5日第1页，共56页，星期日，2025年，2月5日本章主要内容2.1EXCEL的迭代用法2.2逐步逼近法解方程2.3Newton-Raphson法解方程2.4单变量求解解方程2.5求微分方程数值解注：每一种方法讲完，会找同学上台操作。第2页，共56页，星期日，2025年，2月5日2.1EXCEL的迭代用法EXCEL迭代：A2=F(A1)或B1=F(A1)F(A1)为含引用单元格A1的公式，地址为相对引用。向下或向右拖曳填充柄时，EXCEL将自动更新每一步。即：A3=F(A2),A4=F(A3),……C1=F(B1),D1=F(C1),……一般是向下拖曳。第3页，共56页，星期日，2025年，2月5日2.1EXCEL的迭代用法2.循环引用公式引用自己所在的单元格：A1=F(A1)无论是直接引用还是间接引用。只要打开的工作簿中有一个包含循环引用，Excel都将无法自动计算所有打开的工作簿。第4页，共56页，星期日，2025年，2月5日2.1EXCEL的迭代用法2.循环引用但有些化工计算需要用到循环引用，计算此类公式时，EXCEL必须使用前一次迭代的结果计算循环引用中的每个单