数值分析第二章非线性方程解法.pptVIP

* 定义1(初等变换) 第二章 非线性方程的解法 第一节 引言 对于方程f ( )=0,如果存在 使f ( )=0,则 是f ( )=0的根, 或称函数f ( )的零点.本章都把方程的精确解记为 . 一般说来非线性方程没有求根公式,所以本章我们采用的都是迭代法,即用某种方法构造一个近似解序列 ,并使之收敛于精确解 . 当然我们还要考虑它的收敛速度,在同样的精度要求下,收敛速度越快意味着计算量越小. * 第二节 二分法 矩阵的等价 在[a,b]上连续; 且f (a)f (b) 0 如果f (x)满足: 那么方程f (x)=0在[a,b]上至少有一个实根 把区间逐次对分,使包含 的区间长度逐次缩小直到 算法: 令 判别: 再令: 由于 就是满足精度要求的近似解 * 二分法的精度控制 矩阵的等价 二分法的特点: 1.条件宽松,计算简单; 2.全局收敛; 3.在x*的附近收敛慢. 事后估计:每计算一次判别是否成立 一旦满足,计算即告结束. 事先估计:由 二分法的一种改进称0.618法: 这样可以事先确定计算次数k. * 第三节 简单迭代法 基本思想:把方程f (x)=0转化为一种等价的形式x = g (x) 注意：由f (x)=0形式转化为x=g (x)的形式是不唯一的. 显然:这个序列的产生取决于函数g (x)的选择,g (x)称为迭代函数。 序列的收敛与否以及收敛的速度也取决于迭代函数的性质. 行最简形,标准形,等价类 定理(收敛基本定理):对于方程x=g (x),如果满足: * 第三章 矩阵的初等变换与线性方程组 一,矩阵的初等变换 。 几点注意:1.这个定理的条件是收敛的充分条件而非充要条件; 2.收敛不仅与迭代函数g(x)有关,且L越小收敛越快,另外它还与所 给的初值有关,即它是局部收敛的. 3.由于 的条件判别困难,所以常用条件 来代替。 4.精度控制一般采用事后估计的方法,以 作为停机准则. 收敛速度是衡量算法优劣的一个重要标准,我们引进刻画收敛速度的一个量. 有时就直接称该迭代函数或该迭代法是p阶收敛. 当P=1,0c1时,称线性收敛，这时称c是收敛率。 * 把一阶收敛和二阶收敛做一个比较： 设有两个收敛于同一个x*的序列，它们有相同的初始误差1和收敛常数c=0.75,p分别等于1和2。 注意:二阶收敛对初值有更高的要求. * 二,初等矩阵 定义4(初等矩阵) 收敛阶的判别 定理:对于方程x = g (x),如果满足: (1) g (x)在x*附近有连续的m阶导数; (2) 则该迭代是m阶收敛. 迭代的加速收敛方法(Atiken方法) 对于线性收敛的序列,可以通过简单的处理达到加速收敛的效果 * 　 把上述过程写为迭代形式就称为Aitken方法: 在线性收敛的情况下,若收敛率接近于1时,上述方法加速收敛的效果很好. * 第四节 牛顿法 基本思想:把非线性方程线性化,用线性方程的解来逐步逼近线性方程的解. 对方程f (x)=0,若f (x)的二阶导数连续,则对f (x)作二阶Tayor展开 如此重复上述线性化过程,可得: 这个算法称牛顿法. 牛顿法也是一种迭代法,其迭代函数为: * Th1 下面的定理给出了牛顿法的收敛特性 定理 设f (x)有二阶连续的导数,则牛顿法在单根附近是二阶收敛. 所以牛顿法在单根附近是二阶收敛的. * Th5及推论 与二分法和简单迭代法比较,牛顿法有如下特点: 关于牛顿法的几点说明; 1.二分法只要求函数连续,简单迭代法要求迭代函数的一阶导 数连续,而牛顿法则要求函数的二阶导数连续; 2.每次迭代的计算,二分法只需要判别函数值的正负,简单迭代法要计算一次函数值,牛顿法除了函数值以外还需计算导数值; 3.牛顿法可求复数根(简单迭代法也可求复根),且在单根附近是二阶收敛(但对初值的要求也相应提高); 4.基于牛顿法的基本原理和特点,容易作相应的改进和推广. 注意:牛顿法在重根的情况下依然可以使用,但收敛是 一阶的,且重数越大收敛越慢. * 利用初等变换求逆矩阵 关于重根的处理 设x*是方程的m重根,可以用以下两种改进牛顿法来计算 显然两个方程 有相同的解 这