数值计算学习报告.doc.docVIP

《高斯消元法》研究与思考 第一章《高斯消元法》描述 1、《高斯消元法》的相关概念 在自然科学研究和工程技术中有许多问题可归结为求解线性代数方程组的问题,线性方程组求解是科学计算中最常遇到的问题。如在应力分析、电路分析、分子结构、测量学中都会遇到解线性方程组问题。在很多广泛应用的数学问题的数值方法中，如三次样条、最小二乘法、微分方程边值问题的差分法与有限元法也都涉及到求解线性方程组。直接法是在没有舍入误差的情况下，通过有限步四则运算来求的方程组精确解的方法。直接法基本方法是高斯消元法，其改进方法包括高斯列主元消去法，三角分解法，追赶法等的基本思想和原理。 高斯消去法（Gauss Elimination Method）是一种规则化的加减消元法。基本思想是通过逐次消元计算把需要求解的线性方程组转化为上三角形方程组，即把线性方程组的系数矩阵转化为上三角矩阵，从而使一般线性方程组的求解转化为等价(同解)的上三角形方程组的求解。 1.1 Gauss消去法的计算过程 通过一系列的加减消元运算，也就是代数中的加减消去法，以使A对角线以下的元素化为零，将方程组化为上三角矩阵；然后，再逐一回代求解出x向量。现举例说明如下： 1.1.1消元过程 第一步：将(1)/3使x1的系数化为1,再将(2)、(3)式中x1的系数都化为零，即由(2)-2×(1)(1)得 由(3)-4×(1)(1)得 第二步：将(2)(1)除以2/3，使x2系数化为1，得 再将(3)(1)式中x2系数化为零，由(3)(1)-(-14/3)*(2)(2) ，得 第三步：将(3)(2)除以18/3，使x3系数化为1，得 经消元后，得到如下三角代数方程组： 1.1.2回代过程 由(3)(3)得 x3=1,将x3代入(2)(2)得x2=-2,将x2 、x3代入(1)(1)得x2=1, 所以，本题解为[x]=[1,2,-1]T 1.1.3 用矩阵演示进行消元过程 第一步： 先将方程写成增广矩阵的形式 第二步：然后对矩阵进行初等行变换 第三步：将增广矩阵变换成上三角矩阵，主对角线全为1，左下三角矩阵全为0.即原方程组被等价转化成为上三角方程组，然后，逐步回代得原方程组的解即可。 1.1.4高斯消元的公式 综合以上讨论，不难看出，高斯消元法解方程组的公式为 第一步，消元 令 对k=1到n-1，若akk(k)≠0，进行 第二步，回代 若 2、高斯消元法的运算量 由公式 ，可得出 消去过程的第步共含有除法运算次，乘法和减法运算各次，所以消去过程共含有乘除法次数为 含加减法次数为 而回代过程含乘除法次数为，加减法次数为，所以Gauss消去法总的乘除法次数为，加减法次数为 《高斯消元法》的相关问题 为什么说高斯消元法是中国古法？ 2. 哪种线性方程组可用平方根法求解？为什么说平方根法计算稳定？ 3. 什么样的线性方程组可用追赶法求解并能保证计算稳定？ 4. 何为向量范数？给出三种常用的向量范数。 5. 何为矩阵范数? 6.高斯消去法与LU分解有什么关系？用它们解线性方程组Ax?b有何不同？A要满足什么条件？ 《高斯消元法》的相关理论 3.1.特征值和特征向量 设A是一个n?n阶实矩阵，若对于数?，存在非零向量x，使得Ax??x成立。则称?是A的特征值(Characteristic Value)，x为A的对应于?的特征向量(Characteristic Vector)。 3.2 向量和矩阵 用表示全部实矩阵的向量空间,表示全部复矩阵的向量空间. （称为m行n列矩阵）. （称为n维列向量） ,其中 为A的第列. 同理 ,其中 为A的第行. 矩阵基本运算: (1)矩阵加法 . (2)矩阵与标量的乘法 . (3)矩阵与矩阵的乘法 . （4）单位矩阵 ,其中 3.3 特殊矩阵 设,则有 A为: (1)对角矩阵 如果当时,; (2)三对角矩阵 如果当; (3)上三角矩阵 如果当; (4)对称矩阵 如果; (5)正定矩阵 如果设A是n阶实系数对称矩阵, 如果对任何非零向量 　　 都有,就称A正定. 四、《高斯消元法》国外研究进展 十几年来直接法在求解具有较大型稀疏矩阵方程组方面取得了较大进展。关于三对角线性方程组的直接求解已经有大量并行算法, 其中Wang 的分裂法是最早针对实际硬件环境, 基于分治策略提出的并行算法。它不仅通信结构简单, 容易推广到一般带状线性方程组的并行求解, 而且为相继出现的许多其它并行算法提