数值研究课件解线性方程组的直接法.docVIP

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