zucc科学计算实验八实验报告.docVIP

浙江大学城市学院实验报告 课程名称 科学计算 实验项目名称 数值积分 实验成绩 指导老师（签名 ） 日期 2013/12/12 一. 实验目的和要求 用Matlab软件掌握梯形公式、辛普森公式和高斯公式计算数值积分； 通过实例学习用数值积分解决实际问题。 二. 实验内容和原理 分析应用题2-1，2-2，2-3要求将问题的分析过程、Matlab源程序和运行结果和结果的解释、算法的分析写在实验报告上。 分析应用题 用复合左(右)矩形、复合梯形、复合Simpson三种积分公式计算由下表数据给出的积分。 1 2 3 4 5 6 7 0.3 0.5 0.7 0.9 1.1 1.3 1.5 其中函数为，将计算值与精确值作比较。 2-2用以下几种数值积分公式计算积分，对区间进行2等分，使用Simpson积分公式，并在该区间上使用两点和三点Gauss积分公式，输出它们的数值结果，并与精确值作比较，说明理由。 提示：两点Gauss积分公式为 三点Gauss积分公式为 使用两点guass公式计算出的结果 2-3分析应用题 一辆汽车用84s的时间驶完跑道的一圈。汽车的速度每隔6s用雷达测速仪确定，并从开始以ft / s在下表中的各项给出。 求跑道有多长？ 时间 0 6 12 18 24 30 36 42 48 54 60 66 72 78 84 速度 124 134 148 156 147 133 121 109 99 85 78 89 104 116 123 注：可从数值积分的角度来求解。 2-1分析应用题 用复合左(右)矩形、复合梯形、复合Simpson三种积分公式计算由下表数据给出的积分。 1 2 3 4 5 6 7 0.3 0.5 0.7 0.9 1.1 1.3 1.5 其中函数为，将计算值与精确值作比较。 关于这道题目，我的理解是求出从到1.5的积分，因此下面都是 function y=jingq(a) y=0.5*1.5.^2-3*cos(1.5/3)-(0.5*a.^2-3*cos(a/3)); end xk 0.3 0.5 0.7 0.9 1.1 1.3 1.5 精确 1.1683 1.2600 1.1135 0.7291 0.4768 0.1720 0 simpson：利用函数quad() xk 0.3 0.5 0.7 0.9 1.1 1.3 1.5 simpson 1.4323 1.3257 1.1660 0.9533 0.6878 0.3700 0 注：以下步长都取为0.1 对于此题的复合左矩形编的程序 function z=jux(a) x=a; y=x+sin(x/3); n=(1.5-a)/0.1; z=0; for i=1:n z=z+y* 0.1; x=x+0.1; y=x+sin(x/3); end end xk 0.3 0.5 0.7 0.9 1.1 1.3 1.5 jux 1.1683 1.2600 1.1135 0.7291 0.4768 0.1720 0 复合矩形： function z=tix(a) x1=a; x2=x1+0.1; n=(1.5-a)/0.1; z=0; for i=2:n z=z+(x1+sin(x1)+x2+sin(x2))/2; x1=x2; x2=x1+0.1; end z=z*0.1; end xk 0.3 0.5 0.7 0.9 1.1 1.3 1.5 tix 1.4873 1.5620 1.3294 0.7938 0.4259 0.4761 0 由上面的表格可以看出，在小数点位数较少的情况下利用复合simpson公式求出的积分结果与精确值一模一样 分析应用题 用以下几种数值积分公式计算积分，对区间进行2等分，使用Simpson积分公式，并在该区间上使用两点和三点Gauss积分公式，输出它们的数值结果，并与精确值作比较，说明理由。 提示：两点Gauss积分公式为 三点Gauss积分公式为 使用两点guass公式计算出的结果 G2=pi/8*(sin(pi/8-pi/24*sqrt(3))+sin(pi/8+pi/24*sqrt(3)))+pi/8*(sin(3*pi/8-pi/24*sqrt(3))+sin(pi*3/8+pi/24*sqrt(3))) G2 = 0.99991016676