计算函数零点和极值点的迭代法.pptVIP

定理4.1-5设Φ在闭区域上满足条件(1)存在q，0q1，使?x，y?，有||Φ(x)–Φ(y)||?q||x–y||(4.1-9)(2)?x??Φ(x)?则(1)存在唯一不动点x*?，且?x(0)?，迭代向量列收敛于x*。(2)注：证明与定理2.4-3及定理4.1-1相似*第31页，共77页，星期日，2025年，2月5日注：(1)保证序列(2)若?i(x)在上可微，Φ(x)=(?1(x)，?2(x)，…，?n(x))T记则压缩条件可用下式替代其中DΦ(x)是Φ(x)在点x处的Jacobi矩阵*第32页，共77页，星期日，2025年，2月5日例5用不动点迭代法求非线性方程在闭域内的根。解：(1)取x(0)=(1，-1.5)T，按迭代公式：k=0，1，2，…产生的序列收敛（见p339）k=1,x0=[1,-1.5]x1=[log(1-x0(2)),-sqrt(4-x0(1)^2)]whileabs(x1-x0)0.0001k=k+1,x0=x1;x1=[log(1-x0(2)),-sqrt(4-x0(1)^2)]end*第33页，共77页，星期日，2025年，2月5日例5用不动点迭代法求非线性方程在闭域内的根。解：(2)按迭代公式：k=0，1，2，…产生的序列未必收敛（见p340）k=1,x0=[1,-1.5]x1=[sqrt(4-x0(2)^2),1-exp(x0(1))]whileabs(x1-x0)0.0001k10k=k+1,x0=x1;x1=[sqrt(4-x0(2)^2),1-exp(x0(1))]end*第34页，共77页，星期日，2025年，2月5日4.2.1一元非线性方程f(x)=01.牛顿迭代方法线性化：设f(x)在点x(k)处可微，则展开式用线性部分近似表示f(x)?f(x(k))+f(x(k))(x–x(k))即考虑方程f(x(k))+f(x(k))(x–x(k))=04.2 Newton迭代法及其变形*第35页，共77页，星期日，2025年，2月5日若f(x(k))≠0，则有令：k=0，1，2，…(4.2-4)称为Newton迭代公式，其迭代函数为(4.2-5)*第36页，共77页，星期日，2025年，2月5日2.收敛性(1)若x*是f(x)的单重根：f(x*)=0，f(x*)≠0因为而一般不为0，所以，Newton法是局部二阶收敛的——由定理4.1-3*第37页，共77页，星期日，2025年，2月5日例1用Newton法求非线性方程f(x)=xex–1=0在(0，1)内的根，取x(0)=0.5。解：因为f(x)=(1+x)ex，故其Newton迭代公式为，k=0，1，2，…从x(0)=0.5出发，计算结果见p342表4.2-1。k=0,x0=0.5x1=x0-(x0*exp(x0)-1)/((1+x0)*exp(x0))whileabs(x1-x0)0.00001k10k=