数值分析报告某实验——数值积分

《数值分析报告某实验——数值积分》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、实用标准文案桂林电子科技大学数学与计算科学学院实验报告实验室：06406实验日期：2014年11月21日院（系）数学与应用数学年级、专业、班姓名成绩课程名称数值分析实验实验项目名称实验积分指导教师李光云一、实验目的通过实验掌握利用Matlab进行数值积分的操作，掌握Matlab中的几种内置求积分函数，进一步理解复化梯形，复化辛普生公式，并编程实现求数值积分二、实验原理Matlab中，有内置函数计算积分：>>z=trapz(x,y)其中，输入x，y分别为已知数据的自变量和因变量构成的向量，输出为积分值。>>z=quad(

2、fun,a,b)这个命令是使用自适应求积的方法计算积分的命令。其中，fun为被积函数，a，b为积分区间。我们还可以利用复化梯形公式三、使用仪器，材料电脑MATLAB四、实验内容与步骤1.编写复化辛普生公式的Matlab的程序。2.利用复化梯形法程序计算，记录下计算结果随着n增加的变化情况，画图与复化梯形公式的情况比较收敛速度。3.积分的原函数无法用初等函数表达，结合Matlab复化梯形程序，用描点法绘制其原函数在区间的图形。五、实验过程原始记录（数据，图表，计算等）精彩文档实用标准文案一、复化Simpson公式程序：f

3、unctions=Simpson(a,b,n)%输出s为积分的数值解，输入(a,b)为积分区间，n为等分区间的个数.h=(b-a)/(n*2);s1=0;s2=0;s=h*(f(a)+f(b))/3;%先计算特殊两点相加.fork=1:nx1=a+h*(2*k-1);%利用循环计算其他点的相加.s1=s1+f(x1);endfork=1:(n-1)x2=a+h*2*k;%利用循环计算其他点的相加.s2=s2+f(x2);ends=s+h*(4*s1+2*s2)/3;画图程序formatlong;%k为等分区间个数，t存

4、储积分值.k=2:1:40;fori=1:length(k)t(i)=Simpson(0,1,k(i));disp([k(i),t(i)]);endplot(k,t,'.','MarkerSize',20)二、Simpson(0,1,26)ans=3.14159265358779>>Untitled62.000000000000003.141568627450983.000000000000003.141591780936044.000000000000003.141592502458715.0000000000000

5、03.141592613939226.000000000000003.141592640305387.000000000000003.141592648320658.000000000000003.14159265122482精彩文档实用标准文案9.000000000000003.1415926524231710.000000000000003.1415926529697911.000000000000003.1415926532398112.000000000000003.1415926533821513.00000

6、0000000003.1415926534613414.000000000000003.1415926535074515.000000000000003.1415926535353616.000000000000003.1415926535528417.000000000000003.1415926535641118.000000000000003.1415926535715619.000000000000003.1415926535766120.000000000000003.1415926535801121.000

7、000000000003.1415926535825622.000000000000003.1415926535843223.000000000000003.1415926535856024.000000000000003.1415926535865525.000000000000003.1415926535872526.000000000000003.1415926535877927.000000000000003.1415926535881928.000000000000003.1415926535885129.0

8、00000000000003.1415926535887530.000000000000003.14159265358894精彩文档实用标准文案31.000000000000003.1415926535890932.000000000000003.1415926535892233.000000000000003.141592653