计算方法2010-2.ppt

第二章 非线性方程的数值解法 非线性方程的求解过程 判断是否有根，如果有，有几个根，根的存在范围； n次的代数方程最多有n个根，可能有一对共轭复根 根的隔离，求隔根区间 试值法 作图法 扫描法 根的精确化 对分法 迭代法 牛顿法 弦割法 本章重点：对分法、迭代法和牛顿法 2.2 根的隔离-试值法 试值法 2.2 根的隔离-作图法 y=x3-3x-1 2.2 根的隔离-扫描法 扫描法就是将有根区间等分为若干个子区间，然后逐个小区间检验是不是隔根区间，检验的办法就是判断区间端点的函数值是否异号。 扫描法 扫描法的算法 代数方程的实根的上、下界 2.3对分法 对分法的算法 对分法的程序 Private Sub Form_Click() Dim a!, b!, x!, c! a = 3: b = 4 While Abs(b - a) 0.00001 c = (a + b) / 2 If f(c) = 0 Then Print c End Else If f(a) * f(c) 0 Then b = c Else a = c End If End If Wend Print c End Sub 对分法程序 float f(float x) {return x*x*x-2*x*x-4*x-7;} main () { float a,b,c,eps; scanf(%f%f,a,b); scanf(%f,eps); while (b-a=eps) { c=(a+b)/2; if(f(c)==0) break; else if (f(a)*f(c)0) b=c; else a=c; } printf(%f,c); } 2.4 迭代法 迭代法举例 迭代法举例 例2.3 迭代法的收敛性 思考 对于同一个方程，不同的做法收敛性是不一样的，那么，收敛还是发散受什么条件的影响？ 两种做法主要的区别在什么地方？ Φ(x)满足什么条件会收敛？ 迭代收敛性定理，也称不动点原理 迭代收敛定理 定理分析 已知条件 结论 两个不等式代表的意义 定理的应用 局部收敛 迭代法的算法 迭代法的加速（选学） 如何加快迭代收敛的速度？ 根据不等式 迭代法的加速 1 松弛法 在x=φ(x) 的两端同时加上λx，得 迭代法的加速 迭代法的加速 2 埃特金法（Altken） 松弛法中计算导数不方便，用差商代替导数 迭代法的几何意义 2.5 牛顿法 收敛速度 简单迭代的收敛速度 牛顿法的收敛速度 牛顿法举例 2.6 弦割法 第一次作业 构造满足定理条件的等价形式一般难于做到。 要构造收敛迭代格式有两个要素： 1、等价形式 2、初值选取 迭代过程往往就在根的附近进行，只要假定 在 附近连续，且满足 ，则根据连续函数的性质，一定存在 的某邻域S： ，使得 在S上满足定理2.1的条件，故在S中任取初始值x0，迭代公式 必将收敛于方程 的根 这种收敛称为局部收敛。 (1)输入初始近似值x0，精度要求eps，控制最大迭代次数M； (3)当k M且 eps时做循环 ， ， (4)如果 eps，则输出x1；否则输出迭代失败信息。 (2) ② 知L越小，收敛速度越快 加快收敛速度的办法就是设法让L变小 变形为： 相当于 根据迭代收敛定理，导数要尽可能地小，所以取 这样，迭代公式变为： 令 则 此即松弛法，不收敛的迭代函数，经加速后一般也能获得收敛。 ， 用 代替 得 得迭代公式 埃特金公式有时可能使一个本来不收敛的迭代格式获得收敛。 x y y = x x* y=g(x) x0 p0 x1 p1 ? x y y = x x* y=g(x) x0 p0 x1 p1 ? 将f(x)在初值处作Taylor展开 取线性部分作为f(x)的近似，有： 若 ，则有 记为 类似，我们可以得到 x y x* x0 这样一直下去，我们可以得到迭代序列 Newton迭代的等价方程为： 所以 若f(x)在a处为单根，则 所以，迭代格式收敛 定义2.1 设由迭代公式 产生的迭代序列 (k=0,1,2,…)收敛于方程 的根 ，记 ，若存在实数 及非零常数c，使得 则称迭代过程是p阶收敛的。 当p=1时，称为线性收敛； 当p1时，称为超线性收敛；当p=2时，称为平方收敛。 显然，p越大收敛速度越快。 由微分中值定理可知，必存在一点 介于xk与 之间，使得 于是 由此可知，简单迭代法至少是线性收敛的。 收敛速度 函数在a处作Taylor展开 若a为p重根，取迭代格式为： 即 牛顿迭代收敛速度快，格式简单，应用广泛