高斯数值计算方法.docVIP

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