大学数学实验报告----怎样计算∏
圆周率的实验报告
圆周率的实验报告圆周率的实验报告引言：圆周率（π）是数学中一个重要的常数，它表示圆的周长与直径的比值。
圆周率的数值约等于3.14159，是一个无限不循环的小数。
在本次实验中，我们将通过不同的方法来计算圆周率，并探讨其性质和应用。
实验一：测量圆的周长和直径首先，我们需要测量一个圆的周长和直径，以便计算圆周率。
选择一个圆形物体，如一个硬币或者一个圆盘，使用一个软尺或者卷尺测量其周长和直径。
将测量结果记录下来，并计算周长与直径的比值。
实验二：使用几何方法计算圆周率在几何学中，我们可以通过正多边形的外接圆和内接圆来近似计算圆周率。
选择一个正多边形，如正六边形或正十二边形，测量其边长和内切圆的半径。
然后，计算正多边形的周长与内切圆的周长的比值。
随着正多边形的边数增加，这个比值会越来越接近圆周率。
实验三：使用概率方法计算圆周率概率方法是一种基于随机事件的方法来计算圆周率。
我们可以在一个正方形内随机撒点，并计算落在正方形内的点中，落在内切圆内的点的比例。
根据概率理论，这个比例会接近于圆的面积与正方形的面积之比，即π/4。
通过将这个比例乘以4，我们可以得到一个近似的圆周率值。
实验四：使用级数方法计算圆周率在数学中，圆周率可以通过级数来计算。
其中一个著名的级数是莱布尼茨级数：π/4 = 1 - 1/3 + 1/5 - 1/7 + 1/9 - ...通过不断计算级数的和，我们可以逼近圆周率的数值。
在实验中，我们可以计算不同级数的和，并观察其逼近圆周率的速度。
实验五：使用计算机模拟计算圆周率计算机的出现为计算圆周率提供了更加精确和高效的方法。
我们可以使用计算机编写程序，通过数值方法来计算圆周率。
例如，可以使用蒙特卡洛方法，在一个正方形内随机生成大量点，并计算落在内切圆内的点的比例。
根据概率理论，这个比例会逼近圆周率的数值。
结论：通过以上实验，我们可以发现不同方法计算的圆周率值会有一定的误差，但随着方法的改进和精确度的提高，这个误差可以被不断减小。
计算pi
一、实验目的探索精确计算π值的方法，并且比较不同方法之间的不同之处和优缺点。
掌握数值积分的辛普森公式。
二、问题描述1. 任务11) 用反正切函数的幂级数展开式结合有关公式求π，若要精确到40位、50位数字，试比较简单公式和Machin 公式所用的项数。
2) 验证公式111=arctan arctan arctan 4258π++ 试试此公式右端做幂级数展开完成任务1所需要的步数。
2. 任务2用数值积分计算π，分别用梯形法和Simpson 法精确到10位数字，用Simpson 法精确到15位数字。
3. 任务3用Monte Carlo 法计算π，除了加大随机数，在随机数一定时可重复算若干次后求平均值，看能否求得5位精确数字？设计方案用计算机模拟Buffon 实验4. 任务4利用积分20(1)!!sin !!2n n xdx n ππ-=⎰ ，n 为奇数 推导公式224422213352121n n n n π=-+ ……… 用此公式计算π的近似值，效果如何？5. 任务5利用学过的知识（或查阅资料），提出其他计算π的方法（先用你学过的知识证明），然后实践这种方法。
对你在实验中应用的计算π的方法进行比较分析。
6. 任务6e 是一个重要的超越数1e lim 1)n n n→∞=+（ 1111...2!!(1)!e e n n θ=++++++ 试用上述公式或其他方法近似计算e 。
三、问题解法1. 任务11) 根据幂级数展开的相关知识，易知：24122211(1)1n n x x x x--=-+-+-++……… 因为21(arctan )'1x x =+，故可以求得arctan x 的幂级数展开式为： 35211arctan (1)3521n n x x x x x n --=-+-+-+-……… 当x=1时，-11111--(-1)4352-1n n π=+⋯⋯++⋯ 当叠加了十万次以后得到结果π=3.141582654…只有五位有效数字，可见其精度与效率极低。
连乘符号∏的运算法则
连乘符号∏的运算法则
连乘符号∏的运算法则：
连乘符号是“∏”，代表“求乘积”。
用法：上下添加的为求乘积的初始值和终止值，例如：符号下面可写“i=1”，上面写“n”，就代表后面的求积式子中的i从1开始一直加到n。
即（1+D1/P1)(1+D2/P2)……（1+Dn/Pn)。
∏是希腊字母，即π的大写形式，在数学中表示求积运算或直积运算，形式上类似于Σ。
∏，这个符号就是连续求积的意思，把满足∏这个符号下面条件的所有项,都乘起来，求积；
∏是各项连乘的运算符号, 读大写的π（pai）。
∏i=1（符号下面）n（符号上面）ai(符号右面)表示a1×a2....×an。
符号下面表示右面式子可变参量的下限（或初值）。
符号上面表示右面式子可变参量的上限（或终值）。
实验二 怎样计算Pi
数学实验实验报告学院：数学与统计学院班级：数学与应用数学3班学号：0314姓名：康萍时间：实验二怎样计算一、实验目的分别用下列三种方法计算π的近似值，并比较三种方法的精确度： 数值积分法：通过使用编写梯形公式和辛普森公式的程序语言计算π。
泰勒级数法：利用反正切函数泰勒级数计算π。
蒙特卡罗（Monte Carlo ）法：通过使用编写蒙特卡罗公式的程序语言来计算π。
二、实验环境基于Windows 环境下的软件。
三、实验的基本理论和方法1、数值积分法以单位圆的圆心为原点建立直角坐标系，则单位圆在第一象限内的部分G 是一个扇形，由曲线])1,0[(12∈-=x x y 及两条坐标轴围成，它的面积4π=S 。
算出了S 的近似值，它的4倍就是π的近似值。
而扇形面积S 实际上就是定积分4112π=-⎰dx x 。
与π有关的定积分有很多，比如211x +的定积分411102π=+⎰dx x 就比21x -的定积分更容易计算，更适合于用来计算π。
一般地，要计算定积分()dx x f ba ⎰，也就是计算曲线()x f y =与直线b x a x y ===,,0所围成的曲边梯形G 的面积S 。
为此，用一组平行于y 轴的直线()b x x x x x a n i x x n n i =<<<<<=-≤≤=-1210,11 将曲边梯形T 分成n 个小曲边梯形，总面积S 分成这些小曲边梯形的面积之和。
如果取n 很大，使每个小曲边梯形的宽度都很小，可以将它上方的边界()()i i x x x x f ≤≤-1近似的看作直线段，将每个小曲边梯形近似的看作梯形来求面积，就得到梯形公式。
如果更准确些，将每个小曲边梯形的上边界近似的看作抛物线段，就得到辛普森公式。
具体公式如下：梯形公式 设分点11,,-n x x 将积分区间],[b a 分成n 等份，即()n i n a b i a x i ≤≤-+=0,/。
Pi值的计算(mathematica数学实验报告)
arctan x x x3 x5 (1)k1 x 2k1
35
2k 1
来计算 。 从反正切函数的泰勒级数，进行如下编程来计算 ，实验运行如下：
从实验过程可以看出，这种方法花费的时间很长。原因是当 x=1 时得到的 arctan1的
展开式收敛太慢。要使泰勒级数收敛得快，容易想到，应当使 x 的绝对值小于 1，最好
实验基本理论和方法：
1、Mathematica中常用绘图函数Plot在绘制高次函数时的方法；
2、计算圆周率 的数值积分法、泰勒级数法、蒙特卡罗法，并且利用特定的公式来
计算圆周率 。
实验内容和步骤：
（1）数值积分法计算
半径为 1 的圆称为单位圆，它的面积等于 。只要计算出单位圆的面积，就算出了 。 在坐标轴上画出以圆点为圆心，以 1 为半径的单位圆（如下图），则这个单位圆在第一 象限的部分是一个扇形，而且面积是单位圆的 1/4，于是，我们只要算出此扇形的面积， 便可以计算出 。
0
4
利用 Mathematics 编程计算上式，过程如下：
从而得到 的近似值为 3.14159265358979323846264338328，可以看出，用这种方法 计算所得到的 值是相当精确的。n 越大，计算出来的扇形面积的近似值就越接近 的 准确值。
（2）泰勒级数法计算 利用反正切函数的泰勒级数
只要计算出单位圆的面积就算出了为半径的单位圆如下图则这个单位圆在第一象限的部分是一个扇形而且面积是单位圆的14于是我们只要算出此扇形的面积便可以计算出在计算扇形面积时很容易想到使用数学分析中积分的方法第一象限中的扇形由曲线及两条坐标轴围成实际操作中我们不能准确地计算它的面积于是就通过分割的方法将其划分为许多小的梯形通过利用梯形的面积近似于扇形面积来计算利用mathematics编程计算上式过程如下
用蒙特卡罗方法计算π值实验报告
本科生实验报告实验课程蒙特卡罗模拟学院名称核技术与自动化工程学院专业名称核技术及应用学生姓名王明学生学号**********指导教师邮箱****************实验成绩二〇一七年九月二〇一八年一月实验一、选择一种编程语言模拟出π的值一、实验目的1、理解并掌握蒙特卡罗模拟的基本原理；2、运用蒙特卡洛思想解决实际问题；3、分析总结蒙特卡洛解决问题的优缺点。
二、实验原理用蒙特卡洛思想计算π的值分为如下几部：第一步构建几何原理：构建单位圆外切正方形的几何图形。
单位圆的面积为S0=π，正方形的面积S1=4；第二步产生随机数进行打把：这里用MATLAB产生均匀随机数。
分别生产均匀随机数(x,y)二维坐标。
X,y的范围为-1到1.总共生成N个坐标(x,y).统计随机生成的坐标(x,y)在单位圆内的个数M。
第三步打把结构处理：根据S0/S1=M/N计算出π的值。
因此π=4*M/N。
第四步改变N的值分析π的收敛性：总数1000开始打把，依次增长10倍到1百万个计数。
三、实验内容1、用matlab编写的实验代码，总计数率为1000。
zfx_x=[1,-1,-1,1,1];zfx_y=[1,1,-1,-1,1];plot(zfx_x,zfx_y)axis([-3 3 -3 3]);hold on;r=1; theta=0:pi/100:2*pi;x=r*cos(theta); y=r*sin(theta);rho=r*sin(theta);figure(1)plot(x,y,'-')N=1000;mcnp_x=zeros(1,N);mcnp_y=zeros(1,N);M=0;for i=1:Nx=2*(rand(1,1)-0.5);y=2*(rand(1,1)-0.5);if((x^2+y^2)<1)M=M+1;mcnp_x(i)=x;mcnp_y(i)=y;endendplot(mcnp_x,mcnp_y,'.')PI1=4*M/N;2、用matlab绘制的图形四、实验结果1.当模拟总计数为1000时，某次计算结果： PI=3.128。
怎样计算Pi
学 实 验
（2） ） （3） ）
使用上面的两个公式（ ）或者（ ）， ），并应用级展 使用上面的两个公式（2）或者（3），并应用级展 ），则可得计算速度更快的公式 式（1），则可得计算速度更快的公式，例如： ），则可得计算速度更快的公式，例如：
上页 下页 退出
In[20]:=
哈 工 程 大 学 数 学 实
上页 下页 退出
一、祖冲之的计算方法
哈 工 程 大 学 数 学 实 验
据史料记载，某年祖冲之家新建了一所房屋，他将最大的 据史料记载，某年祖冲之家新建了一所房屋， 一个房间铺成地板，选择房间的中心点作为圆心， 一个房间铺成地板，选择房间的中心点作为圆心，根据房间的 大小，画出一个最大的圆，首先在圆周上作出一个三角形， 大小，画出一个最大的圆，首先在圆周上作出一个三角形，将 圆周三等分，然后将每段弧两等分，即将圆周六等分， 圆周三等分，然后将每段弧两等分，即将圆周六等分，依此进 行下去，一直进行到不能够分下去为止。 行下去，一直进行到不能够分下去为止。 则第N次 则第 次，圆周被分 等分， 成3N等分，用3N等分后 等分 等分后 正3N边形的边长来近似 边形的边长来近似 计算3N等分后的每段弧 计算 等分后的每段弧 长，祖冲之就是用这种 方法， 方法，将圆周率的精确 度计算到3.15926与 度计算到 与 3.15927之间， 之间， 之间
上页 下页 退出
二、数值积分方法
哈 工 程 大 学 数 学 实 验
4 ∫0 1 + x 2 dx = π
1
因此，只要计算出来上面的积分近似值，就得到了 的 因此，只要计算出来上面的积分近似值，就得到了π的 近似值。 近似值。
上页 下页 退出
哈 工 程 大 学 数 学 实 验 电 子 教 案
pi的计算 实验报告
Ramanujan公式1914年，印度数学家Srinivasa Ramanujan在他的论文里发表了一系列共14条圆周率的计算公式，这是其中之一。
这个公式每计算一项可以得到8位的十进制精度。
1985年Gosper用这个公式计算到了圆周率的17,500,000位。
3、AGM(Arithmetic-Geometric Mean)算法Gauss-Legendre公式：初值：重复计算：最后计算：这个公式每迭代一次将得到双倍的十进制精度，比如要计算100万位，迭代20次就够了。
1999年9月Takahashi和Kanada用这个算法计算到了圆周率的206,158,430,000位，创出新的世界纪录。
4、Borwein四次迭代式：初值：重复计算：最后计算：这个公式由Jonathan Borwein 和Peter Borwein 于1985年发表，它四次收敛于圆周率。
5、Bailey-Borwein-Plouffe 算法014211()1681848586n n n n n n π∞==---++++∑这个公式简称BBP 公式，由David Bailey, Peter Borwein 和Simon Plouffe 于1995年共同发表。
它打破了传统的圆周率的算法，可以计算圆周率的任意第n 位，而不用计算前面的n-1位。
这为圆周率的分布式计算提供了可行性。
1997年，Fabrice Bellard 找到了一个比BBP 快40％的公式：第三部分：对于π的几种计算的研究和讨论： 1、数值积分法（I ）利用积分公式⎰-=10214dx x π计算πn=10 ans =; n=20 ans =; n=50 ans =; n=100 ans =; n=200 ans =; n=500 ans =; n=1000 ans =; n=2000 ans =;半径为1的圆称为单位圆，它的面积等于π。
只要计算出单位圆的面积，就算出了π。
实验2 怎样计算圆周率
这就是梯形公式。
Simpson公式：仍用分点 xi = a + i(b − a) / n, 0 ≤ i ≤ n 公式： 公式 将积分区间 [a, b] 分成n等分，再作每个小区间 1 x = a + ( i − )( b − a ) / n [ x i −1 , x i ] 的中点 ，将第i个小曲边 2 梯形的上边界近似地看作经过三点 ( x, f ( x)) ( x = x , x , x ) 的抛物线段，则可求得
π
1 β = arctan 239
π
1 1 = 4α − β = 4 arctan − arctan 4 5 239
于是得到
π = 16 arctan − 4 arctan
1 5 1 239
这称为Maqin公式，利用 arctan x 的Taylor展开式 1 1 arctan arctan ， 的近似值，就可以由Maqin 求出 239 5 公式求出 π 的近似值。 Taylor级数是无穷级数，实际计算时必然 只能取它的前n项，导致截断误差
数 学 实 验
π的计算
南通大学信息与计算科学教研室
实验2 实验 怎样计算圆周率 π 一、实验目的 本实验主要涉及数值积分法、 Taylor级数法和Monte Carlo法， 要求尝试利用所知道的数学知 识来进行近似计算。
二、实际问题 众所周知，圆周率 π 是平面上圆的周长与 直径之比，它等于3.141 592 6…。古代人 把3作为它的近似值。古希腊阿基米德 10 1 （Archimedes）曾得到 3 71 < π < 3 7 。我国宋 22 代的祖冲之得到 π 的近似值 7 （约率）和 355 （密率），后者化为小数后等于3.141 113 π 的准确值的误差在 10−6 以下。 592…，与
派的计算实验报告
一、实验目的1. 了解π的定义及其重要性。
2. 掌握使用不同方法计算π的原理和步骤。
3. 比较不同方法计算π的精度和效率。
二、实验原理π（派）是一个数学常数，表示圆的周长与直径的比值。
在数学、物理、工程等领域中，π具有广泛的应用。
π的近似值通常取3.14159，但实际上π是一个无理数，其小数位数无限不循环。
本实验通过以下几种方法计算π的近似值：1. 牛顿迭代法2. 阿基米德法3. 蒙特卡洛法三、实验步骤1. 牛顿迭代法（1）选择初始值x0，通常取x0=3。
（2）根据牛顿迭代公式x1 = x0 - f(x0)/f'(x0)计算下一个近似值。
（3）重复步骤（2），直到满足精度要求。
2. 阿基米德法（1）在坐标轴上画一个半径为1的圆。
（2）画一个内接正六边形，计算其面积S1。
（3）画一个外切正六边形，计算其面积S2。
（4）计算π的近似值：π ≈ 6(S2 - S1)。
3. 蒙特卡洛法（1）在坐标轴上画一个半径为1的圆。
（2）随机生成N个点，计算其中落在圆内的点的数量M。
（3）计算π的近似值：π ≈ 4M/N。
四、实验结果与分析1. 牛顿迭代法选择初始值x0=3，精度要求为10^-6。
经过迭代，计算得到π的近似值为3.1415926535。
2. 阿基米德法选取内接正六边形的边长为1，外切正六边形的边长为2。
计算得到π的近似值为3.1415926535。
3. 蒙特卡洛法选取N=10000，计算得到π的近似值为3.1415926535。
三种方法计算得到的π近似值相差不大，但牛顿迭代法和阿基米德法在计算过程中具有较高的精度。
蒙特卡洛法虽然精度较低，但计算简单，适合大规模计算。
五、实验结论1. 本实验通过三种方法计算π的近似值，结果表明，牛顿迭代法和阿基米德法具有较高的精度。
2. 蒙特卡洛法虽然精度较低，但计算简单，适用于大规模计算。
3. π在数学、物理、工程等领域具有广泛的应用，掌握计算π的方法具有重要意义。
π的计算公式简单方法
π的计算公式简单方法π是一个著名的数学常数，它代表圆的周长与直径的比值，也被称为圆周率。
在数学中，π是一个重要的数，它出现在许多公式中，如圆的面积公式、弧长公式、三角函数公式等。
因此，计算π的方法一直备受关注。
在过去，人们采用的是几何方法和机械方法来计算π，这些方法比较繁琐且耗时。
随着计算机技术的发展，人们开始使用计算机来计算π，但这仍然需要很高的计算能力和时间。
如今，有许多简单的方法可以用来计算π，这篇文章将介绍其中一种简单的方法。
该方法是由印度数学家拉马努金提出的，其基本思想是使用无穷级数来计算π。
这个级数是著名的莱布尼茨级数，它的形式如下：π/4 = 1 - 1/3 + 1/5 - 1/7 + 1/9 - 1/11 + ...在这个级数中，每个项都是相继加或减的奇数分数，随着项数的增加，级数的和越来越接近π/4。
因此，我们可以通过计算这个级数来得到π的近似值。
下面是通过计算前几项级数得到的π的近似值：1项：π/4 ≈ 12项：π/4 ≈ 1 - 1/3 = 2/33项：π/4 ≈ 1 - 1/3 + 1/5 = 8/154项：π/4 ≈ 1 - 1/3 + 1/5 - 1/7 = 16/355项：π/4 ≈ 1 - 1/3 + 1/5 - 1/7 + 1/9 = 128/315 可以看出，随着项数的增加，计算得到的π的近似值越来越接近π的真实值。
当项数无限增加时，级数的和将趋近于π/4，因此，我们可以通过计算足够多的项数来得到更精确的π的近似值。
这个方法虽然简单易行，但也存在一些缺点。
首先，计算出足够多的项数需要一定的时间和计算能力。
其次，级数的收敛速度相对较慢，因此需要计算大量的项数才能达到较高的精度。
总之，计算π的方法虽然多种多样，但采用简单的无穷级数法仍然是一种非常实用的方法。
在实际应用中，我们可以根据具体情况选择不同的计算方法来计算π，以达到更好的效果。
计算π的方法
计算π的方法
计算π的方法有多种，其中比较著名的方法有以下几种：
1. 随机抽样法：将一个正方形内切圆的比例定义为π/4，随机产生大量的点，统计落在圆内的点数和总点数，通过比例计算出π的近似值。
2. 蒙特卡罗法：通过随机模拟实验，计算出嵌入在一定区域内的曲线长度，从而得到π的值。
3. 差分法：通过求解圆的面积和正方形面积之差，推导出π的近似值。
4. 马青公式：利用级数展开式计算π的值，具体公式为
π/4=1-1/3+1/5-1/7+1/9-1/11+...
这些方法都可以大致计算出π的值，但精确程度可能有差异，具体取决于计算过程中的误差。
大学数学实验报告----怎样计算∏
验
过
程
及
结
果
1、数值积分法
n=1000
n=10000
n=1000
n=10000
2、泰勒级数法
练习1、
练习2、
3、蒙特卡罗法
n=1000
n=10000
n=50000
实验结果分析
1、数值积分法中s1使用梯形公式计算出的积分的近似值，s2使用辛普森公式计算出的积分的近似值，计算出的积分值与π的值十分接近；
（5）同时按Shift和Enter键运行。
（6）n=10000时，只需将上面两个程序中的第一行“n=1000”改为“n=10000”.即可。
2、泰勒积分法计算Π
练习1、将x=1代入反正切函数的泰勒级数得到 ，在此级数中取n=20000计算 的近似值。观察所得的结果和所花的时间。
(1)在计算机程序中打开Mathematica4.0系统；
取n=1000,10000,用梯形公式和辛普森公式计算 和 的近似值（取二十位有效数字）。将所得的结果与 的准确值相比较。
（1）在计算机程序中打开Mathematica4.0系统；
（2）点击鼠标进入工作区后，输入以下语句：
（3）同时按Shift和Enter键运行。
（4）点击鼠标进入工作区后，输入以下语句：
2、泰勒级数法可以较为精确的计算出π的值；
3、蒙特卡罗法可以近似的计算出π的值，但是精度很低，随着n的增大精度有所提高，但是仍比数值积分法与泰勒级数法计算出的π的值精度低；
以上三种方法中将n的值取得更大，算出的积分的值将更加接近π的值。
（3）同时按Shift和Enter键运行。
3、蒙特卡罗法计算Π
取n=1000，10000，50000，用蒙特卡罗法计算 的近似值。
派的实验报告
一、实验目的1. 了解派（π）的概念及其在数学中的重要性。
2. 掌握派（π）的计算方法，包括近似值和精确值。
3. 探索派（π）在几何、物理等领域的应用。
4. 通过实验，提高学生的数学思维能力和实践操作能力。
二、实验原理派（π）是数学中一个重要的无理数，表示圆的周长与直径的比值。
派（π）的值约等于3.14159，但它的精确值是无限不循环小数。
在数学、物理、工程等领域，派（π）都有广泛的应用。
三、实验内容1. 派（π）的计算方法2. 派（π）在几何、物理等领域的应用3. 派（π）的近似值与精确值四、实验步骤1. 派（π）的计算方法（1）近似值计算① 用圆的周长除以直径，得到派（π）的近似值。
② 用圆的面积除以半径的平方，得到派（π）的近似值。
（2）精确值计算① 利用公式π=4×（1-1/3+1/5-1/7+1/9-1/11+...）计算派（π）的精确值。
② 利用公式π=6×（1-1/5+1/7-1/9+1/11-1/13+...）计算派（π）的精确值。
2. 派（π）在几何、物理等领域的应用（1）几何领域① 圆的周长和面积计算。
② 椭圆、双曲线、抛物线等曲线的方程。
（2）物理领域① 圆周运动的速度、加速度计算。
② 力学中的转动惯量、扭矩等计算。
3. 派（π）的近似值与精确值（1）近似值派（π）的近似值有3.14、22/7、355/113等。
（2）精确值派（π）的精确值是无限不循环小数，可以用计算器或计算机软件得到。
五、实验结果与分析1. 派（π）的计算方法（1）近似值计算通过实验，我们得到了派（π）的近似值，如3.14、22/7、355/113等。
（2）精确值计算通过公式计算，我们得到了派（π）的精确值，如3.14159、3.1415926535等。
2. 派（π）在几何、物理等领域的应用在几何领域，我们利用派（π）计算了圆的周长和面积；在物理领域，我们利用派（π）计算了圆周运动的速度、加速度等。
π的计算公式的推导
π的计算公式的推导π（pi）是一个著名的数学常数，它表示圆周率，也就是将圆的周长除以圆的直径所得的值。
π的值大约是3.14159，但实际上π的值是无限小数，所以我们通常使用近似值来表示π。
现在我们来看一看π的计算公式的推导过程。
首先，我们可以用一个圆和圆周上的正多边形来求π的值。
在这种情况下，我们假设正多边形的边数为n，正多边形的边长为a，圆的周长为C，圆的直径为d。
根据圆的定义，我们可以得到如下公式：C=πd根据正多边形的定义，我们可以得到如下公式：C=na将这两个公式带入可得：πd=na将圆的直径d带入可得：π=n(a/d)根据正多边形的定义，我们可以得到如下公式：a/d=2sin(180/n)将这个公式带入可得：π=2n sin(180/n)当n趋近于无限大时，sin(180/n)趋近于1，所以π可以近似地表示为：π≈2n这就是π的计算公式的推导过程。
当然，π的计算公式并不止这一个，还有许多其他的计算公式。
比如有比如有一个经典的计算公式，叫做“Leibniz公式”，它可以用来求π的值。
这个公式的形式如下：π/4=1-1/3+1/5-1/7+1/9-1/11+…根据这个公式，我们可以用计算机来求π的值，只需要计算出后面的无限级数的值，再乘以4就可以得到π的值。
当然，这只是其中的一个计算公式，还有很多其他的计算公式，比如“Madhava公式”、“Euler 公式”等等。
这些公式都可以用来求π的值，但是它们的计算方法不尽相同。
总的来说，π是一个著名的数学常数，它的值表示圆的周长与圆的直径的比值。
π的值是无限小数，所以我们通常使用近似值来表示它。
π的计算公式有很多种，比如“Leibniz公式”、“Madhava公式”、“Euler公式”等等，这些公式都可以用来求π的值。
实验03 Π的计算
实验3 π的计算实验目的：1．掌握数学实验的方法和过程，学会撰写数学实验报告；2．掌握π的几种计算方法和思想，并能用其中的一些思想方法计算e ； 实验内容：1．描述刘徽割圆术计算π的原理、方法和计算步骤，并编写实现计算的函数式M 文件。
采取不同的分割计算π的近似值，并将计算的结果与较准确的π值进行比较，对算法进行分析。
2．编写采用级数展开式1114(1)21n i n π∞-==--∑ 和 1212111114(1)2123n n n i n π∞---=⎛⎫=-+ ⎪-⎝⎭∑ 计算π的函数式M 文件。
调用编写的程序，尝试n 取不同整数以求π的近似值，并将计算的结果与较准确的π值进行比较，对算法进行分析和比较。
3．描述Monte Carlo 方法计算π的原理、方法和计算步骤，并编写实现计算的函数式M 文件。
调用编写的程序计算π的近似值，并将计算的结果与较准确的π值进行比较，并对算法进行分析。
4．给出两种计算超越数e 的方法，并通过数值计算实验进行算法分析。
实验仪器与软件：1．CPU 主频在2GHz 以上，内存在512Mb 以上的PC ；2．Matlab 2010a 及以上版本。
实验讲评：实验成绩：评阅教师：20 年 月 日实验3 π的计算一、计算π的刘徽割圆术function capi(n)a(1)=1;for i=i:n-1a(i+1)=sqrt(2-sqrt(4-a(i)^2));endS=3*2^(n-1)*a(n)运行结果改进后function calpi1(n)a(1)=sym(1);%ÉèΪ·ûºÅ¶ÔÏófor i=1:n-1a(i+1)=sym(sqrt(2-sqrt(4-a(i)^2)));%ÉèΪ·ûºÅ¶ÔÏóendS=3*2^(n-1)*a(n);vpa(S,50)运行结果二、function calpi2(n) s=0;for i=1:nif mod(i,2)==0s=s-1/(2*i-1);elses=s+1/(2*i-1);endends=4*s运行结果改进后function calpi3(n)s=0;for i=1:nif mod(i,2)==0s=s-1/(2*i-1)*(1/(2^(2*i-1))+1/(3^(2*i-1))); elses=s+1/(2*i-1)*(1/(2^(2*i-1))+1/(3^(2*i-1))); endends=vpa(4*s,30)运行结果三、function calpi4(n)m=0;for n=1:nif rand(1)^2+rand(1)^2<=1 m=m+1;endends=4*m/n运行结果四、（1）使用泰勒级数的方法进行计算。
实验2_统计的方法求pi
#include <iostream> #include <iomanip> int main(intargc,char*argv[]) { // 声明所用的变量并正确初始化 // 编写求 pi 的循环 // 在你的循环体内 { //调用 rand（）函数，产生随机二维点（两个随机数：x 坐标和 y 坐标） // 将随机数变换到正方形所在的区间，利用系统宏 RAND_MAX // 计算点是否在圆内 // 更新统计数据，求 pi 值 } // 在控制台输出 pi 值 return0; }
3. 练习单步跟踪调试程序。 4. 注意变量的可取值范围。 提示：库函数 rand()可以产生在[0,RAND_MAX]范围内的随机整数，可在程序 中直接使用，其中 RAND_MAX 是宏，是 rand()函数可以产生的最大的整数。 通过强制类型转换和缩放可以把生成的整数变换到任意的范围。例如：
intv1 = rand() % 100; // v1 是 0 到 99 之间的整数 float v2 = float(rand())/float(RAND_MAX); // v2 是 0 到 1 之间的浮点数
如果在正方形区域内随机产生大量的均匀分布的点那么落入内切圆和正方形中的随机点个数的比值等于它们的面积之比
实验 2 统计的方法求π
实验目的
熟练运用循环和数值计算
实验内容
有一个正方形和它的内切圆。圆的半径为 R，则正方形的边长为 2R，内切圆 与正方形的面积比：
S圆 S 圆正方形
=
π R2 4R 2
=
π 4
如果在正方形区域内随机产生大量的均匀分布的点， 那么落入内切圆和正方 形中的随机点个数的比值等于它们的面积之比。该比值乘以 4，即为π 值。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
3、蒙特卡罗法可以近似的计算出π的值，但是精度很低，随着n的增大精度有所提高，但是仍比数值积分法与泰勒级数法计算出的π的值精度低；
以上三种方法中将n的值取得更大，算出的积分的值将更加接近π的值。
取n=1000,10000,用梯形公式和辛普森公式计算 和 的近似值（取二十位有效数字）。将所得的结果与 的准确值相比较。
（1）在计算机程序中打开Mathematica4.0系统；
（2）点ft和Enter键运行。
（4）点击鼠标进入工作区后，输入以下语句：
实
验
过
程
及
结
果
1、数值积分法
n=1000
n=10000
n=1000
n=10000
2、泰勒级数法
练习1、
练习2、
3、蒙特卡罗法
n=1000
n=10000
n=50000
实验结果分析
1、数值积分法中s1使用梯形公式计算出的积分的近似值，s2使用辛普森公式计算出的积分的近似值，计算出的积分值与π的值十分接近；
（5）同时按Shift和Enter键运行。
（6）n=10000时，只需将上面两个程序中的第一行“n=1000”改为“n=10000”.即可。
2、泰勒积分法计算Π
练习1、将x=1代入反正切函数的泰勒级数得到 ，在此级数中取n=20000计算 的近似值。观察所得的结果和所花的时间。
(1)在计算机程序中打开Mathematica4.0系统；
数学实验报告
实验一怎样计算∏
学院：数学与信息科学学院
班级：09级数学（4）班
姓名：***
学号：***
实验二：怎样计算∏
实验名称
利用计算机来计算圆周率∏的近似值。
实验目的
学习使用在Mathematica4.0中用数值积分法、泰勒级数法以及蒙特卡罗法计算Π的值,并观察用不同方法计算出来的Π有何差异。
实验环境
(2)点击鼠标进入工作区后，输入以下语句：
（3）同时按Shift和Enter键运行。
练习2、利用arctanx的泰勒展开示及公式 ， 计算 的近似值，精确到前150位有效数字，并将计算得到的结果与 的精确值得前150位有效数字相比较。
(1)在计算机程序中打开Mathematica4.0系统；
(2)点击鼠标进入工作区后，输入以下语句：
Mathematica4.0系统。
实
验
理
论
1、数值积分法是用求出单位圆的面积来求得Π的值；
2、泰勒级数法是利用反正切函数的泰勒级数 计算Π的值；
3、蒙特卡罗法是用随机数落在单位圆的1/4来计算Π的值。
实
验
内
容
及
步
骤
一、内容：
1、数值积分法计算Π
2、泰勒积分法计算Π
3、蒙特卡罗法计算Π
二、步骤：
1、数值积分法计算Π
（3）同时按Shift和Enter键运行。
3、蒙特卡罗法计算Π
取n=1000，10000，50000，用蒙特卡罗法计算 的近似值。
(1)在计算机程序中打开Mathematica4.0系统；
(2)点击鼠标进入工作区后，输入以下语句：
（3）同时按Shift和Enter键运行。
（4）n=10000，50000时，只需将上面程序中的第一行“n=1000”改为“n=10000”“n=50000”.即可。