《数值分析》课件-第八章 非线性方程组的求解.pptVIP

设将原方程改写成下列形式 据此建立迭代公式 但若采用方程另一种等价形式 建立迭代公式 第二种方法 仍取迭代初值 ，则有 结果会越来越大，不可能趋于某个极限. 迭代法的收敛条件 定理 注 条件(2)可用更强更便于应用的条件代替: 由误差估计可以得到迭代终止的条件 由误差估计可以知道L越小，收敛越快以及迭代最 少次数 局部收敛性 前述定理条件为充分条件，非必要条件。在实际 应用当中 条件并不易检验。 定义 （局部收敛性） 设 有不动点 ， 如果存在 的某个邻域 对任意 ， 迭代产生的序列 且收敛到 ， 则称迭代法局部收敛. 定理 设s为? (x)的不动点, ? (x)在s的某个邻域内连续, 且|? (x*)| 1, 则迭代法是局部收敛的. 本定理给出的就是局部收敛性条件。具体解题时，虽然无法判别隔根区间是否为以x*为中心的邻域，但只要它足够小，且在邻域中满足： x y y = x s y=g(x) x0 p0 x1 p1 ? x y y = x s y=g(x) x0 p0 x1 p1 ? x y y = x s y=g(x) x0 p0 x1 p1 ? x y y = x s y=g(x) x0 p0 x1 p1 ? 收敛性图像说明 用不同方法求方程 的根 例 这里 其不动点为 由此构造不同的迭代法： 可改写为各种不同的等价形式 从计算结果看到迭代法（1）及（2）均不收敛，且它 们均不满足定理3中的局部收敛条件. 迭代法（3）和（4）均满足局部收敛条件，且迭代法（4）比（3）收敛快，因在迭代法（4）中 . function y=iterate(x) x1=g(x) ; n=1 ; while (abs(x1–x)= 1.0e –6 ) ( n=1000 ) x = x1 ; x1=g(x); n = n+1 ; end x1 n matlab program  收敛阶(描述收敛速度)  收敛速度（收敛速度的阶） 设迭代过程 收敛于方程 的根 ， 如果迭代误差 当 时成立下列 渐近关系式 则称{xn}是p阶收敛 ?若 p = 1 , 称{xk}为线性收敛, 这时 0 C ≤1。 ? p1, 称{xk}为超线性收敛; p=2, 称其为平方收敛. 收敛阶定理  数p的大小反映了迭代法的收敛速度的快慢，P越大，收 敛越快，所以说收敛阶是对迭代法收敛速度的一种度量。 （收敛阶定理） 对于迭代过程 ，如果 在所求根 的邻近连续，并且: 则该迭代过程在点 邻近是P阶收敛的。 上述定理说明，迭代过程的收敛速度依赖于迭代函数 的选取. Newton迭代法 Newton法的基本思想 将非线性方程线性化，以线性方程的解逐步逼近非线性方程的解 具体而言： 设xk是非线性方程 f(x)=0的一个近似根，把 f(x)在xk处作一阶泰勒展开，即用前两项近似代替 则近似方程转化为 设 ，上式解为 Newton迭代法 于是方程 f(x)=0的新的近似根xk+1，可得牛顿迭代公式 牛顿迭代公式为特殊的不动点迭代。 其迭代函数为 Newton迭代法几何解释 设 是根 的某个近似值，过曲线 上横坐标为 的点 引切线，并将该切线与 轴的交点的横坐标 作为 的新的近似值. 注意到切线方程为 牛顿法的收敛性  定理 设f(x*)=0, ,且在 x* 的邻域 上 存在, 连续,则可得 (1)Newton迭代公式在单根情况下至少2阶局部收敛. (2) 定理表明，当初值x0充分接近x*时，Newton法的 收敛速度较快，但当初值不够好时，可能会不收敛 或收敛于别的根，这可从Newton法的几何意义看到： 牛顿法的计算步骤： 步骤1 准备 选定初始近似值 ， 计算 步骤