1.1三角形方程组与三角分解.pptVIP

第一章 线性方程组的直接解法;§1.1 三角形方程组和三角分解 1.1.1 三角形方程组的解法;这种解方程组(1.1.1)的方法称为前代法。;第i个式子的计算可先计算分子，算法为;第一步: b(1)=b(1)/L(1,1) b(2)=b(2)-L(2,1)*b(1) b(3)=b(3)-L(3,1)*b(1) …… b(n)=b(n)-L(n,1)*b(1);写成算法的形式：;上三角形方程组的解法 考虑下三角形方程组 Ux = y (1.1.2) y∈Rn已知, x∈Rn 未知, U∈Rn×n已知非奇异的上三角阵;类似算法1.1.1可得算法1.1.2(回代法). 运算量亦为n2.;1.1.2 Gauss变换;Gauss变换的性质 设 x∈Rn, xk≠0, 则存在Gauss变换 Lk 使得 Lk x = (x1 ··· xk 0 ··· 0)T . 事实上，;Gauss变换的逆( I－lkekT )－1 = I＋lkekT. 事实上，;1.1.3 三角分解的计算;使用Gauss变换实现矩阵的LU分解: 一个例子;解: (2);使用Gauss变换实现矩阵的LU分解: 一般步骤;(2) 对A(1), 计算Gauss变换L2使得;(3) 假定已求出L1, ·??· , Lk－1, 使得;0 ┇ 0;(4) 重复上述步骤 n-1 次, 得到A(n－1)为 n 阶上三角阵, 即 Ln－1···L2 L1A = A(n－1).;(6) 当Lk作用于A(k-1)后, A(k-1)的元素发生了什么变化？;(7) Lk和A(k)的元素怎样存储？;以4×4阶矩阵为例，;(8) 算法1.1.3 (计算三角分解: Gauss消去法) for k=1:n-1 A(k+1:n, k)=A(k+1:n, k)/A(k, k) A(k+1:n, k+1:n)=A(k+1:n, k+1:n) -A(k+1:n, k)A(k, k+1:n) end;(9) 两个定理 定义: 称Gauss消去法过程中的akk(k-1)为主元. 算法1.1.3要求akk(k-1) (k = 1, 2,···, n-1) 均不为0. 矩阵A满足什么条件才能保证所有主元均不为0？ 定理1.1.1 主元aii(i-1) (i = 1, 2,···, k) 均不为0;作业： 思考：举例说明一个非奇异矩阵不一定存在LU分解. 利用Doolittle分解法求解 L为下三角阵，U为单位上三角阵，则A=LU称为Crout分解. 利用Crout分解法求解上题. 预习算法1.2.1(全主元Gauss消去法).