数值分析实验数值积分和数值微分.docVIP

黑河学院数学系 实 验 报 告 课 程 名 称： 数值分析 实 验 项 目： 数值积分和数值微分 实验室名称： 数学系实验室 实 验 台 号： 02 指 导 教 师： 王国贤 学 生 学 号： 2011562003 学 生 姓 名： 柳雪梦 实验日期：__2013___年__11___月___28__日 1.实验目的及要求： 要求学生掌握牛顿-科特斯求积公式、 复化辛甫生求积法、 高斯型求积分方法，初步掌握用Matlab、C语言等软件求解问题。 2.实验内容： 应用自适应复化辛甫生求积法实习、应用软件编写程序，求给定问题的数值积分，并与真实值作比较，进行验证，同时达到熟练掌握和应用的程度。 3.实验过程: 用复化的辛生公式计算积分给定下列表格值X 50 55 60 65 Y 1.6990 1.7404 1.7782 1.8129 复化辛甫生公式 利用C语言编程，程序源代码如下： #include stdio.h #include math.h float f(float x) { return(sqrt(x)); } void main() { float a,b,h,x[10],sum=0,sum1=0,sum2=0,sum3=0; int i,m; printf(input a:); scanf(%f,a); printf(input b:); scanf(%f,b); printf(input m:); scanf(%d,m); h=(b-a)/m;x[0]=a;x[m]=b; for(i=1;im;i++) x[i]=a+h*i; sum1=f(x[0])+f(x[m]); for(i=1;i=m-1;i++,i++) sum2=sum2+f(x[i]); for(i=2;i=m-2;i++,i++) sum3=sum3+f(x[i]); sum=(sum1+4*sum2+2*sum3)*h/3; printf(the result is %f,sum); getch(); } 运行后，得到如下结果： input a:0 input b:9 input m:4 the result is 17.726210_ 所以： 计算积分：的真实值 解： 4.实验结果与分析 本次试验是使用复化的辛甫生公式求积。首先，应用C语言进行程序的编写；其次，代入所给数值得出结果。最终得到的值为17.726210比较接近真实值18。但是此方法也有不足之处，当h太小的时候，计算量较大，所以我们需要对该方法进行改进，其思想为让程序自动改变步长，从而减少计算次数，使计算量减小。 对于本次实验，我们可以更熟练的掌握自适应复化辛甫生求积法，可以更合理的应用它，从而实验目的。 5.指导教师评语 1.完成所有规定的实验内容，实验步骤正确，结果正确，分析透彻，实验严谨、认真。 优秀 2.完成绝大部分规定的实验内容，实验步骤正确，结果正确，分析透彻，实验严谨、认真。 良好 3.基本完成规定实验内容，实验步骤正确，结果正确，分析合理，实验比较严谨、认真。 合格 4.未能很好地完成规定实验内容或实验步骤不正确或结果不正确或实验态度不端正。 不及格 教师签字： 日期： 年 月 日 2