php网站建设素材,网站开发方案及报价,网盘建网站,白云手机网站建设用蒙特卡洛方法算pi#xff0d;基于python和R语言 最近follow了MOOC上一门python课#xff0c;开始学Python。同时#xff0c;买来了概率论与数理统计#xff0c;准备自学一下统计。#xff08;因为被鄙视过不是统计专业却想搞数据分析#xff09; 有趣的是书里面有一块讲…用蒙特卡洛方法算pi基于python和R语言 最近follow了MOOC上一门python课开始学Python。同时买来了概率论与数理统计准备自学一下统计。因为被鄙视过不是统计专业却想搞数据分析 有趣的是书里面有一块讲蒲丰投针计算Pi这是一种随机模拟法也就是蒙特卡洛法。蒲丰投针之于我太难暂时没想到怎么用计算机模拟这一过程。 python课中老师也提到用随机模拟法也就是蒙特卡洛法MonteCarlo用计算机模拟几千次实验计算pi的近似值。好巧。 就拿python课中的方法来近似计算pi分别用python和R实现一下。 至于实验是怎样的截图老师的PPT。 我用的是python 3 python代码 from random import random
from math import sqrt
from time import clock
darts2**22
hist0
clock()
for i in range(1,darts):x,yrandom(),random()distsqrt(x**2y**2)if dist1.0:histhist1
pi4*(hist/darts)
print(pi is %s%pi)
print(elaspe is %ss%clock())python运行结果 pi is 3.143444061279297
elaspe is 85.991785 R 代码 #蒙特卡洛方法求pi
hist - 0
darts - 2^22
start -proc.time()
for (i in 1:darts){x - runif(1);y - runif(1)if(sqrt(x^2y^2)1) nexthisthist1
}
pi - 4*hist/darts
proc.time()-start
print(paste0(pi is ,pi))R运行结果 proc.time()-start用户 系统 流逝
31.537 2.477 34.153 print(paste0(pi is ,pi))
[1] pi is 3.14076137542725总结R和python都挺好用的下一步准备试着用python写点小爬虫程序。转载于:https://www.cnblogs.com/litao1105/p/4903079.html
