线性代数方程组的解法.ppt

湘潭大学数学与计算科学学院 * 湘潭大学数学与计算科学学院 * 湘潭大学数学与计算科学学院 * 设 总可由消元过程得到系数矩阵为上三角阵的线性代数 方程组，其第 步的结果为 湘潭大学数学与计算科学学院 * 其中 这里取 2. 回代过程 若通过消元过程原方程组已化为等价的三角形 方程组 湘潭大学数学与计算科学学院 * 且 , 则逐步回代可得原方程组的解 湘潭大学数学与计算科学学院 * Gauss逐步消去法有如下的缺点： 任一主元 ，就无法做下去 任一 绝对值很小时,也不行（舍入误差的影响大） 2.2 Gauss主元素消去法 下面我们讨论列主元消去法. 假设Gauss消去法的消元过程进行到 第 步， 设 湘潭大学数学与计算科学学院 * 并令 为达到最大值 的最小行标 ， 则交换 和 中的第 行和第 行， 再进行消元过程的第 步. 这时每个乘子 都满足 可以防止有效数字大量丢失而产生误差. 湘潭大学数学与计算科学学院 * 例 用列主元消去法解如下方程组 解 对增广矩阵按列选主元再进行高斯消元 湘潭大学数学与计算科学学院 * 湘潭大学数学与计算科学学院 * 回代求解得 湘潭大学数学与计算科学学院 * %magauss2.m function x=magauss2(A,b,flag) %用途：列主元Gauss消去法解线性方程组Ax=b %格式：x=magauss(A,b,flag), A为系数矩阵, b为右端项, 若flag=0, % 则不显示中间过程，否则显示中间过程, 默认为0, x为解向量 if nargin3,flag=0;end n=length(b); for k=1:(n-1) %选主元 [ap,p]=max(abs(A(k:n,k))); p=p+k-1; if pk t=A(k,:); A(k,:)=A(p,:); A(p,:)=t; t=b(k); b(k)=b(p); b(p)=t; end 湘潭大学数学与计算科学学院 * %消元 m=A(k+1:n,k)/A(k,k); A(k+1:n,k+1:n)=A(k+1:n,k+1:n)-m*A(k,k+1:n); b(k+1:n)=b(k+1:n)-m*b(k); A(k+1:n,k)=zeros(n-k,1); if flag~=0, Ab=[A,b], end end %回代 x=zeros(n,1); x(n)=b(n)/A(n,n); for k=n-1:-1:1 x(k)=(b(k)-A(k,k+1:n)*x(k+1:n))/A(k,k); end 湘潭大学数学与计算科学学院 * 全主元消去法 定义 此时交换 和 的行及A的列，使主元位置的元素 的绝对值具有给出的最大值 ， 然后进行第 步消元过程 注意：因为有列的交换，因此未知量的次序有改变，待消元过程结束时必须还原.多使用列主元消去法. 湘潭大学数学与计算科学学院 * Gauss消去法的实质是将矩阵 分解为 其中 --单位下三角矩阵， --上三角矩阵. 事实上，线性方程组 经过 步消元过程后，有等价方程组 其中： ，而 和 的形式为： 2.3 矩阵的三角分解与Gauss消去法的变形 湘潭大学数学与计算科学学院 * （1） 可以直接验证 ， 湘潭大学数学与计算科学学院 * 湘潭大学数学与计算科学学院 * 其中 湘潭大学数学与计算科学学院 * 乘积 是下三角阵，且对角元全部等于1 则 也是对角元等于1的下三角阵 用矩阵 依次左乘原给方程组 两边，得等价方程组 则 其中 湘潭大学数学与计算科学学院 * 由于 为上三角阵，记 ,于是得到 (2) 湘潭大学数学与计算科学学院 * Gauss逐步消去法等价于下述过程： 2. 求解三角形方程组 （回代过程）. (注意上面的全部讨论中要求 ) 湘潭大学数学与计算科学学院 * 比较等式两边对应元素算出 Doolittle