数值分析课件非线性方程组的求解.pptVIP

functiony=iterate(x)x1=g(x);n=1;while(abs(x1–x)=1.0e–6)(n=1000)x=x1;x1=g(x);n=n+1;endx1nmatlabprogram

第28页，共64页，星期日，2025年，2月5日收敛阶(描述收敛速度)

收敛速度（收敛速度的阶）设迭代过程收敛于方程的根，如果迭代误差当时成立下列渐近关系式则称{xn}是p阶收敛?若p=1,称{xk}为线性收敛,这时0C≤1。?p1,称{xk}为超线性收敛;p=2,称其为平方收敛.第29页，共64页，星期日，2025年，2月5日收敛阶定理

数p的大小反映了迭代法的收敛速度的快慢，P越大，收敛越快，所以说收敛阶是对迭代法收敛速度的一种度量。（收敛阶定理）对于迭代过程，如果在所求根的邻近连续，并且:则该迭代过程在点邻近是P阶收敛的。上述定理说明，迭代过程的收敛速度依赖于迭代函数的选取.第30页，共64页，星期日，2025年，2月5日Newton迭代法

Newton法的基本思想将非线性方程线性化，以线性方程的解逐步逼近非线性方程的解具体而言：设xk是非线性方程f(x)=0的一个近似根，把f(x)在xk处作一阶泰勒展开，即用前两项近似代替则近似方程转化为设，上式解为第31页，共64页，星期日，2025年，2月5日Newton迭代法

于是方程f(x)=0的新的近似根xk+1，可得牛顿迭代公式牛顿迭代公式为特殊的不动点迭代。其迭代函数为第32页，共64页，星期日，2025年，2月5日Newton迭代法几何解释

设是根的某个近似值，过曲线上横坐标为的点引切线，并将该切线与轴的交点的横坐标作为的新的近似值.注意到切线方程为第33页，共64页，星期日，2025年，2月5日牛顿法的收敛性

定理设f(x*)=0,,且在x*的邻域上存在,连续,则可得(1)Newton迭代公式在单根情况下至少2阶局部收敛.(2)定理表明，当初值x0充分接近x*时，Newton法的收敛速度较快，但当初值不够好时，可能会不收敛或收敛于别的根，这可从Newton法的几何意义看到：第34页，共64页，星期日，2025年，2月5日牛顿法的计算步骤：

步骤1准备选定初始近似值，计算步骤2迭代按公式迭代一次，得新的近似值，计算步骤3控制如果满足或，则终止迭代，以作为所求的根；否则转步骤4.此处是允许误差，而第35页，共64页，星期日，2025年，2月5日牛顿法的计算步骤：

其中是取绝对误差或相对误差的控制常数，步骤4修改或者，则方法失败；否则以代替转步骤2继续迭代.如果迭代次数达到预先指定的次数，第36页，共64页，星期日，2025年，2月5日functiony=newton(x0)x1=x0–fc(x0)/df(x0);n=1;while(abs(x1–x0)=1.0e–6)(n=100000000)x0=x1;x1=x0–fc(x0)/df(x0);n=n+1;endx1nmatlabprogram

第37页，共64页，星期日，2025年，2月5日计算重根的牛顿迭代法Newton法在重根情形下的收敛阶为线性阶若是方程的重根，即此时有第38页，共64页，星期日，2025年，2月5日牛顿法应用举例可导出求开方值的计算程序：开方公式对于任意给定的正初值均为平方收敛。求无理数的近似值