randtest

import random
import time
import matplotlib.pyplot as plt
import numpy as np
from hist import build_hist


class MyPRNG:
    b = 5.65465463546133
    def __init__(self, seed=None):
        # timestamp
        s = seed if seed is not None else time.time()
        self.z = s-int(s)

    def foo(self):
        x = 3**(self.z+MyPRNG.b)
        return x-int(x)

    def __call__(self):
        self.z = self.foo()
        return self.z


def sample(rand, n):
    return [rand() for i in range(n)]

def expval(X):
    return sum(X) / len(X)

def disp(X, M):
    return sum((x-M)**2 for x in X) / (len(X)-1)


if __name__ == "__main__":
    real_M = 0.5
    real_D = 1/12
    size = 10**6

    bars = 20
    random.seed()
    S1 = sample(random.random, size)
    G1 = build_hist(S1, bars)
    M = expval(S1)
    D = np.var(S1)

    print(f'random.random| M: {M}, D: {D}')
    print(f'random.random| errors: M: {abs(real_M-M)/real_M}, D: {abs(real_D-D)/real_D}')

    R = MyPRNG()
    S2 = sample(R, size)
    G2 = build_hist(S2, bars)
    M = expval(S2)
    D = np.var(S2)

    print(f'My| M: {M}, D: {D}')
    print(f'My| errors: M: {abs(real_M-M)/real_M}, D: {abs(real_D-D)/real_D}')

    # plt.style.use('dark_background')
    plt.style.use('seaborn-dark')
    fig, ax = plt.subplots()
    plt.ylim(ymax = 0.06, ymin = 0)
    x_axis = [f'{i/bars:.2f}' for i in range(bars)]

    x = np.arange(bars)
    width = 0.37
    rects1 = ax.bar(x - width/2, G1, width, label='random.random')
    rects2 = ax.bar(x + width/2, G2, width, label='My' )
    ax.set_title('Сравнение функций распределения')
    ax.set_xticks(x)

    ax.set_ylabel('p')
    ax.set_xticklabels(x_axis)
    ax.legend()
    plt.show()