第2章非线性方程(组)的数值解法1解说.ppt

通常将方程f(x)=0化为与它同解的方程 的方法不止一种,有的收敛,有的不收敛,这取决于 的性态. 方程 的求根问题在几何上就是确定曲线y= 与直线y=x的交点P的横坐标(如图所示) (a) (b) 2、迭代法的几何意义 数值分析——解非线性方程（组） 迭代法的几何意义 数值分析——解非线性方程（组） 例4 三种迭代结果见表2.为什么（1）（2）发散，（3）收敛 取 正根为 数值分析——解非线性方程（组） 4、构造迭代函数 的方法 （情形1） （情形2） （情形3） 0 1 2 3 2.0 3.0 9.0 87.0 2.0 1.5 2.0 1.5 2.0 1.75 1.732143 1.732051 表2-3 的迭代例子 问题：如何构造 ，才能使迭代序列 一定收敛于不动点？近似解的误差怎样估计？ 设迭代公式为 即 而且序列 二、不动点迭代法收敛性 满足微分中值定理条件时，有 当 显然只要 时，式（*）成立。 收敛于不动点 因而有 x *， （*） 分析： 定理3（映内压缩性定理）（收敛的充分性条件） 设方程 在 上存在唯一解, 是迭代函数，如果 1、映内性： 对任何 称正数 则迭代公式 收敛于方程 在 上的唯一根。 2、压缩性： 为压缩因子 且 对任意的初值 并有误差估计式 (2.10) 证明 收敛性是显然的.下面证明误差估计式，因为 据此递推，可得 于是对任意正整数,有 在上式中令 即得(2.10)式，定理证毕. 注意到 对任意正整数，有 令,则 有 【注】 （1） 事前误差估计：对于收敛的迭代序列， （2） 控制迭代次数：根据事前误差估计公式（2.10）， 公式（2.10）右端项可用于误差上限的估计. 其中， 为取定的迭代次数，压缩因子取为 为使近似解达到精度要求，由 （2.11） 即用前后两次迭代结果之差的绝对值是否小于允许误差 来判断迭代是否终止；也可用于误差上限的事后估计. （3） 事后误差估计：由公式（2.11）， 可作为迭代算法的终止条件， 可得所需要的最少迭代次数 k ，即 迭代法的算法框图 例5 求方程 的正根 因为是求正根，所以不考虑 所以在 解： 解得驻点 令 （1）确定有根区间 内有正根 2） 选取恰当的迭代函数，判别根的收敛性. 取 由于 所以在 内有正根 映内性 压缩性 因此迭代公式 （3）迭代计算.迭代公式： 必收敛 取初值为 由表2-4可知，x7已达到6位有效数字，可取1.32472作为所求根的近似值. 数值结果如表2-4. 对于任意的值 1.32472 1.32472 1.32473 1.32476 1.32494 1.32588 1.33086 1.35721 xk 8 7 6 5 4 3 2 1 k 表2-4 例5迭代值 数值分析——解非线性方程（组） 若将原方程改写为 ，用迭代公式 计算，则迭代过程是发散的. 四、局部收敛性与收敛阶 1、全局收敛（大范围收敛）若对 中的任意一点 作初始值，迭代均收敛，这种形式的收敛称为全局收敛或 大范围收敛. 2、局部收敛：若在 的某个邻域 内的每一点，迭代均收敛，称这种形式的收敛为局部收敛. 【定理4】 （局部收敛性）设 在 的邻近连续， 且 则迭代过程 在 邻近具有局部收敛性. 数值分析——解非线性方程（组） 证明 因为 连续，所以存在 的一个 邻域 ： ，使得对任意 ,有 成立，并有 即对任意 ,有 ,因此 满足定理3的映内压缩性条件,从而,对任意 ,迭代过程局部收敛. 第 2 章 非线性方程（组）的数值解法 基础教学部数学教研室 彭 晓 华 “方程是很多工程和科学工作的发动机”.非线性现象广泛存在于物质世界与社会生活中.在工程和科学计算中，常涉及到非线性方程或非线性方程组的求解问题.例如： 2.1 引 言 (1) 在光的衍射理论（the theory of diffraction of light)中,我们需要求 x – tan x = 0 的根 (2) 在行星轨道( planetary orbits)的计算中,对任意的a和b，我们需要求 x – a sin x = b 的根 *数值求解方程组的必要性 (3) 在数学中，需要求n次多项式 （4）在天体力学中，有如下开普勒（Kepler）方程 其中 表示时间， 表示弧度，行星运动的轨道 是 的函数.也就是说，对每个时刻 ，上述方程 （