如何计算π的值

1、蒙特卡罗(Monte Carlo )法

思想:

取一正方形A ,以A 的一个顶点为圆心,A 的边长为半径画圆,取四分之一圆(正方形内的四分之一圆)为扇形B 。已知A 的面积,只要求出B 的面积与A 的面积之比B A

S k S =,就能得出B S ,再由B 的面积为圆面积的四分之一,利用公式2=S R π圆即可求出π的值。因此,我

们的目的就是要找出k 的值。

可以把A 和B 看成是由无限多个点组成,而B 内的所有点都在A 内。随机产生n 个点,若落在B 内的有m 个点(假定A 的边长为1,以扇形圆心为坐标系原点。则只要使随机产生横纵坐标x 、y 满足

221x y +≤的点,就是落在B 内的点)

,则可近似得出k 的值,即m k n =,由此就可以求出π的值。

程序(1):

i=1;m=0;n=1000;

for i=1:n

a=rand(1,2);

if a(1)^2+a(2)^2<=1

m=m+1;

end

end

p=vpa(4*m/n,30)

matlab计算π,如何计算π的值(MATLAB)相关推荐

  1. matlab 计算指北角,用MATLAB如何计算出矩形区域内的最大值?

    在MATLAB图像处理工具箱(Image Processing Toolbox)的GUI函数中,有一些用来在图像或坐标轴上选择矩形或椭圆区域的函数,使用示例如下: 第一行代码:显示matlab内置的一 ...

  2. 用Matlab计算多项式的值

    本篇文章中,将会介绍Matlab计算一元多项式的基本原理和仿真过程.我们假设多项式P(x)有如下的形式:P(x)=anxn+an-1xn-1+-+a1x+a0 我们计算这个多项式的方法是synthet ...

  3. matlab 求复数模值,matlab计算带有复数的函数,最后求复数函数的模,结果里面却有...

    该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 matlab计算带有复数的函数,最后求复数函数的模,结果里面却有复数标志i,输入自变量,得到的结果确实是实数,不知道怎么回事; 代码如下: syms n2 ...

  4. matlab二元多项式求值,matlab多项式代入求值

    Matlab 多项式运算与方程求根 ? Matlab多项式运算无论是在线性代数中,还是信号处理.自动控制等理论 中,多项式运算都有着十分重要的地位,因此,MATLAB 为多项式的操作提供了相应的函数库 ...

  5. spwm matlab 计算波表,利用MATLAB计算SPWM脉冲宽度与并应用STM32输出

    先上一张spwm波形生成原理图: 首先利用MATLAB产生三角波与正弦波叠加: %% 产生正弦波与三角波叠加 y1=abs(sawtooth(a*2*pi*m,0.5));%三角波 y2=0.8*si ...

  6. matlab计算重心的函数,MATLAB计算出图像中多个物体的个数及重心.doc

    MATLAB计算出图像中多个物体的个数及重心 计算出原始图像中物体的个数及每个物体的中心,并显示出来 I=imread('keshe17.bmp'); //主函数 figure(1),imshow(I ...

  7. matlab计算膜孔径分布,应用Matlab绘制活性炭吸附等温线和孔径分布图

    第 31 卷 第 9 期 2012 年 9 月 实 验 室 研 究 与 探 索 RESEARCH AND EXPLORATION IN LABORATORY Vol. 31 No. 9 Sep. 20 ...

  8. Matlab计算多项式的值(数值)

    MATLAB 中,多项式用一个行向量表示,行向量的元素值为多项式系数按幂次的降序排列: 例如多项式, P(x) = 2*x^4 + 3*x^3 - 2*x^2 + 7*x + 11 可表示为, p = ...

  9. matlab计算的幅值裕度准确吗,单选(2分) 以下matlab函数中,可用于幅值裕度和相位裕度的计算的是( )。...

    单选(2分) 以下matlab函数中,可用于幅值裕度和相位裕度的计算的是( ). 更多相关问题 [多选题]预紧力的控制方法包括 [单选题]一般油藏的供油面积是圆形的,形状系数CA等于( ). MCGS ...

最新文章

  1. TCP/IP协议的编写《转载》
  2. Java虚拟机的研究与实现
  3. Hadoop书籍和网络资源介绍
  4. 23种设计模式C++源码与UML实现--适配器模式
  5. 二十不惑,三十而已|网易互联网人的“焦虑必杀技”
  6. 扬州大学广陵学院计算机网络试卷,扬州大学广陵学院控制工程考试样卷B.doc
  7. 避免CRM_IB 013 error message
  8. Android 基础性能数据获取(/proc/)
  9. php装箱,php兑现装箱算法
  10. Caffe-SSD相关源码说明和调试记录
  11. npumpy.insert(arr, obj, values[, axis])
  12. Maven学习总结(48)——利用wagon-maven-plugin插件自动化部署Jar包/war包到Linux远程服务器
  13. (转)关于做android+J2ee系统集成开发的一点心得
  14. 继电器设备改造需要的设备清单
  15. 基于微信小程序的疫情防控平台设计与实现-计算机毕业设计源码+LW文档
  16. local function definitions are illegal
  17. 100代码搞定C语言游戏开发,编程原来如此简单
  18. 【丁丁历险记】ESP8266/ESP32发送数据到微信
  19. led显示屏十大品牌
  20. 推特营销引流入门指南

热门文章

  1. C#时间格式转换问题(12小时制和24小时制)
  2. 我的世界服务器清道夫不显示,我的世界清道夫指令 | 手游网游页游攻略大全
  3. 入驻宣言:凡所见,皆可杀
  4. Unity3D插件 AnyPortrait 2D骨骼动画制作
  5. 小学教育专业有计算机课程吗,小学教育专科是属于师范专业吗
  6. 注册邮箱接收验证码小结
  7. 我的世界服务器无限漏斗,我的世界漏斗bug 我的世界怎么刷东西
  8. HTML 自定义实现emoji
  9. android安全漏洞(八)伪造短信(利用原生android4.0漏洞)
  10. 跨站点请求伪造(CSRF)