解線性方程組的直接方法實(shí)驗(yàn)報(bào)告

1、解線性方程組的直接方法實(shí)驗(yàn)報(bào)告1. 實(shí)驗(yàn)?zāi)康模?、通過(guò)該課題的實(shí)驗(yàn)，體會(huì)模塊化結(jié)構(gòu)程序設(shè)計(jì)方法的優(yōu)點(diǎn)；2、運(yùn)用所學(xué)的計(jì)算方法，解決各類線性方程組的直接算法；3、提高分析和解決問(wèn)題的能力，做到學(xué)以致用；4、通過(guò)三對(duì)角形線性方程組的解法，體會(huì)稀疏線性方程組解法的特點(diǎn)。2. 實(shí)驗(yàn)過(guò)程：實(shí)驗(yàn)代碼：#include"stdio.h"#include"math.h"#includeusingnamespacestd;/Gauss法voidlzy(double*a,double*b,intn)inti,j,k;doublel,x10,temp;for(k=0;k<

2、;n-1;k+)for(j=k,i=k;j<n;j+)if(j=k)temp=fabs(ajk);elseif(temp<fabs(ajk)temp=fabs(ajk);i=j;if(temp=0)cout<<"無(wú)解"return;elsefor(j=k;j<n;j+)temp=akj;akj=aij;aij=temp;temp=bk;bk=bi;bi=temp;for(i=k+1;i<n;i+)l=aik/akk;for(j=k;j<n;j+)aij=aij-l*akj;bi=bi-l*bk;if(an-1n-1=0)cout&l

3、t;<"無(wú)解"return;xn-1=bn-1/an-1n-1;for(i=n-2;i>=0;i-)temp=0;for(j=i+1;j<n;j+)temp=temp+aij*xj;xi=(bi-temp)/aii;for(i=0;i<n;i+)printf("x%d=%lf",i+1,xi);printf("");/平方根法voidpfg(double*a,double*b,intn)inti,k,m;doublex8,y8,temp;for(k=0;k<n;k+)temp=0;for(m=0;m<

4、;k;m+)temp=temp+pow(akm,2);if(akk<temp)return;akk=pow(akk-temp),1.0/2.0);for(i=k+1;i<n;i+)temp=0;for(m=0;m<k;m+)temp=temp+aim*akm;aik=(aik-temp)/akk;temp=0;for(m=0;m<k;m+)temp=temp+akm*ym;yk=(bk-temp)/akk;xn-1=yn-1/an-1n-1;for(k=n-2;k>=0;k-)temp=0;for(m=k+1;m<n;m+)temp=temp+amk*xm;

5、xk=(yk-temp)/akk;for(i=0;i<n;i+)printf("x%d=%lf",i+1,xi);printf("");/追趕法voidzgf(double*a,double*b,intn)inti;doublea010,c10,d10,a110,b110,x10,y10;for(i=0;i<n;i+)a0i=aii;if(i<n-1)ci=aii+1;if(i>0)di-1=aii-1;a10=a00;for(i=0;i<n-1;i+)b1i=ci/a1i;a1i+1=a0i+1-di+1*b1i;y0=b

6、0/a10;for(i=1;i<n;i+)yi=(bi-di*yi-1)/a1i;xn-1=yn-1;for(i=n-2;i>=0;i-)xi=yi-b1i*xi+1;for(i=0;i<n;i+)printf("x%d=%lf",i+1,xi);printf("");intmainintn,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*siz

7、eof(double);B1=(double*)malloc(n*sizeof(double);B2=(double*)malloc(n*sizeof(double);B3=(double*)malloc(n*sizeof(double);for(i=0;i<n;i+)Ai=(double*)malloc(n)*sizeof(double);Bi=(double*)malloc(n)*sizeof(double);Ci=(double*)malloc(n)*sizeof(double);cout<<"第一題(Gauss列主元消去法):"<<endl<<endl;cout<<&