《计算机数值计算方法》由会员分享，可在线阅读，更多相关《计算机数值计算方法（27页珍藏版）》请在金锄头文库上搜索。
1、本科实验报告课程名称：计算机数值方法实验项目：二分法牛顿法迭代法实验地点：学院楼606专业班级：学号：学生姓名：指导教师：王丽娟年05月8日太原理工大学学生实验报告学院名称计算机科学与技术专业班级学号学生姓名实验日期5月8日成绩课程名称数值计算方法实验题目二分法迭代法求解方程、实验目的和要求(必填)熟悉使用二分法、迭代法、牛顿法、割线法等方法对给定的方程进行根的求解。选择上述方法中的两种方法求方程：f(x)=x3+4x2-10=0在1,2内的一个实根，且要求满足精度lx*-xlvO.5x1O-5n二、实验内容和原理(必填)f(x)在区间(x，y)上连续先找到a、b属于区间(x，y)，使f(a)
2、，f(b)异号，说明在区间(a,b)内一定有零点，然后求f(a+b)/2,现在假设f(a)O,ab如果f(a+b)/2=0，该点就是零点,如果f(a+b)/2a，从开始继续使用中点函数值判断。如果f(a+b)/20，则在区间(a,(a+b)/2)内有零点，(a+b)/2二b，从开始继续使用中点函数值判断。这样就可以不断接近零点。通过每次把f(x)的零点所在小区间收缩一半的方法，使区间的两个端点逐步迫近函数的零点，以求得零点的近似值三、主要仪器设备(必填)PCQEV-C四、操作方法与实验步骤(可选)1.迭代法#includevstdio.hintmain()inti=0;doublea;doub
3、leb100;printf(”请输入任意一个在1与2之间的数n);scanf(%lf,&b0);dobi+1=bi-(bi*bi*bi+4*bi*bi-10)/(3*bi*bi+8*bi);i+;while(bi-bi+l0.000005);printf(%.81f,bi);getch();LC:Usersbos5Deskto方法实迭代法心亡请输入任意一个在1与2之间的数1.2L36S230B12牛顿法#includestdio.h#includemath.hmain()doublex=1.5,a=0;while(pow(pow(x-a),2),0.5)/20.5*0.00001)a=x;x=
4、0.5*pow(10-pow(x,3),0.5);printf(%fn,x);getchar();28695440254134545837517&360094367S4736388736591736487836541&365138365277136520613652421365224k.365233*G:de叭DE/-CPPEx吕nnpl亡或H已IIclexe六、实验结果与分析（必填）结果与书上的答案一致，编译过程中对do-while循环中条件的存在一定的不合适。实验地点综合楼八层606室指导教师王丽娟本科实验报告课程名称：计算机数值方法实验项目：高斯消元法，LU分解法，追赶法实验地点：学院楼
5、606专业班级：_学号：_学生姓名：指导教师：王丽娟年5月8日太原理工大学学生实验报告学院名称计算机科学与技术专业班级学号学生姓名实验日期5月8日成绩课程名称计算机数值方法实验题目高斯消兀法，LU分解法，追赶法、实验目的和要求（必填）2）线性方程组的直接解法合理利用Gauss消元法、LU分解法、追赶法求解下列方程组:123x114012x28241x1-3J130.3x10-1559.1431_x1x59.17_5.291-6.130-1246.7811.2952乙x3x1-4112112_4215_x1-2_872101x2-74836x3-71261120_x1-4_-3_(n=5,10,
6、100,)21实验内容和原理（必填）高斯分解法:将原方程组化为三角形方阵的方程组：l=a/aikikkka=a-l*ak=l,2,，n-1ijijikkji=k+1,k+2,，nj=k+1,k+2,，n+1由回代过程求得原方程组的解：x=a/annn+1nn(k=n-1,n_2,2,1)x=(a工ax)/akkn+1kjjkkLU分解法：将系数矩阵A转化为A=L*U,L为单位下三角矩阵，U为普通上三角矩阵，然后通过解方程组l*y=b,u*x=y,来求解x.追赶法：用来求对角方程组；将系数矩阵A转化为A=L*U,L为普通下n-1对角矩阵，U为单位上n-1对角矩阵，然后通过解方程组l*y=b,u*
7、x=y,来求解x.三、主要仪器设备(必填)PC,DEV-C+四、操作方法与实验步骤(可选)1高斯消元法#includevstdio.hintmain()inti,j;floatz,b,c,d,xl,x2,x3;floata34;printf(请输入矩阵n);for(i=0;iv3;i+)for(j=0;jv4;j+)scanf(%f,&aij);z=a10/a00;if(z!=0);for(j=0;j4;j+)a1j=a1j-z*a0j;b=a20/a00;if(b!=0);for(j=0;j4;j+)a2j=a2j-b*a0j;c=a21/a11;if(c!=0);for(j=0;jintm
8、ain()inti,j,n;floatz,b,c,k,xl,x2,x3,x4;floata45;printf(请输入矩阵n);for(i=0;i4;i+)for(j=0;jv5;j+)scanf(%f,&aij);for(i=l;i4;i+)if(a00vai0)for(j=0;jv5;j+)k=a0j;a0j=aij;aij=k;for(n=l;n4;n+)z=an0/a00;if(z!=0);for(j=0;jv5;j+)anj=anj-z*a0j;for(n=2;nv4;n+)c=anl/all;if(c!=O)for(j=0;jv5;j+)anj=anj-c*a1j;b=a32/a22
9、;if(b!=0)for(j=0;jv5;j+)a3j=a3j-b*a2j;for(i=0;iintmain()intn,i,j,k,l,m;floata1010,b1010,c1010;floaty10,w10;floatz;printf(”请输入为几阶矩阵n”)；scanf(%d,&n);printf(请输入矩阵n);for(i=0;ivn;i+)for(j=0;jvn;j+)scanf(%f,&aij);for(l=0;ln-1;l+)for(k=l+l;kvn;k+)bkl=akl/all;z=akl/all;printf(b%d%d=%fnz=%fn,k,l,bkl,z);if(bk
10、l!=0);for(j=0;jn;j+)akj=akj-bkl*alj;for(i=0;in;i+)bii=1;for(i=0;in;i+)for(j=0;jn;j+)if(ij)bij=0;for(m=0;mvn;m+)for(i=0;in;i+)wm=wm+bnO*yn;for(i=0;in;i+)for(j=0;jn;j+)printf(”%f,aij);printf(n);printf(n);for(i=0;in;i+)for(j=0;jn;j+)printf(”%f,bij);printf(n);getch();M国笛谢国303国莹30Q03国莹30萤黑麻刖站色1站席商耳口3霰助常商商M霰助常商商验霰商商商商色5站商商口3站商商口商附1站商商口商附1站商商嚴口附1smQ葩色腳smQ葩色腳1smQ葩色腳3s尉国葩腳算m詣尬捕CKn国1詣尬捕CKn国M詣尬捕CKn国m詣战捕CKn国騎口捕腳口国騎口捕腳口国1騎口捕腳口国4騎口捕腳口国騎肪尉口口国騎肪尉口口国騎肪尉口口国1騎肪尉口口国4丽mMinclude-saio.hmaino(doublea15L
