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

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

ID:29446925

大小:60.00 KB

页数:6页

时间:2018-12-19

数值分析报告某实验——数值积分_第1页
数值分析报告某实验——数值积分_第2页
数值分析报告某实验——数值积分_第3页
数值分析报告某实验——数值积分_第4页
数值分析报告某实验——数值积分_第5页
资源描述:

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

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

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。