在mathematica的 NMinimise 函数中,我进行了以下设置:
NMinimize[{1/rij, x1^2 + y1^2 <= 25, x2^2 + y2^2 <= 25,
x3^2 + y3^2 <= 25, x4^2 + y4^2 <= 25, x5^2 + y5^2 <= 25,
x6^2 + y6^2 <= 25, x7^2 + y7^2 <= 25, x8^2 + y8^2 <= 25,
x9^2 + y9^2 <= 25, x10^2 + y10^2 <= 25},
Join[Take[xi, number], Take[yi, number]]]
其中 xi 和 yi 表示生成变量的列表,x1,x2,x3,x4 等等,直到 x100。而不是设置很多约束,如 x1^2 + y1^2 <= 25、x2^2 + y2^2 <= 25 等等,就像上面 id 一样,为变量类设置一个约束,类似于 x#^2 +y#^2<=25 所以这个问题可以推广到非常大的 n。
我试过输入一个不等式列表,但是,mathematica 似乎不接受它作为输入。
谢谢。