造船用钢优化问题.docVIP

一、问题的描述 有一家钢铁公司收到一份500吨造船用钢的订单。对这些造船用钢有如下要求（品质）： 表格6.1：造船钢品质要求 化学元素 最低含量 最高含量碳（C） 2 3 钢（Cu） 0.4 0.6 锰（Mn） 1.2 1.65 此公司储存有七种不同的原材料，都可以用于制造这种钢。表6.1列出了这些原材料的品质、库存量以及价格。 原材料 C% Cu% Mn% 可用库存（吨） 单价（欧元/吨）铁碳合金1 2.5 0 1.3 400 200 铁碳合金1 3 0 0.8 300 250 铁碳合金2 0 0.3 0 600 150 钢碳合金1 0 90 0 500 220 钢碳合金2 0 96 4 200 240 吕碳合金1 0 0.4 1.2 300 100 吕碳合金2 0 0.6 0 250 165我们的目标是求出各种原材料各取多少才能使成本最低。 二、问题假设 设xi（i=1、2、3......7）分别为这七种原材料的用量，pi为第i种原材料的单价，ci为第i种原材料的碳含量，Cui为第i种原材料的钢含量，Mni为第i种原材料的锰含量。Yi为第i种原材料的可用库存。 三、模型建立与求解 目标函数：min 约束条件：  xi=yi xi=yi ∑xi=500  xi=yi  Lingo代码为： model: sets:kind/1..7/:x,p,c,cu,mn,y; j/1,2/:t,g,m; endsets data: p=200 250 150 220 240 100 165; y=400 300 600 500 200 300 250; c=2.5 3 0 0 0 0 0; cu=0 0 0.3 90 96 0.4 0.6; mn=1.3 0.8 0 0 4 1.2 0; t=2 3; g=0.4 0.6; m=1.2 1.65; endsets min=@sum(kind(i):(p(i)*x(i))); @for(kind(i):@sum(ci*xi)500*t1); @for(kind(i):@sum(ci*xi)500*t2); @for(kind(i):@sum(cui*xi)500*g1); @for(kind(i):@sum(cui*xi)500*g2); @for(kind(i):@sum(mni*xi)500*m1); @for(kind(i):@sum(mni*xi)500*m2); @for(kind(i):@sum(xi)=500); @for(kind(i):@sum(x(i))=500); @for(kind(i):x(i)y(i)); 输入lingo运行的程序为： min=200*x1+250*x2+150*x3+220*x4+240*x5+100*x6+165*x7; 2.5*x1+3*x2=2*500; 2.5*x1+3*x2=3*500; 0.3*x3+90*x4+96*x5+0.4*x6+0.6*x7=0.4*500; 0.3*x3+90*x4+96*x5+0.4*x6+0.6*x7=0.6*500; 1.3*x1+0.8*x2+4*x5+1.2*x6=1.2*500; 1.3*x1+0.8*x2+4*x5+1.2*x6=1.65*500; X