I did a sampling on a first interval and I would like to add more points. I have a condition is to take as arguments the precedent points to complete the first sampling. I don't really know/found how is possible to di it?
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D
from smt.sampling_methods import LHS
%matplotlib notebook
#first interval
# a=[1.0E-06, 1.0E-07]
# z=[1.0E3, 1.0E4]
# e=[1.0e5, 1.0e6]
#first sampling
lambda1 = np.array([[7.075e-07, 1.225e+03, 7.075e+05],
[1.225e-07, 6.175e+03, 7.975e+05],
[2.125e-07, 6.625e+03, 9.775e+05],
[9.775e-07, 5.725e+03, 5.275e+05],
[3.475e-07, 4.825e+03, 4.825e+05],
[8.425e-07, 7.525e+03, 2.125e+05],
[9.325e-07, 7.975e+03, 6.175e+05],
[4.825e-07, 8.425e+03, 7.525e+05],
[4.375e-07, 5.275e+03, 2.575e+05],
[7.525e-07, 2.125e+03, 6.625e+05],
[3.925e-07, 3.925e+03, 4.375e+05],
[2.575e-07, 3.475e+03, 1.675e+05],
[6.175e-07, 2.575e+03, 3.025e+05],
[3.025e-07, 1.675e+03, 3.925e+05],
[7.975e-07, 9.325e+03, 3.475e+05],
[5.275e-07, 8.875e+03, 8.875e+05],
[1.675e-07, 3.025e+03, 5.725e+05],
[8.875e-07, 9.775e+03, 1.225e+05],
[5.725e-07, 7.075e+03, 8.425e+05],
[6.625e-07, 4.375e+03, 9.325e+05]])
#new interval
a=[10**-5.0, 10**-8.0]
z=[10**2.0, 10**5.0]
e=[10**4.0, 10**7.0]
xlimits = np.array([a, z,e])
sampling = LHS(xlimits=xlimits,criterion='ese')
num = 20
x = sampling(num)
Aucun commentaire:
Enregistrer un commentaire