请问lingo程序可以运行,怎么才可以灵敏度分析呢?下面是题目和输入的内容(已激活灵敏度计算)model:sets:sell/1 2 3 4 5 6 7/:p,s,x,m;need/1 2 3 4 5 6 7 8 9 10 11 12 13 14 15/:a,b,d;link(sell,need):c,y;endsetsdata:p=16

来源:学生作业学帮网 编辑:学帮网 时间:2024/05/07 00:49:38

请问lingo程序可以运行,怎么才可以灵敏度分析呢?下面是题目和输入的内容(已激活灵敏度计算)

model:
sets:
sell/1 2 3 4 5 6 7/:p,s,x,m;
need/1 2 3 4 5 6 7 8 9 10 11 12 13 14 15/:a,b,d;
link(sell,need):c,y;
endsets
data:
p=160 155 155 160 155 150 152;
s=800 800 1000 2000 2000 2000 3000;
d=104,301,750,606,194,205,201,680,480,300,220,210,420,500,0;
c=170.7 160.3 140.2 98.6 38 20.5 3.1 21.2\x0964.2 92 96 106 121.2 128 142
  215.7 205.3 190.2 171.6 111\x0964.6 86 71.2 114.2 142\x09146 156 171.2 178\x09192
  230.7 220.3 200.2 181.6 121\x09105.5\x0996 86.2 48.2 82 86 96 111.2 118 132
  260.7 250.3 235.2 216.6 156\x09139.6\x09131 116.2 84.2 62\x0951 61\x0976.2 83 97
  255.7 245.3 225.2 206.6 146\x09130.5\x09121 111.2 79.2 57\x0933 51\x0971.2 73 87
  265.7 255.3 235.2 216.6 156\x09140.5\x09131 121.2 84.2 62\x0951 45\x0926.2 11 28
  275.7 265.3 245.2 226.6 166\x09150.5\x09141 131.2 99.2 76\x0966 56\x0938.2 26 2;
enddata
min=@sum(sell(i):p(i)*x(i))+0.1*@sum(need(j):a(j)*(a(j)+1)/2+b(j)*(b(j)+1)/2)+@sum(link(i,j):c(i,j)*y(i,j));
@for(sell(i):x(i)>=500*m(i));
@for(sell(i):x(i)<=s(i)*m(i));
@for(sell(i):@sum(need(j):y(i,j))=x(i));
@for(need(j):@sum(sell(i):y(i,j))=a(j)+b(j));
@for(need(j)|j#ne#15:a(j)+b(j+1)=d(j));
a(1)=0;b(15)=0;
@for(sell(i):@bin(m(i)));
end

非线性整数规划没法做灵敏度分析