数值分析第二章数值试验.docx

一、实验名称：解线性方程组的列主元素高斯消去法和LU分解法二、实验目的：通过数值实验，从中体会解线性方程组选主元的必要性和LU分解法的优点，以及方程组系数矩阵和右端向量的微小变化对解向量的影响。三、实验内容：解下列两个线性方程组（1）（2）列主元素高斯消去法：function [RA,RB,n,X]=liezhu(A,b)B=[A b]; n=length(b); RA=rank(A); RB=rank(B);zhica=RB-RA;if zhica>0,disp('请注意：因为RA~=RB，所以此方程组无解.')returnendif RA==RB if RA==ndisp('请注意：因为RA=RB=n，所以此方程组有唯一解.') X=zeros(n,1); C=zeros(1,n+1); for p= 1:n-1[Y,j]=max(abs(B(p:n,p))); C=B(p,:);B(p,:)= B(j+p-1,:); B(j+p-1,:)=C;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);endend 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); endelse disp('请注意：因为RA=RB<n，所以此方程组有无穷多解.')endend对于LU分解法：function hl=zhjLU(A)[n n] =size(A); RA=rank(A); if RA~=ndisp('请注意：因为A的n阶行列式hl等于零，所以A不能进行LU分解.A的秩RA如下:'), RA,hl=det(A);returnendif RA==n for p=1:nh(p)=det(A(1:p, 1:p));endhl=h(1:n);for i=1:nif h(1,i)==0disp('请注意：因为A的r阶主子式等于零，所以A不能进行LU分解.A的秩RA和各阶顺序主子式值hl依次如下：'), hl;RAreturnendend if h(1,i)~=0 disp('请注意：因为A的各阶主子式都不等于零，所以A能进行LU分解.A的秩RA和各阶顺序主子式值hl依次如下：')for j=1:nU(1,j)=A(1,j);endfor k=2:nfor i=2:n for j=2:n L(1,1)=1;L(i,i)=1; if i>jL(1,1)=1;L(2,1)=A(2,1)/U(1,1); L(i,1)=A(i,1)/U(1,1);L(i,k)=(A(i,k)- L(i,1:k-1)*U(1:k-1,k))/U(k,k);elseU(k,j)=A(k,j)-L(k,1:k-1)*U(1:k-1,j);endendendendhl;RA,U,Lendend四、实验结果：1.（1）对于第一题，运用列主元素高斯消去法输入：A=[3.016.03 1.99;1.27 4.16 -1.23;0.987 -4.81 9.34];b=[1;1;1];[RA,RB,n,x]=liezhu(A,b)得：请注意：因为RA=RB=n，所有此方程组有唯一解RA = 3RB = 3n = 3x = 1.0e+003 * 1.592599624841381 -0.631911376202549 -0.493617724759390（2）对于第一题运用：LU分解输入：A=[3.01 6.03 1.99;1.27 4.16 -1.23;0.987 -4.81 9.34];h1=zhjLU(A)得：请注意：因为A的各阶主子式都不等于零，所以A能进行LU分解.A的秩RA和各阶顺序主子式值h1依次如下：RA = 3U = 3.010000000000000 6.030000000000000 1.990000000000000 0 4.160000000000000 -2.069634551495017 0 0 5.310737406320019L = 1.000000000000000 0 0 0.421926910299