高斯 数值计算方法

《高斯 数值计算方法》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、数学与软件科学学院实验报告学期：_____至___第_学期年月日课程名称:___数值计算方法与算法专业:级__班实验编号：05实验项目Gauss列主元素消元法求解线性方程组指导教师_张莉_姓名：林海学号：2009060619实验成绩：一、实验目的及要求实验目的：掌握使用列主元素法求解线性方程组的算法，比较Gauss列主元素消元法与Gauss消元的区别和算法实现时的特点。实验要求：使用C语言实现Gauss列主元素消元法。二、实验内容用Gauss列主元素消元法求解方程组：要求：（1）使用P108的算法，存储结构按照上课讲的要求：用A[4]

2、[5]（或A[3][4]）来存储增广阵，消元最后得到的上三角阵为数组的上三角部分。（2）最后的输出结果输出增广阵的最后一列作为X的近似解，并且输出数组的上三角部分查看结果与P116页例5.2的U矩阵作比较。三、实验步骤(该部分不够填写.请填写附页)程序代码如下：#include#includevoidmain(){/*----------step1-------------------*/intn,i,j,m,u,v,k,s;double**Arr,*Brr,*Xrr,temp;clrscr();pr

3、intf("Pleaseinputtheordernis::");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("PleaseinputArr[][]::");for(i=1;i

4、f("%lf",&Arr[i][j]);}printf("");}printf("PleaseinputBrr[]::");for(i=1;is)m=u;}for(v=k;v

5、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;i=1;i--){temp=0;for(j=i+1;j

6、emp=temp+Arr[i][j]*Xrr[j];}Xrr[i]=(Brr[i]-temp)/Arr[i][i];printf("TheX[%d]=%lf",i,Xrr[i]);}getch();}1，头文件。2，定义所要使用的变量。n为阶数，其余整形变量为控制循环；然后用指针定义数组；temp为中间变量。3，输入阶数。4，动态的分配空间大小。5，输入系数矩阵和常数矩阵。6，将系数矩阵Arr[][]中列向量中系数最大的化到对角矩阵中。7，化矩阵为斜三角形矩阵。8，通过三角矩阵求解Xrr[][]的值。9，运行结果。实验报告附页四

7、、实验结果分析与评价(该部分不够填写.请填写附页)注:实验成绩等级分为(90-100分)优,(80-89分)良,(70-79分)中,(60-69分)及格,(59分)不及格