数值分析课件 第5 解线性方程组的直接法.docVIP

第五章 线性代数方程组的数值解法 线性方程求解问题是科学研究和工程计算中最常见的问题。如电学中的网络问题、工程力学中求解连续力学体（微分方程）问题的差分方法、有限元法、边界元法及函数的样条插值、最小二乘拟合等，都包含了解线性方程组问题。因此，线性方程组的解法在数值计算中占有极其重要的地位。 对于阶线性方程组，若，则方程组有惟一解。由克莱姆（Cramer）法则，其解为,其中为用向量代替中第列向量所得矩阵。每个阶行列式共有项，每项都有个因子，所以计算一个阶行列式需做次乘法，我们共需要计算个行列式，要计算出，还要做次除法，因此用Cramer法则求解要做次乘除法（不计加减法），计算量十分惊人。如时，就需作约次乘法。可见Cramer法则在理论上是绝对正确的，但当较大时，在实际计算中却是不可行的。因此寻求有效的数值计算方法就成为非常必要的课题。 线性方程组的类型很多，若按其系数矩阵阶数的高低和含零元素多少，大致可分为两类：一类是低阶稠密线性方程组，即系数矩阵阶数不高，含零元素很少。另一类是高阶稀疏线性方程组，即系数矩阵阶数高，零元素占绝对优势（比如占以上）。 线性方程组的数值解法也可分为两大类：直接法和迭代法。 直接法是在没有舍入误差的情况下，通过有限步运算可以得到方程组精确解的方法。但是，在实际计算时，由于初始数据变为机器数而产生的误差以及计算过程中所产生的舍入误差等都要对解的精确度产生