Monte Carlo method

1 蒙特卡洛法(Monte Carlo method)计算\(\pi\)的值 用随机函数runif()随机生成N个在1个单位正方形内的散点,根据三角公式筛选出在圆内的点,数据如下: ## 生成基础数据 set.seed(2021) N <- 5000 df01 <- tibble( x = runif(N, -0.5, 0.5), y = runif(N, -0.5, 0.5) ) %>% mutate( z = if_else( x^2 + y ^2 < 0.25, 1, 0), p = runif(N,0,2 * pi), x2 = cos(p) * 0.5, y2 = sin(p) * 0.5 ) knitr::kable(head(df01), caption = "随机数据", col.names = c("散点X", "散点Y", "是否在圆内", "随机生成0-180度", "圆X", "圆Y"), digits = 3, align = "c" ) 表1....

January 7, 2022 · 1 min · luofei