数值分析英文版课件.pptVIP

数值方法 主讲：郭同彤 哈尔滨工业大学深圳研究生院 基础科学学科部 今日主题 第三章：非线性方程的数值解法 3.1 引言 3.2 二分法和试位法 3.3 不动点迭代法 3.4 迭代加速收敛的方法 3.5 Newton 迭代法 今日主题 第三章：非线性方程的数值解法 3.1 引言 3.2 二分法和试位法 3.3 不动点迭代法 3.4 迭代加速收敛的方法 3.5 Newton 迭代法 3.1 引言 (1) 本章主要讨论数值求解方程 f(x)=0 （4.1.1） 其中x?R, f ?C[a,b]。 根和零点的定义。 如果函数f(x)是多项式函数，即 3.1 引言 (2) 设f(x)可分解为 今日主题 第三章：非线性方程的数值解法 3.1 引言 3.2 二分法和试位法 3.3 不动点迭代法 3.4 迭代加速收敛的方法 3.5 Newton 迭代法 3.2.1 二分法 (1) 假设 已找到方程 f(x)=0 的一个有根区间 [a,b]； f(a)f(b)0； 方程在区间 [a,b] 只有一个根。 二分法步骤： 令 [a1,b1]=[a,b] , 执行以下迭代步骤： 对于区间 [an,bn] , 其中点为 xn=1/2(an+bn)； 若 f(an)f(xn)0 , 则将 [an+1,bn+1] 替换为 [an,xn]； 若 f(an)f(xn) 0 , 则将 [an+1,bn+1] 替换为 [xn,bn]。 3.2.1 二分法 (2) 3.2.1 二分法 (3) 3.2.1 二分法 (4) 例 4.2.2 用二分法求解 f(x)=sinx-x2/4，要求误差小于 10-2，需要用二分法计算多少次？ 解法1 f(x)=0等价于求两条曲线的交点 取区间[a,b]=[1.5,2] 应用公式 计算对分5次即可，x5≈1.93。 3.2.1 二分法 (5) 解法2 f(p/2)0, f(p)0 f(x)=cosx- p/20, x ?[p/2, p] 因此，f(x) 在 [p/2, p] 内为单调下降函数，有唯一的实根 可取[a,b]= [p/2, p] 类似解法1进行求解。 3.2.2 试位法 今日主题 第三章：非线性方程的数值解法 3.1 引言 3.2 二分法和试位法 3.3 不动点迭代法 3.4 迭代加速收敛的方法 3.5 Newton 迭代法 3.3.1 不动点和不动点的迭代 (1) 为解方程 f(x)=0 (4.3.1) 将其变换为等价的方程 x = f (x) (4.3.2) 其中?是连续函数。构造迭代公式： x k+1 = f (xk) (4.3.3) 如果 则称?为迭代函数， x*是函数 ?的一个不动点，也就是方程的根。此迭代法称为不动点迭代法。 3.3.1 不动点和不动点的迭代 (2) 可以通过不同的途径将方程 f(x)=0 变换成方程 x = f (x) 的形式。 例如， ?（x)＝x - f(x) ?（x) ＝ x - Af(x), 其中A为常数， 也可以用其他的方法。 3.3.1 不动点和不动点的迭代 (3) 设序列 {xk} 收敛到 x*，记误差 ek=xk-x*. 若存在实数p≥1及非零常数C，使： 则称 {xk} 为 p 阶收敛，C称为渐进误差常数。 当 p=1 时，称 {xk} 线性收敛 当 p1 时，称 {xk} 超线性收敛 当 p=2 时，称 {xk} 平方收敛 今日主题 第三章：非线性方程的数值解法 3.1 引言 3.2 二分法和试位法 3.3 不动点迭代法 3.4 迭代加速收敛的方法 3.5 Newton 迭代法 3.4.1 Aitken加速方法 (1) 设 {xk} 线性收敛到 x*，记 ek = xk – x* ，有 当ｋ充分大时有 其中 |c| =C ，由 可解出 3.4.1 Aitken加速方法 (2) 定义差分记号， 写成 它是x*的一个新的近似值 从序列｛xk}, 用上式得到序列 { } 的方法，称为Ａitken加速方法。 3.4.1 Aitken加速方法 (3) 可以证明，只要{xk}满足 　　　　　　　 且 则序列 { }是完全确定的，而且有 说明 ? 今日主题 第三章：非线性方程和方程组的数值解法 3.1 引言 3.2 二分法和试位法 3.3 不动点迭代法 3.4 迭代加速收敛的方法 3.5 Newton 迭代法 3.5.1 Newton迭代法公式和收敛性 (1) 为求解方程 f(x)=0 的根 x＊，假设 有一个近似值　xk ≈ x＊ f ’’存在且连续 因 f (x*)=0, 则： 若