计算方法非线性方程求根.docVIP

PAGE PAGE 1 计算方法第7章 非线性方程求根 本章主要内容： 1.区间二分法. 2切线法. 3.弦位法. 4.一般迭代法. 重点、难点 一、区间二分法 区间二分法是求方程f(x)=0根的近似值的常用方法。 基本思想：利用有根区间的判别方法确定方程根的区间[a,b] ,将有根区间平分为二；再利用有根区间的判别方法判断那一个区间是有根区间；重复上述步骤，直到小区间端点差的绝对值小于等于精度要求的数值，则用将上一区间的分半值作为方程的根的近似值。 区间二分法的计算步骤如下： 计算区间端点的函数值f(a) , f(b)（不妨设f(a)＜0,f(b)＞0）； 确定初始有根区间[a,b]. 2.二分有根区间[a,b]，并计算 取 3.判断： 若，则方程的根为； 若 ，则有根区间为 ；令 若 ，则有根区间为；令 4. 如果│b-a│ε（ε为误差限），则方程的根为；否则转向步骤2，继续二分有根区间[a1,b1]，并计算中点值，继续有根区间的判断，直到满足精度要求为止，即│bn-an│ε 二分次数的确定：如果给定误差限ε，则需要二分的次数可由公式 确定应二分的次数。 例1 用区间二分法求方程在某区间内实根的近似值（精确到0.001） 【思路】参见上述区间二分法的计算步骤 解 ∵f(1.8)=-0.168＜0, f(1.9)=0.359＞0 ∴f(x)在区间[1.8 ，1.9]内有一个根。 由公式 取n=6， 计算结果列表如下： n an bn xn f(xn) 1 1.8 1.9 1.85 + 2 1.8 1.85 1.825 - 3 1.825 1.85 1.8375 + 4 1.825 1.8375 1.83125 - 5 1.83125 1.8375 1.834375 + 6 1.83125 1.834375 1.8328125 则方程在区间[1.8，1.9]内所求近似值为x*≈ x = 1.8328125 区间二分法的优点是计算程序简单，只要f(x)在区间[a,b]上连续，区间二分法就可使用，但区间二分法不能用来求偶次重根，由于区间二分法收敛比较慢，在实际计算中，区间二分法常用来求比较好的含根区间和初始近似值，以便进一步使用收敛更快的迭代法求出更精确的近似值。 迭代序列收敛阶的概念 设迭代序列收敛于，如果存在实数与正常数c，使得 ，则称序列是阶收敛于。 特别地，当时，称序列为线性（一次）收敛； 为线性收敛时，必须要求。 当时，称序列为平方（二次）收敛； 当时，称序列为超线性收敛； 收敛阶越大，则序列与的误差缩减越快，也就是序列收敛越快。 二、切线法(牛顿法) 1. 切线法的基本思想：假设方程f(x)=0在区间[a,b]上有唯一根x*,过曲线y= f(x)上的一点（x0,f(x0)），作曲线的切线，用此切线与x轴的交点的横坐标x1作为方程的根x*的新的近似值, 再过点（x1,f(x1)），作曲线的切线，则又得到新的近似值,按此方法进行迭代计算，直到满足精度要求为止。 切线法(牛顿法)的迭代公式为 2.切线法的收敛性 我们利用定理（7.1）来判断切线法的收敛性。定理（7.1）还给出了一个初始值x0的选择方法， 定理7.1. 设f(x)在[a,b]上存在二阶连续导数，且满足条件 ⑴ f(a)f(b) 0； ⑵f /(x) 在[a,b]上不等于零 （3）f //(x) 在[a,b]上不变号 则对任意初值x0∈[a,b] ，只要满足f(x0) f //(x)≥0. 则由切线法迭代公式得到的近似根序列平方收敛于方程f(x)=0在区间[a,b]的唯一根x*。 切线法的计算步骤：先判断有根区间[a,b]，然后选择初始值x0（一般地，若f//(x)0，则选择区间的右端点；若f//(x) 0，则选择区间的左端点），再建立迭代公式进行计算（列表计算）。 例2 用切线法求例1中方程在[1，2]内根的近似值，精确到0.000001 【思路】根据f(x0)f//(x)0在有根区间上选择初始值x0，代入迭代公式进行计算 解 计算得 n 0 1 2 3 4 xn 1 1.857142857 11.834243504 1.834243185 例3 证明 计算的切线法迭代公式为 （n=0,1,…） 解 因为计算等同于求方程的根，