1、第三章 非线性方程的数值解法,科学技术研究与工程实践中，经常会遇到求解非线性方 程的问题，一般归纳为求解方程其中 是一元非线性实函数。根据 是代数多项式还 是超越函数（指数函数、对数函数、三角函数等），方程分 别称为代数方程或超越方程。 次代数方程： 超越方程：,第三章 非线性方程的数值解法,对于代数方程，根的数目与方程的次数相同，不高于4 次的代数方程已有求根公式，高于4次的代数方程没有精确 的求根公式；超越方程则复杂得多，如果有解，可能是一个 或多个，或无穷多个，没有精确的求解公式。因此，需要研 究如何采用一定的数值算法求得满足一定精度要求的近似根。 内容提要：隔根区间的确定、二分法、迭代
2、法、牛顿迭代法、弦截法 重点内容：迭代法、牛顿迭代法,第三章 非线性方程的数值解法,数值方法求方程 根的近似值，要解决三个问题： 1、根的存在性：方程有没有根，如果有根，有几个； 2、根的隔离：找出有根区间，把有根区间分成较小的子区间，每个子区间只有一个根（隔根区间）； 3、根的精确化：确定了隔根区间后，可以用各种方法将某一近似根 逐步精确化。按照一定的方法产生一个序列，此序列在一定的条件下收敛于方程的根 。产生序列 的不同方法就构成了不同的方程求根方法。,第三章 非线性方程的数值解法,应用举例：本征声线求解 本征声线：从声源出发经过一定的传播路径到达接收点的声 线，接收点处的声场是所有本征声
3、线能量叠加的结果。 折射定律(Snell)：,第三章 非线性方程的数值解法,对于给定的海洋环境，每条声线的轨迹由声线的起始掠射角 唯一决定。设声源位于 处，接收点位于 处，声 线轨迹方程：需要导出声线以掠射角 从声源出发，经过一定水平距离 到达深度 的计算公式，有改变 ，给定精度要求，进行求解。,第三章 非线性方程的数值解法,1 根的搜索与二分法,根的搜索 定义（隔根区间）：如果在区间 内只有方程 的一个根，则称区间 为隔根区间。隔根区间的确定 描图法 逐步搜索法,介值定理： 定理1：设函数 在区间 上连续，且 ， 则方程 在 上至少有一个根。 定理2：设函数 在区间 上是单调连续函数，且，则
4、方程 在 有且仅有一个根。,1 根的搜索与二分法,1 根的搜索与二分法,描图法：画出 的简图，由曲线与 轴的交点位置 确定出隔根区间。或者将方程等价变形为 ，画 出函数 和 的简图，从两条曲线交点的横坐 标位置确定隔根区间。 例：求方程 的隔根区间。 解：由图可知，方程仅有一个实根， 隔根区间为 。,1 根的搜索与二分法,逐步搜索法：首先确定方程 的实根所在区间 再按照选定的步长 （n为正整数），逐点计算处的函数值 ，当 与 异号时，则 即为方程的一个隔根区间。 对于 次代数方程其根的绝对值的上下界有如下结论： （1）如果 ，则方程的根的绝对值小于； （2）如果 ，则方程根的绝对值大于 。,1
5、 根的搜索与二分法,例：求方程 的隔根区间。 解：利用逐步搜索法设方程的根为 ，则 即根的所在区间为 和 。 取 计算 ，隔根区间为,1 根的搜索与二分法,二分法 基本思想：通过计算隔根区间的中点，逐步缩小隔根区间， 从而得到方程的近似根序列 。 过程：设 为连续函数，方程 的隔根区间为 ，则函数值 异号。 （1）将区间 二分得到中点 ，计算函数值 ； （2）根据计算得到的函数值进行判断，如果 ，则 方程的根为 ；否则判断 的符号，如果其与异号，则隔根区间变为 ；如果其与 异号，则 隔根区间变为 ；,1 根的搜索与二分法,（3）把新的隔根区间记为 ，继续按照上述规则对隔根 区间 进行二分。满足
6、精度要求则停止。通过不断的二分，就得到一系列隔根区间有 。由于后一区间的长度都是前一 区间长度的一半，所以 的长度为 当 时，区间 的长度趋于零，即区间最终会收敛于 方程 的解 。实际计算中，只要二分的次数足够多 ，可取最后区间中点 作为方程根的近似值，即,1 根的搜索与二分法,误差：如果事先给定精度要求为 ，则满足时，停止计算。,1 根的搜索与二分法,例：用二分法求方程 在 内的根的近似 值，要求绝对误差不超过 。 解： 即 严格单调增加，又 ，所以方程在 上有 唯一实根。 令 ，得到 ，取 ，即至少二分7次 。计算过程如下：,1 根的搜索与二分法,function f=f3(x) f=x3
7、+4*x*x-10;,a0=1;b0=2;k=0; while 1 k=k+1;x=(a0+b0)/2; if f3(x)*f3(a0)0b0=x; else a0=x; end if (b0-a0)/20.005break;end end x=(b0+a0)/2;,1 根的搜索与二分法,所以最后得,1 根的搜索与二分法,例：用二分法求方程 在 内的根，要求精确到小数点后第三位小数，需要二分多少次？ 解：设 ，由于 ， 所以在区间 内方程 有唯一实根。令 ，求得所需对分次数至少是10次。特点：运算简单，方法可靠，对函数只要求在区间上连续 ；但收敛速度慢，不能用来求复数根及偶数重根。常用于为 其
8、它求根方法提供较好的近似初始值。,2 迭代法及其迭代收敛的加速方法,迭代法（逐次逼近） 基本思想：利用某种递推算式，使某个预知的近似根（初 值）逐次精确化，直到得到满足精度要求的近似根。 算法：给定方程 ，其中 在有根区间 上连 续。设 是方程的一个近似根，将方程 改写为 ，代入初值 构造迭代序列,2 迭代法及其迭代收敛的加速方法,这种求方程近似根的方法称为迭代法（或简单迭代法）， 称为迭代函数。如果由迭代法产生的迭代序列 极限存在 ，即 ，则称 收敛，否则称 发散。 称为迭 代函数 的不动点，或是方程 的解。由 转化为 时，迭代函数 不是唯一的，不同，会产生不同的序列 ，从而收敛情况也不 一
9、样。,2 迭代法及其迭代收敛的加速方法,几何意义： 求方程 的根 ，在几何上就是求直线 与曲线交点 的横坐标，如图所示。从图中可以看出， 当迭代函数 的导数 在根 处满足不同条件时，迭 代过程的收敛情况也有所不同。所以迭代过程的收敛依赖 于迭代函数 的构造，为使迭代法有效，必须保证其收 敛性。,2 迭代法及其迭代收敛的加速方法,2 迭代法及其迭代收敛的加速方法,2 迭代法及其迭代收敛的加速方法,例：已知方程 在 内有一个根，用两种不同的迭代公式 （1） ；（2） 进行迭代，观察 所得序列的收敛性。 解：计算结果如表所示,2 迭代法及其迭代收敛的加速方法,x0=0.3;k=0;while 1k=
10、k+1;xk=10x0-2; %公式1%xk=log10(x0+2) %公式2if abs(xk-x0)=0.0005break;endx0=xk;end,2 迭代法及其迭代收敛的加速方法,例：对于方程 ，试采用如下迭代法求方程的根 (1) ；(2) ；(3) 选取初值 。 解：计算结果见下表,x0=1.5;k=0; while 1 k=k+1; xk=x06-1; %公式1，发散%xk=(x0+1)(1/6);%公式2，收敛%xk=1/(x05-1);%公式3，发散if abs(xk-x0)=0.0005break;endx0=xk; end,2 迭代法及其迭代收敛的加速方法,2 迭代法及其
11、迭代收敛的加速方法,迭代法收敛定理：设方程 ，若迭代函数 在有根 区间 上满足： （1）当 时， ； （2） 在 上可导，且有 ， ； 则有 （1）方程 在 上有唯一的根 ； （2）对任意初值 ，迭代公式 产生的数列 收敛于方程的唯一根 ，即 ； （3）误差估计,2 迭代法及其迭代收敛的加速方法,拉格朗日中值定理：如果函数 在闭区间 上 连续，在开区间 内可导，则至少存在一点 使得,2 迭代法及其迭代收敛的加速方法,证明： （1）根的存在性：令 ，则 在 上连续 ，由条件1有 由连续函数的介值定理，必存在 使得 ， 即为方程的根。 根的唯一性：假设另有 也是方程的根， 则有两式相减并应用微分中
12、值定理得,2 迭代法及其迭代收敛的加速方法,式中 在 和 之间，所以 ，由条件（2）有上式不成立，所以有 。 （2）收敛性：取 时，有 ，由微分中值定理可知，在 和 之间存在 使得： 反复应用上述关系可得由条件（2） ，所以 即,2 迭代法及其迭代收敛的加速方法,（3）误差估计：由微分中值定理，在 和 之间存在 使得 同样得到 同时即 因为 所以 反复应用 ，最后得到,2 迭代法及其迭代收敛的加速方法,迭代次数估计： 给定精度要求 ，即 ， 只要从而迭代次数满足,2 迭代法及其迭代收敛的加速方法,注意： 由式 可以看出，当 越小时，序列 收敛越快。只要相邻两次迭代的偏差 足够小，就可以保证近似
13、解 有足够的精度。所以，常采用 条件 来控制迭代次数。,2 迭代法及其迭代收敛的加速方法,例：求方程 的最大根，要求精度 。 解：（1）求隔根区间 方程变为 ， 做 和 的图形，由图可知方程 的最大根在区间 内。,2 迭代法及其迭代收敛的加速方法,2）建立迭代公式，判断收敛性。 方程变形为 ，迭代公式为在区间 内可导，并且 单调递增。 又 ，所以，当 时， 由于 ，所以 迭代法收敛。,2 迭代法及其迭代收敛的加速方法,（3）迭代计算过程如下表所示。所以方程的最大根,x0=3.5;k=0; while 1k=k+1;xk=(log10(x0)+7)/2;if abs(xk-x0)=0.0001b
14、reak;endx0=xk; end,2 迭代法及其迭代收敛的加速方法,2 迭代法及其迭代收敛的加速方法,定理（局部收敛性定理）： 设 是方程 的根， 在 的某一邻域连续，且，则必存在 的一个邻域 ，对任 意选取的初值 ，迭代公式 产生 的数列 收敛于方程的根 ，称迭代法在 的邻域 具有 局部收敛性。,注意：在实际应用中 事先不知道，因此条件无法验证。但是如果已知根的初始值 在 根 的附近，又根据 的连续性，可采用条件来代替 。,2 迭代法及其迭代收敛的加速方法,2 迭代法及其迭代收敛的加速方法,例：用迭代法求方程 在隔根区间 内的根 ，要求精确到小数点后第4位。 解： (1)构造迭代公式：，
15、迭代公式为 (2)判断收敛性： 由局部收敛性定理 ， 所以迭代法收敛,2 迭代法及其迭代收敛的加速方法,(3)迭代计算过程如下表所示：,2 迭代法及其迭代收敛的加速方法,x0=1.5; k=0; while 1k=k+1;xk=(x0*x0+1)(1/3);if abs(xk-x0)=0.00005break;endx0=xk; end,2 迭代法及其迭代收敛的加速方法,用迭代法求方程根的近似值的计算步骤如下： (1)准备：选定初值 ，确定方程 的等价方程 (2)迭代：按迭代公式 计算出 (3)判别：如果 ，则终止迭代计算，取 作为方 程的近似根。否则，转到第(2)步继续迭代计算。,2 迭代法
16、及其迭代收敛的加速方法,迭代收敛的加速方法 迭代-加速公式:记 ，则由微分中值定理有其中 在 和 之间。 设 在根 附近变化不大，又设 ，由迭代收敛条 件 ，有 上式整理为,2 迭代法及其迭代收敛的加速方法,上式说明 作为根的近似值时绝对误差大致为 如果把该误差作为一种补偿，可以得到更好的近似值记得到迭代-加速公式,2 迭代法及其迭代收敛的加速方法,例:用迭代-加速公式求 在区间 内的根。 解: ，加速-迭代公式为,2 迭代法及其迭代收敛的加速方法,x0=1.5; k=0; while 1k=k+1;xk1=(x0*x0+1)(1/3);xk=20*xk1/11-9*x0/11;if abs(
17、xk-x0)=0.00005break;endx0=xk; end,3 牛顿（Newton）迭代法,应用范围： 用于求解高次代数方程、超越方程等非线性方程；可求方程 的实根、复根，也可求单根和重根。 基本思想： 将非线性方程 逐步转化为线性方程来进行求解。 牛顿迭代法的最大优点是在方程的单根附近具有较高的收敛 速度，是一种将近似根精确化的相当有效的迭代法。,3 牛顿（Newton）迭代法,对于非线性方程 ，设 连续可微，将 在 处 泰勒（Taylor）展开如果 ，则可在点 附近取线性部分近似替代 ， 得到 的近似方程进行变换，得,3 牛顿（Newton）迭代法,由迭代法思想，上式写为迭代公式的
18、一般形式称上式为牛顿（Newton）迭代公式，其迭代函数为,3 牛顿（Newton）迭代法,几何意义：,3 牛顿（Newton）迭代法,如图所示，假设 是非线性方程 在隔根区间 内 的根， 在 内可导，且对于 有 。任取初值 ，过曲线 上的点 作切线以其作为曲线 的近似表达式。切线交 轴于 ，则以 作为方程根的第一次近似值。同理，过曲线 上点做切线,3 牛顿（Newton）迭代法,切线交 轴于 ，则 然后以 作为方程根的第二次近似值。 反复作切线，第k+1条切线方程为它与 轴的交点横坐标为由此得到方程的近似根数列 。 几何意义即依次用切线代替曲线，用线性函数 的零点作 为函数零点的近似值，牛顿
19、迭代法又称为切线法。,3 牛顿（Newton）迭代法,初始值 的选取： 定理（局部收敛性定理）：设 是方程 的根，若 （1）函数 在 的邻域内具有连续的二阶导数； （2）在 的邻域内 ， 则存在 的某个邻域 ，对于任意的初值 ，由牛顿迭代法产生的数列收敛于根 。,3 牛顿（Newton）迭代法,证明：迭代函数 的导数为由条件(1)和(2)可知， 在 的邻域内可导。又 ， 所以在 的某个邻域 内恒有,3 牛顿（Newton）迭代法,定理（非局部收敛定理）：设 是方程 在隔根区间 内的根，如果 （1）对于 连续且不变号； （2）选取初始值 ，使 ， 则由牛顿迭代公式产生的数列收敛于根 。如果在 上
20、 的符号不容易判断，则利用下式 确定初值即如果在 处， 满足上式，并且 ，则大多数情 况下可以保证牛顿迭代法是收敛的。,3 牛顿（Newton）迭代法,计算步骤： 1. 选定初值 ，构造迭代公式 ； 2. 迭代计算，得到 ； 3. 迭代终止判断，由迭代终止控制条件 控制迭代计算是否终止； 4. 迭代计算终止，得到方程根的近似值。,3 牛顿（Newton）迭代法,例：用牛顿迭代法求方程 在隔根区间 内的根，要求精确到小数点后第4位小数。 解: (1)令 ，牛顿迭代公式为(2)判断收敛性所以取 时满足 ，牛顿迭代法收敛。,3 牛顿（Newton）迭代法,(3)计算过程如下表,x0=1.5; k=0
21、; while 1k=k+1;xk=(2*x03-x02+1)/(3*x02-2*x0);if abs(xk-x0)=0.00005break;endx0=xk; end,3 牛顿（Newton）迭代法,3 牛顿（Newton）迭代法,例：用牛顿迭代法求方程 在 附近的根，精确到小数点后第4位。 解：所以由此 成立，所以可取 为初始值。 迭代公式为：,3 牛顿（Newton）迭代法,迭代计算过程如下表所示,x0=-1; k=0; while 1k=k+1;xk=(40*x041+2*x03-1)/(41*x040+3*x02);if abs(xk-x0)=0.00005break;endx0=
22、xk; end,3 牛顿（Newton）迭代法,3 牛顿（Newton）迭代法,例：用牛顿迭代法建立计算 近似值的迭代公式。 解：令 ，则以上问题化为求方程 的正根，牛顿迭代公式为继续证明对于任意初值 ，上面迭代公式恒收敛。 证明：对于任意 有即序列 有下界；,3 牛顿（Newton）迭代法,又即 单调下降，所以迭代公式收敛。计算 的近似值。的值在1011之间，取初值 ，迭代公式为,3 牛顿（Newton）迭代法,x0=10; C=115; k=0; while 1k=k+1;xk=(x0+C/x0)/2err=abs(xk-x0)if err=0.000005break;endx0=xk;
23、end,3 牛顿（Newton）迭代法,迭代计算过程如下（ ）：,迭代过程的收敛速度,定义：如果从任何可取得的初始值出发，由迭代法产生的迭代 序列都收敛于方程的根，则称该迭代法是大范围收敛；如果选取的初始值充分接近于方程的根，或者说存在根 的某一邻域，使得在其中任意选取初始值，迭代序列都收敛 于方程的根，则称该迭代法为局部收敛。一种迭代法具有使用价值，不但需要是收敛的，还要求 收敛的比较快。所谓迭代过程的收敛速度，是指在接近收敛 时迭代误差的下降速度。,迭代过程的收敛速度,定义：设迭代法 产生的数列 收敛于 的根 ，令误差 ，如果存在某个实数 以及 正常数 ，使则称数列 是 阶收敛的，相应的迭
24、代法是 阶方法，常 数 称为渐进误差常数。,迭代过程的收敛速度,可以看出， 和 为同阶无穷小量，即以 为基本 无穷小量时， 为 阶无穷小量，阶数 越高，收敛速度 越快。收敛速度是误差的收缩率，阶数越高，误差下降的越 快。当 时，称数列 线性收敛；当 时 ，称数列 为平方收敛（或二阶收敛）；当 时，称数 列为超线性收敛。显然， 越大，数列收敛的越快，所以迭 代法的收敛阶是对迭代法收敛速度的一种度量。,迭代过程的收敛速度,证明简单迭代法是线性收敛： 设 是方程 的根， 在 的某一邻域连续，在邻 域内有 ，且 ， 则简单迭代法线性收敛。 证明：设迭代过程 收敛于 ，则式中 在 和 之间。令 ，则 于
25、是迭代过程一阶收敛（线性收敛）。,迭代过程的收敛速度,证明牛顿迭代法平方收敛： 设 是 在隔根区间 内的根，如果（1）对于连续且不变号；（2）选取初始值，使 ，则牛顿迭代法平方收敛。 证明：将 在 处做泰勒展开，有式中 介于 之间。整理上式得,迭代过程的收敛速度,即两边取极限，并利用得所以牛顿迭代法平方收敛。,4 弦截法,如果函数 比较复杂，求导可能比较困难，有时甚至导函 数不存在。尤其是当 很小时，计算中容易产生很大的 误差。这时可将牛顿迭代公式中的 用差商来代， 即就得到求解非线性方程的弦截法。,4 弦截法,4 弦截法,设方程 的一个隔根区间为 ，如图所示。连结曲 线 上的两点 得到弦AB
26、，令，则弦AB的方程为则弦AB与 轴交点的横坐标 为以 作为方程根 的一个近似值。,4 弦截法,然后过曲线 上的两点 做弦，得到 其与 轴交点的横坐标 为于是把 作为方程根 的一个新的近似值。 反复做弦，得到迭代公式的一般形式为利用上述公式求方程 根的近似值的方法就叫做 弦截法。,4 弦截法,几何意义：依次用弦来代替曲线，用线性函数的零点作为函数 零点的近似值。 定理（弦截法收敛定理）：如果 在根 的某个邻域内有二阶连续导数，且对任意 ，有 ，则当 邻域充分小时，对 邻域内任意的 ，由弦截法迭代公式得到的近似值序列 收敛到方程 的根 。并可证明弦截法是按阶收敛的。,4 弦截法,计算步骤： （1
27、）准备：选定初始近似值 ，并计算 ； （2）迭代计算：依迭代公式计算得到 ，并计算 ； （3）控制：如果 或 ，则终止迭代。否则 以 分别代替 ，转（2）继续迭代计算。 弦截法和牛顿法比较：都是先把 线性化，但方式不同， 牛顿法用切线方程来近似代替非线性方程，而弦截法是用弦 线来近似代替非线性方程。,4 弦截法,例：用弦截法求方程 在区间 内根的近似值。精确到小数点后第4位。 解：令 ，取 ，弦截法迭代公式为,4 弦截法,迭代计算过程如下表所示：,4 弦截法,x0=0.5; x1=1; k=0; while 1k=k+1;xk=x1-(exp(2*x1)+x1-4)*(x1-x0)/(exp(
28、2*x1)+x1-exp(2*x0)-x0);if abs(xk-x0)=0.00005break;endx0=x1;x1=xk; end,课后题,1、二分法求方程 在区间 内的根，要求误 差不超过 。 2、已知方程 在区间 内有一根，试问用二分 法求根，使其具有5位有效数字至少应二分多少次？ 3、判断下列方程由几个根，求出隔根区间，并写出收敛的 迭代公式。 （1） （2） 5、用迭代法求的 正根，要求准确到小数点 后第5位。,课后题,4、方程 在1.5附近有根，把方程写成4种不同的等价形式，并建立相应的迭代公式：判断每种迭代公式产生的数列在1.5附近的收敛性，并用（1） 的迭代公式求方程的根
29、，要求误差不超过 。,课后题,6、用迭代法求方程 的根，要求准确到小数点 后第4位。 7、用迭代-加速公式求方程 在 附近的根，要 求准确到小数点后第4位。 9、应用牛顿法于方程 和 ， 分别导出求 的迭代公式，并求 10、用牛顿法求方程 在 附近的根，要 求准确到小数点后第3位。,课后题,11、证明计算 的牛顿迭代公式为取 ，计算 ，要求 。 12、用弦截法求 的根，取 ，计 算到 为止。,1、解：令 ，则 所以 单调增； 所以 在区间 内单调减。又 ，所以方程有唯一实根。 有,课后题,2、解：根的区间为 ，近似值具有5位有效数字，则精确 到小数点后4为小数，即取 。 4、解： （1） ，所
30、以迭代公式收敛。,课后题,（2） ，收敛（3）（4） 利用公式（1）计算方程的根，列表如下,课后题,5、解：方程变形为 ，分别画出 和 的图形，可知交点在 的范围；取 ，则 ，可知交点在 的范围，所以隔根区间可取为 。 令 ， ，从而 单调减； 单调增； ， 满足条件1； ，满足条件2， 所以迭代公式收敛。,课后题,迭代公式为： 计算过程如下：,课后题,6、解：（1）求隔根区间 方程变形为 ，做 和 的图形。由 图可知，方程的根在区间 内。 （2）建立迭代公式，判断其收敛性 将方程等价变形为 ，迭代函数 ，迭代公式为 。,课后题,因为 ，当 时， 所以 单调递减； 又 ，所以 ； ，所以 ；
31、由定理可知，迭代公式收敛。 （3）计算：取 ，列表计算,课后题,7、解： 取 ，迭代-加速公式为：计算过程如下：,课后题,9、证明： （1）对于方程 ， ，牛顿迭代公 式为所以,课后题,（2）对于方程 ， ，牛顿迭 代公式为所以,课后题,10、解：令 ， ， 牛顿迭代公式为 牛顿迭代公式为满足局部收敛性定理，计算过程如下表所示,课后题,11、证明：由 ， ， 牛顿迭代公式为取 ，计算过程如下：,课后题,课后题,12、解：令 ， ，满足收敛 条件，取 ，弦截法迭代公式为计算过程如下表所示,第三章 非线性方程的数值解法,主要内容： 1.确定隔根区间； 2.二分法； 3.简单迭代法； 4.迭代收敛的
32、加速方法； 5.牛顿迭代法； 6.弦截法； 7.迭代法的收敛阶。,第三章 非线性方程的数值解法,二分法：,迭代法收敛定理：设方程 ，若迭代函数 在有根 区间 上满足： （1）当 时， ； （2） 在 上可导，且有 ， ； 则有 （1）方程 在 上有唯一的根 ； （2）对任意初值 ，迭代公式 产生的数列 收敛于方程的唯一根 ，即 ； （3）误差估计,第三章 非线性方程的数值解法,第三章 非线性方程的数值解法,定理（局部收敛性定理）： 设 是方程 的根， 在 的某一邻域连续，且，则必存在 的一个邻域 ，对任 意选取的初值 ，迭代公式 产生 的数列 收敛于方程的根 ，称迭代法在 的邻域 具有 局部收
33、敛性。,第三章 非线性方程的数值解法,迭代-加速公式：牛顿迭代法,第三章 非线性方程的数值解法,初始值 的选取： 定理（局部收敛性定理）：设 是方程 的根，若 （1）函数 在 的邻域内具有连续的二阶导数； （2）在 的邻域内 ， 则存在 的某个邻域 ，对于任意的初值 ，由牛顿迭代法产生的数列收敛于根 。,第三章 非线性方程的数值解法,定理（非局部收敛定理）：设 是方程 在隔根区间 内的根，如果 （1）对于 连续且不变号； （2）选取初始值 ，使 ， 则由牛顿迭代公式产生的数列收敛于根 。如果在 上 的符号不容易判断，则利用下式 确定初值即如果在 处， 满足上式，并且 ，则大多数情 况下可以保证
34、牛顿迭代法是收敛的。,第三章 非线性方程的数值解法,二分法： 由公式有取 。即采用二分法求解，得到满足精度要求的根需 要二分13次。,第三章 非线性方程的数值解法,简单迭代法： 方程 变形为而 即迭代函数单调减，又 ；所以 单调减， ，所以 所以在区间 上迭代法收敛，取初始值 ；,第三章 非线性方程的数值解法,迭代计算过程如下表所示：,第三章 非线性方程的数值解法,迭代-加速公式：取 ，公式为取初始值 ，,第三章 非线性方程的数值解法,迭代计算过程如下表所示：,第三章 非线性方程的数值解法,牛顿迭代法：当 时， 不变号； 又 ，所以取 ； 迭代公式为,第三章 非线性方程的数值解法,迭代计算过程如下表所示：,总结： 1、对非线性方程 求根，首选要对 的性态及根的 近似位置有一个大致了解，使用较小的有根区间把方程的根 分离出来； 2、简单迭代法是一种逐次逼近的方法，需要判断收敛性， 了解收敛速度，可以采用加速公式； 3、牛顿迭代法是最常用的一种迭代方法，具有二阶收敛速 度，但是对初值的选取比较苛刻，需要靠近方程的根，否则 可能不收敛； 4、弦截法是牛顿迭代法的一种修改，虽然收敛速度较牛顿 法慢，但不需要计算导数，使用更加方便。同样要求初始值 的选取比较靠近方程的根。,第三章 非线性方程的数值解法,
