第2章解线性方程组直接法本章讨论n元线性方程组若矩阵A非奇异,即det(A)≠0,则方程组(2.1)有唯一解。消去后两个方程中x1得上述求解消元过程可用矩阵表示为:第一步.设,依次用其中如此继续消元下去,第n-1步结束后得到矩阵:次序Gauss消去法求解n元线性方程组乘除运算量是：§1.2列主元Gauss消去法用次序Gauss消去法,消元得然后进行第二步消元得增广矩阵B(3)=(A(3),b(3)).方程组含有四位有效数字准确解为x1*=17.46,x2*=-45.76,x3*=5.5462.用列主元Gauss消去法求解,消元过程为回代得:x3=5.545,x2=-45.77,x1=17.46§2矩阵三角分解法则有A(3)=也就是:所以有:式A=LU称为矩阵A三角分解.下面介绍矩阵三角分解Doolittle分解方法,对k=2,3,…,n,计算由利用三角分解方法解线性方程组先解由为了提升数值稳定性,可考虑列主元三角分解法,设已完成A=LUk-1步分解计算,矩阵分解成比如,用列主元三角分解解例3中方程组.则有设A为对称正定矩阵,则有唯一分解A=LU,且ukk>0.分解A=GGT称为对称正定矩阵Cholesky分解.解三角方程Gy=b,GTx=y可得平方根法是求对称正定系数线性方程组三角分解法,对称正定矩阵Cholesky分解计算量和存贮量均约为普通矩阵LU分解二分之一.且Cholesky分解含有数值稳定性.追赶法是求三对角线性方程组三角分解法.即方程其中解当满足条件|a1|>|c1|>0;|an|>|dn|>0;|ai||ci|+|di|,cidi0,i=2,3,…,n-1.时,追赶法是数值稳定,追赶法含有计算程序简单,存贮量少,计算量小优点.§3向量和矩阵范数记x=(x1,x2,…,xn)T,惯用向量范数有:惯用三种向量范数满足以下等价关系‖x‖‖x‖1n‖x‖,xRn§3.2矩阵范数矩阵-范数:‖A‖设‖‖是一个向量范数,则定义矩阵范数与矩阵特征值之间也有亲密联络.任何两种矩阵范数也含有等价性m‖A‖‖A‖M‖A‖,ARnn若把线性方程组变为设线性方程组Ax=b再设b是准确,A有误差A,此时记解为x+x,则(A+A)(x+x)=b经常使用条件数有Condp(A)=‖A‖p‖A-1‖pp=1，2，。(1)矩阵元素间数量级差很大,且无一定规律;称之为预条件方程组,显然与原方程组等价.可逆矩阵C称为预条件矩阵.矩阵C应满足条件练习题练习题练习题课堂练习课间