第2章非线性方程求根.pptVIP

第30页，共54页，星期日，2025年，2月5日§4迭代法的收敛性一、收敛性定义2(p26)设序列{Xk}收敛到X*，记误差εk=Xk-X*(1)?若存在实数P≥1及非零常数C，使limP=Ck-∞?则称{Xk}为P阶收敛，C称为渐进误差常数，当P=1，2时分别称{Xk}为线形收敛及平方收敛第31页，共54页，星期日，2025年，2月5日例证明xk+1=xk(xk2+3a)/(3xk2+a)是计算的三阶方法要证明|xk+1-√a|/|xk-√a|3为常数第32页，共54页，星期日，2025年，2月5日定理3：设x*为Φ的不动点，整数p1,Φ(p)在x*某邻域存在且连续,且满足Φ（k)(x*)=0(k=1,2,….p-1)Φ(p)(x*)≠0则由Xk+1=Φ(Xk)产生的序列{Xk}在x*的邻域是p阶收敛的特别地，（p26)(1)当Φ’(x*)≠0，则为线形收敛（2）当Φ’(x*)=0，Φ”(x*)≠0，则为平方收敛第33页，共54页，星期日，2025年，2月5日例考查取Φ(X)=1+1/x+1/x2时，用Xk+1=Φ(Xk）求解x3-x2-x-1=0在区间[1，2]根的收敛性和收敛阶（x*=1.839)解：由Φ(X)=1+1/x+1/x2得知Φ’(X)=-（1/x2+2/x3)对任意[1，2]区间上的x不满足|Φ’(X)|1,也即不满足全局收敛定理但可验证：当x∈[1.7,2]时|Φ’(X)|1若取x0=2可算出x1=1.91293….x13=1.83929对Φ’(X)=-（1/x2+2/x3)，显然Φ’(X*)≠0一阶收敛第34页，共54页，星期日，2025年，2月5日二、加速收敛迭代加速公式设xk+1*为近似值xk迭代后的结果，x*为x=Φ(x)的根则x*-xk+1*=Φ(x*)-Φ(xk)=Φ’(ζ)(x*-xk)设Φ’(x)在求根范围内变化不大即Φ’(x)≈L，|L|〈1所以x*-xk+1*≈L(x*-xk)整理得：x*-Lx*-xk+1*+Lxk+1*=Lxk+1*-Lxk（1-L）x*-（1-L）xk+1*=L（xk+1*-xk）x*-xk+1*≈(xk+1*-xK)第35页，共54页，星期日，2025年，2月5日取第k+1次的近似值为:Xk+1=xk+1*+(xk+1*-xK)第36页，共54页，星期日，2025年，2月5日得迭代加速公式（1）迭代xk+1*=Φ(xk)（2）改进Xk+1=xk+1*+(xk+1*-xK)（例p28）第37页，共54页，星期日，2025年，2月5日2、Aitken加速公式（1）迭代yk=Φ(xk)(2)校正zk=Φ(yk)(3)改进xk+1=zk-可达到二阶收敛第38页，共54页，星期日，2025年，2月5日第2章非线性方程求根《计算方法》第2章非线性方程求根第1页，共54页，星期日，2025年，2月5日第二章方程求根我们已经熟悉求解一元一次方程、一元二次方程以及某些特殊类型的高次代数方程或非线性方程的方法。这些方法都是代数解法,求出的根是方程的准确根。但是在许多实际问题中遇到的方程,例如代数方程x3-x-1=0?或超越方程§1引言第2页，共54页，星期日，2025年，2月5日方程的形式很多