生成指数分布随机数———python代码

如下为生成指数分布随机数python代码，并将生成的数保存到csv文件内

# 生成指数分布随机变量
import matplotlib.pyplot as plt
import numpy as np
import pandas as pd
import csv


class Random(object):

    def __init__(self, x=None):  # 初始函数，x为初始种子，缺省为时间值
        if x is None:
            import time
            self.seed = int(time.time())
        else:
            self.seed = int(x)
        self._seed = self.seed

    def random_rand(self):
        '''
        线性同余产生0~1之间的随机数

        X(k) = [a * X(k-1) + c] mod m
        '''
        # 以下参数值参照GCC编译器
        m = 2 ** 32
        a = 1103515245
        c = 12345
        self._seed = (a * self._seed + c) % m
        return self._seed / float(m - 1)

JT_random = Random(1)


def exponential_rand(lam):
    if lam <= 0:                      # 检验参数λ的合法性
        return -1
    u = JT_random.random_rand()       # 生成随机数
    return (-1.0 / lam) * np.log(u)   # 根据实现原理的公式 计算指数分布随机变量


a = []
b={}
for i in range(0,1000):
    c = int(JT_random.random_rand()+1)
    k= exponential_rand(c)*10
    a.append(int(k))
for i in a:
    if i in b.keys():
        b[i] += 1
        continue
    b[i]=1n
for i,j in b.items():
    print(i,"  ",j)


with open("kkk.csv", "a", newline='') as f:
    writer = csv.writer(f)
    writer.writerow(b.keys())
    writer.writerow(b.values())