如何使用AMPL中的正态分布生成一个随机数矩阵？

Question

如何使用AMPL中的正态分布生成一个随机数矩阵，其平均值和方差如下所示？

param mean := 
1 45 
2 35 
3 40; 

param variance : 
    1  2  3 := 
1 1  -2  -1 
2 -2  36  -8 
3 -1  -8  9; 

Answer 1

AMPL为您的参数支持一些随机数函数。 你可以尝试这样的事情：

param matrix {x in dimx, y in dimy} 
      = Normal(mean[x], variance[x, y]); 

如果你想确保只得到正数，你应该用最大功能包的正常功能。

max(Normal(mean[x], variance[x, y]), 0);