请问lingo程序可以运行,怎么才可以灵敏度分析呢?下面是题目和输入的内容(已激活灵敏度计算)
![](https://img.yulucn.com/upload/9/4f/94f04636d69b20692cc2f9e5d783bfa1_thumb.jpg)
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.2x0964.2 92 96 106 121.2 128 142
215.7 205.3 190.2 171.6 111x0964.6 86 71.2 114.2 142x09146 156 171.2 178x09192
230.7 220.3 200.2 181.6 121x09105.5x0996 86.2 48.2 82 86 96 111.2 118 132
260.7 250.3 235.2 216.6 156x09139.6x09131 116.2 84.2 62x0951 61x0976.2 83 97
255.7 245.3 225.2 206.6 146x09130.5x09121 111.2 79.2 57x0933 51x0971.2 73 87
265.7 255.3 235.2 216.6 156x09140.5x09131 121.2 84.2 62x0951 45x0926.2 11 28
275.7 265.3 245.2 226.6 166x09150.5x09141 131.2 99.2 76x0966 56x0938.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