def stratified_sampling4(M,nsims):
L=nsims/M
i=0
noises=np.empty((M,nsims))
for i in range(M):
noises[i,:]=norm.ppf(np.random.uniform(i/M,(i+1)/M,nsims))
return noises
微信
公众号
扫码关注公众号