數(shù)值計算實驗報告1

《數(shù)值計算基礎》實驗報告專業(yè)：班級：學號：姓名：指導老師：

實驗名稱全選主元消去法，Doolittle分解實驗日期4月實驗目的（10分）實驗內(nèi)容一(35分)實驗內(nèi)容二（35分）實驗小結(jié)（20分）總分【實驗目的】熟悉線性代數(shù)方程組高斯消去法，高斯列主元消去法，高斯全主元消去【實驗內(nèi)容一】一、問題用Gauss消去法求解方程組二、算法描述1.輸入方程組維數(shù)n，矩陣A，右端項b和控制精度eps2.對于k=1：n=1， （1）|A（u，k）|=max|A（i，k）|；%選主元 （2）如果|A（u，k）|<eps則停止；%控制小主元 （3）如果u=k則轉(zhuǎn)（4），否則 A（k，k：n+1）A（u，k：n+1）；%換行 （4）A（k+1：n,k）:=A(k+1:n,k)/A(k,k), A（k+1：n,k+1:n）:=A(k+1:n,k+1:n)-A(k+1:n,k)*A(k,k+1:n),%消元按行進行 b-A(k+1:n,k)*b(k)3.如果A（n，n）=0則停止4.b（n）:=b(n)/A(n,n), 對于i=n-1：-1:1， b（i）:=[b(i)-A(I,i+1:n)*b(i+1:n)]/A(i,i)%回代求解5.輸出解b（1：n)三、關鍵源代碼voidallselect(doubleA[][100],double*b,int*order,intN){//全選主元 //變量申明 for(intk=0;k<N;k++){ p1=p2=k; for(inti=k;i<N;i++)//找最大 for(intj=k;j<N;j++){ if(fabs(A[i][j])>m){ m=fabs(A[i][j]); p1=i;p2=j; } } if(p1!=k){//換行 for(i=0;i<N;i++){ var=A[p1][i];A[p1][i]=A[k][i];A[k][i]=var; } bm=b[p1];b[p1]=b[k];b[k]=bm; } if(p2!=k){//換列 for(i=0;i<N;i++){ var=A[i][p2];A[i][p2]=A[i][k];A[i][k]=var; } v=order[k];order[k]=order[p2];order[p2]=v;} }}voidtobeonlyone(doubleA[][100],double*b,intN){//歸一化for(inti=0;i<N;i++){ b[i]=b[i]/A[i][i]; for(intj=N-1;j>=i;j--) A[i][j]=A[i][j]/A[i][i]; for(intk=i+1;k<N;k++) { b[k]=b[k]-A[k][i]*b[i]; for(j=N-1;j>=i;j--) A[k][j]=A[k][j]-A[k][i]*A[i][j];} }}voidcomeback(doubleA[][100],double*b,double*x,intN){//回代 doubletotal; x[N-1]=b[N-1]; for(inti=N-2;i>=0;i--) { total=0; for(intj=i+1;j<=N-1;j++) total+=x[j]*A[i][j]; x[i]=b[i]-total;}}

四、運行結(jié)果【實驗內(nèi)容二】一、問題用Doolittle法解方程組：二、算法描述、 第一步： 第二步：三、關鍵源代碼voiddoolite(double*p,intn){//Doolittle分解inti,j,k; doublesum1,sum2,*L,*U;L=creat(n);U=creat(n);for(i=0;i<n;i++){ *(U+0*n+i)=*(p+0*n+i); *(L+i*n+0)=*(p+i*n+0)/(*(U+0*n+0)); } for(i=1;i<n;i++){ for(j=i;j<n;j++){ sum1=0; for(k=0;k<=i-1;k++){ sum1+=*(L+i*n+k)*(*(U+k*n+j)); } *(U+i*n+j)=*(p+i*n+j)-sum1; } for(j=i+1;j<n;j++){ sum2=0; for(k=0;k<=i-1;k++){sum2+=*(L+j*n+k)*(*(U+k*n+i)); } *(L+j*n+i)=(*(p+j*n+i)-sum2)/(*(U+i*n+i)); } }for(i=0;i<n;i++) for(j=0;j<n;j++){ if(i>j)*(p+i*n+j)=*(L+i*n+j); else*(p+i*n+j)=*(U+i*n+j); } for(i=0;i<n;i++){for(j=0;j<n;j++)printf("%.4lf\t",*(p+i