实验10 使用矩形、梯形、辛普森和高斯公式求积分.docVIP

PAGE PAGE １ 数学实验报告 实验序号： 日期： 班级 姓名 学号 实验名称 使用矩形、梯形、辛普森和高斯公式求积分 问题背景与实验目的： 1、利用MATLAB掌握梯形公式、辛普森公式和随机模拟方法计算数值积分； 2、通过实例学习用数值积分和数值微分解决实际问题 实验内容： 用矩形、梯形(trapz)和辛普森(simp)三种公式计算由下表数据给出的积分. k 1 2 3 4 5 6 7 0.3 0.5 0.7 0.9 1.1 1.3 1.5 0.3895 0.6598 0.9147 1.1611 1.3971 1.6212 1.8325 已知该表数据为函数 所产生，将计算值与精确值作比较。 选择一些函数用梯形(trapz)、辛普森(quad)两种方法计算积分. 改变步长(对梯形公式)，改变精度要求(对辛普森公式)，进行比较、分析。如下函数供选择参考： d. , 实验原理与数学模型： 实验过程记录（含：基本步骤、主要程序清单及异常情况记录等）： 第一题 h=0.2; x=0.3:h:1.5; y=x+sin(x/3); format long z1=sum(y(1:6))*h z1 = 1.27417511985366 z2=sum(y(2:7))*h z2 = 1.59009354424514 z3=trapz(y)*h z3 = 1.43213433204940 z4=quad(x+sin(x/3),0.3,1.5) z4 = 1.43226481015968 z5=quadl(x+sin(x/3),0.3,1.5) z5 = 1.43226481016296 format short y=x+sin(x/3) y = 0.3998 0.6659 0.9312 1.1955 1.4585 1.7199 1.9794 y1=[0.3895 0.6598 0.9147 1.1611 1.3971 1.6212 1.8325]; y-y1 ans = 0.0103 0.0061 0.0165 0.0344 0.0614 0.0987 0.1469 第二题 步长为0.2 h=0.2; x=-2:h:2； y=(1./sqrt(2*pi))*exp(-x.^2./2) z1=sum(y(1:20))*h z1 = 0.9538 z2=sum(y(2:21))*h z2 = 0.9538 z3= quad((1./sqrt(2*pi))*exp(-x.^2./2),-2,2) z3 = 0.9545 z4=quadl((1./sqrt(2*pi))*exp(-x.^2./2),-2,2) z4 = 0.9545 步长为0.15 h=0.15; x=-2:h:2； y=(1./sqrt(2*pi))*exp(-x.^2./2)； zz1=sum(y(1:26))*h zz1 = 0.9472 zz2=sum(y(2:27))*h zz2 = 0.9490 zz3=quad((1./sqrt(2*pi))*exp(-x.^2./2),-2,2) zz3 = 0.9545 zz4=quadl((1./sqrt(2*pi))*exp(-x.^2./2),-2,2) zz4 = 0.9545 实验结果报告： 第一题 z1=sum(y(1:6))*h z1 = 1.27417511985366 z2=sum(y(2:7))*h z2 = 1.59009354424514 z3=trapz(y)*h z3 = 1.43213433204940 z4=quad(x+sin(x/3),0.3,1.5) z4 = 1.43226481015968 z5=quadl(x+sin(x/3),0.3,1.5) z5 = 1.43226481016296 format short y=x+sin(x/3) y = 0.3998 0.6659 0.9312 1.1955 1.4585 1.7199 1.9794 y1=[0.3895 0.6598 0.9147 1.1611 1.3971 1.6212 1.83