simulationexam/gibbs_sampler.ipynb

In [1]:

%pylab inline
from network.network import (is_valid_network,
                             get_number_of_activated_antenna,
                             can_antenna_be_turned_on,
                             pp_network,
                             turn_off_antenna,
                             turn_on_antenna)
from network.gibbs import get_random_networks, estimate_number_of_activated_antenna
from numpy import mean
from random import Random

Populating the interactive namespace from numpy and matplotlib

In [5]:

random = Random()
antennas = []
box_size = 8
N = 500
for _, network in zip(range(N), get_random_networks(box_size, random)):
    antennas.append(get_number_of_activated_antenna(network))
    
    

In [6]:

hist(antennas, bins=100)
title("Estimation de la distribution du nombre \n d'antennes allumées dans un réseau")
savefig("/home/didier/Downloads/gibbs.png")

In [ ]: