[理学]线性方程组的直接解法.pptVIP

第二章 线性方程组的直接解法 第一节 三角形方程组和三角分解 一、三角形方程组的解法 三角形方程组简单易于求解，而且它又是用分解方法解一般 线性方程组的基础，所以首先考虑这种特殊类型的线性方程 组的解法 求解方程组的Gauss算法 for(k=1:n-1) |A(p,q)|=max{|A(i,j)|:i=k:n,j=k:n}%确定p,q使得 A(k,1:n)?A(p,1:n) %交换p,k行 A(1:n,k)?A(1:n,q)%交换k,q列 u(k)=p %记录置换矩阵Pk v(k)=q %记录置换矩阵Qk If(A(k,k≠0) 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,k)-A(k+1:n,k)A(k,k+1,n)); else stop(矩阵奇异) end 算法1.3.2 改进平方根法 for(j=1:n) for(i=1:j-1) v(i)=A(j,i)A(i,i) end A(j,j)=A(j,j)-A(j,1:j-1)v(1:j-1) A(j+1:n,j)=(A(j+1:n,j)-A(j+1:n,1:j-1)v(1:j-1))/A(j,j) end 这种计算三角分解的方法称作Gauss消去法. 根据Ｌ和Ｕ的形状以及它们的形成顺