第4章非線性方程求根.docxVIP

第4章 非线性方程求根问题的引入我们知道，在多项式方程中，求根公式有一、二、三、四次方程，当n大于等于已经证明不能用公式计算，因此需要寻找另一种计算方法；同时在工程和科学技术中许多问题常常归结为求解非线性方程式的问题，非线性方程的解法也需要给出一种方法，本章来讨论这个问题。例1 关于真实气体的状态方程为 如果p与T都已知，则求体积V的方程为 这是一个非线性方程，如何求解呢？通常，非线性方程的根不止一个，对于非线性方程一般用对分法与迭代法求解。在用迭代法时，要给定初始值或求解范围。4.1 实根的对分法设有非线性方程 为[a,b]上的连续函数，且（不妨设方程只有一个实根），二分法叙述如下：第1步：，如果 则根一定在区间。于是我们得到长度缩小一半的含根区间，即 设已经完成了第1，第2，……，第k-1步，得到分半计算的含根区间 且满足： （1） （2），现我们看第k步： 否则有根区间为 且有 一般这样的过程不会一直下去，通常预先给定一个精度ε，使得 ，取对数计算得，来控制二分的次数。 例： 用二分法小数点后第3位（即要求解：显然，由，算得，k=11。计算结果如下表。 1 1.0 2.0 1.5 8.890625 2 1.0 1.5 1.25 1.564697 3 1.0 1.25 1.125 -0.097713 4 1.125 1.25 1.1875 0.616653 5 1.125 1.1875 1.15625 0.233269 6 1.125 1.15625 1.140625 0.0615778 7 1.125 1.140625 1.132813 -0.0195756 8 1.132813 1.140625 1.136719 0.0206190 9 1.132813 1.136719 1.134766 0.0004307 10 1.132813 1.134766 1.133789 -011 1.133789 1.134766 1.134277 -0.0045915注意：若函数在给定区间上有几个零点时，对分法只能算出其中一个，且是实根。对分法求根算法计算的一般步骤如下：输入有根区间[a，b]和误差控制量，定义函数。 计算中点以及的值分情况处理 停止计算，转向步骤4 修正区间 修正区间 得根输出近似根。X2X0X1YX 在算法中，常用代替的判断，以避免数值的溢出。4.2 迭 代 法 对给定的方程，将它转换成等价形式：。给定初始值，由此来构造迭代序列，如果迭代收敛，即有，则就是方程的根。在计算中当小于给定的精度控制量时，取为方程的根。例： 求方程的等价形式解： （a）(b) 注意：在方程转化为不同的等价形式后，会有不同的迭代函数，从而有不同的迭代序列，由此产生的迭代效果也可能不同：例： 对上例中的方程，考察用迭代法求根： （a）(b) 解：取初始值，分别代入迭代公式 计算结果如图。0 1.0 1.01 1.341471 0.5235992 1.473820 0.0236013 1.495301 -0.4965554 1.497152 -1.4877615 1.4972896 1.4973007 1.497300 -0计算看出，我们选取的两个迭代函数构造的序列的收敛情况不一样，一个收敛到1.497300，另一个的值无意义。问题在于迭代函数的选取上。因此，我们需要解决两个问题：如何选取迭代函数，使迭代过程收敛。若收敛较慢，怎样加速收敛