|
代码如下:
require(DiceDesign)
require(DoE.wrapper)
set.seed(27662)
LHD=lhs.design(200,7,type="dmax",factor.names=list(
Angle1=c(80,90),Angle2=c(90,100),
Angle3=c(32,40),Radial2=c(10,15),
Radial3=c(25,35),Length1=c(120,130),
Length2=c(0,10)),range=0.98,niter_max=2000)
LHDnew=maximinESE_LHS(LHD,T0=0.005*phiP(LHD,p=50),inner_it=200,J=50,it=4)
plot(LHDnew$design$Angle1,LHDnew$design$Angle2)
plot(LHDnew$critValues,type="l")
LHDopt=LHDnew$design
运行上述代码时,每次运行得到的LHDopt的结果都不一样,如何才能确保运行得出的LHDopt不变?
|
|