资源描述:
《运筹学实验二 图论、动态规划求解》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、西华大学能源与环境学院学生上机实验报告西华大学上机实验报告课程名称:运筹学年级/专业:2009/水利水电工程实验成绩:指导教师:施浩然姓名:丁冬冬实验日期:2011年11月实验名称:图论、动态规划求解学号:312009080801417实验学时:3一、实验目的掌握网络图的计算机输入,求解最小树、最短路、最大流问题二、实验内容或设计思想首先将欲求的网络用计算机语言表达,再用lingo计算软件求出模型问题的解。最小树应用破圈方法求解最大流算法用找流量可增链的方法求解网络最短路运用了动态规划,函数迭代、矩阵运算等的基本原理求解三、实验环境与工具计算机,li
2、ngo软件,运筹学软件。四、实验过程或实验数据1、最小数问题给定6个点,各点之间的距离如下程序所示。试用LINGO软件解网络最短路。model:sets:city/1..5/:u;link(city,city):dist,x;endsetsdata:dist=04614999999406876608999999148806999999799999960;enddataN=@size(city);min=@sum(link:dist*x);@for(city(k)
3、k#gt#1:@sum(city(i)
4、i#ne#k:x(i,k))=1;@for(ci
5、ty(j)
6、j#gt#1#and#j#ne#k:U(j)>=U(k)+x(k,j)-(n-2)*(1-x(k,j))+(n-3)*x(j,k);););@sum(city(j)
7、j#gt#1:x(1,j))>=1;@for(link:@bin(x););@for(city(k)
8、k#gt#1:@bnd(1,U(k),999999);U(k)<=n-1-(n-2)*x(1,k););end程序运行结果为:Globaloptimalsolutionfound.Objectivevalue:23.00000Extendedsolversteps:08西华大
9、学能源与环境学院学生上机实验报告Totalsolveriterations:8VariableValueReducedCostN5.0000000.000000U(1)0.0000000.000000U(2)1.0000000.000000U(3)1.0000000.000000U(4)3.0000000.000000U(5)2.0000000.000000DIST(1,1)0.0000000.000000DIST(1,2)4.0000000.000000DIST(1,3)6.0000000.000000DIST(1,4)14.000000.0000
10、00DIST(1,5)999999.00.000000DIST(2,1)4.0000000.000000DIST(2,2)0.0000000.000000DIST(2,3)6.0000000.000000DIST(2,4)8.0000000.000000DIST(2,5)7.0000000.000000DIST(3,1)6.0000000.000000DIST(3,2)6.0000000.000000DIST(3,3)0.0000000.000000DIST(3,4)8.0000000.000000DIST(3,5)999999.00.000000D
11、IST(4,1)14.000000.000000DIST(4,2)8.0000000.000000DIST(4,3)8.0000000.000000DIST(4,4)0.0000000.000000DIST(4,5)6.0000000.000000DIST(5,1)999999.00.000000DIST(5,2)7.0000000.000000DIST(5,3)999999.00.000000DIST(5,4)6.0000000.000000DIST(5,5)0.0000000.000000X(1,1)0.0000000.000000X(1,2)1
12、.0000004.000000X(1,3)1.0000006.000000X(1,4)0.00000014.00000X(1,5)0.000000999999.0X(2,1)0.0000004.000000X(2,2)0.0000000.000000X(2,3)0.0000006.000000X(2,4)0.0000008.000000X(2,5)1.0000007.000000X(3,1)0.0000006.000000X(3,2)0.0000006.000000X(3,3)0.0000000.000000X(3,4)0.0000008.00000
13、0X(3,5)0.000000999999.0X(4,1)0.00000014.00000X(4,2)0.0