第六章(线性方程组直接解法).docVIP

第六章（上） 解线性方程组的直接法 §1、引言 1、　问题的提出 本章讨论个变量个线性方程的方程组求解，其表达式为　　　　　　通常用向量矩阵表示，则上述方程可写成　　　　　(1)为系数矩阵，为未知向量，为常数列（右端向量）） 其中　　　　，　　并记做，分别表示为阶实矩阵，为维实向量. 非奇异，即，方程组(式1)有唯一解，并可用Cramer法则将解用公式表示出来，但由于计算量太大，因此不能用于实际求解.分解法；（方法本身无误差） 　 迭代法：用某种极限过程去逐步逼近精确解，（无穷向量序列解向量）。（逐步精确法），Jacobi、G-S迭代法、SOR。（概念并不陌生，如在迭代法求非线性方程根时，从迭代初始值出发，反复迭代得到这样一个迭代序（数）列，它以方程的根为极限，类似） 3、阵特征值与谱半径　　定义1　设，若存在一个数 (实数或复数)和非零向量，使,则称为的特征值，为对应的特征向量，可使齐次方程 有非零解，故系数行列式，即　　　 则称 特征矩阵， 特征多项式，方程 特征方程，在复数域中有个根（的全体特征值为的谱为的谱半径. 　 的特征值和特征向量还有如下性质：　　与有相同的特征值及相同的特征向量.　　若非奇异，则的特征值为，特征向量为.. 例1　求的.　　解的特征方程为　　　）　　　故的特征值为.谱半径为.设，若，，，称为对称正定矩阵，若对有，则称为半正定矩阵.当为对称时　　() 对称正定的所有特征值；　　( 对称正定的所有顺序主子式　顺序主子阵；　顺序主子式； 性质：对称正定，则的对角元素中左上角前行列的子矩阵为的阶顺序主子阵，记为，称行列式为阶顺序主子式。（矩阵左上角的的元素组成子方阵的行列式值，或顺主子阵的行列式） 　　，对称：，， ， 注：对称正定，则非奇异，且也对称正定；定义 若满足,（*）（主对角元素的绝对值严格大于其所在行的其它元素的绝对值之和） 则称A为严格对角占优矩阵为严格对角占优矩阵。 则：，， 性质：若为严格对角占优， 变换类型 方程组同解性 行列式值 倍乘:, 同解 倍加: 同解 行交换: 同解 列交换: 不同解 　　备注：初中学过的消元法解方程二元一次方程，即利用一个方程将另一方程的某个未知量消去，变成一元一次方程后直接求解，再代回求另一个未知量；。§2 Gauss消去法 一、顺序消去法（简单消去法） 基本思想（通过消元将方程组变一个等价易于求解的同价上三角方程组，再回代求解） （的情况）： ，() （消元） （回代求解） 注：（1）步骤：消元过程：，回代求解： （2）称为主元素，（简主元，实现第次消元的关键元素，这里为顺序主元） 消元时要求不为0，回代求解：均不为0 ，方可用顺序消去法求解。 （3）(主元素的累乘积为系数矩阵行式值) Gauss消去法就是将方程组(式1)通过()步消元，将(式1)转化为上三角方程组　　　再回代求此方程组的解.　　下面记增广矩阵，即 　　第1步设，计算－×，可消去下方元素均为0) 公式：，　　一般地假定已完成了()步消元，即已将转化为以下形式：　　　　　第步，假定，计算　　　　－×　　　其中　　　　　当则可得到，即方程组　　直接回代解得，　　　　　并且有　在整个消去法消元过程中，共需乘除法运算次数为/*++)*/回代过程除法次数为，于是Gauss消去法的乘除法总次数为 　次乘除;注： 当：高斯法为430次，Cramer法则（计算量事前就可估计，所以为直接法） 可行性分析 定理　各阶顺序主子式 .可用顺序消去法求解。 证明“”（必要性）用归纳法，当，现假设()成立，即，对，故能进行步消元，即（顺序消元不会改变矩阵行列式值）由前提　从而有均成立，证毕.”（充分性）: 均不为0,消元同上 注： （1）当均不为0 a)　 b），（，反例见例3）， （2）不能保证一定能用顺序消去法求解。（仅满足时） 说明：在顺序消去（不引入行列交换）的求解时，（通过消元将化为上三角阵，然后将对角元素相乘即得行行列式值）；若系数矩阵不满足各阶顺序主子式不为0则不能用顺序高斯消去法求解。 推论1：若系数矩阵为对称正定 推论2：若系数矩阵为严格对角占优 例用Gauss消去法解方程组　　　　　　 解 　　，　　回代，得解 例3,为非奇异阵解。阵　解：显然为非奇异 ，，,由定理1知,不能用顺序消去法求解. 方法: 引入行交换,则（对于解方程组属于等价变换） ,，，, 由定理1知,能用顺序消去法求解 补充：（1）方程组的三种初等变换（所得的方程组与原方程组均为同解方程组）：交换方程组中任意两个方程的位置