数值分析第五章解线性方程组直接法.pdf

郑州大学研究生课程 （2011-2012学年第一学期） 数值分析 Numerical Analysis 第五章 解线性代数方程组的直接法 第五章 解线性代数方程组的直接法 §5.1 引言 §5.2 高斯消去法 §5.3 主元素高斯消去法 §5.4 三角分解法 §5.5 对称正定矩阵的平方根法 §5.6 三对角方程组的追赶法 §5.7 范数与误差分析 2/130 郑州大学研究生2011-2012学年课程 数值分析 Numerical Analysis §5.1 引言 自然科学和工程计算中的很多问题的解决常常归 结为求解线性方程组。如三次样条插值函数问题 、用最小二乘远离确定拟合曲线、求解微分方程 的数值解等，最终都要转化为求解线性方程组。 3/130 郑州大学研究生2011-2012学年课程 数值分析 Numerical Analysis §5.1 引言 从微观的薛定谔方程、分子动力学方程到宏观 的结构计算、工程力学中弹塑性方程、热弹性 方程组，数值求解方法包括差分法及有限元方 法等。这些离散方法最终的结果是常化为线性 方程组的求解。 4/130 郑州大学研究生2011-2012学年课程 数值分析 Numerical Analysis §5.1 引言 快速、高效地求解线性方程组是数值线性代数研究中 的核心问题，也是目前科学计算中的重大研究课题之一。 各种各样的科学和工程问题，往往最终都要归结为求 解一个线性方程组。 5/130 郑州大学研究生2011-2012学年课程 数值分析 Numerical Analysis §5.1 引言 系数矩阵的分类 第一类：低阶稠密方程组，即系数矩阵的阶数 不高，含零元素很少，在线性代数等课程学习 中通常见到的，都属这类方程组； 第二类：高阶稀疏方程组，系数矩阵的阶数很 高，如几百阶、甚至成千上万阶，其中零元素 成片分布，数量上绝对占优。 6/130 郑州大学研究生2011-2012学年课程 数值分析 Numerical Analysis §5.1 引言 线性方程组的数值解法有：直接法和迭代法。 直接法：在假定没有舍入误差的情况下，经过有限次 运算可以求得方程组的精确解； 迭代法：从一个初始向量出发，按照一定的迭代格 式，构造出一个趋向于真解的无穷序列。 7/130 郑州大学研究生2011-2012学年课程 数值分析 Numerical Analysis §5.1 引言 如果线性方程组的系数行列式不为零，即det(A) ≠0, 则该方程组有唯一解。由克莱姆(cramer)法则，其解为 det(A ) xi i (i 1, 2, , n) det(A) n 1 n n 这种方法需要计算 + 个 阶行列式并作 次除法，而每个 n阶行列式计算需作(n 1) n !次乘法，计算量十分惊人。 − × 如n 30,需2.38=×1035 次乘法。可见其在理论上是绝对正确， 但在n较大时，在实际计算中确实不可行的。 8/130 郑州大学研究生2011-2012学年课程 数值分析 Numerical Analysis §5.1 引言 例：直接法解线性方程组 ⎧⎪x1 −2x 2 +2x 3 =−2 2x 3x 3x 4