c语言如何精确求得pi值,[updated]C语言和python实现蒙特卡里算法求pi值

逃离我推掉我的手 2023-10-04 14:20 11阅读 0赞

高中数学课本上曾经写过用计算机算得pi的近似值，用的是蒙特卡里方法。但是那时的我还不知道编程为何物，今天突然又想到这个问题，尝试了一下，发现并不是那么难，用C语言和python各实现了一遍。

不多说上代码。

c语言实现

\#include

\#include

\#include

\#include

\#define MAXN 1200000//总随机点数，这个数值越大，pi值计算的越精确，相应的耗费更多的时间

int main(void)

\{

int i=0;

double count=0.0;

double x=0.0,y=0.0,dist=0.0;

double my\_pi=0;

srand((unsigned)time(NULL));

//随机数种子被初始化

clock\_t start = clock();//时间起点

for(i=0;i

\{

x = ((double)rand()/RAND\_MAX);//生成0-1的随机数

y = ((double)rand()/RAND\_MAX);

dist = sqrt(x\*x + y\*y);//计算(x，y)到原点的距离

if(dist <= 1.0)

count++;//记录落在半圆内的随机点的个数

\}

my\_pi = 4\*(count/MAXN);

clock\_t finish = clock();//函数运行结束时间记录

printf("%lf\\n",my\_pi);

double duration = (double)(finish - start) / CLOCKS\_PER\_SEC;//计算出函数所耗费的时间

printf("It costs %lf s",duration);

return 0;

\}

python版 python版本为3.5X

\# -\*- coding:utf-8 -\*-

from random import random

from math import sqrt

from time import clock

MAXN = 120000

count =0

clock()

for i in range (1,MAXN):

x, y = random(),random()

dist = sqrt(x\*\*2+y\*\*2)

if dist <= 1.0:

count = count + 1

pi = 4\*(count/MAXN)

print("PI的值为") % pi

print("程序的运行时间是%-5.5ss") % clock()

![074a880891d6][]

pi.PNG

看的出来计算出的pi值离精确的pi值还有一定距离，并且数值会波动，毕竟这是用概率去模拟的，样本点足够多，才可能精确不是？

第一次使用Markdown写东西，将就看吧o(╯□╰)o

┴┬┴┬／￣＼＿／￣＼

┬┴┬┴▏　　▏▔▔▔▔＼

┴┬┴／＼　／　　　　　　﹨

┬┴∕　　　　　　　／　　　)

┴┬▏　　　　　　　　●　　▏

┬┴▏　　　　　　　　　　　▔█

┴◢██◣　　　　　＼＿＿＿／

┬█████◣　　　　　　　／

┴█████████████◣

◢██████████████▆▄

█◤◢██◣◥█████████◤＼

◥◢████　████████◤　　 ＼

┴█████　██████◤　　　　　 ﹨

┬│　　　│█████◤　　　　　　　　▏

┴│　　　│　　　　　　　　　　　　　　▏

┬ ∕　　　 ∕　　　　／▔▔▔＼　　　　 ∕

┴/＿＿＿／﹨　　　∕　　　　　﹨　　／＼

┬┴┬┴┬┴＼ 　　 ＼ 　　　　　﹨／　　 ﹨

┴┬┴┬┴┬┴ ＼＿＿＿＼　　　　 ﹨／▔＼﹨ ▔＼

▲△▲▲╓╥╥╥╥╥╥╥╥＼　　 ∕　 ／▔﹨／▔﹨

＊＊╠╬╬╬╬╬╬╬╬＊﹨　　／　　／／

[074a880891d6]: https://www.jianshu.com/p/074a880891d6

发表评论取消回复

表情：

评论列表（有 0 条评论，11人围观）

还没有评论，来说两句吧...

相关阅读

相关 Raspberry Pi Pico C/C++语言开发环境搭建 Raspberry Pi Pico C/C++ SDK

以下在Windows 10系统下测试 Raspberry Pi Pico介绍关于Pico板的介绍请参考Raspberry官网：https://www.raspberry

刺骨的言语ヽ痛彻心扉/ 2024年03月27日 17:08/ 0 赞/ 120 阅读

相关 C语言表达式求值

表达式求值一、整型提升 1.什么是整型提升 2.整型提升的意义 3.如何进行整型提升 4.证明整型提升确实存

灰太狼/ 2024年03月22日 16:30/ 0 赞/ 168 阅读

相关 c语言如何精确求得pi值,[updated]C语言和python实现蒙特卡里算法求pi值

高中数学课本上曾经写过用计算机算得pi的近似值，用的是蒙特卡里方法。但是那时的我还不知道编程为何物，今天突然又想到这个问题，尝试了一下，发现并不是那么难，用C语言和python

逃离我推掉我的手/ 2023年10月04日 14:20/ 0 赞/ 12 阅读

相关蒙特卡罗模拟（python实现）

博主专业是电气，电气科学实验中，有时用到蒙特卡罗模拟，今天讲解下。先附上个人电气博文目录传送门。见下。点击即可 [学好电气全靠它，个人电气博文

心已赠人/ 2022年10月14日 00:44/ 0 赞/ 203 阅读

相关泰勒公式求e值——c语言算法实现

题目：通过泰勒公式的变形：e = 1 + (1/1!) + (1/2!) + (1/3!) + … + (1/n!) (1/n! >= 10^-7)来求e的值

ゞ浴缸里的玫瑰/ 2022年09月18日 12:44/ 0 赞/ 307 阅读

相关用蒙特卡洛方法求取PI值

include <cstdio> include <cstring> include <cmath> int x[100001],y[10

骑猪看日落/ 2022年08月21日 04:23/ 0 赞/ 217 阅读

相关【C语言】优先级和求值顺序

![1][] ![2][] ![3][] ![4][] ![5][] > 文章摘至《C和指针》 > > Kenneth A.Reek 著 > 徐波

桃扇骨/ 2022年06月17日 00:00/ 0 赞/ 203 阅读

相关 java蒙特卡洛算法求PI值（可视化）及三门问题

蒙特卡洛算法求PI值工具类 package com.egeo.jframetext; import java.awt.BasicStro

╰+哭是因爲堅強的太久メ/ 2022年05月24日 02:29/ 0 赞/ 210 阅读

相关 MPI_计算pi值程序

include "mpi.h" include <iostream> include <cmath> using namespace

- 日理万妓/ 2021年11月19日 14:08/ 0 赞/ 422 阅读

相关蒙特罗卡π算法（C++语言描述）

圆的面积计算公式为：S=π\r\r 将圆放到一个直角坐标系中，如图黄色部分的面积是S/4=(π\r\r)/4;如果我们将取一个单位圆，则S/4=π/4. ![SouthEa

小咪咪/ 2021年09月29日 08:20/ 0 赞/ 375 阅读