计算方法电子教案 第六章解线性代数方程组及直接法.pptVIP

改进平方根法解对称正定方程组的算法 令LTX = y，先解下三角形方程组LDY = b得 （二）矩阵范数 其中 解上三角形方程组LTX = Y得 6.5 误差分析 （一）向量范数 向量范数是n维Euclid空间中长度概念的推广，其定义如下： 定义 对任一向量 ,按照一定规则确定一个实数与对应，该实数记为 满足下面三个性质。 （1） ，当且仅当X=0; （2）对任意实数a， ； （3）对任意向量 , 那么称该实数 为向量的范数。 在 中，常用的几种范数有 其中 分别是X的n个分量。 定义 设A是n × n阶的矩阵，定义‖A‖为矩阵的范数。 矩阵范数具有下面的性质： （1）对任意非零矩阵A，‖A‖恒为正数，当且仅当矩阵为零时，其范数为零。 （2）对任意实数 （3）对任意两个阶相同的矩阵A，B有 （4）对于矩阵A具有相同维数的向量X，恒有 （5）对同阶矩阵A,B有 前三种矩阵范数是： （ 是A*A的最大特征值） 定义 称Rn中的向量序列 收敛于Rn中的向量X，如果 定义 称nxn阶矩阵序列 收敛于nxn阶矩阵A，如果 其中 是nxn阶矩阵范数 范数的等价性可以推出，在某一种范数意义下向量序列收敛，则在任何一种范数意义下该向量序列收敛。 把向量序列 收敛于向量X，记为 把矩阵序序 收敛到矩阵A，记为 定理6.2 Rn中的向量序列 收敛于Rn中的向量X的充分必要条件是 其中 分别是 中的第j个分量。 证明:由范数的等价性，结知，只要对某一种范数证明定理即可，即 ，结论是明显的。 定理说明，Rn空间中向量序列的收敛可以归结为各坐标分量元素的收敛。 定理 nxn阶矩阵序列 收敛于nxn阶矩阵A的必要充分条件是 其中 分别是 在（i,j）位置上元素。 谱半径：设nxn阶矩阵A的特征值为 则称 为矩阵的谱半径。 谱半径和矩阵范数有关系： 证明：对矩阵A的任一特征值 ，与其对应的特征向量x1，有关系 ，两端取范数，再利用矩阵范数性质有 由于 ，则 ，所以 ，故 定理6.4设A是对任意nxn阶矩阵，由A的各次幂所组成的矩阵列 收敛于零，即 的必要充分条件是 * * 第六章 解线性代数方程组的直接法 本章研究的对象是 阶线性代数方程组： 用矩阵和向量的记号来表示： 其中 若 非奇异，方程组有唯一解 6.1 高斯消去法 一、过程（消元过程） 二、回代过程（求解过程） 计算过程： 利用 去乘以第1个方程后，第2个方程减去第1个方程， 消去第一项，变成： 取 乘以第1个方程， 利用第3个方程减去第1个方程消去 如此下去得： 用矩阵符号表示为： 令 则系数矩阵 和右端项 中的元素计算公式 可以写成： 类似地，若 ， 取 乘以第2个方程，再用第3个方程 减去第2个方程得： 如此下去有： 用矩阵 令 ,则系数矩阵 重复 次后，即可得到等价的三角形方程组 用矩阵表示： 消元过程中元素的计算公式可以归纳为: 回代过程： 即可以求得方程组的解 三、高斯消去法的运算量 只计算乘除法的运算量， 在消去第一列的 个系数， 需作乘 法运算 次 ,消去第二列需作 次 乘法运算，总起来， 把 化成 所需要的乘法运算的量为： 要做的除法运算量是： 回代需要乘除法总运算量是： 高斯消去法解 阶线性方程 组总共需要的乘除法运算是： （五）高斯消去法的矩阵解释与 分解， 高斯消去法的消元过程 ,从矩阵运算的角度看， 其中 即 所以 步以后有： 因此 令 我们有 其中 有形式： 6.2 主元素消去法 先看一个例子： 准确解为： 现在用高斯消去法求解， 取五位有效数字， 把第2个方程中的 x1消去，得到 从第二个方程解出x2=0.6667, 再代入第一个方程得 这时求得x1=0 为了控制舍入误差， 我们采用另一种消去过程， 先将方程组中 的两个方程交换一下，变为： 消去第二个方程中的x1，得 从而求得x2=0.6667, x1=0.3333 结果与准确解十分相近， 这个例子告诉我们控制舍入误差是十 分必要的。 主元素消去法的思想： 在消去过程中， 作为除数的有 , ,…， 。 我们在系数矩阵中选取绝对值最大的数作为 , ,…， 。 例2 解方程组 取四位有效数字，它的准确