python实现统计模拟【第3章】【利用随机数求积分】

糖逗

1. #利用随机数求积分
   
2. #对pi进行估计
   
3. import numpy as np
   
4. import math
   
5. 
   
6. class comput:
   
7.     def __init__(self, n = 7000000):
   
8.         self.n = n
   
9.    
   
10.     def random_number2(self, seed = 1, m = pow(2, 31)-1, c = 19,a = pow(7, 5)):   
    
11.         res = []
    
12.         xi_1 = seed
    
13.         for i in range(self.n):
    
14.             xi = (a * xi_1 + c) % m
    
15.             res.append(xi / m)
    
16.             xi_1 = xi
    
17.         return np.array(res)
    
18.    
    
19.     def comput_pi(self):
    
20.         x = self.random_number2(seed = 23)
    
21.         y = self.random_number2(seed = 61)
    
22.         res_temp = pow((2*x-1), 2) + pow((2*y-1), 2)
    
23.         res = np.sum(res_temp <= 1) / self.n
    
24.         return res*4
    
25. 
    
26. res = comput().comput_pi()
    
27. print(res)

复制代码

参考：
1.Sheldon M.Ross《统计模拟》