2017-2018年度第一学期 00151102 00151103.PDF

2017-2018年度第一学期00151103 计算方法（B） 童 伟 华 管理科研楼1205 室 E-mail: tongwh@ 中国科学技术大学 数学科学学院 / 第五章 解线性方程组的直接法 2 线性方程组的求解  在科学研究和工程应用中，求解线性方程组是非常基 础的问题  一般的线性方程组 a x ++a x b  11 1 1n n 1    ⇔ Ax = b  a x ++a x b  n 1 1 nn n n  （Gram公式）当且仅当 时，方程组有唯一解 det(A) ≠0 D xi i , i 1,2,,n D a  a i b a i  a n  11 1 −1 1 1 +1 1   D det(A), Di det        a  a b a  a   n 1 ni−1 n ni+1 nn  3 线性方程组的求解  求解线性方程组的方法可分为  直接解法：采用消元（初等变换）、矩阵分解等技巧；从理论上来说， 直接法经过有限次四则运算（假设运算无舍入误差） ，可以得到线性方 程组的精确解  迭代解法：采用迭代、分块、预条件处理等技巧；将线性方程的解视为 某种极限过程的向量序列，近似解  直接法与迭代法的使用建议  一般来说，对同等规模的线性方程组，直接法对计算机的要求高于迭代 法  对于中等规模的线性方程组，考虑到直接法的准确性与可靠性，建议使 用直接法求解  对于高阶稀疏线性方程组（非零元素较少），建议使用迭代法求解  要充分考虑线性方程的特性，譬如对称性、正定性、稀疏性等，来选用 合适的算法 4 线性方程组的求解  求解线性方程组的常用软件  Matlab  LAPACK/EISPACK/LINPACK （一般的线性方程组求解）  TACUS （大型稀疏线性方程组）  SuperLU （大型稀疏线性方程组）  Eigen （开源，C++ ，线性代数模板库）  MKL （商业软件）  …  推荐参考书  G. H. Golub, C. F. V. Loan. Matrix Computations.