高斯 数值计算方法

数学与软件科学学院实验报告学期： 至__第__学期 年月曰课程名称:___数值计算方法与算法 专业:—级-班实验编号：0乞实验项目Gauss列主元素消元法求解线性方程组指导教师_张莉—姓名：林海学号：2009060619 实验成绩： 一、实验目的及要求实验目的：掌握使用列主元素法求解线性方程组的算法，比较Gauss列主元素消元法与Gauss消元的区别和算法实现时的特点。实验要求：使用C语言实现Gauss列主元素消元法。二、实验内容用Gauss列主元素消元法求解方程组：2x+x+x=41 2 3<x+3x+2x=6123x+2x+2x=5123要求：使用P108的算法，存储结构按照上课讲的要求：用A[4][5](或A[3][4])来存储增广阵，消元最后得到的上三角阵为数组的上三角部分。最后的输出结果输出增广阵的最后一列作为X的近似解，并且输出数组的上三角部分查看结果与P116页例5.2的U矩阵作比较。三、实验步骤(该部分不够填写.请填写附页)程序代码如下：#include<stdio.h>#include<math.h>voidmain(){/* step1 */intn,i,j,m,u,v,k,s;double**Arr,*Brr,*Xrr,temp;clrscr();printf("Pleaseinputtheordernis::\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("PleaseinputArr[][]::\n");for(i=1;i<n;i++){for(j=1;j<n;j++){printf("A[%d][%d]=\t",i,j);scanf("%lf",&Arr[i][j]);}printf("\n");}printf("\nPleaseinputBrr[]::\n");for(i=1;i<n;i++){printf("Brr[%d]=\t",i);scanf("%lf",&Brr[i]);}/* step2 */for(k=1;k<n;k++){s=fabs(Arr[k][k]);m=k;for(u=k+1;u<n;u++){if(fabs(Arr[u][k])>s)m=u;}for(v=k;v<n;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<n;i++){temp=Arr[i][k]/Arr[k][k];for(j=k+1;j<n;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;j<n;j++)temp=temp+Arr[i][j]*Xrr[j];}Xrr[i]=(Brr[i]-temp)/Arr[i][i];printf("\nTheX[%d]=%lf",i,Xrr[i]);}getch();}

8,通过三角矩阵求解Xrr[]□的值。8,通过三角矩阵求解Xrr[]□的值。9,运行结果。■ A!indow-Vsystem37\cind.e5cefor(i=n-l;i>=l;i—){temp=6;for(j=i<：j<n：j—){temp='teinp-,-Arr[i][j]*Xrr[j];1Hrr[i]=(Brr[il-temp)/firr[i][i];printf("\nTheX[阀]=Klf",i,Xrr[i]);P1总日weinputtheor'dernis;;3TOC\o"1-5"\h\zPIeas&inpuiflrrlIII::A[1][L1= 2nmr?i- ifl[1II3]= 1A[2J[11=AE21E3]-A[2J[11=AE21E3]-132122fHHIIl]=122A[3J[21=A[3J[3]=PleaseinputBrr[1:;HrrllI-kKrrl2|-6BrrIO]=5IheXI曰丨--乩刖:圮储TheXI2]=1.060060*TI®Hill-