我该如何在Ampl cplex中进行蒙特卡洛模拟

Question

我正在最小化目标函数。

i有一个参数> param preq {i in 1..n：=统一(0.002，pmax/n(pmax为.0095，n为12。

如何使用此统一参数在模型文件中引入蒙特卡洛模拟？

Answer 1

ampl语言具有循环。下一个示例应适用于任何求解器：

param n := 3;
set I := 1..n;
# 
# random numbers
#
param p{i in I} := Uniform(0,1);
display p;
#
# minimal model
# uses a scalar q
#
param q;
var x;
minimize z:x;
e: x = 2*q;
#
# solve in loop
# assign q before solving
#
param results{I};
for {i in I} {
   let q := p[i];
   solve;
   let results[i] := z;
} 
display results;
end;

AMPL网站有更多信息：例如https://ampl.com/new/loop1.html