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

第2章 解线性方程组的直接法 2.1 GAUSS消元法 2.2 改进平方根法 2.3 追赶法 2.4 LU分解法 2.5 直接法的稳定性分析 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列主元法 方法要点 在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消元法 2.1.5 GAUSS列主元法的计算步骤 1. 输入 2. 输出 3. 计算：(1)消元，(2)回代 2.1.6 GAUSS列主元法的计算实例 [例] 求解方程组 2.1 GAUSS消元法 2.1.7 GAUSS全主元素法 [全主元素] 在GAUSS消元法的第i轮消元前，从子矩阵Ain中按绝对值最大的原则所找到的元素称为全主元素。 [全主元消元法] 找全主元素 ? 作列和行的交换 ? 消元 ? 回代 2.4 LU分解法 从GAUSS消元法可知，系数矩阵A消元后，形成了一个上三角矩阵和一个下三角矩阵，因此可重点研究系数矩阵A的分解方法。 LU分解法： 2.4.1 LU分解法的算法推导 对线性方程组Ax=b，令A=LU，其中 2.4 LU分解法 依据的原理：矩阵的乘积法则 2.4 LU分解法 LU分解公式(计算顺序)的推导： 2.4 LU分解法 2.4.2 LU分解法的求解公式 由 Ax=b 及 A=LU，有 LUx=b ; 令Ux=y， 有 Ly=b , 可解得： 再由 Ux=y 可解得： 2.4 LU分解法 2.4.3 LU分解法的计算步骤 1.输入; 2.输出; 3.计算过程：作LU分解，计算yi，计算xi 2.4.4 LU分解法的计算实例 [例1] 对矩阵A作LU分解 [例2] 用LU分解法求线性方程组 2.2 改进平方根法 适用情况：系数矩阵A非奇异、是正定矩阵 2.2.1 正定矩阵的定义和性质 正定矩阵的定义 设A∈Rn×n，若有： (1) AT=A, (2) 对任意x∈Rn(x≠0)都有xTAx0 则称A是正定对称矩阵。 正定矩阵的性质 (简记: 二个“正定”，三个“0”) A的各阶主子矩阵也是正定对称矩阵 (正定) A非奇异，且A-1也是正定对称矩阵 (正定) A的对角线元素aii0 (0) A的所有特征值0 (0) A的各阶主子式0 (0) 2.2 改进平方根法 2.2.2 改进平方根法的算式推导 方法要点: 对矩阵A作A=LDLT分解，再求解 算法组成: 1. 三角分解：A=LDLT 2.2 改进平方根法 [计算顺序的推导] 2.2 改进平方根法 [计算公式的推导] 2.2 改进平方根法 [计算公式的推导