3非线性方程matlab计算方法.pptVIP

定义. :设序列{xk}收敛于x*，若存在p≥1和正数c， 使得成立 则称{xk}为 p 阶收敛的 特别, p = 1，要求c1, 称线性收敛; 1p2,称超线性收敛 p=2,称平方收敛。 迭代法的收敛阶(收敛速度) 几何意义: x y y = x x* y=g(x) x0 p0 x1 p1 ? x y y = x x* y=g(x) x0 p0 x1 p1 ? x y y = x x y y = x x y y = x x y y = x x* x* x* x* y=g(x) y=g(x) y=g(x) y=φ(x) x0 p0 x1 p1 ? x0 p0 x1 p1 ? x0 p0 x1 p1 ? x0 p0 x1 p1 ? 取x0作为初始近似值，将f(x)在x0做Taylor展开: 重复上述过程 ? 作为第一次近似值 Newton 迭代公式 基本思路：将非线性方程f(x)=0 线性化 §2.3 Newton-Raphson迭代法 (1) 选定初值x0 ,计算f (x0) , f ?(x0) 计算步骤 (2) 按公式 迭代 得新的近似值xk+1 (3) 对于给定的允许精度?，如果 则终止迭代，取 ;否则k=k+1，再转 步骤(2)计算 允许精度 最大迭代次数 迭代信息 用迭代法求方程在０，1.5处的解 例 用Newton迭代法求方程xex-1=0在0.5附近的根,精度要求?=10-5. 解 Newton迭代格式为 k xk ?(xk) |xk-xk-1| 0 1 2 3 4 0.5 0000-0000.0000000003 0.0000000003 0000牛顿法的几何意义 x y x* x0 x 1 x 2 牛顿法也称为切线法 这样一直下去，我们可以得到迭代序列 Newton迭代的等价方程为： 所以 若f(x)在a处为单根，则 所以，迭代格式收敛 注：Newton’s Method 收敛性依赖于x0 的选取。 x* x0 ? x0 ? x0 1、优点：牛顿迭代法具有平方收敛的速度，所以在迭代过程中只要迭代几次就会得到很精确的解。这是牛顿迭代法比简单迭代法优越的地方。 牛顿迭代法的优缺点 2、缺点：选定的初值要接近方程的解，否则有可能得不到收敛的结果。再者，牛顿迭代法计算量比较大。因每次迭代除计算函数值外还要计算微商值。 牛顿法主要有两个缺点：局部收敛，计算量大。 2.4 弦截法 将Newton迭代中的导数，用差商代替，有格式 是2步格式。收敛速度比Newton迭代慢 x0 x1 切线 割线 习题: 要求:习题二:1,2,3,4,6,9,10 选作：数值实验题 复习：分别用切法和割线法解下面方程 nIter=1,x=3.000000 nIter=2,x=2.333333 nIter=3,x=2.238095 nIter=4,x=2.236069 nIter=5,x=2.236068 输入初始值x0,x1,误差eps，最大迭代次数 1 2 0.00005 10 nIter=2,x=2.333333 nIter=3,x=2.230769 nIter=4,x=2.235955 nIter=5,x=2.236068 nIter=6,x=2.236068 * 方程是在科学研究中不可缺少的工具.方程求解是科学计算中一个重要的研究对象. 几百年前就已经找到了代数方程中二次至五次方程的求解公式.但是,对于更高次数的代数方程目前仍无有效的精确解法.对于无规律的非代数方程的求解也无精确解法.因此,研究非线性方程的数值解法成为必然. 计 算 方 法 第三章 非线性方程的数值解法  史晓非 大连海事大学信息科学技术学院 信号与图像处理研究所 简介（Introduction） 我们知道在实际应用中有许多非线性方程的例子，例如 （1）在光的衍射理论（the theory of diffraction of light)中,我们需要求x-tanx=0的根 （2）在行星轨道（ planetary orbits）的计算中,对任意的a和b，我们需要求x-as