原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
《数值分析》实验报告
标题:用复化梯形公式和复化辛普森公式求积分
姓名:刘先俊班级:信科142学号:2014035031
一、实验目的
1.掌握复合梯形公式及复化辛普森公式的基本思想。2.编程实现用复合梯形公式求积分。
3.掌握利用复化辛普森公式计算积分。
4.熟悉matlab软件的使用。.
二、实验题目
1、根据matlab编写复合梯形公式求如下积分问题
1,,,,,
0
2、根据matlab编写复合辛普森公式求如下积分问题
14,,,1+,^20
三、实验步骤
1、复化梯形公式
1基础理论?
b?a将区间[ab]划分为n等分，分点为,=,+,?，i=0,1,?,n,h=,在每个,n
小区间[,,,]上应用梯形公式,,+1
,,?,,,,,,,,,,,,,?,,+,,,i=1,2,?,n,?1,2,,?1
则得
,,,,,?,,,,,,,,,,,,,=,,,(,),,?,,,+,,,?1,2,,,?1,=1,=1
若记上式右端?,并注意到,=,,,=,则上式可写成,0,
《数值分析》实验报告
,?,?1,,,,,?=[,,+2,,+,(,)](1),(,),,?,,,=1,,,2
2基础算法?
构造trap.m文件
functiont=trap(fun,a,b,n)
%用途:用复化梯形公式求积分
%格式:s=trap(fun,a,b,n).fun为被积函数，a,b为积分区间的左右端点
%n为区间的等分数，S返回数值积分值
h=(b-a)/n;
sum=0;
fork=1:n-1
x=a+k*h;
sum=sum+feval(f,x);
end
t=h/2*(feval(f,a)+feval(f,b)+2*sum);2、复化辛普森公式
1基础理论?
1将区间[ab]划分为n等分，记区间[,,,]的中点为,，则得如下的复化,?1,,?2
辛普森公式
,??,,,,1,,,,,(,),,?,=,,,+4,,,,+,,,=[,,+,,,=1,?1,,?,662
,,?1,1,,,4,,,,+2,,+,(,)](2),?1,=1,,?2
2基础算法?
构造simp.m文件
functions=simp(fun,a,b,n)%用途:用复化辛普森公式求积分
%格式:s=simp(fun,a,b,n).fun为被积函数，a,b为积分区间的左右端点
%n为区间的等分数，S返回数值积分值
h=(b-a)/n;s1=0;s2=0;
fork=1:(n-1)
x=a+k*h;
s1=s1+feval(fun,x);
end
fork=0:(n-1)
x=a+h*(k+1/2);
s2=s2+feval(fun,x);
end
s=h/6*(feval(fun,a)+feval(fun,b)+2*s1+4*s2);
《数值分析》实验报告
四、运行结果
1、复化梯形公式运行结果
调用MATLAB程序trap.m，在MATLAB命令窗口输入:formatlong
fun=inline(exp(x));trap(fun,0,1,5)
计算结果为:ans=1.724005619782788
实验数据及结果记录:
实验数据1、trap(fun,0,1,5)ans=1.724005619782788实验数据2、trap(fun,0,1,10)ans=1.719713491389315实验数据3、trap(fun,0,1,50)ans=1.718339104138158实验数据4、trap(fun,0,1,100)ans=1.718296147450418实验数据5、trap(fun,0,1,1000)ans=1.7182819716491962、复化辛普森公式运行结果
调用MATLAB程序simp.m，在命令窗口输入:
formatlong
fun=inline(‘4./(1+x.^2)’);
s=simp(fun,0,1,10)
计算结果为:ans=3.141592652969785
实验数据及结果记录:
simp(fun,0,1,10)实验数据1、
ans=3.141592652969785实验数据2、simp(fun,0,1,50)ans=3.141592653589753实验数据3、simp(fun,0,1,100)ans=3.141592653589793实验数据4、simp(fun,0,1,500)ans=3.141592653589796实验数据5、simp(fun,0,1,1000)ans=3.14
专注地铁、铁路、市政领域安全管理资料的定制、修改及润色，本人已有7年专业领域工作经验，可承接安全方案、安全培训、安全交底、贯标外审、公路一级达标审核及安全生产许可证延期资料编制等工作，欢迎大家咨询~
文档评论（0）