数值计算课件第二章非线性方程的数值解法.ppt

牛顿法每一步要计算f和f’，相当于2个函数值，且有些导数难求。为了避开导数的计算，用差商代替导数。x0切线割线切线斜率?割线斜率2.4弦截法：x2x1第45页,共54页，星期六，2024年，5月用割线斜率（差商）替换切线斜率，代入牛顿法迭代公式：上式中，固定弦的一个端点（x0,f(x0)），而另一端点变动，称为单点弦法。2.4.1单点弦法：第46页,共54页，星期六，2024年，5月单点弦法几何意义第47页,共54页，星期六，2024年，5月因为f(x*)=0，故求导数得 所以0??’(x*)?1，所以单点弦法仅为线性收敛。单点弦法收敛速度：迭代函数：当初值x0充分接近时很接近f’(x*)第48页,共54页，星期六，2024年，5月2.4.2双点弦法：为了加快收敛速度，弦的两个端点都在变动，称为双点弦法或称快速弦截法。迭代时需要2个初值xk和xk－1。双点弦法迭代公式：。第49页,共54页，星期六，2024年，5月双点弦法几何意义P1P2第50页,共54页，星期六，2024年，5月双点弦法收敛速度：双点弦截法的收敛性与牛顿迭代法一样，即在根的某个邻域内，f(x)有直至二阶的连续导数，且f’(x*)?0，具有局部收敛性，同时在邻域任取初值x0、x1，迭代均收敛。可以证明，双点弦截法具有超线性敛速度，收敛的阶为：第51页,共54页，星期六，2024年，5月例4.4用Newton法和弦截法解下面方程的根，并比较解:由Newton法由弦截法第52页,共54页，星期六，2024年，5月x0=0.5;x1=0.4;x2=0.3430962343x3=0.3473897274x4=0.3472965093x5=0.3472963553x6=0.3472963553Newton法由弦截法要达到精度10-8弦截法迭代5次Newton迭代法迭代4次x0=0.5;x1=0.3333333333x2=0.3472222222x3=0.3472963532x4=0.3472963553第53页,共54页，星期六，2024年，5月2.5各种迭代方法的误差特性总结方法的误差特性即方法的收敛速度。二分法：线性收敛牛顿法：单根：平方收敛复根：线性收敛改进牛顿法：线性收敛或平方收敛单点弦截法：双点弦截法：超线性收敛线性收敛简单迭代法：视迭代函数不同而定第54页,共54页，星期六，2024年，5月***例题例2.2.1试用迭代法求方程在区间(1，2)内的实根。解：由建立迭代关系k=10,1,2,3…….计算结果如下:第13页,共54页，星期六，2024年，5月例题精确到小数点后五位第14页,共54页，星期六，2024年，5月例题但如果由建立迭代公式仍取，则有，显然结果越来越大，是发散序列第15页,共54页，星期六，2024年，5月第16页,共54页，星期六，2024年，5月xyy=xxyy=xxyy=xxyy=xx*x*x*x*y=φ(x)y=φ(x)y=φ(x)y=φ(x)x0p0x1p1?x0p0x1p1?x0p0x1p1?x0p0x1p1?简单迭代法第17页,共54页，星期六，2024年，5月收敛定理考虑方程x=φ(x),φ(x)在[a,b]上连续,若(I)对所有x?[a,b]，有φ(x)?[a,b]；(II)存在0?L1，使所有x?[a,b]有|φ’(x)|?L1。则：1）方程x=φ(x)在[a,b]上的解x*存在且唯一。2）任取x0?[a,b]，由迭代过程xk+1=φ(xk)收敛于x*?????简单迭代法推论验后误差估计：??误差估计式：验前误差估计：第18页,共54页，星期六，2024年，5月证明：①φ(x)在[a,b]上有根？令有根②根唯一？反证：若不然，设还有