高斯 数值计算方法.doc

数学与软件科学学院 实验报告 学期：_至_ 第_ 学期 年 月 日 课程名称:_数值计算方法与算法 专业: 级_班实验编号：05 实验项目Gauss列主元素消元法求解线性方程组 指导教师_张莉_姓名： 林海 学号： 2009060619 实验成绩： 一、实验目的及要求实验目的：掌握使用列主元素法求解线性方程组的算法，比较Gauss列主元素消元法与Gauss消元的区别和算法实现时的特点。实验要求：使用C语言实现Gauss列主元素消元法。二、实验内容用Gauss列主元素消元法求解方程组：要求：（1）使用P108的算法，存储结构按照上课讲的要求：用A45（或A34）来存储增广阵，消元最后得到的上三角阵为数组的上三角部分。（2）最后的输出结果输出增广阵的最后一列作为X的近似解，并且输出数组的上三角部分查看结果与P116页例5.2的U矩阵作比较。三、实验步骤(该部分不够填写.请填写附页)程序代码如下：#include#includevoid main() /*-step 1-*/ int n,i,j,m,u,v,k,s; double *Arr,*Brr,*Xrr,temp; clrscr(); printf(Please input the order n is:n); scanf(%d,&i); n=i+1; Arr=(double *)malloc(sizeof(double)*n*n); Brr=(double *)malloc(sizeof(double)*n); Xrr=(double *)malloc(sizeof(double)*n); printf(Please input Arr:n); for(i=1;in;i+) for(j=1;jn;j+) printf(A%d%d=t,i,j);scanf(%lf,&Arrij); printf(n); printf(nPlease input Brr:n); for(i=1;in;i+) printf(Brr%d=t,i); scanf(%lf,&Brri); /*-step 2-*/ for(k=1;kn;k+) s=fabs(Arrkk); m=k; for(u=k+1;us) m=u; for(v=k;vn;v+) temp=Arrkv;Arrkv=Arrmv;Arrmv=temp;temp=Brrk;Brrk=Brrm;Brrm=temp; for(i=k+1;in;i+) temp=Arrik/Arrkk;for(j=k+1;j=1;i-) temp=0; for(j=i+1;jn;j+) temp=temp+Arrij*Xrrj; Xrri=(Brri-temp)/Arrii; printf(nThe X%d = %lf,i,Xrri); getch();1，头文件。2，定义所要使用的变量。n为阶数，其余整形变量为控制循环；然后用指针定义数组；temp为中间变量。3，输入阶数。4，动态的分配空间大小。5，输入系数矩阵和常数矩阵。6，将系数矩阵Arr中列向量中系数最大的化到对角矩阵中。7，化矩阵为斜三角形矩阵。8，通过三角矩阵求解Xrr的值。9，运行结果。实 验 报 告 附 页四、实验结果分析与评价 (该部分不