寻找最小的六峰驼峰背函数所在的区域-2.1 <= X（I）<= 2.1。此功能与目标函数值的两个全球最低-1.0316284……以及四个目标函数值更高的局部极小值。

rng默认％用于重现objconstr = @ (x) (4 * x (: 1)。^ 2 - 2.1 * x (: 1)。^ 4 + x (: 1)。^ 6/3...+ X（：，1）* X（：，2） -  4 * X（：，2）^ 2 + 4 * X（：，2）。。。^ 4）;磅= [-2.1,-2.1];乌兰巴托=磅;x = surrogateopt (objconstr磅,乌兰巴托)

Surrogateopt停止，因为它超过了“options.MaxFunctionEvaluations”功能评价极限集。

X =1×20.0895 - -0.7130

求Rosenbrock函数的最小值

受非线性约束，该溶液在于半径1/3的周围的点[1 / 3,1 / 3]的盘：

$$(x(1)-1/3.{)}^2+(x(2)-1/3.{)}^2\le (1/3.{)}^2$$。

要做到这一点，写一个函数objconstr(x)返回结构字段中Rosenbrock函数的值Fval，并在表单中返回非线性约束值 $$c(x)\le 0$$在结构领域Ineq。

类型objconstr

函数f = objconstr（x）的f.Fval = 100 *（X（2） -  X（1）^ 2）^ 2 +（1  -  X（1））^ 2;f.Ineq =（X（1）-1/3）^ 2 +（X（2）-1/3）^ 2  - （1/3）^ 2;

呼叫surrogateopt使用0下界和每个组件上的2/3的上限。

磅= (0,0);乌兰巴托= (2/3,2/3);[x, fval exitflag] = surrogateopt (@objconstr磅,乌兰巴托)

Surrogateopt停止，因为它超过了“options.MaxFunctionEvaluations”功能评价极限集。

X =1×20.6541 - 0.4277

FVAL = 0.1197

exitflag = 0

检查解的非线性约束的值。

DISP（objconstr（x）的.Ineq）

6.7195e-04

约束函数值是接近零，表示该约束是在溶液活性。

找出最小的ps_example二维变量的函数x它的第一个组件被限制为整数值，所有的组件都在-5和5之间。

INTCON = 1;rng默认％用于重现objconstr = @ps_example;磅= [-5，-5];UB = [5,5];X = surrogateopt（objconstr，LB，UB，INTCON）

Surrogateopt停止，因为它超过了“options.MaxFunctionEvaluations”功能评价极限集。

X =1×2-5.0000 -0.0005

尽量减少在该地区六驼峰的骆驼回功能-2.1 <= X（I）<= 2.1。此功能与目标函数值的两个全球最低-1.0316284……以及四个目标函数值更高的局部极小值。

为了系统地搜索区域，使用出发点的规则网格。设置120函数计算的最大数量。使用'surrogateoptplot'绘图功能。要了解'surrogateoptplot'情节，看解读surrogateoptplot。

rng默认％用于重现objconstr = @ (x) (4 * x (: 1)。^ 2 - 2.1 * x (: 1)。^ 4 + x (: 1)。^ 6/3...+ X（：，1）* X（：，2） -  4 * X（：，2）^ 2 + 4 * X（：，2）。。。^ 4）;磅= [-2.1,-2.1];乌兰巴托=磅;[xpt材料,Ypts] = meshgrid(三3);startpts = [xpt材料(:),Ypts (:));选择= optimoptions (“surrogateopt”,“PlotFcn”,'surrogateoptplot',...“InitialPoints”，startpts，“MaxFunctionEvaluations”,120);x = surrogateopt (objconstr、磅、乌兰巴托、期权)

Surrogateopt停止，因为它超过了“options.MaxFunctionEvaluations”功能评价极限集。

X =1×2-0.0904 0.7127

创建表示该地区的六驼峰的骆驼背功能的结构问题-2.1 <= X（I）<= 2.1。设置120函数计算的最大数量。

rng默认％用于重现objconstr = @ (x) (4 * x (: 1)。^ 2 - 2.1 * x (: 1)。^ 4 + x (: 1)。^ 6/3...+ X（：，1）* X（：，2） -  4 * X（：，2）^ 2 + 4 * X（：，2）。。。^ 4）;选择= optimoptions (“surrogateopt”,“MaxFunctionEvaluations”,120);问题=结构('目的'objconstr,...'磅'[-2.1,-2.1],...乌兰巴托的[2.1,2.1],...“选项”，期权，...“求解”,“surrogateopt”);x = surrogateopt(问题)

Surrogateopt停止，因为它超过了“options.MaxFunctionEvaluations”功能评价极限集。

X =1×20.0895 - -0.7130

将六峰骆驼背函数最小化，同时返回最小值点和目标函数值。设置选项以禁止所有其他显示。

rng默认％用于重现objconstr = @ (x) (4 * x (: 1)。^ 2 - 2.1 * x (: 1)。^ 4 + x (: 1)。^ 6/3...+ X（：，1）* X（：，2） -  4 * X（：，2）^ 2 + 4 * X（：，2）。。。^ 4）;磅= [-2.1,-2.1];乌兰巴托=磅;选择= optimoptions (“surrogateopt”,'显示',“关”,“PlotFcn”，[]）;[X，FVAL] = surrogateopt（objconstr，LB，UB，选项）

X =1×20.0895 - -0.7130

FVAL = -1.0316

通过请求监视代理优化过程surrogateopt返回输出。使用'surrogateoptplot'绘图功能。要了解'surrogateoptplot'情节，看解读surrogateoptplot。

rng默认％用于重现objconstr = @ (x) (4 * x (: 1)。^ 2 - 2.1 * x (: 1)。^ 4 + x (: 1)。^ 6/3...+ X（：，1）* X（：，2） -  4 * X（：，2）^ 2 + 4 * X（：，2）。。。^ 4）;磅= [-2.1,-2.1];乌兰巴托=磅;选择= optimoptions (“surrogateopt”,“PlotFcn”,'surrogateoptplot');[X，FVAL，exitflag，输出] = surrogateopt（objconstr，LB，UB，选项）

Surrogateopt停止，因为它超过了“options.MaxFunctionEvaluations”功能评价极限集。

X =1×20.0895 - -0.7130

FVAL = -1.0316

exitflag = 0

输出=结构体字段:elapsedtime：46.9309 funccount：200 constrviolation：0 ineq：[1X0双] rngstate：1x1的结构]消息： 'Surrogateopt停止，因为它由超过功能评价极限集......'

通过设置最小的最大函数计算次数来快速完成代理优化。为了准备重新启动优化的可能性，请求所有解算器输出。

rng默认％用于重现objconstr = @ (x) (4 * x (: 1)。^ 2 - 2.1 * x (: 1)。^ 4 + x (: 1)。^ 6/3...+ X（：，1）* X（：，2） -  4 * X（：，2）^ 2 + 4 * X（：，2）。。。^ 4）;磅= [-2.1,-2.1];乌兰巴托=磅;选择= optimoptions (“surrogateopt”,“MaxFunctionEvaluations”,20);[x, fval exitflag、输出试验]= surrogateopt (objconstr、磅、乌兰巴托、期权);

Surrogateopt停止，因为它超过了“options.MaxFunctionEvaluations”功能评价极限集。

从之前评估的点开始，优化另外20个功能评估。

options.InitialPoints =试验;[x, fval exitflag、输出试验]= surrogateopt (objconstr、磅、乌兰巴托、期权);

Surrogateopt停止，因为它超过了“options.MaxFunctionEvaluations”功能评价极限集。

通过比较这40个函数的求值图和搜索全球最低，你会看到重新启动代理优化是不一样的具有连续运行的求解器。

为了使因碰撞或任何其他原因重新启动替代优化，设置检查点的文件名。

选择= optimoptions (“surrogateopt”,'CheckpointFile','checkfile.mat');

创建一个优化问题并设置少量的函数计算。

rng默认％用于重现objconstr = @ (x) (4 * x (: 1)。^ 2 - 2.1 * x (: 1)。^ 4 + x (: 1)。^ 6/3...+ X（：，1）* X（：，2） -  4 * X（：，2）^ 2 + 4 * X（：，2）。。。^ 4）;磅= [-2.1,-2.1];乌兰巴托=磅;opts.MaxFunctionEvaluations = 30;[X，FVAL，exitflag，输出] = surrogateopt（objconstr，LB，UB，OPTS）

Surrogateopt停止，因为它超过了“options.MaxFunctionEvaluations”功能评价极限集。

X =1×20.0067 - -0.7343

FVAL = -0.9986

exitflag = 0

输出=结构体字段:elapsedtime：11.7528 funccount：30 constrviolation：0 ineq：[1×0双] rngstate：[1×1结构]消息： 'Surrogateopt停止，因为它超过了↵'options.MaxFunctionEvaluations功能评价极限集'”

设置选项使用100分功能的评价（这意味着70多家已经完成），然后重新启动优化。

opts.MaxFunctionEvaluations = 100;[X2，fval2，exitflag2，输出2] = surrogateopt（'checkfile.mat'，OPTS）

Surrogateopt停止，因为它超过了“options.MaxFunctionEvaluations”功能评价极限集。

x2 =1×20.0895 - -0.7130

fval2 = -1.0316

exitflag2 = 0

输出2 =结构体字段:elapsedtime: 64.8767 funccount: 100 constrviolation: 0 ineq:[1×0双]rngstate: [1×1 struct]消息:“Surrogateopt停止,因为它超过了函数评估限制由↵options.MaxFunctionEvaluations’。”