非线性方程数值解法.pptVIP

例：给出下列迭代过程的收敛速度. 改进： 方法一：如果 的根 的重根数已知，那么可 将迭代函数改写为 ,由上面的推导过程不 难验证 ,从而下面迭代公式至少平方收敛 方法二：如果 的根 的重根数m 未知，那么可将方程 改写为 由于 是 的 重根， 是 的 重根。 从而 是函数 的单根。对于 使用牛顿迭代法，则至少具有平方收敛速度。 第四章 非线性方程的数值解法 设f(x)为一元连续函数，称方程f(x)=0为函数方程。 当f(x)不是x的线性函数时，称对应的方程为非线性方程。 定义 分类 当f(x)是多项式时称为代数方程，否则称为超越方程。 方程的根也称为函数的零点．若可表示为: 则称x* 为 f (x) = 0 的 m 重根，当m=1时称为单根， 此时有: 定义 §4.2 二分法和试位法 零点定理：若 f ?C[a, b]，且 f (a) · f (b) 0，则 f 在 (a, b) 上必有一根。 二分法是方程求根中最常用且最简单的方法．其基本思想是：先利用零点定理确定根x* 的存在区间，然后将含根x* 的区间对分，通过判别对分点函数值的符号，将有根区间缩小一半．重复以上过程，将根的存在区间缩到充分小，从而求出满足精度要求的根的近似值． 一、基本原理及做法 a1 b1 x1 a2 b2 x* 此时的[a2, b2]必是根的存在区间，重复此过程就得到一个包含根x* 的区间套，满足 二、误差分析 第1步产生的 有误差 对于给定的精度 ? ,可估计二分法所需的步数 k ： 例1 证明方程 在[1,2]上有一个实根 ，并用二分法求这个根，要求 。若要求 ，需二分区间[1,2]多少次？ 第 k 步产生的 xk 有误差 三、算法 1．对 做到第5步 2． 3． 5．若 则 6．输出超过最大迭代次数 N 的信息，停机。 输入 ，区间 ，精度 ，最大的迭代次数N 4．若 ，则输出 停机 ①简单; ② 对f (x) 要求不高(只要连续即可) . 收敛慢 注：用二分法求根，最好先给出 f (x) 草图以确定根的大概位置 , 或用搜索程序，将[a, b]分为若干小区间，对每一个满足 f (ak)·f (bk) 0 的区间调用二分法程序，可找出区间[a, b]内的多个根，且不必要求 f (a)·f (b) 0 。 四、二分法的特点 注：在应用中常用二分法来判别根的存在区间，如区间 较大可用二分法适当收缩区间，并选择初值为该区间中 点．再用收敛速度快的迭代法，迭代求根． 例2 求x3-3x+1=0的实根分布情况，要求区间长度不大于1， 并求出最小正根的近似值。 53 19 3 -1 1 3 -1 -17 -51 f(x) 4 3 2 1 0 -1 -2 -3 -4 x 可见，在(-2,-1)区间、(0,1)区间、(1,2)区间各有一实根， 下面求(0,1)区间上的实根，列表： f(xk) xk bk ak k 0.000017 0.347290 0.347412 0.347168 13 -0.000305 0.34741 0.34766 0.34717 12 … … … … … -0.07227 0.375 0.5 0.25 3 0.265625 0.25 0.5 0 2 -3.75 0.5 1 0 1 可见x?[0.347167968,0.347412109] 若取?=0.0005,由表5.2可知,当k=13时, bk-ak=0.0002440.0005,此时过程结束,取 x*=(0.347168+0.347412)/2=0.34729?0.347 §4.3 不动点迭代法 f (x) = 0 等价变换 f (x) 的根 一、 迭代法的一般形式 迭代格式 迭代函数 的不动点 x = 从x0 出发，计算 x1 = , x2 = ,…, xk+1 = ,… 思路： 若 收敛，即存在x* 使得 ，且 连续，则由