chap02-解线性方程组及直接法.ppt

第2章 解线性方程组的直接法 2.1 GAUSS消元法 2.2 矩阵三角分解法 2.3 范数简介 2.4 直接法的稳定性分析 解线性方程组方法---概述 解Ax=b的数值算法： (1)直接法 (2)迭代法 解Ax=b的直接法 用于求解中小规模的线性方程组 包括三小类方法： Gauss消元法 三角分解法 正交变换法 2.1 GAUSS消元法 特点: 矩阵的初等变换 (是对中学代数中加减消元法、代入消元法的综合利用和升华) 2.1.1 基本GAUSS消元法 问题：求解n阶线性方程组Ax=b的方法 ? 解法：消元法(分2步) 1. 消元 (1) 目的：得到系数矩阵为上三角矩阵的方程组 (2) 方法：加减消元法 (3) 过程：利用aii消aji(j=i+1,…,n)，其余对应元素和bi作相应变动。 2.1 GAUSS消元法 (4) 计算公式推导： 假设已进行了i-1轮消元，方程组Ax=b已变成： 2.1 GAUSS消元法 相应要计算和变动的矩阵元素为： 2.1 GAUSS消元法 2. 回代求解xi： 2.1 GAUSS消元法 2.1.2 基本GAUSS消元法的计算实例 [例1-例4] 用GAUSS消元法求解方程组 2.1 GAUSS消元法 2.1.3 GAUSS消元法的计算量 1、消元计算量 除法次数： 乘法次数： 2、回代计算量 乘法次数： 2.1 GAUSS消元法 2.1.3 GAUSS列主元法 方法要点 在GAUSS消元法基础上，第i轮消元前，先选出列主元素，并将列主元素所在行与第i行进行交换。 [列主元素] 在GAUSS消元法第i轮消元前， 中绝对值最大的元素称为列主元素，记为： 2.1 GAUSS消元法 [求解过程] 1.消元（找列主元素?交换?消元） (1) 确定列主元素所在行 (记行号=p)。 方法--按定义，若列主元素为0，则选不出。 (2) 将列主元素所在的第p行与第i行元素交换 (3) 消元 2.1 GAUSS消元法 2. 回代求解xi： 2.1 GAUSS消元法 2.1.4 GAUSS列主元法的计算量 1. 消元过程的计算量 (1) 找列主元素：比较和绝对值计算 (2) 交换2行的次数 (3) 除法的次数： (4) 乘法和减法的次数： 2. 回代过程的计算量 除法次数 乘法、加法、减法的次数： 2.1 GAUSS消元法 GAUSS列主元法的计算实例 [例] 求解方程组 2.1 GAUSS消元法 2.1.5 GAUSS全主元素法 [全主元素] 在GAUSS消元法的第i轮消元前，从子矩阵Ain中按绝对值最大的原则所找到的元素称为全主元素。 [全主元消元法] 找全主元素 ? 作列和行的交换 ? 消元 ? 回代 2.1 GAUSS消元法 2.1.6 GAUSS列主元法和全主元素法的计算量 2.1.7 GAUSS全主元素法计算程序 2.1 GAUSS消元法 2.1.8 消元法的适用范围 GAUSS消元法能需满足的基本条件： [定义3] 矩阵A中前i行i列构成的矩阵称为i阶主子矩阵，对应行列式称为主子行列式。 [定理1] GAUSS消元法能进行到底的充要条件是各阶主子行列式都不为零。 2.2 矩阵三角分解法 从GAUSS消元法可知，系数矩阵A消元后，形成了一个上三角矩阵和一个下三角矩阵，可重点研究系数矩阵A的分解方法。 2.2.1 LU分解法 对线性方程组Ax=b，令：A=LU，其中 L---下三角矩阵 U---上三角矩阵 2.2.1 LU分解法 若L为单位下三角矩阵，则称为Doolittle分解；若U为单位上三角矩阵，则称为Crout分解。 Doolittle分解： [定理1] 若A各阶主子矩阵非奇异，则可进行LU分解。 2.2.1 LU分解法 依据的原理：矩阵的乘积法则 2.2.1 LU分解法 LU分解公式(计算顺序)的推导： 2.2.1 LU分解法 2.2.2 LU分解算例 [例1] 求下面矩阵的LU分解 2.2.1 LU分解法 2.2.3 利用LU分解法解线性方程组 由 Ax=b 及 A=LU，有 LUx=b ; 令Ux=y， 有 Ly=b , 可解得： 再由 Ux=y 可解得： 2.2.1 LU分解法 2.2.4 LU分解法的计算实例 [例2] 用LU分解法求线性方程组 2.2.5 平方根法和改进平方根法 1、现实背景 在许多工程物理问题所构成的线性方程组中，其系数矩阵A通常具有以下特点： 是正定的(对称) 对角严格占优 系数非零元素少，零元素多，且非零元素