python实现统计模拟【第3章】【伪随机数的产生】

糖逗

1. #伪随机数的产生
   
2. #生成（0,1）上的均匀分布
   
3. #方法1：乘同余法
   
4. def random_number1(seed = 1233223, m = pow(2, 31)-1, a = pow(7, 5), number = 10):
   
5.     res = []
   
6.     xi_1 = seed
   
7.     for i in range(number):
   
8.         xi = (a * xi_1)%m
   
9.         res.append(xi/m)
   
10.         xi_1 = xi
    
11.     return res
    
12. temp = random_number1(number = 12)
    
13. print(temp)
    
14. 
    
15. #方法2:混合同余法
    
16. def random_number2(seed = 1233223, m = pow(2, 31)-1, c = 1,a = pow(7, 5), number = 10):
    
17.     res = []
    
18.     xi_1 = seed
    
19.     for i in range(number):
    
20.         xi = (a * xi_1 + c) % m
    
21.         res.append(xi / m)
    
22.         xi_1 = xi
    
23.     return res
    
24. temp = random_number2(number = 10)
    
25. print(temp)
    

复制代码

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