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

§1 引言 学习线性方程组数值解法的必要性 科学计算中经常遇到线性方程组求解问题 如电路分析、分子结构、测量学、运筹学、流体力学、数值逼近及微分方程的数值解法等 当线性方程组的阶数较大时，人工求解已不可能。 当求解方法不得当时，即使计算机求解都很难实现。 对20阶的线性方程组，用Cramer法则求解，乘除法的运算次数达到9.7×1020，若用每秒钟一亿次的计算机计算也要30万年；而用Gauss消去法求解，则只需乘除法次数为3060次，不需1秒钟就可计算出来。 线性方程组求解方法的分类 直接法、间接法 矩阵的特征值和谱半径 特征值的性质及特征多项式 特征值的有关性质 AT和A具有相同的特征向量及特征值； 若A非奇异，则A－1与A的特征值互为倒数，特征向量相同； 相似矩阵具有相同的特征值。 特征多项式 对称正定矩阵 定义 性质 非奇异，且其逆矩阵也对称正定； 所有特征值大于零； 所有对角元也大于零； 所有顺序主子式都大于零。 初等矩阵 定义 性质 初等置换矩阵 初等下三角阵（Gauss变换矩阵） Household变换（初等反射矩阵） 定义 性质 例题：用Gauss消去法解方程组并求其系数矩阵行列式的值。 消去法与三角分解 Guass消元的过程从矩阵变换角度出发，实质上是进行了(n-1)次Gauss变换，即 三. 矩阵的直接三角分解法 当矩阵A非奇异且顺序主子式 ，则A可做LU分解，即A=LU. 这时方程组Ax=b可以转化为求解LUx=b，等价于以下两个三角形方程组 Ly＝b，及 Ux＝y. 现在通过矩阵乘法直接求得L及U的元素。 对比两边元素，直接可得 若己经求得U的第r-1行及L的第r-1列，则由矩阵乘法可得 即 从而可以求得U的第r行元素。L的第r列元素可由 若 ,可得 求解三角形方程组的算法： 对下三角方程组Ly=b,有 对上三角方程组Ux=y,有 例5 利用直接三角分解法求下列方程组的解 四、改进的平方根法 对称正定矩阵A可分解为 A=LDLT 其中 改进平方根法求解对称正定方程组 现设定理对n-1阶满足条件（5.2)的三对角阵成立，求证对满足条件(5.2)的n阶三对角阵定理亦成立。 因 ， 消去法执行一步后得 显然 其中 于是由归纳法假设，则有 定理4 设 ，其中A为满足条件（5.2）的三 对角阵，则A的所有顺序主子式都不为零，即 故 证明 由于A是满足（5.2）的n阵三对角阵，因此 ，A的任一个顺序主子阵亦是满足（5.2） 的三对角阵，由定理3，则有： 即 于是，由矩阵的三角分解定理，则有 由矩阵乘法，可得计算待定系数 的计算公式，即 于是，得到解（5.1)的追赶法公式： (1) 分解计算公式 A=LU (2) 求解 的递推公式 (3) 求解 的递推公式 将计算 的过程称为追的过程,计算方程组解 的过程称为赶的过程. 追赶法解 仅需要5n-4次乘除运算。只需要用3个一维数组分别存贮A的系数 且还需要用三个一维数组保存计算的中间结果 例5 用追赶法解方程组 解 (1) 计算 (2) 计算 (3) 求解计算 §6 解对称正定方程组的平方根法 一、引言 在工程技术问题中,例如用有限元方法解结构力学中问题时,常常需要求解具有对称正定矩阵方程组,对于这种具有特殊性质系数矩阵,利用矩阵的三角分解法求解就得到解对称正定矩阵方程组的平方根法，这是一种解对称正定矩阵方程组的有效方法。 二、对称正定矩阵 设有方程组 (6.1) 其中