数值计算方法 孙延风版.pptVIP

第一章 解线性方程组的直接法 线性方程组 对于n个变元m个方程所组成的线性方程组 (1) 当右端常数项b1=b2=…=bm=0时，称为n元齐次线性方程组，否则称为n元非齐次线性方程组。 　 对于一般的线性方程组来说，所谓方程组(1)的一个解就是指由n个数k1，k2，…，kn组成的一个有序数组( k1，k2，…，kn )，当x1，x2，…，xn分别用k1，k2，…，kn代入后，使(1)中的每个等式都变成恒等式。方程组(1)的解的全体称为它的解集合。如果两个方程组有相同的解集合，我们就称它们是同解的。 写为矩阵形式 简记为Ax＝b 求解线性代数方程组的方法分类 直接法 假设计算过程中不产生舍入误差，经过有限次运算可求得方程组的精确解的方法。 低阶稠密矩阵方程组 迭代法 从解的某个近似值出发，通过构造一个无穷序列去逼近精确解的方法。 大型稀疏矩阵方程组 直接法 Gauss 消元法 矩阵的LU分解 预备知识 　分别表示n维实和复向量空间，用　 表示n×n阶实矩阵的向量空间. 转置矩阵 单位矩阵 　其中 非奇异矩阵 　设　　　　　　　　　。如果AB＝BA＝I，则称B是A的逆矩阵，记为A-1，且 。如果A-1存在，则称A为非奇异矩阵。如果 均为非奇异矩阵，则 矩阵的行列式 设 ，则A的行列式可按任一行（或列）展开，即 其中Aij 为aij的代数余子式， Aij=(-1)i+j Mij ,Mij为元素aij的余子式。 行列式性质 (a)det(AB)=det(A)det(B), (b) det(AT)=det(A), (c)det(cA)=cndet(A) (d) 是非奇异矩阵. 设 对角矩阵 如果当 时，aij=0 对称矩阵 如果 AT=A 对称正定矩阵 如果(a) AT=A,(b)对任意非零向量 ， 正交矩阵 如果A-1=AT 上三角矩阵及单位上三角矩阵 下三角矩阵及单位下三角矩阵 设A是一个m×n矩阵，对A施行一次初等行变换，相当于在A的左边乘以相应的m阶初等方阵；对A施行一次初等列变换，相当于在A的右边乘以相应的n阶初等方阵。 定理 设 ，则下述命题等价： (1)对任何 ,方程组Ax=b有唯一解； (2)齐次方程组Ax=0只有唯一解x=0； (3) (4)A-1存在 (5)A的秩rank(A)=n. 1.1 Gauss消元法 举例说明Gauss消元法求解过程 Gauss 消元法分为两个过程 消元过程 把原方程组化为阶梯形方程组 回代过程 求解方程组的解 Gauss消元法的基本思想 用消元法解方程组实际上就是反复地对方程组进行以下三种变换将其化为阶梯形式求解： (1) 用一个非零的数乘某一个方程； (2) 把一个方程的倍数加到另一个方程上； (3) 互换两个方程的位置。 我们称这样的三种变换为方程组的初等变换。 可以证明，初等变换总是把方程组变成同解的方程组。 例1：解线性方程组 (2) 把第二个方程的-3，2倍分别加到第三、四个方程上， 使得在第三、四个方程中消去未知量x2: 消去第四个方程中的x3: (3) 由此我们容易求出方程组的解为(-16，5，3，1)。 例 2 解方程组 解： 解： 用初等变换消去第二、三个方程中的x1: 再施行一次初等变换，得 上述方程组无解。 解： 用初等变换消去第二、三个方程中的x1: 再施行一次初等变换，得 改写成 最后得到 其中x2是自由未知量。