列主高斯消去法

《列主高斯消去法》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、武夷学院实验报告课程名称：数值分析项目名称：线性方程组的直接解法姓名：专业：班级：（1）学号：20104161003同组成员一、实验准备注：1、实验准备部分包括实验环境准备和实验所需知识点准备。2、若是单人单组实验，同组成员填无。：线性方程组的直接解法用列主高斯消去法解下面线性方程组输出解向量x一、实验过程记录注：实验过程记录要包含实验步骤，页码不够可自行添加。：%列主元Gauss校区法解线性方程组A=[4111;13-11;1-120;1102;];%系数矩阵b=[0.650.0500.5]';%n维向量y=inv(A)*b%matlab的计算结果n=length(b)

2、;%方程个数nx=zeros(n,1);%未知向量%-------------消去-----------fork=1:n-1Auk=A(k:n,k);[m,u]=max(abs(Auk));u=u+k-1%u为最大元所在的列%------交换最大的行和当前行的值-------forj=k:ntemp=A(u,j);A(u,j)=A(k,j);A(k,j)=temp;endtemp=b(k);b(k)=b(u);b(u)=temp;%ifA(k,k)==0;%error('Error');%endfori=k+1:n%A(i,k)=A(i,k)/A(k,k);Aik=A(i,k)/A

3、(k,k)forj=k:nA(i,j)=A(i,j)-Aik*A(k,j);endAb(i)=b(i)-Aik*b(k)endend%-------------回代-----------x(n)=b(n)/A(n,n)fork=n-1:-1:1S=b(k);forj=k+1:nS=S-A(k,j)*x(j);endx(k)=S/A(k,k)endx%程序的计算结果error=abs(x-ones(n,1))%误差%误差小于直接进行高斯消去的计算误差实验操作成绩（百分制）__________实验指导教师签字：__________一、实验小结：注：1、实验小结应包含实验所需知识点和

4、实验方法的总结，实验心得体会等。2、分组实验需包含同组讨论内容。y=0.2000-0.2000-0.20000.2500u=1Aik=0.2500A=4.00001.00001.00001.000002.7500-1.25000.75001.0000-1.00002.000001.00001.000002.0000b=0.6500-0.112500.5000Aik=0.2500A=4.00001.00001.00001.000002.7500-1.25000.75000-1.25001.7500-0.25001.00001.000002.0000b=0.6500-0.1125-0.

5、16250.5000Aik=0.2500A=4.00001.00001.00001.000002.7500-1.25000.75000-1.25001.7500-0.250000.7500-0.25001.7500b=0.6500-0.1125-0.16250.3375u=2Aik=-0.4545A=4.00001.00001.00001.000002.7500-1.25000.7500001.18180.090900.7500-0.25001.7500b=0.6500-0.1125-0.21360.3375Aik=0.2727A=4.00001.00001.00001.00000

6、2.7500-1.25000.7500001.18180.0909000.09091.5455b=0.6500-0.1125-0.21360.3682u=3Aik=0.0769A=4.00001.00001.00001.000002.7500-1.25000.7500001.18180.09090001.5385b=0.6500-0.1125-0.21360.3846x=0000.2500x=00-0.20000.2500x=0-0.2000-0.20000.2500x=0.2000-0.2000-0.20000.2500x=0.2000-0.2000-0.20000.2500er

7、ror=0.80001.20001.20000.7500实验报告成绩（百分制）__________实验指导教师签字：__________