VBA实现蒙特卡罗方法Excel示例-计算圆周率PI

6关注
36粉丝

版主

已卖：5810份资源

学科带头人

64%

还不是VIP/贵宾

-

TA的文库 其他...

私募

信托行业研究

0%

威望: 0 级
论坛币: 2066 个
通用积分: 4061.6661
学术水平: 122 点
热心指数: 146 点
信用等级: 111 点
经验: 55803 点
帖子: 1836
精华: 0
在线时间: 2523 小时
注册时间: 2009-2-16
最后登录: 2025-3-9

楼主

peng3409 发表于 2010-6-3 20:11:40 |AI写论文

是否 +2 论坛币

k人参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

这是个比较经典的例子，就是采用蒙特卡罗方法计算圆周率PI，计算原理就是采用圆的面积与相应外切正方形的面积比的4倍就是PI值，
因此采用随机模拟，统计落入圆内的数据点的数目与总数相比，使用VBA实现的例子源代码如下：
01.Option Explicit
02.'================================
03.' 蒙特卡罗方法计算圆周率示例
04.' Monte_Carlo_PI
05.'
06.'================================
07.Sub Monte_Carlo_PI()
08. Dim i As Long
09. Dim x As Double
10. Dim y As Double
11. Dim dist As Double
12.

13. Dim cntInCircle As Long
14. Dim iterMAX As Long
15. iterMAX = 100000
16. cntInCircle = 0
17.

18. Randomize
19. For i = 1 To iterMAX
20.       x = Rnd()
21.       y = Rnd()
22.       dist = Sqr((x - 0.5) ^ 2 + (y - 0.5) ^ 2)
23.       If dist <= 0.5 Then
24.          cntInCircle = cntInCircle + 1
25.       End If
26. Next
27. Debug.Print 4 * cntInCircle / iterMAX, _
28.    Application.WorksheetFunction.Pi(), _
29.    Abs(Application.WorksheetFunction.Pi() _
30.          - 4 * cntInCircle / iterMAX)
31.End Sub

注：增加iterMAX的值，可以更加逼近PI值