非线性方程选读.ppt

第二章 非线性方程的数值解法  简介（Introduction） 我们知道在实际应用中有许多非线性方程的例子，例如 （1）在光的衍射理论（the theory of diffraction of light)中,我们需要求x-tanx=0的根 （2）在行星轨道（ planetary orbits）的计算中,对任意的a和b，我们需要求x-asinx=b的根 （3） 在数学中，需要求n次多项式xn + a1 xn-1+...+an-1 x + an ＝0的根 方程根的数值计算步骤 判断根的存在 确定根的分布范围 根的精确化 迭代法流图 几何意义: 收敛性分析P32 定义 若存在常数?(0≤ ?1),使得对一切x1,x2∈［a,b］, 成立不等式 |g(x1)-g(x2)|≤ ? |x1-x2|， 则称g(x)是［a,b］上的一个压缩映射， ?称为压缩系数 重要： （迭代法的局部收敛定理） 定理条件非必要条件，而且定理中的压缩条件不好验证，一般来讲， 若知道迭代函数g(x)∈C1『a,b］(1),并且满足|g′(x)|≤ ? ≤1(2),对任意的x∈［a,b］(1),则g（x)是［a,b］上的压缩映射 例题 已知方程2x-7-lgx＝0，求方程的含根区间，考查用迭代法解此方程的收敛性。 在这里我们考查在区间[3.5,4]的迭代法的收敛性 很容易验证：f(3.5)0,f(4)0 将方程变形成等价形式：x＝(lgx+7)/2 举例 用一般迭代法求x3-x-1=0的正实根x* 例题 用一般迭代法求方程x-lnx＝2在区间（2，?）内的根，要求|xk-xk-1|/|xk|=10-8 将方程化为等价方程：x＝2＋lnx k xi 0 3.000000000 1 3.098612289 2 3.130954362 3 3.141337866 4 3.144648781 5 3.145702209 6 3.146037143 另一种迭代格式： 计算步骤 用迭代法求方程在０，1.5处的解 2.4 弦截法 习题: 要求:习题二:2,3,4,5,6,7 选做:编程实现: 1:给定方程的隔根区间 2:二分法 3:简单迭代法 4:牛顿迭代法 5:割线法 Sipi.dlmu.edu.cn 复习：线性代数中的“行列式,矩阵运算 　和线性方程”知识 复习：分别用切法和割线法解下面方程 容易得到：g′(x)在包含x*的某邻域U(x*) 内 连续，且|g′(x*)|1 解：令f(x)=x-lnx-2 f(2)0,f(4)0,故方程在（2,4）内至少有一个根 因此，? x0?（2，?），xk+1＝2＋lnxk产生的序列? xk ?收敛于X* 取初值x0＝3.0,计算结果如下： 7 3.146143611 8 3.146177452 9 3.146188209 10 3.146191628 11 3.146192714 12 3.146193060 13 3.146193169 14 3.146193204 0 3.000000000 1 3.147918433 2 3.146193441 3 3.146193221 定义. :设序列{xk}收敛于x*，若存在p≥1和正数c， 使得成立 则称{xk}为 p 阶收敛的 特别, p = 1，要求c1, 称线性收敛; 1p2,称超线性收敛 p=2,称平方收敛。 迭代法的收敛阶(收敛速度) 取x0作为初始近似值，将f(x)在x0做Taylor展开: 重复上述过程 ? 作为第一次近似值 Newton 迭代公式 基本思路：将非线性方程f(x)=0 线性化 §2.3 Newton-Raphson迭代法 牛顿法的几何意义 x y x* x0 x 1 x 2 牛顿法也称为切线法 这样一直下去，我们可以得到迭代序列 Newton迭代的等价方程为： 所以 若f(x)在a处为单根，则 所以，迭代格式收敛 注：Newton’s Method 收敛性依赖于x0 的选取。 x* x0 ? x0 ? x0 (1) 选定初值x0 ,计算f (x0) , f ?(x0) (2) 按公式 迭代 得新的近似值xk+1 (3) 对于给定的允许精度?，如果 则终止迭代，取 ;否则k=k+1，再转 步骤(2)计算 允许精度 最大迭代次数 迭代信息 例 用Newton迭代法求方程xex-1=0在0.5