西安交大 数学实验五 matlab计算圆周率pi-课件·PPT
用计算机计算圆周率课件
这里是一个示例课件的简要概述，用于介绍如何使用计算机来计算圆周率：
1. 引言：介绍圆周率的概念和重要性，以及计算圆周率的方法。
2. 数学公式：讲解数学公式中圆周率的含义和计算方法，如利用圆的周长和直径的比值。
3. 近似方法：介绍一些近似计算圆周率的方法，如阿基米德法、蒙特卡洛方法和无穷级数法。
4. 计算机的角色：讨论计算机在计算圆周率中的作用，包括高精度计算、并行计算和分布式计算。
5. 编程语言和软件工具：介绍一些常用的编程语言和软件工具，如Python、MATLAB和Mathematica，用于编写计算圆周率的程序。
6. 编程实例：提供一个简单的编程实例，演示如何使用编程语言计算圆周率的近似值。
7. 程序优化：讨论如何优化计算圆周率的程序，提高计算速度和精度，如使用多线程、并行计算和算法优化。
8. 应用领域：讨论计算圆周率在实际应用中的重要性，如在科学计算、密码学和通信领域的应用。
9. 结论和总结：总结课程内容，强调计算圆周率的重要性和应用前景。
这只是一个简要的概述，您可以根据具体需要和时间安排，进一步展开和深入讲解各个部分。
在课件中可以添加图表、示意图和实例代码等辅助材料，以帮助学生更好地理解和应用计算圆周率的方法。
matlab圆周率
Matlab圆周率1. 引言圆周率（π）是数学中的一个重要常数，表示圆的周长与直径之间的比值。
圆周率的精确值是一个无限不循环的小数，被广泛应用于数学、物理和工程等领域。
在本文中，我们将使用Matlab来计算圆周率，并探讨一些与圆周率相关的数值计算方法和应用。
2. 圆周率的定义与计算方法2.1 定义圆周率可以通过圆的周长与直径的比值来定义，即π = C/d，其中C表示圆的周长，d表示圆的直径。
2.2 常用计算方法2.2.1 几何法几何法是最直观的计算圆周率的方法之一。
它基于圆的定义，将一个正方形内切于一个圆，并通过计算正方形的边长与圆的直径的比值来逼近圆周率。
几何法的思想简单易懂，但精度较低，随着划分的正方形数量增加，计算结果会逼近真实值。
2.2.2 随机法随机法是一种通过生成随机点，并统计这些点落在圆内的比例来估计圆周率的方法。
随机法利用概率统计的思想，通过大量的随机模拟来逼近圆周率。
使用Matlab可以方便地生成随机数，并进行统计计算，从而得到较为准确的圆周率估计值。
2.2.3 数列法数列法是通过一些特定公式或数列来逼近圆周率的方法。
其中最著名的数列是莱布尼兹级数和无穷乘积法。
莱布尼兹级数是一个交替级数，通过逐项相加可以逼近圆周率的值。
无穷乘积法则是使用无穷乘积的形式来逼近圆周率，其中最著名的是狄利克雷（Dirichlet）级数。
这些数列法适用于在计算机中进行快速近似计算，但对于精确值的计算效果较差。
3. Matlab中的圆周率计算函数Matlab提供了多种计算圆周率的函数，方便了用户进行圆周率的计算和应用。
下面介绍几个常用的Matlab函数。
### 3.1 pi pi是Matlab内置的常量，表示圆周率的值。
使用pi函数可以直接获得圆周率的近似值，具有较高的精度。
3.2 montecarlo_pimontecarlo_pi是一个基于随机法的圆周率计算函数。
它通过生成随机点，并统计这些点落在圆内的比例来估计圆周率。
西安交大 数学实验五 matlab计算圆周率pi
并在1994年计算到了4044000000位．它的另一 种形式是
426880 10005 . (6n)!(545140134 13591409) n (n! )3 (3n)!(640320)3n n 0
1995 年 ， 由 David Bailey,Peter Borwein 和 Simon Plouffe 共同发表了下面的圆周率计算公式 （简称BBP公式）
x5
b x
现对每个子区间再二等分，得到2n个子区间
b a
n 1 n x i 1 x i h f ( x )dx T2 n ( f (a ) 2 f ( x i ) 2 f ( ) f (b)) , 4 2 i 1 i 1 ba h , xk a kh k 0,1,2, , n n n n 1 h x i 1 x i 1 (Tn h f ( )) (Tn h f (a ih )) 2 2 2 2 i 1 i 1
2.圆周率的数值积分计算方法 1 1 1 1 0 1 x 2 dx 4 40 1 x 2 dx
表1给出的是不同等分区间 数计算出的 圆周率的近似值.
等分区间数n 10 20 50 100 500 1000 5000
圆周率的近似值
3.14242598500110 3.14180098689309 3.14162598692300 3.14160098692312 3.14159298692312 3.14159273692313 3.14159265692313
1
0
e
x2
dx
抛物线方法 y=inline('exp(-x.^2)'); quad(y,0,1) ans = 0.746826
MATLAB经典教程(全)PPT课件
MATLAB的优势
易于学习、使用灵活、高效的数值计 算和可视化功能、强大的工具箱支持。
发展历程
从最初的数值计算工具，逐渐发展成 为一款功能强大的科学计算软件，广 泛应用于工程、科学、经济等领域。
MATLAB工作环境与界面
MATLAB工作环境
包括命令窗口、工作空间、命令历史窗口、当 前文件夹窗口等。
界面介绍
详细讲解MATLAB界面的各个组成部分，如菜 单栏、工具栏、编辑器窗口等。
基本操作
介绍如何在MATLAB环境中创建、保存、运行脚本和函数，以及如何进行基本 的文件操作。
基本数据类型与运算
矩阵大小
使用`size`函数获取矩阵的行数 和列数。
矩阵元素访问
通过下标访问矩阵元素，如 `A(i,j)`表示访问矩阵A的第i行第j 列元素。
矩阵基本操作
包括矩阵的加、减、数乘、转置 等操作。
矩阵运算及性质
矩阵乘法 满足乘法交换律和结合律，但不满足 乘法交换律。
矩阵的逆
对于方阵，若存在一矩阵B，使得 AB=BA=I（I为单位矩阵），则称B 为A的逆矩阵。
Hale Waihona Puke 03 数据分析与可视化数据导入、导出及预处理
数据导入
介绍如何使用MATLAB导入各种格式的数据文件， 如.csv、.txt、.xlsx等。
数据导出
讲解如何将MATLAB中的数据导出为常见的数据文件格式，以 便于数据共享和交换。
数据预处理
阐述数据清洗、数据变换、数据规约等预处理技术，为后续的数 据分析和可视化奠定基础。
01
02
matlab 圆周率
本实验涉及的方法
3、蒙特卡罗法——蒲丰(Buffon)掷针法
x表示针的中点与最近的一条平行线间的距离, 为针与平行线 的夹角.则0 x d 2 ,0 , 确定x 平面内一个矩形 .如图：
平面内一个矩形确定的夹角为针与平行线一条平行线间的距离表示针的中点与最近的针与平行线相交本实验涉及的方法3蒙特卡罗法蒲丰buffon掷针法针与平行线相交时针在平行线矩形区域内所组成的区域面积与矩形区域面积比即为针与直线相交的概率
实验三
圆周率的求法
问
题
• 设半径为1的圆的周长与圆的直径比为A; • 设半径为2的圆的周长与圆的直径比为B. • 请问A与B哪个大?
计算历程的几个阶段
3.分析方法时期 17世纪出现了数学分析，这锐利的工具 使得许多初等数学束手无策的问题迎刃 而解。 π 的计算历史也随之进入了一 个新的阶段。
计算历程的几个阶段
1593年，韦达给出
1671年，苏格兰数学家詹姆斯.格雷戈里 公开了他发现的公式： 3 5 7
arctan x x x3 x5 x7 (1 x 1)
故
2 nl md
本实验涉及的方法
3、蒙特卡罗法——蒲丰(Buffon)掷针法
• 随机生成一个针的中点与最近一条平行直线间的 距离x,满足0≤x≤d/2; • 随机生成针与平行线的夹角a,满足0≤a≤ ； • 判断x≤L/2*sin[a],若成立，则针与直线相交， 否则不相交。
本实验涉及的方法
计算历程的几个阶段
圆周长大于内接正四边形而小于外切正四边形
计算历程的几个阶段
在我国，首先是由数学家刘徽得出较精确的 圆周率。他采用割圆术。
西安交大数学实验五matlab计算圆周率pi-精选文档
3
5
7
2 n 1
并利用这个公式计算到了圆周率的100位.
1914年，印度数学家Srinivasa Ramanujan 表了下面的公式：
9801 ( 4 n )! ( 1103 26390 n ) 22 4 n 4 4 n 4 ( n ! ) 99 n 0
发
在1985年，Gosper用这个公式计算到了圆周率 的17500000位．
在中国
• 祖冲之: 在刘徽研究的基础上，进一步地发展， 经过既漫长又烦琐的计算，一直算到圆内接正 24576边形，而得到一个结论： • 3.1415926 ＜ π ＜ 3.1415927 同时得到π 的两个近似分数：约率为22／7； 密率为355／113。
• 他算出的 π 的8位可靠数字，不但在当时是最精 密的圆周率，而且保持世界记录九百多年。以致 于有数学史家提议将这一结果命名为“祖率”。
1 .圆周率 的幂级数计算方法
例1
利用 arctan x 的 Maclaurin 展开式 ， 计算 的近似值 .
解
1 2 4 n 12 n 1 x x (1 ) x . 2 1 x
xx ( 1 )x arctan x x . 35 2 n 1
1989年，David 和 Gregory Chudnovsky 发表 了下面的公式
n 1 ( 1 ) ( 6 n )! 13591409 5451401 n 12 , 3 3 3 n n ( 3 n )! ( n ! ) 2 0 640320
并在1994年计算到了4044000000位．它的另一 种形式是
在中国
• 刘徽：公元263年前后，刘徽提出著名的 “割圆 术”求出了比较精确的圆周率。他发现：当圆内 接正多边形的边数不断增加后，多边形的周长会 越来越逼近圆周长，而多边形的面积也会越来越 逼近圆面积。于是，刘徽利用正多边形面积和圆 面积之间的关系，从正六边形开始，逐步把边数 加倍：正十二边形、正二十四边形，正四十八边 形……，一直到正三○七二边形，算出圆周率等 于三点一四一六，将圆周率的精度提高到小数点 后第四位。
MATLAB求圆周率
编写M 文件,用公式 +-+-≈7
1513114π
求π的近似值,计算到第n 项,n 由用户输入 【代码】：n=input('n=');
sum=0;
for i=1:n;
sum=sum+(-1)^(i+1)/(2*i-1);
end
※ input 有两种用法（1）. 在命令窗口中输入Val=input('请输入一个整数')，这样在命
令窗口中便会显示“请输入一个整数”提示用户进行输入操作，当用户输入一个整数后，便会被赋给Val 。
（2）. 当你在命令窗口输入str=input('Please input', 's')然后
从键盘输入：1+2+3，这样str 实际得到的是：'1+2+3'而不是6。
值得注意的是：需要注意的是，如果执行本函数时，用户敲了回车而不是输入了一个数，则该函数返回一个空矩阵。
可以用MATLAB 中的isempty 函数判断输入的是否为空。
※ for i=a:b; end MATLAB 中的循环语句，i=a 开始第一次计算从for 的下一行计算到end
截止，之后开始计算a+1，a+2…….b
※ a=a+b 在MATLAB 中表示将a+b 的值赋予a 通常在循环语句中出现
※ a=b 表示将b 的值赋给a
※ a= =b 表示判断a b 的值是否相等。
数学实验matlab ppt课件
MATLAB已经发展成为多学科、多种工作平台 的功能强大的大型软件;
MATLAB已经成为线性代数、自动控制理论、 概率论及数理统计、数字信号处理、时间序列 分析、动态系统仿真等高级课程的基本教学工 具，是攻读学位的大学生、硕士生、博士生必 须掌握的基本技能.
数学实验matlab
❖语言简洁紧凑，使用方便灵活。MATLAB的基本数据单元是既不需要指
主窗口除了嵌入一些子窗口外，还主要包括 菜单栏和工具栏。
1．菜单栏
在MATLAB 6.5主窗口的菜单栏，共包含File、 Edit、View、Web、Window和Help 6个菜单项。
(1) File菜单项：File菜单项实现有关文件的操作。
(2) Edit菜单项：Edit菜单项用于命令窗口的编辑操作。
(3) View菜单项：View菜单项用于设置MATLAB集成环 境的显示方式。 (4) Web菜单项：Web菜单项用于设置MATLAB的Web 操作。 (5) Window菜单项：主窗口菜单栏上的Window菜单， 只包含一个子菜单Close all，用于关闭所有打开的编辑 器窗口，包括M-file、Figure、Model和GUI窗口。
数学实验matlab
Matlab简介与入门
数学实验matlab
什么是数学实验 ？
简单讲就是利用计算机和数学软件平台. 一方面，对学习知识过程中的某些问题 进行实验探究、发现规律； 另一方面，结合已掌握的数学（微积分、 代数与几何等）知识，去探究、解决一些简单 实际问题，从而熟悉从数学建模、解法研究到 实验分析的科学研究的方法。
在MATLAB里，有很多的控制键和方向键可用于命令行 的编辑。
工作空间窗口
工作空间是MATLAB用于存储各种变量和结 果的内存空间。在该窗口中显示工作空间中所有 变量的名称、大小、字节数和变量类型说明，可 对变量进行观察、编辑、保存和删除。