高斯 數(shù)值計算方法

數(shù)學(xué)與軟件科學(xué)學(xué)院實驗報告學(xué)期： 至__第__學(xué)期 年月曰課程名稱:___數(shù)值計算方法與算法 專業(yè):—級-班實驗編號：0乞?qū)嶒烅椖縂auss列主元素消元法求解線性方程組指導(dǎo)教師_張莉—姓名：林海學(xué)號：2009060619 實驗成績： 一、實驗?zāi)康募耙髮嶒災(zāi)康模赫莆帐褂昧兄髟胤ㄇ蠼饩€性方程組的算法，比較Gauss列主元素消元法與Gauss消元的區(qū)別和算法實現(xiàn)時的特點。實驗要求：使用C語言實現(xiàn)Gauss列主元素消元法。二、實驗內(nèi)容用Gauss列主元素消元法求解方程組：2x+x+x=41 2 3<x+3x+2x=6123x+2x+2x=5123要求：使用P108的算法，存儲結(jié)構(gòu)按照上課講的要求：用A[4][5](或A[3][4])來存儲增廣陣，消元最后得到的上三角陣為數(shù)組的上三角部分。最后的輸出結(jié)果輸出增廣陣的最后一列作為X的近似解，并且輸出數(shù)組的上三角部分查看結(jié)果與P116頁例5.2的U矩陣作比較。三、實驗步驟(該部分不夠填寫.請?zhí)顚懜巾?程序代碼如下：#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,運行結(jié)果?！?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總?cè)誻einputtheor'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-