第四部分解非线性方程的迭代法学习课件.pptVIP

* 第4章 解非线性方程的迭代法本章讨论求非线性方程x04.1 的根的问题.其中?x是高次多项式函数或超越函数.如x3x5-2x4+8x2-7x+1xe2x+1-xlnsinx-2 等等. §1二 分 法设?x在区间[a,b]上连续且?a?b0,根据连续函数的介值定理,区间[a,b]上必有方程?x0的根,称[a,b]为方程?x0的有根区间.,得到新的有根区间[a1,b1],设?x在区间[a,b]上连续且?a?b0 . 0 a b ? y x y?x记a0a,b0b,计算 若|?x0|? ,则取??x0 ;否则,若?a0?x00,取a1a0,b1x0 ;若?a0?x00,取a1x0,b1b0 而且有根区间[a1,b1]长度是有根区间[a0,b0]长度的一半, x0 再对有根区间[a1,b1]重复上面运算, 即: 计算 若|?x1|?, 则取??x1; 否则,若?a1?x10,取a2a1 , b2x1 ;若?a1?x10, 取a2x1 ,b2b1,得到新的有根区间[a2,b2]. x1而且有根区间[a2,b2]长度是有根区间[a1,b1]长度的一半.一直进行下去,直到求出有根区间[ak,bk]. 此时,再计算或者有|?xk|? ,或者有 可见,k趋向无穷大时,xk收敛于? . 而且,若要|xk-?|? ,只要 此时可取近似根??xk .在计算过程中,若出现|?xk|?1,或bk-ak?2 .则可取xk作为方程?x0的近似根,终止运算. 例1用二分法求x3+4x-100在区间[1,2]内根的近似值,并估计误差.解 这里?xx3+4x-7， ?1?2-180,而且 ??x3x2+40,所以?x0在[1,2]区间有唯一根.取x01.5,由于?x02.375,得新有根区间[1,1.5], x11.25,由于?x1-0.0468,得新有根区间[1.25,1.5], x21.375,由于?x21.0996,得新有根区间[1.25,1.375], x31.3125,由于?x30.511,得新有根区间[1.25,1.3125],…………………………………………………. x91.座机电话号码3,得有根区间[1.座机电话号码3,1.座机电话号码5], x101.座机电话号码4, ?x10-0.000105285 取??x101.座机电话号码4作为方程根的近似值,且有 只需k5ln210-1?15.61.即需取??x16.如果取精度?10-5,则要使二分法要求函数在区间[a，b]上连续，且在区间两端点函数值符号相反，二分法运算简便、可靠、易于在计算机上实现。但是，若方程?x0在区间[a，b]上根多于1个时，也只能求出其中的一个根。另外，若方程?x0在区间[a，b]有重根时，也未必满足?a?b0. 而且由于二分法收敛的速度不是很快,一般不单独使用,而多用于为其他方法提供一个比较好的初始近似值.§2.1 简单迭代法的一般形式 §2 简 单 迭 代 法首先把方程?x0改写成等价同解形式x?x4.2 得到迭代序列xk , 如果xk?? ,则有???, 即?是方程?x0的根. 取一个合适的初始值x0,然后作迭代xk+1?xk ,k0,1,2,…4.3这种求方程根的方法称为简单迭代法,或逐次逼近法.其中?x 称为迭代函数 ,式4.3称为迭代格式. 若迭代序列xk 收敛 , 则称简单迭代法是收敛的.解 改写原方程为等价方程求方程x3-2x-30在[1,2]内的根. 例2,建立迭代格式 如果取初值x01.9 ,计算得 1.座机电话号码 1.座机电话号码 1.座机电话号码 1.座机电话号码 1.座机电话号码 …… 6 7 8 9 10 … 1.9 1.座机电话号码 1.座机电话号码 1.座机电话号码 1.座机电话号码 1.座机电话号码 0 1 2 3 4 5 xk k xk k由计算结果有，x10x9，因此可取??x101.座机电话号码.定义4.1 设?x为定义在区间I上的函数, 且对任何x?I,均有?x?I,则称?x为I到自身上的映射.方程也可改写成xx3-3/2, 建立迭代格式xk+1xk3-3/2 , k0,1,2,…仍取初值x01.9, 则有x11.9295, x22.0917, x33.0760, x413.0529 可见,xk??,此迭代格式是发散的. §2.2 简单迭代法的收敛条件定义4.2 设?x为I到自身上的映射,且存在0L1,使对任何x1,x2?I,有|?x2-?x1|? L|x2-x1|,则称?x为I上的压缩映射, L称为Lipschitz常数.若?x为I上的压缩映射,则?x在I上连续.定理4.2 若?x为I到自身上的映射,且?x?C1I, |??x|? L1,则?x为I上的压缩映射.证 对任意x1