数值分析实验报告

1 / 15数值分析实验报告数值分析实验报告12345课程实验报告课程实验报告数值分析实验报告课题一：解线性方程组的直接方法2 / 151.实验目的：1、通过该课题的实验，体会模块化结构程序设计方法的优点；2、运用所学的计算方法，解决各类线性方程组的直接算法；3、提高分析和解决问题的能力，做到学以致用；4、通过三对角形线性方程组的解法，体会稀疏线性方程组解法的特点。2.实验过程：实验代码：#include “#include “#include3 / 15using namespace std;/Gauss 法void lzy(double *a,double *b,int n)int i,j,k;double l,x10,temp;for(k=0;k for(j=k,i=k;j if(j=k)temp=fabs(ajk);else if(temp 4 / 15temp=fabs(ajk);i=j;if(temp=0)cout return;elsefor(j=k;j temp=akj;5 / 15akj=aij;aij=temp;temp=bk;bk=bi;bi=temp;for(i=k+1;i l=aik/akk;for(j=k;j aij=aij-l*akj;bi=bi-l*bk;6 / 15if(an-1n-1=0)cout return;xn-1=bn-1/an-1n-1;for(i=n-2;i=0;i-)temp=0;for(j=i+1;j temp=temp+aij*xj;xi=(bi-temp)/aii;for(i=0;i 7 / 15printf(“x%d=%lft”,i+1,xi);printf(“n”);/平方根法void pfg(double *a,double *b,int n)int i,k,m;double x8,y8,temp;for(k=0;k temp=0;for(m=0;m temp=temp+pow(akm,2);8 / 15if(akk return;akk=pow(akk-temp),/);for(i=k+1;i temp=0;for(m=0;m temp=temp+aim*akm; aik=(aik-temp)/akk;temp=0;for(m=0;m temp=temp+akm*ym;yk=(bk-temp)/akk;xn-1=yn-1/an-1n-1;for(k=n-2;k=0;k-)9 / 15temp=0;for(m=k+1;m temp=temp+amk*xm;xk=(yk-temp)/akk;(来自: 海达范文网:数值分析实验报告) for(i=0;i printf(“x%d=%lft”,i+1,xi);printf(“n”);/追赶法void zgf(double *a,double *b,int n)10 / 15int i;double a010,c10,d10,a110,b110,x10,y10; for(i=0;i a0i=aii;if(i ci=aii+1;if(i0)di-1=aii-1;a10=a00;for(i=0;i b1i=ci/a1i;11 / 15a1i+1=a0i+1-di+1*b1i;y0=b0/a10;for(i=1;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 printf(“x%d=%lft”,i+1,xi);printf(“n”);12 / 15int 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;i Ai=(double *)malloc(n)*sizeof(double); 13 / 15Bi=(double *)malloc(n)*sizeof(double); Ci=(double *)malloc(n)*sizeof(double); cout cinn;cout for(i=0;i for(j=0;j 武汉理工大学学 生 实 验 报 告 书实验课程名称： 数 值 分 析 开 课 名 学 生 姓 名：20011 2016 学年 第 二 学期第一次试验二分法计算流程图：简单迭代法算法流程图：牛顿迭代法流程图：14 / 15弦截法算法程