内江师范学院数值分析实验三资料.docVIP

内江师范学院 数 值 分 析 实 验 报 告 册 编制 张莉 审定 牟廉明 专业： 班级： 级 班 学号： 姓名： 数学与信息科学学院 2013年9月 说 明 一、学生在做实验之前必须要准备实验，主要包括预习与本次实验相关的理论知识，熟练与本次实验相关的软件操作，收集整理相关的实验参考资料，要求学生在做实验时能带上充足的参考资料；若准备不充分，则学生不得参加本次实验，不得书写实验报告； 二、要求学生要认真做实验，主要是指不得迟到、早退和旷课，在做实验过程中要严格遵守实验室规章制度，认真完成实验内容，极积主动地向实验教师提问等； 三、各个实验按照学生水平分别设置了A、B、C、D四个等级，其中对应的难度系数为1、0.8、0.7、0.6，也可根据实际完成情况制定相应地的难度系数，但总体保证难度排序为A级难度最大，B级次之，C级较易，D级最简单。 四、学生可以根据自己对各个实验涉及到的知识点掌握的程度自由选取A、B、C、D等级的实验题目。 五、学生要认真工整地书写实验报告，实验报告的内容要紧扣实验的要求和目的，不得抄袭他人的实验报告； 四、根据实验准备、实验态度、实验报告的书写、实验报告的内容进行综合评定，并给出实验成绩评定分。 实验名称： 实验三 数值积分 指导教师： 吴开腾 张莉 实验时数： 2 实验设备：安装了Matlab、C++、VF软件的计算机 实验日期：2013年 10 月 29 日 实验地点： 第五教学楼北902 实验目的： 1. 掌握数值积分的基本思想和基本步骤。 2. 理解各类数值积分方法的优缺点，并能自行编程求解。 3. 比较各类数值积分的代数精度，体会复化、变步长及其加速的思想和实现步骤。 实验准备： 1. 在开始本实验之前，请回顾教科书的相关内容； 2. 需要一台准备安装Windows XP Professional操作系统和装有数学软件的计算机。 实验内容及要求 A题 计算积分 。 B题 用Romberg法求函数积分，精度为。 C题 用复化梯形公式计算积分 并分析剖分区间数对误差的影响，取，积分精确值。 D题 计算。 说明：实验过程应包括对问题的简要分析、求解方法、求解步骤、程序及其必要的图表等内容。 实验过程： A题的实验过程 实验中问题的重述 计算积分 实验分析 由于积分为是从负无穷到正无穷积分，用常规的求积分的int函数求解时，会显示Warning:Explicit integral could not be found。说明函数y=exp(sin(x)-x.^2/100)不是初等原函数，可用数值积分进行分析。为方便计算，将积分区间分为两个区间，以0为分界点，分为和，则有： 所以在两个区间分别计算，然后将结果相加。 实验求解过程 首先利用辛普生法，从结果得到：在积分区间积分上限增加，积分值趋于12.3599；在积分区间积分下限减小，积分值趋于10.0804； 然后利用牛顿－柯特斯法，观察结果得到，的积分区间，积分上限增加，积分值趋于12.3599；在的积分区间，积分下限减小，积分值趋于10.0804； 实验求解结果 从上述两种方法的实验结果表明，辛普生法和牛顿－柯特斯法在的积分区间，积分上限增加，积分值均趋于12.3599；在的积分区间，积分下限减小，积分值均趋于10.0804，故在两个积分区间都是可积的，所以 故I=10.0804+12.3599=22.4403 B题的实验过程 实验中问题的重述 用Romberg法求函数积分，精度为。 实验求解程序 M文件 function [quad,R] = Romberg(f,a,b,eps) %f表示被积函数 %a，b表示积分区间的端点 %eps表示精度 %quad是用Romberg加速算法求得的积分值 %R为Romberg表 %err表示误差的估计 h=b-a; R(1,1)=h*(feval(f,a)+feval(f,b))/2; M=1;J=0;err=1; while erreps J=J+1; h=h/2; S=0; for p=1:M x=a+h*(2*p-1); S=S+feval(f,x); end R(