线性方程组的数值解法讲课文档.pptVIP

3.4 高斯消去法和选主元（续12） 一个线性方程组称为是病态的，如果其系数矩阵接近奇异且它的行列式接近0 矩阵条件数 cond(A)=||A||||A-1|| 现在二十七页，总共五十二页。 3.5 三角分解法 A=LU:下三角矩阵L的主对角线为1，上三角矩阵U的对角线元素非零 定义3.4 如果非奇异矩阵A可表示为下三角矩阵L和上三角矩阵U的乘积：A=LU，则A存在一个三角分解 A非奇异蕴含着对所有的k有ukk≠0，k=1,2,3,4. 现在二十八页，总共五十二页。 矩阵的LU分解 是否所有的非奇异矩阵A都能作LU分解呢？ 一个例子： N阶方阵A有唯一LU分解的充要条件是A的各阶顺序主子式均不为零 现在二十九页，总共五十二页。 3.5 三角分解法（续1） 利用前代/回代算法求解形如Lx=b或Ux=b的线性方程组是容易的 如果对一个给定的矩阵A，能够找到一个下三角矩阵L和一个上三角矩阵U，使A=LU 则求解线性方程组Ax=b的问题可以分解成两个简单的问题： Ly=b Ux=y 易见：Ax=(LU)x=L(Ux)=Ly=b 现在三十页，总共五十二页。 3.5 三角分解法（续2） 假设已有矩阵A: 对A作LU分解: 检验分解结果: 现在三十一页，总共五十二页。 3.5 三角分解法（续3） 构造一系列乘数矩阵M1, M2