随机模拟算法求解圆周率

圆周率(π)这个东西是从小学开始一直陪伴我们的,这里使用使用蒙特卡洛算法来产生大量的随机数求解π的近似值。

计算方式

首先我们知道 正方形的面积公式是S1 = a * a,圆形的面积S2 = π * r * r;
所以以圆的直径为正方形边长,可以得出π的表达式。

π = 4 * S2 / S1 

这样一来,重点就是求解正方形和圆形的面积,这里使用在一正方形区域内圆内产生相应的随机点,
为了便于可视化分析,在圆内和圆外的点分别用不同的颜色来表示,最后圆内产生的点数就近似记作圆的面积,区域内的点数记作正方形的面积。

可视化分析-小数据集

这里用Java swing来求解。

首先看看Circle的model类

public class Circle {
   

    private int x, y, r;

    public Circle(int x, int y, int r){
        this.x = x;
        this.y = y;
        this.r = r;
    }

    public int getX(){ return x; }
    public int getY(){ return y; }
    public int getR(){ return r; }
  // 判断点是否包含在圆内
    public boolean contain(Point p){
        return Math.pow(p.x - x, 2) + Math.pow(p.y - y, 2) <= r*r;
    }
}

这里定义了圆的坐标,以及圆的半径。
然后看看点的model类

public class MonteCarloPiData {
   

    private Circle circle;
    private LinkedList<Point> points;
    private int insideCircle = 0;

    public MonteCarloPiData(Circle circle){
        this.circle = circle;
        points = new LinkedList<Point>();
    }

    public Circle getCircle(){
        return circle;
    }
    // 得到点的数量
    public int getPointsNumber(){
        return points.size();
    }

    
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
编写思路: 求解圆周率的问题,从数学角度来看是一个无限极限的求解问题,这是一个非常典型的数值计算问题。而在计算机编程中,可以通过数值积分、蒙特卡罗方法等多种算法来进行求解。 我们采用的是蒙特卡罗方法。蒙特卡罗方法是一种基于随机采样的数值计算方法,通过对随机事件的大量模拟,来近似计算数学问题的解。 在本问题中,我们可以利用蒙特卡罗方法来对圆周率进行求解。具体来说,我们可以通过在一个正方形内随机生成大量的点,然后根据这些点落在正方形内部和圆内部的数量比例,来计算圆周率的近似值。 设正方形边长为2,圆的直径为2,半径为1。在正方形内随机生成n个点,落在圆内的点的数量为m。根据圆的面积公式和正方形的面积公式,可以得到圆周率的近似值为4m/n。 程序特点: 1.简单易懂:该程序采用了蒙特卡罗方法,是一种非常基础的数值计算方法。代码简单易懂,适合初学者学习。 2.高效稳定:该程序使用了Python内置的random库生成随机数,计算速度较快,而且结果稳定。 3.可扩展性强:该程序的计算方法可以很容易地扩展到其他形状的求解问题上,如球体积、曲线长度等。 程序代码: import random def calc_pi(n): count = 0 for i in range(n): x = random.uniform(-1, 1) y = random.uniform(-1, 1) if x ** 2 + y ** 2 <= 1: count += 1 return 4 * count / n if __name__ == '__main__': print(calc_pi(1000000))

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值