计算方法2012-2.pptVIP

解决的问题 主要内容 非线性方程的求解过程 2.2 根的隔离-试值法 试值法 的函数曲线 2.2 根的隔离-作图法 y=x3-3x-1的曲线图 2.2 根的隔离——扫描法 扫描法就是将有根区间等分为若干个子区间，然后逐个小区间检验是不是隔根区间，检验的办法就是判断区间端点的函数值是否异号。 扫描法 扫描法的算法 扫描法的算法 代数方程的实根的上、下界 对分法 2.3 对分法 对分法的算法 对分法的VB程序 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 对分法C程序 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 牛顿法 收敛速度 简单迭代的收敛速度 牛顿法的收敛速度 牛顿法举例 牛顿法求重根 将重根变成另外一个方程的单根 f(x)/f (x)=0 2.6 弦割法 双点弦割 单点弦割 第一次作业 证： 先证根的存在性， 由条件(1)知，在区间[a,b]上有 当F(a)=0或F(b)=0时，a或b就是方程的根，否则有 因 连续，所以F(x)也连续，由连续函数性质可知， 存在 ，使 ，即 ，于是 是方程的根。 作辅助函数 可知，F(x)在区间[a,b]上严格单调上升，所以F(x)在此区间上至多有一个根，根的唯一性得证。 再证根的唯一性。 最后证明迭代序列的极限就是方程的根。由 根据微分中值定理，必存在ξ介于xk与α之间，及 介于xk与xk -1之间，使得 由条件(2)知 (2.8) 于是 整理可得 此即(1)式。 将上式代入(2.6)式后，即得 此即(2.7)式。 再反复利用(2.8)式的第2式，可得 由 可知必有 即 故迭代法收敛于方程的根得证。 构造满足定理条件的等价形式一般难于做到。 要构造收敛迭代格式有两个要素： 1、等价形式Φ(x) 2、初值选取x0 迭代过程往往就在根的附近进行，只要假定 在 附近连续，且满足 ，则根据连续函数的性质，一定存在 的某邻域S： ，使得 在S上满足定理2.1的条件，故在S中任取初始值x0，迭代公式 必将收敛于方程 的根 这种收敛称为局部收敛。 (1)输入初始近似值x0，精度要求eps，控制最大迭代次数M； (3)当k M且 eps时做循环 ， ， (4)如果 eps，则输出x1；否则输出迭代失败信息。 (2) x y y = x x* y=g(x) x0 p0 x1 p1 ? x y y = x x* y=g(x) x0 p0 x1 p1 ? ② 知L越小，收敛速度越快 加快收敛速度的办法就是设法让L变小 变形为： 相当于 根据迭代收敛定理，导数要尽可能