第6章 非线性方程的计算方法.pptVIP

第1页，共20页，星期日，2025年，2月5日§2BisectionMethodabx1x2abWhentostop?或不能保证x的精度x*?2xx*第2页，共20页，星期日，2025年，2月5日6.1.1BisectionMethod误差分析：第1步产生的有误差第k步产生的xk有误差对于给定的精度?,可估计二分法所需的步数k：①简单;②对f(x)要求不高(只要连续即可).①无法求复根及偶重根②收敛慢注：用二分法求根，最好先给出f(x)草图以确定根的大概位置。或用搜索程序，将[a,b]分为若干小区间，对每一个满足f(ak)·f(bk)0的区间调用二分法程序，可找出区间[a,b]内的多个根，且不必要求f(a)·f(b)0。第3页，共20页，星期日，2025年，2月5日6.1.1BisectionMethodLab02.BisectionMethodUsetheBisectionmethodtofindonmgivenintervalsthemrealrootsofagivenpolynomialofdegree5?n?m,.InputThereareseveralsetsofinputs.Foreachset: The1stlinecontainsanintegernwhichisthedegreeofapolynomial.n=?1signalstheendoffile. The2ndlinecontainsn+1realnumberswhicharethecoefficientsofthepolynomial.Thenumbersareseparatedbyspaces. The3rdlinecontainsanintegerMaxandtworealnumberseps1andeps2,whereMaxisthemaximumnumberofiterations,eps1istheaccuracyboundforxandeps2istheaccuracyboundforp(x). The4thlinecontainsanintegerm(n?m?0),followedbympairsofrealnumbersa1b1…ambmwhicharetheendpointsoftheintervals[a1,b1]…[am,bm].第4页，共20页，星期日，2025年，2月5日Output?EachrootistobeprintedasintheCfprintf:fprintf(outfile,%12.7f?,root);/*here?representsaspace*/Ifthereisnorootfoundinaninterval,simplyprint“no?root?”.Theoutputofthenextsetmustbeprintedinanewline.SampleInput(?representsaspace)21?0??11000?002??2??0.5?0.5?231?0?0??11000?002??1?0?0?2?1SampleOutput(?representsaspace)???1.0000000????1.0000000?no?root????1.0000000?6.1.1BisectionMethod第5页，共20页，星期日，2025年，2月5日6.1.2迭代法/*Fixed-PointIteration*/f(x)=0x=φ(x)等价变换f(x)的根φ(x)的不动点思路从一个初值x0出发，计算x1=φ(x0),x2=φ(x1),…,xk+1=φ(xk)