第三章线性代数方程组的数值解探究.pptVIP

一般地说，差分格式被写成如下的形式： 方程组（3.1）可被进一步写成矩阵的形式 : 解线性代数方程组的直接法。 大家知道，线性方程组 显然，按Cramer法则求解方程组（3.3）需要计算n+1个n阶行列式。每个n阶行列式按直接展开办法来算，需作（n-1）×n！次乘法和n!次加法运算。当n=30时，共约需完成 次乘法和加法运算，这是一个十分惊人的数字，即使在一台每秒作一亿次运算的计算机上完成这一计算也是不可能的。所以，尽管这种办法也是一种直接法，并且理论上可行，但实际上是无法进行求解的。即使采用其它办法来计算行列式，按Cramer法则求解的工作量也比通常的直接法大得多。因而，Cramer法则对于数值计算来说是没有什么用处的，仅在一些特殊场合才有用。 另外，矩阵求逆也是人们熟悉的求解线性方程的一种方法。 现在计算机上常用的直接解法大多数是以系数矩阵的三角形化为基础的。就是说，先对方程组进行变换，使其化为等价的（即具有相同解的）三角形方程组。由于三角形方程组的求解十分容易，原方程的求解问题即告解决。 §2 三角形方程组的解法 所谓三角形方程组是指下面两种形式的方程组： 若用矩阵符号可分别写为： 三角形方程组的求解是很简单的。方程组(3.6)的计算公式可归结为： 对于方程组(3.7)，其计算公式可归结为： §3 Gauss消去法 Ga