非线性方程求解.pptVIP

第1页，共32页，星期日，2025年，2月5日4.1化工实际问题的提出求解非线性方程是化工设计及模拟计算中必须解决的一个问题。与线性方程相比，非线性方程问题无论是从理论上还是从计算公式上，都要比线性方程复杂的多。对于一般的非线性f(x)=0，计算方程的根既无一定章程可循也无直接方法可言。例如，求解高次方程组7x6-x3+x-1.5=0的根，求解含有指数和正弦函数的超越方程ex-sin(x)=0的零点。解非线性方程或非线性方程组也是计算方法中的一个主题。一般地，我们用符号f(x)来表示方程左端的函数，方程的一般形式表示为f(x)=0，方程的解称为方程的根或函数的零点。通常，非线性方程的根不止一个，而任何一种方法只能算出一个根。因此，在求解非线性方程时，要给定初始值或求解范围。而对于具体的化工问题，初值和求解范围常常可根据具体的化工知识来决定。常见的雷诺数和摩擦系数关系方程在雷诺数低于4000时有以下关系式：(4-1)4.14.84.74.54.34.2总目录4.94.64.4第2页，共32页，星期日，2025年，2月5日4.1化工实际问题的提出这是一个典型的非线性方程。我们在管路设计中经常碰到。当我们已知雷诺数Re，如何根据公式（4-1）求出摩擦系数λ，这是我们在管路设计中必须首先解决的问题。对于方程（4-1）而言，无法用解析的方法求出摩擦系数，只能用数值求解的方法。如用在下面即将介绍的松弛迭代法，假设：则利用松弛迭代公式可得：经11次迭代可得摩擦系数为0.07593。同样，在n个组分的等温闪蒸计算中，通过物料和相平衡计算，我们可得到如下非线性方程：4.14.84.74.54.34.2总目录4.94.64.4第3页，共32页，星期日，2025年，2月5日4.1化工实际问题的提出在方程（4-3）中只有α是未知数，ki为相平衡常数，zi为进料组分的摩尔浓度，均为已知数。和上面的情况一样，方程（4-3）也无法直接解析求解，必须利用数值的方法，借助于计算机方可精确的计算。对于这个问题的求解，可利用我们下面介绍的牛顿迭代法进行计算，也可利用其他迭代公式进行计算，如采用牛顿迭代公式，则可以得到如下的具体迭代公式：（4-4）饱和蒸气压是我们经常要用到的数据，虽然我们可以通过实验测量来获取饱和蒸气压的数据，但我们通常利用前人已经测量得到的数据或回归的公式来获取，这可以减轻我们大量的基础实验工作。公式（4-5）是一种常用的饱和蒸气压计算公式：(4-5)其中p为饱和蒸气压，单位为mmHg，T为温度，单位为K，A、B、C、D为已知系数。要想得到某一温度下的饱和蒸气压，直接利用公式（4-5）是无法得到的。因为公式（4-5）两边都有未知变量，并且无法用解析的方法求解，必须用数值计算的方法求解。通过上面的一些例子，我们可以发现，如果没有适当的手段和办法来求解非线性方程，那么化学化工中的许多研究、设计等工作将无法展开，这势必影响化学化工的发展，下面我们将介绍一些实用的非线性方程求解方法，并提供计算机程序。4.14.84.74.54.34.2总目录4.94.64.4第4页，共32页，星期日，2025年，2月5日4.2实根的对分法4.2.1使用对分法的条件4.2.2对分法求根算法4.2.3对分法VB程序清单4.14.84.74.54.34.2总目录4.94.64.4第5页，共32页，星期日，2025年，2月5日4.2.1使用对分法的条件对分法或称二分法是求方程近似解