2012数值分析-lec11-非线性方程组的解法.pptVIP

数值分析 朱立永 北京航空航天大学 数学与系统科学学院 Email: numerical_analysis@ Password:beihang2012 答疑时间：星期三下午2:00－5:00 答疑地点：主216 第十讲非线性方程组的迭代解法 第三章非线性方程与非线性方程组的迭代解法 非线性方程组的解法 含n个方程的n元非线性方程组的一般形式是 非线性方程组的一些基本概念 上面的方程组化为：F(X)=0 例子： 将方程组 F(X)=0, 写成与之等价的形式:X=G(X), 然后再利用 X(k+1)=G(X(k)), k=1,2,…,求解原方程的根。 简单迭代法 x=0.0 y=0.0 10 x1=x y1=y x=0.25*(1+y1-0.1*exp(x1)) y=0.25*(x1-0.125*x1*x1) write(10,*) x,y if ((abs(x-x1)+abs(y-y1)).lt.0 then goto 15 endif goto 10 15 end 简单迭代法的收敛性 Newton迭代法 计算步骤（框图）： k x (k) 0 1 2 3 (1.5, 1.0)T (1.5, 0.75)T (1.488095, 0.755952)T (1.488034, 0.755983)T 例：用牛顿法解方程组 取初始值（1，1，1），计算如下 N x y z 0 1.0000000 1.0000000 12.1893260 1.5984751 1.3939006 1.8505896 1.4442514 1.2782240 1.7801611 1.4244359 1.2392924 1.7776747 1.4239609 1.2374738 1.7776719 1.4239605 1.2374711 1.7776719 1.4239605 1.2374711 简化牛顿法。目的是避免计算迭代公式中繁杂的导数，解决方法与一元函数牛顿法类似，即将所有导数取为固定值，如迭代初值的导数值。 与单个方程的情形类似，牛顿法中 f 的导数的元素用合适的差商来近似，如 就可得到拟牛顿法或弦截法。 若用格式 其中下山因子 合适地选取使得 就得到牛顿下山法。 若用格式 ，其中 是 的简单 修正，且满足 则得到Broyden算法。特别，若取 ，其中 u,v 是待定的列向量，使其满足上式，则得到秩一Broyden算法。 比如 小结 1、本章的目的是求解形如 f (x)=0 的方程，而其核心方法是 将所要求解的方程变形为 x = ? (x)，利用 ? (x) 为压缩映射， 通过迭代求出其解。 2、变形中切记要恒等变形！ 3、在恒等变形中，为使变形得到的函数? (x) 为压缩映射，一 个技巧是利用待定参数。 4、恒等变形的一种重要格式是牛顿迭代，证明其迭代收敛阶 的一个常用技巧是泰勒展开。 5、n维空间中代数方程迭代求解的收敛条件是谱半径小于 1 . Newton-Krylov subspace methods 作业 教材P93页习题10、11 * * *