我是一个新手,不太懂如何使用proc optmodel,有些语法还很困惑。
这是我的数据集。
data opt_test;
input ID GRP $ x1 MIN MAX y z;
cards;
2 F 10 9 11 1.5 100
3 F 10 9 11 1.2 50
4 F 11 9 11 .9 20
8 G 5 4 6 1.2 300
9 G 6 4 6 .9 200
1 H 21 18 22 1.2 300
7 H 20 18 22 .8 1000
;
run;
这里有几个要点:
GRP内的ID必须具有相同的x2,这受到MIN和MAX的限制。我现在希望根据y的值进一步限制x2的增加/减少。如果y<1,则不希望x2低于.95*x1。如果y>1,则不希望x2超过1.05*x1。我在网上查找并尝试了一些方法来实现这一点。以下是我的最新尝试,cond_1和cond_2是感兴趣的问题,因为其他所有内容都有效:
proc optmodel;
set<num> ID;
string GRP{ID};
set GRPS = setof{i in ID} GRP[i];
set IDperGRP{gi in GRPS} = {i in ID: GRP[i] = gi};
number x1{ID};
number MIN{ID};
number MAX{ID};
var x2{gi in GRPS} >= max{i in IDperGRP[gi]} MIN[i]
<= min{i in IDperGRP[gi]} MAX[i]
;
impvar x2byID{i in ID} = x2[GRP[i]];
number y{ID};
number z{ID};
read data opt_test into
ID=[ID]
GRP
x1
MIN
MAX
y
z
;
max maximize = sum{gi in GRPS} sum{i in IDperGRP[gi]}
(x2[gi]) * (1-(x2[gi]-x1[i])*y[i]/x1[i]) * z[i];
con cond_1 {i in ID}: x2[i] >=
if y[i]<1 then .95*x1[i] else 0;
con cond_2 {i in ID}: x2[i] <=
if y[i]>=1 then 1.05*x1[i] else 99999999;
solve;
create data results from [ID]={ID} x2=x2byID GRP x1 MIN MAX y z;
print x2 maximize;
quit;
var
,则结果表达式将是非线性的。在这种情况下,条件只涉及到num
,因此约束是线性的。请查看我的答案以获取更多详细信息。 - Leo