计算方法 第六章 非线性方程的数值解法.pptVIP

本章介绍方程的迭代解法，它既可以用来求解代数方程，也可以用来解超越方程，并且仅限于求方程的实根。 运用迭代法求解方程的根应解决以下两个问题： 确定根的初值; 将进一步精确化到所需要的精度。 6.2 二分法 6.1.1确定有根区间的方法 为了确定根的初值，首先必须圈定根所在的范围， 称为圈定根或根的隔离。 在上述基础上，采取适当的数值方法确定具有一定 精度要求的初值。 对于代数方程，其根的个数（实或复的）与其次数 相同。至于超越方程，其根可能是一个、几个或无 解，并没有什么固定的圈根方法 求方程根的问题，就几何上讲,是求曲线 y=f (x)与 x轴交点的横坐标。 由高等数学知识知, 设f (x)为区间[a,b]上的单值连续, 如果f (a)·f (b)0 , 则[a,b]中至少有一个实根。如果f (x)在[a,b]上还是单调地递增或递减，则仅有一个实根。 (1) 画图法 画出y = f (x)的略图，从而看出曲线与x轴交点的 大致位置。 也可将f (x) = 0分解为?1(x)= ?2(x)的形式，?1(x) 与 ?2(x)两曲线交点的横坐标所在的子区间即为含根 区间。 例如 xlogx-1= 0 可以改写为logx=1/x 画出对数曲线y=logx,与双曲线y= 1/x,它们交 点的横坐标位于区间[2,3]内 (1) 画图法 用逐步搜索法进行实根隔离的关键是选取步长h 要选择适当h ，使之既能把根隔离开来，工作量 又不太大。 为获取指定精度要求的初值,可在以上隔离根的 基础上采用对分法继续缩小该含根子区间 二分法可以看作是搜索法的一种改进。 6.3 迭代法 6.4 牛顿迭代法 6.4.3 牛顿迭代法的收敛性 例6.11 用牛顿迭代法求 f(x)=xex-1的根,ε=10-4 解：因 f (x)= x ex –1 , f ′(x)=ex ( x+1) 建立迭代公式 牛顿法例题 牛顿法例题（续） 牛顿法优缺点 由以上图形可以看出，序列{xk}的收敛速度，取决于曲线 在根附近的斜率 ，由拉格朗日定理可知 其中? k在xk和xk-1之间。所以在根x*附近，??’(x) ?恒小于1， 则此迭代序列收敛，若??’(x) ? 1，则此序列发散。 6.3.3 迭代法收敛的条件 对方程f(x)=0可以构造不同的迭代公式, 但迭代公式 并非总是收敛。 那么,当迭代函数 满足什么条件时，相应的迭代公式才收敛呢？即使迭代收敛时，我们也不可能迭代很多次，而是迭代有限次后就停止，这就需要估计迭代值的误差，以便适时终止迭代 定理6.1 设函数 在[a,b]上具有连续的一阶导 数, 且满足 （1）对所有的x∈[a,b] 有 ∈[a,b] （2）存在 0 L 1 ,使所有的x∈[a,b]有 ≤ L 则方程 在[a,b]上的解 存在且唯一 ，对任意的 ∈[a ,b] ，迭代过程 均收敛于 。并有误差估计式 ① ② 证: 构造函数 ,由条件①对任意的x∈[a, b] ∈[a, b]有 由连续函数介值定理知, 必有 ∈[a, b], 使 所以有解存在, 即 假设有两个解 和 , , ∈[a, b]，则 ， 由微分中值定理有其中ξ是介于x*和 之间的点 从而有ξ∈[a,b]，进而有 由条件(2)有 1, 所以 - =0，即 = ，解唯一。 按迭代过程 ,有 由于L1,所以有 ,可见L越小,收敛越快 再证误差估计式 ① ② ∵ ∴ 即 ① 得证。 即 ② 得证。 6.3.4 迭代法的算法框图 例5 对方程 ,构造收敛的迭代格式, 求其最小正根,计算过程保留4位小数。 解 容易判断[1,2]是方程的有根区间, 且在此区间 内 ,所以此方程在区间[1,2]有 且仅有一根.将原方程改写成以下两种等价形式。