蒙特卡洛方法近似计算圆周率

2024-06-05 20:44

本文主要是介绍蒙特卡洛方法近似计算圆周率,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

一、原理图

正方形区域内有1/4圆区域,向区域内随机大量掷点。根据概率论知识,落在每一点的概率相等,当n无限大时,落在1/4区域内的点占总点数量的比例即1/4圆占正方形面积的比例,圆周率PI=4(N1/N)[N1:落在1/4圆中点的数量,N:总点数]。

二、c++试验

#include <iostream>
#include <random>
#include <cmath>
#include <iomanip>  // 为了设置输出精度using namespace std;
typedef long long ll;int main()
{int N = 10; // 点阵的规模N*Nint M = 100;// 掷点数量for (int N = 10; N <= 10000; N *= 10){for (int M = N*N; M <= 100000000; M *= 10){int cnt = 3;while (cnt--){std::random_device rd;  // 真随机数生成器std::mt19937 gen(rd()); // 使用Mersenne Twister算法生成高质量随机数std::uniform_real_distribution<> dis(0.0, N); // 生成0到N之间的均匀分布的浮点数long double PI; // 近似圆周率ll N1 = 0; // 落在半圆内的点数量for (ll i = 0; i < M; i++){long double x = dis(gen);long double y = dis(gen);long double d = sqrt(x * x + y * y);if (d <= N)  N1++;}PI = (N1 / (long double)M) * 4;cout<<"N="<<N<<",M="<<M << setprecision(15) << "时近似计算圆周率π = " << PI << endl; // 设置输出精度}}}
}

改变M、N的值,观察输出圆周率的精度(M最大取了100000000,再大循环太折磨了...)

结果如下:

N=10,M=100时近似计算圆周率π = 3.2
N=10,M=100时近似计算圆周率π = 3.04
N=10,M=100时近似计算圆周率π = 3.28
N=10,M=1000时近似计算圆周率π = 3.076
N=10,M=1000时近似计算圆周率π = 3.264
N=10,M=1000时近似计算圆周率π = 3.12
N=10,M=10000时近似计算圆周率π = 3.1464
N=10,M=10000时近似计算圆周率π = 3.158
N=10,M=10000时近似计算圆周率π = 3.1344
N=10,M=100000时近似计算圆周率π = 3.13672
N=10,M=100000时近似计算圆周率π = 3.13812
N=10,M=100000时近似计算圆周率π = 3.13512
N=10,M=1000000时近似计算圆周率π = 3.14016
N=10,M=1000000时近似计算圆周率π = 3.142876
N=10,M=1000000时近似计算圆周率π = 3.143608
N=10,M=10000000时近似计算圆周率π = 3.1416004
N=10,M=10000000时近似计算圆周率π = 3.1419076
N=10,M=10000000时近似计算圆周率π = 3.1425656
N=10,M=100000000时近似计算圆周率π = 3.14145444
N=10,M=100000000时近似计算圆周率π = 3.14154308
N=10,M=100000000时近似计算圆周率π = 3.14189244
N=100,M=10000时近似计算圆周率π = 3.1204
N=100,M=10000时近似计算圆周率π = 3.1444
N=100,M=10000时近似计算圆周率π = 3.1388
N=100,M=100000时近似计算圆周率π = 3.14372
N=100,M=100000时近似计算圆周率π = 3.14164
N=100,M=100000时近似计算圆周率π = 3.14344
N=100,M=1000000时近似计算圆周率π = 3.13904
N=100,M=1000000时近似计算圆周率π = 3.141212
N=100,M=1000000时近似计算圆周率π = 3.14156
N=100,M=10000000时近似计算圆周率π = 3.1418392
N=100,M=10000000时近似计算圆周率π = 3.1423576
N=100,M=10000000时近似计算圆周率π = 3.1415224
N=100,M=100000000时近似计算圆周率π = 3.14161424
N=100,M=100000000时近似计算圆周率π = 3.14168548
N=100,M=100000000时近似计算圆周率π = 3.14168172
N=1000,M=1000000时近似计算圆周率π = 3.142296
N=1000,M=1000000时近似计算圆周率π = 3.140968
N=1000,M=1000000时近似计算圆周率π = 3.14106
N=1000,M=10000000时近似计算圆周率π = 3.142186
N=1000,M=10000000时近似计算圆周率π = 3.1408512
N=1000,M=10000000时近似计算圆周率π = 3.1420132
N=1000,M=100000000时近似计算圆周率π = 3.14153392
N=1000,M=100000000时近似计算圆周率π = 3.14163044
N=1000,M=100000000时近似计算圆周率π = 3.14182124
N=10000,M=100000000时近似计算圆周率π = 3.14151084
N=10000,M=100000000时近似计算圆周率π = 3.14142316
N=10000,M=100000000时近似计算圆周率π = 3.14153156

从你的实验结果可以看出,随着 M 的增加,计算出的圆周率 π 的精度逐渐提高,结果也更加接近实际值 (  3.14159265 )。同时,增加 N 也可以在一定程度上提高精度。

这篇关于蒙特卡洛方法近似计算圆周率的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/1034129

相关文章

【C++】_list常用方法解析及模拟实现

相信自己的力量,只要对自己始终保持信心,尽自己最大努力去完成任何事,就算事情最终结果是失败了,努力了也不留遗憾。💓💓💓 目录   ✨说在前面 🍋知识点一:什么是list? •🌰1.list的定义 •🌰2.list的基本特性 •🌰3.常用接口介绍 🍋知识点二:list常用接口 •🌰1.默认成员函数 🔥构造函数(⭐) 🔥析构函数 •🌰2.list对象

浅谈主机加固,六种有效的主机加固方法

在数字化时代,数据的价值不言而喻,但随之而来的安全威胁也日益严峻。从勒索病毒到内部泄露,企业的数据安全面临着前所未有的挑战。为了应对这些挑战,一种全新的主机加固解决方案应运而生。 MCK主机加固解决方案,采用先进的安全容器中间件技术,构建起一套内核级的纵深立体防护体系。这一体系突破了传统安全防护的局限,即使在管理员权限被恶意利用的情况下,也能确保服务器的安全稳定运行。 普适主机加固措施:

webm怎么转换成mp4?这几种方法超多人在用!

webm怎么转换成mp4?WebM作为一种新兴的视频编码格式,近年来逐渐进入大众视野,其背后承载着诸多优势,但同时也伴随着不容忽视的局限性,首要挑战在于其兼容性边界,尽管WebM已广泛适应于众多网站与软件平台,但在特定应用环境或老旧设备上,其兼容难题依旧凸显,为用户体验带来不便,再者,WebM格式的非普适性也体现在编辑流程上,由于它并非行业内的通用标准,编辑过程中可能会遭遇格式不兼容的障碍,导致操

透彻!驯服大型语言模型(LLMs)的五种方法,及具体方法选择思路

引言 随着时间的发展,大型语言模型不再停留在演示阶段而是逐步面向生产系统的应用,随着人们期望的不断增加,目标也发生了巨大的变化。在短短的几个月的时间里,人们对大模型的认识已经从对其zero-shot能力感到惊讶,转变为考虑改进模型质量、提高模型可用性。 「大语言模型(LLMs)其实就是利用高容量的模型架构(例如Transformer)对海量的、多种多样的数据分布进行建模得到,它包含了大量的先验

【北交大信息所AI-Max2】使用方法

BJTU信息所集群AI_MAX2使用方法 使用的前提是预约到相应的算力卡,拥有登录权限的账号密码,一般为导师组共用一个。 有浏览器、ssh工具就可以。 1.新建集群Terminal 浏览器登陆10.126.62.75 (如果是1集群把75改成66) 交互式开发 执行器选Terminal 密码随便设一个(需记住) 工作空间:私有数据、全部文件 加速器选GeForce_RTX_2080_Ti

【VUE】跨域问题的概念,以及解决方法。

目录 1.跨域概念 2.解决方法 2.1 配置网络请求代理 2.2 使用@CrossOrigin 注解 2.3 通过配置文件实现跨域 2.4 添加 CorsWebFilter 来解决跨域问题 1.跨域概念 跨域问题是由于浏览器实施了同源策略,该策略要求请求的域名、协议和端口必须与提供资源的服务相同。如果不相同,则需要服务器显式地允许这种跨域请求。一般在springbo

AI(文生语音)-TTS 技术线路探索学习:从拼接式参数化方法到Tacotron端到端输出

AI(文生语音)-TTS 技术线路探索学习:从拼接式参数化方法到Tacotron端到端输出 在数字化时代,文本到语音(Text-to-Speech, TTS)技术已成为人机交互的关键桥梁,无论是为视障人士提供辅助阅读,还是为智能助手注入声音的灵魂,TTS 技术都扮演着至关重要的角色。从最初的拼接式方法到参数化技术,再到现今的深度学习解决方案,TTS 技术经历了一段长足的进步。这篇文章将带您穿越时

模版方法模式template method

学习笔记,原文链接 https://refactoringguru.cn/design-patterns/template-method 超类中定义了一个算法的框架, 允许子类在不修改结构的情况下重写算法的特定步骤。 上层接口有默认实现的方法和子类需要自己实现的方法

使用JS/Jquery获得父窗口的几个方法(笔记)

<pre name="code" class="javascript">取父窗口的元素方法:$(selector, window.parent.document);那么你取父窗口的父窗口的元素就可以用:$(selector, window.parent.parent.document);如题: $(selector, window.top.document);//获得顶级窗口里面的元素 $(

消除安卓SDK更新时的“https://dl-ssl.google.com refused”异常的方法

消除安卓SDK更新时的“https://dl-ssl.google.com refused”异常的方法   消除安卓SDK更新时的“https://dl-ssl.google.com refused”异常的方法 [转载]原地址:http://blog.csdn.net/x605940745/article/details/17911115 消除SDK更新时的“