数值分析第三章解线性方程组直接法.pptVIP

* 第三章 矩阵的初等变换与线性方程组 一,矩阵的初等变换 第三章 解线性方程组的直接法 第一节 引言 直接法的基本原理就是高斯消元法,再根据数值计算的特点做一些适当的处理而得到的一类算法.直接法的特点是没有截断误差,只有计算误差(舍入误差). 解线性方程组的方法可分为两大类:直接法和迭代法. 迭代法是类似于上一章单个方程那样,以某种方式构造一个向量序列,使得这个向量序列收敛到解向量.因此迭代法既有截断误差又有舍入误差. * 第二节 主元素消元法 矩阵的等价 在线性代数中我们通过矩阵的行变换用高斯消元法求解方程组 现在我们称之为简单消元法.其过程如下 * 消元过程完成后,可以通过回代来求出方程组的解 如果主对角线元素的绝对值很小,由第一章的误差分 析可知计算时将会导致很大的计算误差. 行最简形,标准形,等价类 以上是消元过程,其特点是每一步计算都要以主对角元素(称主 元素)作为除数. 回代过程中的每一步计算也都要以主对角元素作为除数. * 求解下列方程组 . 行最简形,标准形,等价类 * 二,初等矩阵 定义4(初等矩阵) 为了避免上述情况的发生,我们可以通过行交换,在需要消元的 这一列中把绝对值最大的元素 换到主对角线的位置 成为主元素 这个过程称为选主元,这种方法称为列主元素消元法. 通过选主元可以使计算结果大为改善. 　 * 对方程组 . 行最简形,标准形,等价类 还可以进一步用标度化的选主元(相对最大) * 第三节 矩阵的三角分解 消元法求解方程组是通过行初等变换把系数矩阵化为对角阵,由 线性代数知识可知,左乘一个初等矩阵,就相当于做一次行变换. 记 这正是第一次消元的结果, 称为消元矩阵,这是一个单位下三角矩阵 　 * 类似地再用消元矩阵 左乘 ,就实现了第二次消元。 　 　 就实现了第二次消元 * Th1 如此反复进行,就完成了消元过程 由于所有的消元矩阵都是单位下三角矩阵,所以其乘积还是单位下三角矩阵,且其逆还是单位下三角的,记之为L,因此A=LU, 这样就把方程组的系数矩阵A分解为下三角矩阵和上三角矩阵的乘积. 这称为矩阵的三角分解. * Th5及推论 通过矩阵的三角分解就可以把原方程组化为两个三角形方程组 求解这两个三角形方程组无需消元,只要回代就行了. 定理:如果方程组Ax=b无须经过行交换就可以用高斯简单消元 法求解( ),则A必可以做三角分解,且分解式唯一. 定义:如果n阶方阵A满足 则称A是严格主对角占优的. 定理:如果A是严格主对角占优的,则必有 用矩阵分解方法求解方程组计算量与消元法相同,但它可以保留 消元结果.这在工程上常常是很有用的. * Th1 由于 仅仅是在计算 时用到,用过以后不需要再保留, 所以我们可以由此节约储存单元.这称为矩阵分解的紧凑格式 由A=LU,根据矩阵乘法 再由的三角特征可知： * Th1 由此得算法： * 利用初等变换求逆矩阵 第四节 平方根法和追赶法 一, 平方根法 若A是正定矩阵(当然是对称),A=LU,现把U再做进一步分解 由于A是正定矩阵，所以D中的元素都大于零，可以把D也再分解 * 利用初等变换求逆矩阵 这就是平方根法的分解，它减少了几乎一半计算量。 还有一类特殊的矩阵称为带状矩阵，它可以分解为两个半带 状的矩阵。带状矩阵中一类更特殊的称为三对角矩阵，它有 一种特别简单的算法，称为追赶法。 * 例 由 类似的： * 利用初等变换求A-1B 这样不仅计算量很小而且只需3个向量的储存空间，分解后 解方程的过程也特别简单