generators.py 文件源码

def stratified_sampling2(nnoises,nsims):
    """On divise l'intervalle [0,1] en plusieurs stratas, m stratas"""
    noises = np.empty((nnoises,nsims))
    m=int(nnoises/nsims)
    i=0
    """Pour chaque strata, on prend l échantillons suivant la loi uniforme
    tel que nsims=m*l et enfin on prend le pdf du cdf."""
    for i in range(m+1,1):
        noises[i,:]=norm.ppf(np.random.uniform(i/m,(i+1)/m,nsims))
    return noises