资源描述：
《高斯列主元消元法解方程组》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。
1、温州大学计算机科学与工程学院 计算方法上机实验报告实验一、高斯列主元消元法解方程组题目：求解：，，，n=4;6;50;100;……求解精度为，并验证：及/随着n的逐渐增大，你会发现什么？一、实验目的： １、学习使用matlab编写数值计算程序。 ２、了解高斯列主元消元法解方程组原理和编程方法。 ３、根据高斯列主元消元法解方程组原理编写matlab程序，并运行得出结果。二、实验原理：如此至多经过n-1步,就得到与之同解的上三角形方程组的增广矩阵,再用回代过程即可得方程组的解.三、实验内容与步骤：物理与电子信息工程学院 计算机科学与技术专业 陈波 20091121 选修温州大学计算机
2、科学与工程学院 计算方法上机实验报告 １、实验内容：依照实验原理编写高斯列主元消元法解方程组的程序。 ２、实验步骤：首先，在电脑上安装matlab，然后，启动matlab，新建一个gauss.m文件。实验程序如下：程序代码程序代码说明functionx=gauss(A,b)[n,n]=size(A);x=zeros(n,1);Aug=[A,b];fork=1:n-1[piv,r]=max(abs(Aug(k:n,k)));r=r+k-1;ifr>ktemp=Aug(k,:);Aug(k,:)=Aug(r,:);Aug(r,:)=temp;endifAug(k,k)==0,error('对
3、角元出现0'),endforp=k+1:nAug(p,:)=Aug(p,:)-Aug(k,:)*Aug(p,k)/Aug(k,k);endend％ usethefunction。％求矩阵A的大小。％ 初始化x为0。％ 增广矩阵％for循环％找列主元所在子矩阵的行r％列主元所在大矩阵的行％矩阵为0，报错％把增广矩阵消元成为上三角物理与电子信息工程学院 计算机科学与技术专业 陈波 20091121 选修温州大学计算机科学与工程学院 计算方法上机实验报告A=Aug(:,1:n);b=Aug(:,n+1);x(n)=b(n)/A(n,n);fork=n-1:-1:1x(k)=b(k);for
4、p=n:-1:k+1x(k)=x(k)-A(k,p)*x(p);endx(k)=x(k)/A(k,k);end％解上三角方程组％求得x值四、实验数据及结果：１、当n为4时,在command窗口输入如下代码：formatlongn=4;fora=1:n;forc=1:n;A(a,c)=1/(a+c-1);endendb=[1:n]';x=gauss(A,b);fprintf('%12.8f',x)实验结果为：-64.00000000900.00000000-2520.000000001820.000000002、当n为6时,在command窗口输入如下代码：formatlongn=6;fora
5、=1:n;forc=1:n;A(a,c)=1/(a+c-1);endend物理与电子信息工程学院 计算机科学与技术专业 陈波 20091121 选修温州大学计算机科学与工程学院 计算方法上机实验报告b=[1:n]';x=gauss(A,b);fprintf('%12.8f',x)实验结果为：-216.000000047350.00000120-57120.00000827166320.00002178-201600.0000242585932.000009623、当n为50时,在command窗口输入如下代码：formatlongn=50;fora=1:n;forc=1:n;A(a,
6、c)=1/(a+c-1);endendb=[1:n]';x=gauss(A,b);fprintf('%12.8f',x)实验结果为：23964.74834020-4122991.37806186176597502.11690632-3260288682.1972771031898227709.80876200-181566275672.69705000619956768219.91504000-1230092760464.934800001193508577123.06670000-71714652666.11485300-587044042244.59863000-454093769462
7、.972900001744941129911.71780000-2715153157131.341300002706193814657.00730000758555481068.28345000-2729174782832.09330000-650678553386.91565000770544482197.08948000661071297466.389400003114329102228.31640000
