程序的循环结构和random库的使用及蒙特卡罗方法求圆周率
创始人
2025-06-01 22:10:33
0

 

第三个参数就是步长 

 

 引入文件时记得指明字符格式,否则读入不了

break结束循环,continue跳出本次循环 

  一个break只能跳出一个循环

一定要注意是是break是退出当前循环而不是条件判断 

我们为什么要用随机数种子呢?

就是因为你编程中给出了随机数种子,那么下次程序再次运行,只要种子相同,那么产生的随机也相同,我们对于随机数的程序我们可以复现或者再现你程序运行的过程,如果不使用随机数种子,那么它默认使用的是当前系统时间,而那个时间是精确到微秒的,而这个时间是我们很难再现的,所以后面产生的程序运行的结果就是完全随机。

因此如果我们再现随机过程,我们就用随机种子

 

 

 

 

 

蒙特卡罗方法是一种通过随机取样来估计数学问题的方法。我们可以用蒙特卡罗方法来估计圆周率。

假设我们有一个边长为 2 的正方形,以及一个内切于该正方形的圆形。圆的半径为 1。

现在我们从正方形中随机取样若干个点,并计算它们与圆心的距离。如果距离小于等于 1,则该点在圆内,否则在圆外。

我们可以用这些采样点的比例来估计圆的面积与正方形面积之比。这个比例等于圆周率除以 4。

因为圆的面积为 πr²,而正方形面积为 4r²,所以它们的比例为 π/4。

因此,我们可以通过采样点的比例乘以 4 来估计圆周率。

具体的做法是,随机生成大量的点,计算在圆内的点的数量,然后用这个数量除以总的采样点数,再乘以 4,就可以得到圆周率的近似值。

需要注意的是,采样点数量越多,估计的圆周率就越精确。

from random import random
from time import perf_counter
DARTS = 1000*1000
hits = 0.0
start = perf_counter()
for i in range(1, DARTS+1):x, y = random(), random()dist = pow(x ** 2 + y ** 2, 0.5)if dist <= 1.0:hits = hits + 1pi = 4 * (hits/DARTS)
print("圆周率值是: {}".format(pi))
print("运行时间是: {:.5f}s".format(perf_counter()-start))

 

这段代码使用了蒙特卡罗方法来估计圆周率的值。下面是对代码的解剖:

导入 random 模块和 perf_counter 函数。

将常量 DARTS 设为 1000*1000,表示采样点数量。

定义变量 hits,初始值为 0.0,用于记录落在圆内的点的数量。

使用 perf_counter 函数记录开始时间。

使用 for 循环生成 DARTS 个采样点,每个点的横坐标和纵坐标都是一个 [0, 1) 之间的随机数。

计算采样点距离圆心的距离,如果距离小于等于 1,则将 hits 加 1。

根据蒙特卡罗方法,用 hits/DARTS 计算落在圆内的点的比例,再乘以 4,得到圆周率的近似值。

输出圆周率的近似值和运行时间。

需要注意的是,蒙特卡罗方法是一种概率方法,估计的圆周率并不是精确的值,而是一个近似值。因此,采样点数量越多,估计的值越接近真实值。在这段代码中,采样点数量为 1000*1000,可以得到较为准确的近似值。

 

相关内容

热门资讯

linux入门---制作进度条 了解缓冲区 我们首先来看看下面的操作: 我们首先创建了一个文件并在这个文件里面添加了...
C++ 机房预约系统(六):学... 8、 学生模块 8.1 学生子菜单、登录和注销 实现步骤: 在Student.cpp的...
JAVA多线程知识整理 Java多线程基础 线程的创建和启动 继承Thread类来创建并启动 自定义Thread类的子类&#...
【洛谷 P1090】[NOIP... [NOIP2004 提高组] 合并果子 / [USACO06NOV] Fence Repair G ...
国民技术LPUART介绍 低功耗通用异步接收器(LPUART) 简介 低功耗通用异步收发器...
城乡供水一体化平台-助力乡村振... 城乡供水一体化管理系统建设方案 城乡供水一体化管理系统是运用云计算、大数据等信息化手段࿰...
程序的循环结构和random库...   第三个参数就是步长     引入文件时记得指明字符格式,否则读入不了 ...
中国版ChatGPT在哪些方面... 目录 一、中国巨大的市场需求 二、中国企业加速创新 三、中国的人工智能发展 四、企业愿景的推进 五、...
报名开启 | 共赴一场 Flu... 2023 年 1 月 25 日,Flutter Forward 大会在肯尼亚首都内罗毕...
汇编00-MASM 和 Vis... Qt源码解析 索引 汇编逆向--- MASM 和 Visual Studio入门 前提知识ÿ...
【简陋Web应用3】实现人脸比... 文章目录🍉 前情提要🌷 效果演示🥝 实现过程1. u...
前缀和与对数器与二分法 1. 前缀和 假设有一个数组,我们想大量频繁的去访问L到R这个区间的和,...
windows安装JDK步骤 一、 下载JDK安装包 下载地址:https://www.oracle.com/jav...
分治法实现合并排序(归并排序)... 🎊【数据结构与算法】专题正在持续更新中,各种数据结构的创建原理与运用✨...
在linux上安装配置node... 目录前言1,关于nodejs2,配置环境变量3,总结 前言...
Linux学习之端口、网络协议... 端口:设备与外界通讯交流的出口 网络协议:   网络协议是指计算机通信网...
Linux内核进程管理并发同步... 并发同步并发 是指在某一时间段内能够处理多个任务的能力,而 并行 是指同一时间能够处理...
opencv学习-HOG LO... 目录1. HOG(Histogram of Oriented Gradients,方向梯度直方图)1...
EEG微状态的功能意义 导读大脑的瞬时全局功能状态反映在其电场结构上。聚类分析方法一致地提取了四种头表面脑电场结构ÿ...
【Unity 手写PBR】Bu... 写在前面 前期积累: GAMES101作业7提高-实现微表面模型你需要了解的知识 【技...