解线性方程组的直接方法实验报告.doc

解线性方程组的直接方法实验报告 解线性方程组的直接方法实验报告 1.实验目的： 1、通过该课题的实验，体会模块化结构程序设计方法的优点； 2、运用所学的计算方法，解决各类线性方程组的直接算法； 3、提高分析和解决问题的能力，做到学以致用； 4、通过三对角形线性方程组的解法，体会稀疏线性方程组解法的特点。 2.实验过程： 实验代码： #include stdio.h #include math.h #include using namespace std; /Gauss法 void lzy(double *a,double *b,int n) int i,j,k; double l,x10,temp; for(k=0;kn-1;k+) for(j=k,i=k;jn;j+) if(j=k) temp=fabs(ajk); else if(tempfabs(ajk) temp=fabs(ajk); i=j; if(temp=0) cout无解 ; return; else for(j=k;jn;j+) temp=akj; akj=aij; aij=temp; temp=bk; bk=bi; bi=temp; for(i=k+1;in;i+) l=aik/akk; for(j=k;jn;j+) aij=aij-l*akj; bi=bi-l*bk; if(an-1n-1=0) cout=0;i-) temp=0; for(j=i+1;jn;j+) temp=temp+aij*xj; xi=(bi-temp)/aii; for(i=0;in;i+) printf(x%d=%lf ,i+1,xi); printf( ); /平方根法 void pfg(double *a,double *b,int n) int i,k,m; double x8,y8,temp; for(k=0;kn;k+) temp=0; for(m=0;mk;m+) temp=temp+pow(akm,2); if(akktemp) return; akk=pow(akk-temp),1.0/2.0); for(i=k+1;in;i+) temp=0; for(m=0;mk;m+) temp=temp+aim*akm; aik=(aik-temp)/akk; temp=0; for(m=0;m=0;k-) temp=0; for(m=k+1;mn;m+) temp=temp+amk*xm; xk=(yk-temp)/akk; for(i=0;in;i+) printf(x%d=%lf ,i+1,xi); printf( ); /追赶法 void zgf(double *a,double *b,int n) int i; double a010,c10,d10,a110,b110,x10,y10; for(i=0;in;i+) a0i=aii; if(i0) di-1=aii-1; a10=a00; for(i=0;in-1;i+) b1i=ci/a1i; a1i+1=a0i+1-di+1*b1i; y0=b0/a10; for(i=1;i=0;i-) xi=yi-b1i*xi+1; for(i=0;in;i+) printf(x%d=%lf ,i+1,xi); printf( ); int main int n,i,j; double *A,*B,*C,*B1,*B2,*B3; A=(double *)malloc(n*sizeof(double); B=(double *)malloc(n*sizeof(double); C=(double *)malloc(n*sizeof(double); B1=(double *)malloc(n*sizeof(double); B2=(double *)malloc(n*sizeof(double); B3=(double *)malloc(n*sizeof(double); for(i=0;in;i+) Ai=(double *)malloc(n)*sizeof(double); Bi=(double *)malloc(n)*sizeof(double); Ci=(double *)malloc(n)*sizeof(double)