[工学]数值分析+-lec.pptVIP

常用的求非线性方程根的方法 二分法（对分法，搜索法） 不动点法 （简单迭代法、压缩映象法）及其加速算法 Newton方法及其变体 二分法（对分法） 理论依据：设函数 f(x)∈C[a,b], f(a)f(b)0，则 f(x)在区间[a,b]上有一实根s使 f(s)=0 二分法工作原理：每次对半分根所在的区间[ak, bk]成两个等长的子区间[ak, 0.5(ak+ bk) ] 和 [0.5(ak+ bk) ，bk]，并剔除根所不在的区间，直至根所在的区间变得充分小！，即|bk-ak|ε。 算法 : Set a0=a, b0=b, given ε, and maximal iteration number N; computing f(a) and f(b). If f(a)f(b)0, stop and print out message of computational failure. If |f(a)| ε stop and print the root, which is a. If |f(b)| εstop and print the root, which is b; Starting the iteration loop, k=0,1,… Computing, x0=0.5(a0+b0) and f(x0). If |f(x0)| ε, stop and print the root, which is x0. Otherwise, if f(a0)f(x0)0, set b0=x0; otherwise, set a0=x0; Set k=k+1; If b0-a0 ε and kN, go back step 3; If kN, print out message of failure and stop. 迭代过程的几何表示 Newton 迭代法 Newton 法的几何解释 作业 课后习题1、2、4、5、7、9 * 割线法的几何表示 x0 X x* x1 x2 x3 Y f(x)0 P0 P2 P1 割线法的收敛性 引理：设f(s)=0，在 s 的某领域 [s- δ, s+δ] 内 f″(x) 连续，f′(x)≠0，又设 xk-1，xk∈[s-δ, s+δ] 且 xk-1, xk, s 互异，记 ek=s-xk, 则有 其中 xk＋1 是由割线法产生，ηk, ξk 在 min(xk-1, xk, s) 与 max(xk-1, xk, s) 之间。 定理：设 f(s)=0,在 s 的某领域 Iδ= [s- δ, s+δ] 内 f″(x) 连续，f′(x)≠0，则存在ε0，当 x-1, x0∈Iδ时，则由割线法产生的序列 {xk} 收敛于 s ，且收敛速度的阶至少为 1.618。 数值分析 朱立永 北京航空航天大学 数学与系统科学学院 Email: numerical_analysis@ Password:beihang 答疑时间：星期三下午2:00－5:00 答疑地点：主216 第十讲非线性方程的迭代解法 第三章非线性方程与非线性方程组的迭代解法 这一部分的主要任务是解 其中f（x）是一个一元非线性函数。 非线性科学是当今科学发展的一个重要研究方向，很多实际工程物理问题都归结为非线性方程（组）的求解。 例如： 无穷组解 无解 一个解 两个解 四个解 非线性方程的求根非常复杂。 求根问题包括下面三个问题： 根的存在性：即f(x)=0有没有根？若有， 有几个根？ 哪儿有根？确定有根区间 根的近似求解 根，重根的定义： 根：如果存在常数s，使得 f(s)=0则称s是f(x)=0的根（零点）； 重根：如果 ， 称s为m重根。特别地对 f(x) 是多项式，则有 简单迭代法及其收敛速度 迭代法的构想 从一个初值x0出发，计算 如果 {xk}收敛，即存在x*, 使得 则由 得 即 x*是 φ(x)的不动点，也就是 f(x) 的根。 O x* x2 x1 x0 x y 问题： {xk} 收敛吗？怎样实现φ(x) 简单迭代法的收敛速度 Steffensen加速收敛方法 非线性问题的最简单解法是线性近似. 将非线性方程线性化，以线性方程的解逐步逼 近非线性方程的解，这就是Newton法的基本思想 Newto