[工学]实验二直接法解线性方程组.ppt

实验二 解线性方程组的直接方法 刘广利 Matlab实验讲解 列主元高斯消去法 A=[-3 2 6; 10 -7 0; 5 -1 5] b=[4 7 6] gauss(A,b) ans= 0.0000 -1.0000 1.000 再输入：[x, det, flag]= gauss(A,b)回车 X=0 -1 1 det=155 flag=OK 用平方根法解线性代数方程组的算法 （1）对矩阵A进行Cholesky分解，即A=LLT，由矩阵乘法： 对于 i = 1, 2,…, n 计算 j=i+1,i+2,…,n （2）求解下三角形方程组 LY=b （3）求解LTX = y 由于此法要将矩阵A作LLT三角分解，且在分解过程中含有开方运算，故称该称为LLT分解法或平方根法。 n = length(A); L=zeros(n); for i=1:n z=0; for k=1:i-1 z=z+L(i,k)^2; end L(i,i)=sqrt(A(i,i)-z); for j=i+1:n z=0; for k=1:i-1 z=z+L(j,k)*L(i,k); end L(j,i)=(A(j,i)-z)/L(i,i); end end j=i+1,i+2,…,n for i=1:n z=0; for k=1:i-1 z=z+L(i,k)*y(k); end y(i)=(b(i)-z)/L(i,i) End for i=n:-1:1 z=0; for k=i+1:n z=z+L(k,i)*x(k); end x(i)=(y(i)-z)/L(i,i); end A=[4 -1 1 ; -1 4.25 2.75; 1 2.75 3.5] [L, flag]=Col_Factor(A) L= 2 0 0 -0.5 2 0 0.5 1.5 1 Flag=OK 用平方根法： A=[9 -36 30; -36 192 30; 30 -180 180] b=ones(3,1) R=chol(A) 回车 再输入：x=R\(R’\b)回车 Jordan消去法 A=[1 2 3; 5 4 10; 3 -0.1 1] b=[1 0 2] Jordan(A,b) ans= 1.2000 2.0000 -1.4000 矩阵的三角形分解 直接三角分解法解AX = b的计算公式 对于r = 2, 3, …, n计算 （2）计算U的第r行元素 （3）计算L的第r 列元素 (r ? n) (1) (4) (5) L = eye(n);U = zeros(n); for i=1:n U(1,i)=A(1,i); L(i,1)=A(i,1)/U(1,1); end for r=2:n for i=r:n z=0; for k=1:r-1 z=z+L(r,k)*U(k,i); end U(r,i)=A(r,i)-z; end for i=r+1:n z=0; for k=1:r-1 z=z+L(i,k)*U(k,r); end L(i,r)=(A(i,r)-z)/U(r,r); end end y(1)=b(1); for i=2:n z=0; for k=1:i-1 z=z+L(i,k)*y(k); end y(i)=b(i)-z; end x(n)=y(n)/U(n,n); for i=n-1:-1:1 z=0; for k=i+1:n z=z+U(i,k)*x(k); end x(i)=(y(i)-z)/U(i,i); end function [L,U,y,x]=LU_Decom1(A,b) 实验5.1 （主元的选取与算法的稳定性）