资源描述:
《数学模型数学建模 第六次作业 图论组合优化实验》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、数学模型第六次作业图论(组合优化)实验6.1实验目的与要求l学会用图论(组合优化)的方法或思想建模l学会用LINGO软件求解组合优化问题l简历相应的数学模型,并对计算结果进行分析和讨论。6.2基本实验1.设备更新问题某公司需要对一台已经使用了2年的机器确定今后4年(n=4)的最优更新策略.公司要求,用了6年的机器必须更新,购买一台新机器的价格是100万元,表6.1给出了该问题的数据,请给出设备的更新策略。解:根据题意可以设A,B,C,D…为判断年度,用数字表示机器已经使用的时间,则有第1年:机器使用的情况可以表示为
2、A2;第2年:若第1年机器不更新则情况表示成B3,若第一年机器更新则42第二年应表示成B1;……使用Lingo软件进行计算并取最优解,编程如下:model:sets:nodes/A2,B3,B1,C4,C2,C1,D5,D3,D2,D1,E6,E4,E3,E2,E1,F/;arcs(nodes,nodes)/A2,B3A2,B1B3,C4B3,C1B1,C2B1,C1C4,D5C4,D1C2,D3C2,D1C1,D2C1,D1D5,E1D5,E6D3,E4D3,E1D2,E3D2,E1D1,E2D1,E1E6,FE
3、4,FE3,FE2,FE1,F/:c,x;endsetsdata:c=17.3-20.215.7-30.218.4-0.213.8-50.217.3-20.218.4-0.212.2-70.215.7-30.217.3-20.218.4-0.2530506080;enddatan=@size(nodes);max=@sum(arcs:c*x);@sum(arcs(i,j)
4、i#eq#1:x(i,j))=1;@for(nodes(i)
5、i#ne#1#and#i#ne#n:@sum(arcs(i,j):x(i,j))-
6、@sum(arcs(j,i):x(j,i))=0);@sum(arcs(j,i)
7、i#eq#n:x(j,i))=1;@for(arcs:@bin(x));End得到结果如下:Globaloptimalsolutionfound.Objectivevalue:139.0000Objectivebound:139.0000Infeasibilities:0.000000Extendedsolversteps:0Totalsolveriterations:042VariableValueReducedCostN16.00
8、0000.000000C(A2,B3)17.300000.000000C(A2,B1)-20.200000.000000C(B3,C4)15.700000.000000C(B3,C1)-30.200000.000000C(B1,C2)18.400000.000000C(B1,C1)-0.20000000.000000C(C4,D5)13.800000.000000C(C4,D1)-50.200000.000000C(C2,D3)17.300000.000000C(C2,D1)-20.200000.000000C(C
9、1,D2)18.400000.000000C(C1,D1)-0.20000000.000000C(D5,E1)12.200000.000000C(D5,E6)-70.200000.000000C(D3,E4)15.700000.000000C(D3,E1)-30.200000.000000C(D2,E3)17.300000.000000C(D2,E1)-20.200000.000000C(D1,E2)18.400000.000000C(D1,E1)-0.20000000.000000C(E6,F)5.0000000
10、.000000C(E4,F)30.000000.000000C(E3,F)50.000000.000000C(E2,F)60.000000.000000C(E1,F)80.000000.000000X(A2,B3)1.000000-17.30000X(A2,B1)0.00000020.20000X(B3,C4)1.000000-15.70000X(B3,C1)0.00000030.20000X(B1,C2)0.000000-18.40000X(B1,C1)0.0000000.2000000X(C4,D5)1.000
11、000-13.80000X(C4,D1)0.00000050.20000X(C2,D3)0.000000-17.30000X(C2,D1)0.00000020.20000X(C1,D2)0.000000-18.40000X(C1,D1)0.0000000.2000000X(D5,E1)1.000000-12.20000X(D5,E6)0.00000070.2