- 1、本文档共12页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
华中科技大学
数值分析实验报告
考生姓名
考生学号
班 级
指导老师 路志宏
2013年4月15日实验
实验目的:复化求积公式计算定积分
试验题目:数值计算下列各式右端定积分的近似值。
(1);
(2);
(3);
(4);
实验要求:
(1)若用复化梯形公式、复化Simpson公式和复化Gauss-Legendre I型公式做计算,要求绝对误差限为,分别利用他们的余项对每种算法做出步长的事前估计。
(2)分别用复化梯形公式、复化Simpson公式和复化Gauss-Legendre I型公式做计算。
(3)将计算结果与精确解做比较,并比较各种算法的计算量。
实验内容:
1.公式介绍
(1)复化梯形公式:
=;
余项:;
(2)复化Simpson公式:
=;
余项:;
(3)复化Gauss-Legendre I型公式:
;
余项:;
2.步长估计
(1);
则可以得到:;;
;
估计步长:;
将上述结果分别带入到复化梯形公式、复化Simpson公式和复化Gauss-Legendre I型公式的余项中可以得到:
复化梯形公式:;;
复化Simpson公式:;;
复化Gauss-Legendre I型公式:;;
(2);
则可以得到:;;
;
估计步长:;
将上述结果分别带入到复化梯形公式、复化Simpson公式和复化Gauss-Legendre I型公式的余项中可以得到:
复化梯形公式:;;
复化Simpson公式:;;
复化Gauss-Legendre I型公式:;;
(3);
则可以得到:;;
;
估计步长:;
将上述结果分别带入到复化梯形公式、复化Simpson公式和复化Gauss-Legendre I型公式的余项中可以得到:
复化梯形公式:;;
复化Simpson公式:;;
复化Gauss-Legendre I型公式:;;
(4);
则可以得到:;;
;
估计步长:;
将上述结果分别带入到复化梯形公式、复化Simpson公式和复化Gauss-Legendre I型公式的余项中可以得到:
复化梯形公式:;;
复化Simpson公式:;;
复化Gauss-Legendre I型公式:;;
++编程计算结果
(1)区间逐次分半求积法:
依据“事后误差法”,将区间逐次分半进行计算,并利用前后两次计算结果来判断误差的大小。在逐次二分进行计算时,可以用与来估计误差,这种直接用计算结果来估计误差的方法通常称作误差的事后估计法,若(为计算结果允许的误差),则停止计算,并取作为积分的近似值;否则将区间再次二分后算出,并检验不等式是否满足。
由于是区间分半,因此区间等分数必定是2的n次方。
结束输出结果选择题号N复化公式函数选择f(x)开始(2)流程图:
结束
输出结果
选择题号N
复化公式函数
选择f(x)
开始
(3)计算结果及误差:
分别对4题作复化Trapezoid、Simpson、Gauss_Legendre计算,并计算计算值与精确值之间的误差,结果如下表:
(1)计算结果表
数据类型
求积类型
(1)
计算值
区间二分相邻误差
绝对误差
区间二分等分
复化Trapezoid公式
-0.
2043509
2^11=2048
复化Simpson公式
-0.
9952908
2^5=32
Gauss_Legendre
I公式
-0.
6878386
2^5=32
(2)计算结果表
数据类型
求积类型
(2)
计算值
区间二分相邻误差
绝对误差
区间二分等分
复化Trapezoid公式
1385336
2^11=2048
复化Simpson公式
5122482
2^3=8
Gauss_Legendre
I公式
5527155
2^3=8
(3)计算结果表
数据类型
求积类型
(3)
计算值
区间二分相邻误差
绝对误差
区间二分等分
复化Trapezoid公式
4969086
2^11=2048
复化Simpson公式
4541318
2^5=32
Gauss_Legendre
I公式
4438881
2^4=16
(4)计算结果表
数据类型
求积类型
(4)
计算值
区间二分相邻误差
绝对误差
区间二分等分
复化Trapezoid公式
1970612
2^13=8192
复化Simpson公式
0756376
2^5=32
Gauss_Legendre
I公式
9317526
2^5=32
由上表中的误差分析可知,利用题目所要求的
文档评论(0)