数值分析 第5章 解线性方程组的直接方法.ppt

《数值分析 第5章 解线性方程组的直接方法.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、5.1引言与预备知识5.2高斯消去法5.3高斯主元素消去法5.4矩阵三角分解法5.5向量和矩阵的范数5.6误差分析Ch5解线性方程组的直接方法在自然科学和工程技术中有很多问题的解决常常归结为解线性代数方程组．如三次样条函数问题，用最小二乘法求实验数据的曲线拟合问题，解非线性方程组问题，用差分法或者有限元方法解常微分方程、偏微分方程的边值问题等都导致求解线性代数方程组，而这些方程组的系数矩阵大致分为两种，一种是低阶稠密矩阵，另一种是大型稀疏矩阵。关于线性方程组的数值解法一般有两类：1.直接法2.迭代法§5.1引言与预备知识5.1.1引言本章讨论n元线性方程组（5.1）的直接解法。方程

2、组(5.1)的矩阵形式为Ax=b其中若矩阵A非奇异,即det(A)≠0,则方程组(2.1)有唯一解。所谓直接解法是指，若不考虑计算过程中的舍入误差，经过有限次算术运算就能求出线性方程组的精确解的方法。但由于实际计算中舍入误差的存在，用直接解法一般也只能求出方程组的近似解。Cramer法则是一种不实用的直接法，本章将介绍几种实用的直接法。5.1.2预备知识M行n列矩阵.n维列向量.矩阵的基本运算:(1)矩阵的加法(7)矩阵的行列式行列式性质:(a)det(AB)=det(A)det(B)(6)非奇异矩阵(5)单位矩阵(4)转置矩阵(3)矩阵与矩阵的乘法(2)矩阵与标量的乘法5.1.3

3、矩阵特征值与谱半径定义1设若存在一个数λ(实数或复数)和非零向量使(1.1)则称λ为A的特征值,x为A对应λ的特征向量,A的全体特征值称为A的谱，记作称为A的谱半径.(1.2)由式(1.1)知,λ可使齐次方程组有非零解,故系数行列式记称为矩阵Ａ的特征多项式，方程(1.3)称为Ａ的特征方程．(1.3)在复数域中有n个根故由行列式(1.3)展开可知：的n个特征值故是它的特征方程(1.3)的几个根,并有(1.4)(1.5)A的迹.A的特征值λ和特征向量x还有以下性质:(1)AT与A有相同的特征值λ及相同的特征向量x.(2)若A非奇异，则A-1的特征值为λ-1,特征向量为x.(3)相似矩阵

4、B=S-1AS有相同的特征多项式.例1求的特征值及谱半径.解:A的特征方程为故A的特征值为A的谱半径为5.1.4特殊矩阵定理1.定理2.定理3.定理4(Jordan标准型)设A为n阶矩阵,则存在一个非奇异矩阵P使得其中:(1)当A的若当标准型中所有若当块Ji均为一阶时,此标准型变成对角矩阵.返回主页求解高斯消去法(逐次消去法)及消去法和矩阵三角分解之间的关系：§5.2高斯消去法例2用消去法解方程组解第1步.将方程(2.2)乘上-2加到方程(2.4)上,消去未知数x1,得到(2.2)(2.3)(2.4)第2步.将方程(2.3)加到方程(2.5)上去,消去方程(2.5)中的x2,得到

5、与原方程组等价的三角形方程组解为:首先举一个简单的例子来说明消去法的基本思想.上述过程相当于思路首先将A化为上三角阵/*upper-triangularmatrix*/，再回代求解/*backwardsubstitution*/。=消元记Step1：设，计算因子将增广矩阵/*augmentedmatrix*/第i行mi1第1行，得到其中Stepk：设，计算因子且计算回代注意1：只要A非奇异，即A1存在，则可通过逐次消元及行交换，将方程组化为三角形方程组，求出唯一解。共进行？步n1注意2:设Ax=b,其中A为非奇异矩阵,如果由于A为非奇异矩阵,所以A的第一列一定有元素不等于零

6、.例如定理5设Ax=b,其中(1)如果则可通过高斯消去法将Ax=b约化为等价的三角形线性方程组(2.10),且计算公式为:①消元计算(k=1,2,…,n-1)②回代计算(2)如果A为非奇异矩阵,则可通过高斯消去法(及交换两行的初等变换)将方程组Ax=b约化为方程组(2.10).定理6约化的主元素aii(i)≠0(i=1,2,…,k)的充要条件是矩阵A的所有顺序主子式/*determinantofleadingprincipalsubmatrices*/Di≠0(i=1,2,…,k).即(2.12)推论如果A的顺序主子式Dk≠0(k=1,2,…,n-1),则§5.2.2三角分解法/*

7、MatrixFactorization*/高斯消元法的矩阵形式/*MatrixFormofG.E.*/：Step1:记L1=，则Stepn1:其中Lk=记为L单位下三角阵/*unitarylower-triangularmatrix*/记U=定理7若A的所有顺序主子式/*determinantofleadingprincipalsubmatrices*/均不为0，则A的LU分解唯一（其中L为单位下三角阵）。证明：由§1中定理可知，LU分解存在。下面证明唯一性。若不