Recipe 576870 revision 1

from numpy.random import normal
from scipy.stats import norm as normStat


CRAZINESS = 0.1
MIN = 50 
MAX = 300
STEP = 5
AVG = 180
DEV = 40 
p = normStat(AVG, DEV).pdf

def noise():
    return 1 + CRAZINESS * normal(1)

for x in range(MIN, MAX, STEP):
    print str(x), "\t",   p(x) * noise()

• « Back to Recipe 576870

History

• revision 1 (14 years ago)