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

问题描述:

请问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
1个回答 分类:数学 2014-11-19

问题解答:

我来补答
非线性整数规划没法做灵敏度分析
再问: 这是2000数学建模的b题,各种优秀论文都用灵敏度分析耶,我们就是命令窗口是暗的。。
再答: 那你看看他怎么写的吧 反正非线性或者整数规划 lingo都不能做灵敏度分析
 
 
展开全文阅读
剩余:2000