弦截法求方程根课件.ppt

第二章 非线性方程求根 ;【本章重点】 1． 不动点迭代法及其收敛性与收敛速度。 2． Newton迭代法;【课前思考】 1．什么是方程f(x)=0求根的二分法？如何估计近似 根xn的误差? 2．什么是不动点迭代法？怎样判断迭代法 的收敛性？ 3．迭代法收敛速度：收敛阶定义，如何加速迭代收 敛？ 4．给出方程f(x)＝0求根的Newton法，它有何优缺点 ？如何用Newton法求方程根？;§1 方程求根与二分法;　 而当n=3,4时方程的根，可在数学手册中查到虽可用公式表示，但表达式太复杂，一般不用；当n≥5已没有直接用公式表达的求根算法。因此对n≥3的代数方程求根方法与一般超越方程(1.1)一样都采用迭代方法求根，设 (表示f在区间 上连续)，若有f(a)f(b)＜0,则f(x)=0在区间 上至少有一个实根， 称为有根区间，通常可用逐次搜索法求得方程(1.1)的有根区间.; 二分法原理 给定方程f(x)=0,设f(x)在区间[a,b]连续,且f(a)f(b)0,则方程f(x)在(a,b)内至少有一根,为便于讨论,不妨设方程f(x)=0在(a,b)内只有一实根x*采取使有根区间逐步缩小,从而得到满足精度要求的实根 x*的近似值xk。 取[a,b]区间二等分的中点x0 =(a+b)/2, 若f(x0)=0,则x0是f(x)=0的实根 若f(a)f(x0)0 成立,则 必在区间(a, x0)内,取a1=a,b1= x0;否则 必在区间(x0,b)内,取a1= x0,b1=b, 这样,得到新区间(a1,b1),其长度为[a,b]的一半,如此继续下去,进行k次等分后,得到一组不断缩小的区间,[a,b],[a1,b1],......[ak,bk].;;a; ; ; 例3:用二分法求方程 在区间(1,2)内的实根,要求误差限为 。;§1 方程求根与二分法;例：利用二分法求方程f(x)= -tanx在区间[0，∏/2] 内的实根，使精度达到10-5。;function [x,k]=deminmethod(a,b,f,emg) fa=feval(f,a); fab=feval(f,(a+b)/2); k=0; while abs(b-a)emg if fab==0 x=(a+b)/2; return; elseif fa*fab0 b=(a+b)/2; else a=(a+b)/2; end fa=feval(f,a); fab=feval(f,(a+b)/2); k=k+1; end x=(a+b)/2;;§2 迭代法;§2 迭代法;;定理;证明：① g(x) 在[a, b]上存在不动点？;④;注：上述定理为全局收敛性，不易检验定理条件。可将[a, b]缩小，定义局部收敛性：若在 x* 的某? 邻域 P:| x ??? x* | ? ? 有 g(x)?C1[a, b] 且 | g’(x*) | 1，则由?x0?P 开始的迭代序列{xk} ?P局部收敛于x* 。即调整初值可得到收敛的结果。;定义：设迭代过程xk+1=g(xk)收敛于方程x=g(x)的根x*, 如果迭代误差ek=xk-x* 当k- ∞时成立下列渐进关系：;分析：g(x)在x*处的p-1阶泰勒展式;发散;2阶收敛;例：利用迭代法求方程 f(x)=x3-x-1在x0=1.5附近 的根x*，使精度达到10-5。;f = 1.50000000000000 1.35720880829745 1.33086095880143 1.32588377423235 1.32493936340188 1.32476001129270 1.32472594522689 1.32471947453436 k = 7;? Aitken（埃特金） 加速：;一般地，记：;? Steffensen （斯蒂芬森）加速：;§3 迭代收敛的加速法;§3 迭代收敛的加速法;程序设计; [f,k]= Steffensen(10^-5,1.5) f = 1.50000000000000 1.41629297458894 1.35565044147664 1.32894877728401 1.32480448904104 1.32471799396881 1.32471795724475 1.32471795724475 k = 7;§4 牛顿法