第五次报告张董强.docVIP

数值实验报告 实验五 解线性方程组的直接法 学院：信 息 学 院 姓名： 张董强 班级：计算机12-1 学号：2011010225 解线性方程组的直接法 解线性方程组的直接法是指经过有限步运算后能求得方程组精确解的方法。 但由于实际计算中舍入误差是客观存在的，因而使用这类方法也只能得到近似 解。目前较实用的直接法是古老的高斯消去法的变形，即主元素消去法及矩阵的 三角分解法。引进选主元的技巧是为了控制计算过程中舍入误差的增长，减少舍 入误差的影响。一般说来，列主元消去法及列主元三角分解法是数值稳定的算法， 它具有精确度较高、计算量不大和算法组织容易等优点，是目前计算机上解中、 小型稠密矩阵方程组可靠而有效的常用方法。 一、 实验目的 1 、了解求线性方程组的直接法的有关理论和方法； 2 、会编制列主元消去法、LU 分解法的程序； 3、 通过实际计算，进一步了解各种方法的优缺点，选择合适的数值方法。 二、 实验任务 1. 用列主元高斯消去法求解方程组。 0.101x1+2.304x2+3.555x3=1.183; -1.347x1+3.172x2+4.623x3=2.137; -2.835x1+1.072x2+5.643x3=3.035; 解：function [RA,RB,n,X]=gaus(A,b) B=[A b]; n=length(b); RA=rank(A); RB=rank(B); zhica=RB-RA; if zhica0, disp(请注意：因为RA~=RB，所以此方程组无解.) return end if RA==RB if RA==n disp(请注意：因为RA=RB=n，所以此方程组有唯一解.) X=zeros(n,1); C=zeros(1,n+1); for p= 1:n-1 for k=p+1:n m= B(k,p)/ B(p,p); B(k,p:n+1)= B(k,p:n+1)-m* B(p,p:n+1); end end b=B(1:n,n+1);A=B(1:n,1:n); X(n)=b(n)/A(n,n); for q=n-1:-1:1 X(q)=(b(q)-sum(A(q,q+1:n)*X(q+1:n)))/A(q,q); end else disp(请注意： 因为RA=RBn， 所以此方程组有无穷多解.) end end 2. 用矩阵直接三角分解法求解方程组Ax =b?，其中 解： A = [0.101 2.304 3.555 -1.347 3.712 4.623 -2.835 1.072 5.643]; b = [1.183; 2.137; 3.035]; % A=[4,-2,0,4;-2,2,-3,1;0,-3,13,-7;4,1,-7,23]; % b=[1 2 3 4]; n=length(b); x=zeros(n,1); A(2:n,1)=A(2:n,1)./A(1,1); for i=2:n-1 A(i,i)=A(i,i)-sum(A(i,1:i-1).*A(1:i-1,i)); for j=i+1:n A(i,j)=A(i,j)-sum(A(i,1:i-1).*A(1:i-1,j)); A(j,i)=(A(j,i)-sum(A(j,1:i-1).*A(1:i-1,i)))/A(i,i); end end A(n,n)=A(n,n)-sum(A(n,1:n-1).*A(1:n-1,n)); A U=A;L=A; for i=1:n L(i,i)=1; end for i=1:n-1 for j=i+1:n L(i,j)=0; end end for i=2:n for j=1:i-1 U(i,j)=0; end end y=zeros(n,1); y(1)=b(1); for i=2:n y(i)=b(i)-sum(L(i,1:i-1).*y(1:i-1)); end y x(n)=y(n)/U(n,n); for i=n-1:-1:1 x(i)=(y(i)-sum(U(i,i+1:n).*x(i+1:n)))/U(i,i); end 3 、设对称正定阵系数阵线方