LINGO求解优化问题.docVIP

首 页  一、Lingo简介 1. 目标函数 一个函数解析式你希望求它的最大或最小值 max=函数解析式 或 min=函数解析式max=3*b+2*c^2; min=b^(1/3)-c*k； Lingo的语句以 ; 号结束min=3*x1+5*x2; !x1,x2是非负变量; 3*x1+2*x2=36; 3*x1+5*x2=45; 5. 变量类型 变量类型 说 明 @bin( 变量名) ; 限制该变量为0或1。 @bnd( a,变量名, b); 限制该变量介于a，b之间。 @free(变量名); 允许该变量为负数。 @gin( 变量名); 限制该变量为整数。  二、Lingo高级 sets语句 连续六个月的产量可以用x1,x2,x3,x4,x5,x6表示, 但十二个月的产量用同样的方法表示就显繁琐Lingo可以通过sets语句设置数组功能使问题变得简单。 例 定义数组x, 有x(1),x(2),x(3),x(4)…x(12)个成员以十二个月的产量sets: r/1..12/:x; !r是组的类型名,x数组名; endsets; sets语句以sets开头endsets结束sets: mat/1..4/: x; !mat是组的类型名,x数组名; endsets min=50*x(1)+20*x(2)+30*x(3)+80*x(4); 400*x(1)+200*x(2)+150*x(3)+500*x(4)=500; 3*x(1)+2*x(2)=6; 2*x(1)+2*x(2)+4*x(3)+4*x(4)=10; 2*x(1)+4*x(2)+x(3)+5*x(4)=8; 有时我们要用到常数数组比如在400*x(1)+200*x(2)+150*x(3)+500*x(4)=500 中x(1), x(2), x(3), x(4)的系数400, 200, 150, 500，此时,可用data语句例 定义数组a, 其中a(1)=a(2)=200，a(3)=150，a(4)=500。 sets: l/1..4/: a,x; endsets data: a=400 200 150 500; enddata data语句是以data开头enddata结尾sets: l/1..4/: x, a; endsets data: a=7 2 3 9; !a(1)=7, a(2)=2, a(3)=3, a(4)=9; enddata max=x(1)*a(3)+x(2)*a(1)+x(3)*a(4)+x(4)*a(2); x(1)+x(4)-x(2)-x(3)a(1); x(4)+2*x(2)a(4); x(1)+x(3)a(1); Lingo含有一些针对数组的命令，方便了数组的使用。 @for循环语句: @for(数组类型名(i): 循环的语句); 示例程序4（@gin语句） sets: r/1..5/: a, b; endsets data: a= 3.3 4.6 2.7 7.1 10.3; enddata max=a(1)*b(1)-a(2)*b(2)+a(3)*b(3)-a(4)*b(4); @for(r(i): b(i) a(i)); !等价于b(1)a(1)；b(2)a(2)；b(3)a(3)；b(4)a(4); @for(r(i): @gin(b(i))); !等价于@gin(b(1))；@gin(b(2))；@gin(b(3))；@gin(b(4)); sets: r/1..5/: a, b; endsets data: a= 3.3 4.6 2.7 7.1 10.3; enddata max=@sum(r(i):b(i))+@sum(r(i):b(i)/a(i)) +@sum(r(i):b(i)*a(i)); !等价于max=b(1)+b(2)+b(3)+b(4)+b(1)/a(1)+b(2)/a(2)+b(3)/a(3)+b(4)/a(4); @for(r(i): b(i) a(i)); @for(r(i): @gin(b(i))); 示例程序6（bnd语句） sets: m/1..4/: x, need, g, y; endsets data: need=4000 2000 3000 10000; enddata min=30000*@sum(m(i): x(i))+30*@sum(m(i): g(i)); g(1)=600+y(1)*(x(2)+x(3)+x(4))-need(1); g(2)=g(1)+y(2)*(@sum(m(i): x(i))-x(2))-need(2); ! @sum(m(i):x(i))-x(2)等价于x(1)+x(3)+x