没有合适的资源?快使用搜索试试~ 我知道了~
MATLAB数学实验.pdf
1.该资源内容由用户上传,如若侵权请联系客服进行举报
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
版权申诉
0 下载量 184 浏览量
2021-09-14
15:39:57
上传
评论
收藏 912KB PDF 举报
温馨提示
试读
24页
。。。
资源推荐
资源详情
资源评论
MATLAB 数学实验
实验三 圆周率的计算
学号: 姓名:XX
一、 实验目的
1. 本实验涉及概率论、定积分、三角函数等有关知识,要求掌握
计算 π 的三种
方法及其原理。
2. 学习和掌握数学软件 MATLAB 的使用方法。
二、 实验内容
圆周率是一个极其驰名的数。从有文字记载的历史开始,这个数
就引起了外行人和学者们的兴趣。作为一个非常重要的常数,圆周率
最早是出于解决有关圆的计算问题。仅凭这一点,求出它的尽量准确
的近似值,就是一个极其迫切的问题了。事实也是如此,几千年来作
为数学家们的奋斗目标,古今中外一代又一代数学家为此献出了自己
的智慧和劳动。回顾历史,人们对 π 的认识过程,反映了数学和计算
技术发展情形的一个侧面。π 的研究,在一定程度上反映这个地区或
时代的数学水平。德国数学家康托说:“历史上一个国家所算的圆周
率的准确程度,可以作为衡量这个国家当时数学发展水平的指标。”
直到 19 世纪初,求圆周率的值还是数学中的头号难题。
1. 圆周率的计算方法
古人计算圆周率,一般是用割圆法。即用圆的内接或外切多边形
来逼近圆的周长。 Archomedes用正 96 边形得到 35 位精度;刘徽用
正 3072 边形得到 5 位精度;Ludolph Van Ceulen用正 2^62 边形得
到了 35 位精度。这种基于几何的算法计算量大,速度慢,吃力不讨好。
随着数学的发展,数学家们在进行数学研究时有意无意得发现了许多
计算圆周率的公式。下面挑选一些经典的常用公式加以介绍。除了这
些经典公式外,还有很多其他公式和由这些经典公式衍生出来的公式,
就不一一列举了。
1) Machin 公式
239
1arctan 451arctan
16-=π ()1
21...753arctan 121753--++-+-=--n x x x x x x n n 这个公式由
英国天文学教授 John Machin 于 1706 年发现。他利用这个公式计算
到 100
位的圆周率。Machin 公式每计算一项可以得到 1.4 位的十进制精
度。因为它的还算过程中
被乘数和被除数都不大于长整数,所以可以很容易地在计算机上
编程实现。
还有很多类似于 Machin 公式的反正切公式。在所有这些公式中,
Machin 公式似乎是最快的了。虽然如此,如果要计算更多的位数,比
如几千万位,Machin 公式就力不从心了。下面介绍的算法,在 PC 机
上计算大约一天时间,就可以得到圆周率的过亿位的精度。这些算法
用程序实现起来比较复杂。因为计算过程中涉及两个大数的乘法运算,
要用 FFT(Fast Fourier Transform)算法。FFT 可以将两个大数的乘法
运算时间由 O (n^2)缩短为 O (nlog(n)).
2) Ramanujan 公式
()()
()∑∞
=+?
=
044499263901103!4!4229801
n n n n n n π 1914 年,印度学家 Srinivasa Ramanujan 在他的
论文里发表了一系列共 14 条圆周率的计
算公式,这是其中之一。这个公式每计算一项可以得到 8 位的十
进 制 精 度 。 1985 年 Gosper 用 这 个 公 式 计 算 得 到 看 圆 周 率 的
17500000 位。
1989 年,David&Gregory Chudnovsky 兄弟将 Ramanujan 公
式改良成为:
()()()
()2
330364032054514013413591409!!3!61121
+∞
=+?-=∑n n n n
n n n π 这个公式被称为 Chudnovsky 公式,每计算一项可以得
到 15 位的十进制精度。1994 年
Chudnovsky 兄弟利用这个公式计算得到了 4044000000 位。
Chudnovsky 公式的另一个更方便于计算机编程的形式是:
()()()()()∑
∞
=-+=
33640320!3!13591409
545140134!610005
426880n n
n n n n π
3) AGM(Arithmetic-Geometric Mean)算法
Gauss-Legendre 公式: 初值:
1==x a 2
1=b 4
1=
c 重复计算:
a y = 2
b a a +=
by b = ()2
y a x c c --= x x 2= 最后计算:
()c
b a 42
+=π
这个公式每迭代一次将得到双倍的十进制精度,比如要计算 100
万位,迭代 20 次就够
了。1999 年 9 月 Takahashi和 Kanada 用这个计算得到了圆周率
的 206 158 430 00 位,创出新的世界纪录。
4) Borwein 四次迭代式 初值:
2460-=a 120-=
y
重复计算: ()()
4144
14
1
1111n
n
n y y y -+--=
+
()()
2
111324
11121++++++++-+=n n n n n n n y y y y a a
最后计算:
n
a 1
=
π 这个公式由 Jonathan Borwein 和 Peter Borwein 于 1985 年
发表,它四次收敛于圆周率。 5) Bailey-Borwein-Plouffe 算法
+-+-+-+=∑
∞
=6815814821841610
n n n n n n
π 这个公式简称 BBP 公式,由 David Bailey,Peter Borwein 和
剩余23页未读,继续阅读
资源评论
苦茶子12138
- 粉丝: 1w+
- 资源: 6万+
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
最新资源
- C++前馈神经网络.zip
- PCB灯光画(哆啦A梦)DIY制作.rar
- 基于c语言-文档+程序 简单潮流计算,这个压缩包包含两种方法
- 基于粒子群优化随机森林的数据分类预测.zip
- 基于matlab实现AR法模拟脉动风场风速.rar
- 基于matlab实现cpp电力系统的牛顿-拉夫逊潮流计算 有程序和PPT.rar
- 基于matlab实现的对电力系统进行交流潮流计算,计算方法为极坐标系下的牛顿拉夫逊法.rar
- 基于matlab实现的蒙特卡洛仿真8psk的蒙特卡洛仿真.rar
- 基于matlab实现的牛顿拉夫逊潮流计算.rar
- 基于matlab实现的牛顿-拉夫逊法的潮流计算程序,并提供9节点、30节点的算例 .rar
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功