第二章代数方程(组)的数值解法.ppt

2.4.3 威格斯坦法 * * * * 2.4.3 威格斯坦法 2.4.4 牛顿-拉夫森法 牛顿-拉夫森法是牛顿法的一种推广，将非线性方程组逐次进行线性化处理，从而构造迭代算法。 * * 2.4.4 牛顿-拉夫森法 * * 2.4.4 牛顿-拉夫森法 令 则化为一组线性方程组 * * 例题 油换热器进出口温度T1、T2。 T1(0）=180℃ ,T2(0）=292 ℃,要求精度 * * 2.4.4 牛顿-拉夫森法 * * * * * * * * 由于每一区间都是前一区间的一半，因此区间[an , bn]的长度为 若每次二分时所取区间中点都不是根，则上述过程将无限进行下去. 当 n→∞ 时，区间必将最终收缩为一点x* ，显然x*就是所求的根. * * 2.3.1 二分法 若取区间[an , bn]的中点 作为x*的近似值，则有下述误差估计式 只要 n 足够大, (即区间二分次数足够多)，误差就可足够小. * * 2.3.1 二分法 　　二分法的计算步骤: 　　步骤1 准备 计算函数f(x)在区间[a, b]端点处的值f(a), f(b). 若f(a)·f((a+b)/2)0, 则以(a+b)/2代替b ，否则以(a+b)/2代替a. 　　步骤2 二分 计算函数f(x)在区间中点(a+b)/2处的值f((a+b)/2). 　　步骤3 判断 若f((a+b)/2)=0，则(a+b)/2即是根，计算过程结束，否则检验. 反复执行步骤2和步骤3，直到区间[a, b]长度小于允许误差ε，此时中点(a+b)/2即为所求近似根. * * 例1 用二分法求方程 f(x)=x3-x-1=0在（1，1.5）内的实根,要求误差不超过0.005. 解： 已知x∈(1, 1.5), 要想满足题意，即： 则要 |x*-xn|≤0.005 由此解得 取n=6, 按二分法计算过程见下表, x6 = 1.3242 为所求之近似根. * * n an bn xn f(xn) 说明 0 1 2 3 4 5 6 1.0 1.25 1.25 1.3125 1.3125 1.3125 1.3203 1.5 1.5 1.375 1.375 1.3438 1.3281 1.3281 1.25 1.375 1.3125 1.3438 1.3281 1.3203 1.3242 - + - + + - - (1) f(a)0, f(b)0 (2) 根据精 度要求， 取到小数 点后四位 即可. 二分法的优点是算法简单，便于理解；缺点是收敛的太慢。故一般不单独将其用于求根，只是用其为根求得一个较好的近似值. * * 2.3.2 迭代法 将方程f(x)=0改写为等价方程形式 x=?(x). 若要求x*满足f(x*)=0，则x*=?(x*)；选择一个初始近似值x0，即可求得 x1=?(x0). xk+1=?(xk) (k=0,1,2,…,n). * * 几何意义： x2 x1 x0 y = x * * 分别按以上三种形式建立迭代公式，并取x0=1进行迭代计算，结果如下： 解 对方程进行如下三种变形： 例2 用迭代法求方程x4+2x2-x-3=0 在区间[1, 1.2]内的实根. * * 准确根 x* = 1.124123029, 可见迭代公式不同, 收敛情况也不同. 第二种公式比第一种公式收敛快得多, 而第三种公式不收敛. * * 收敛 发散 * * 收敛的条件 * * 　　例如，在前面例2中采用的三种迭代公式，在隔根区间(1, 1.2)内，有 故前两个迭代公式收敛，第三个迭代公式不收敛. * * 2.3.3 威格斯坦法 * * 简单迭代法受收敛条件 的约束，且收敛速度较慢，威格斯坦法是对简单迭代法的一种改进。它不仅突破约束，而且加快了收敛速度。 首先给出两点 与曲线交于两点 连接 交 y=x 于一点 x2 由以上两式简化合并可得： 2.3.3 威格斯坦法 找出简单迭代和威