第五章线性方程组的直接解法讲解.pptVIP

第五章 解线性代数方程组的直接方法 第一节 求解线性代数方程组的基本定理 第二节 高斯消元法及其计算机实现 一、三角形方程组的解法 二、高斯消元法 再交换第k,ik两个方程和第k,jk列,使 成为主元. 称这个过程为完全选主元.  不论是哪种方式选出主元,而后再按上面介绍的计 算步骤进行消元的计算,一般都称为选主元的高斯消元法.在实际计算中,常用按列选主元的高斯消元法. 高斯消去法求解线性方程组： 下面导出 L, U 的元素, 因为 例 选主元的三角分解法 例 列选主元三角分解法 三对角方程组与追赶法 若A=LU ,则 对称正定方程组的平方根法 改进的平方根法 二、矩阵的范数 三、谱半径 谱范数 四、方程组的状态,条件数 条件数计算举例 1--范数: 2--范数: --范数: (Euclid欧氏-范数) 例 设 计算三种范数. 解 一、向量的范数 第五节 向量和矩阵的范数 (3) 三角不等式, 则称 为方阵 的范数. (1) ,等号成立当且仅当 (2) 对任意实数 ,有 (4) 定义 对任意n阶方阵A ,若对应一个非负实数 , 满足条件 为方阵 的 范数. (1) 列范数 (2)(欧氏--范数) 定理 设 ,则 (3) 行范数 可以验证该范数满足定义4的四个条件. 定义 设 定义非负实数 例 解 为方阵 的谱半径. 定理 对于 均有 称 为 的谱范数. 若 为实对称矩阵, 定义 设n阶方阵A的特征值为 ， 称 方程组 Ax=b 的系数矩阵和常数项数据,因为观测或实验,或离散化得到,在计算机上会进行舍入,可能会有误差,需要考虑系数矩阵和常数项有微小变化时,解的变化情况. 例 考察方程组 一个解为 病态方程组 上述方程组尽管只是常数项有微小扰动，但解 大不相同： 一个解为 假设求解是在四位浮点十进制数的计算机上进行 0.0001x1+x2=1 x1+x2=2 将两个方程对调，得 x1+x2=2 0.0001x1+x2=1 在四位浮点十进制数的计算机上,上式为 x1+x2=2 即 x1+x2=2 (0.1000×101-0.00001 ×101x2=1 x2=1 （1-0.0001) x2=1 x1+x2=2 消元，得 解得： x1=1， x2=1 现在我们再用列主元法解下列问题 例 用列主元消去法解方程组 解 第一次消元对 因列主元素为a31(1),故先作行交换E1 E3,然后进行消元计算可得 -0.002x1+2x2+2x3 =0.4 x1+0.78125x2 =1.3816 3.996x1+5.5625x2+4x3=7.4178 -0.002 2 2 0.4 [A(1) |b(1) ]= 1 0.78125 0 1.3816 3.996 5.5625 4 7.4178 3.996 5.5625 4 7.4178 [A(2) |b(2) ] = 0 -0.61077 -1.0010 -0.47471 0 2.0029 2.0020 0.40371 由此回代,得x=(1.9272,-0.69841,0.90038)T与精确解 x=(1.9273,-0.69850,0.90042)T相比较是比较准确的. 3.996 5.5625 4