运筹学基础实验报告旆妫.pdf

《运筹学基础》实验报告 5070209136 朱锋 1. 问题描述与抽象 旅行商问题 某旅行商从节点 1 出发，遍历所有节点后回到节点 1。各节点之间的路径长度如下图所 示，利用 OPL进行建模并求解，找出旅行商的最短路径。 2. 整数规划模型建立 矩阵 C 存储上图中各节点之间的边权， 表示从节点 i 到节点 j 的边权， (M 可根据具体问题赋一个相对极大的数 )表示节点到自身距离认为无限大， 此为无向图， 故 得到对称矩阵 C： 同时定义 : 则旅行商问题要求使遍历路径的边权和最小， 用 z 表示遍历所有城市节点的代价和， 则 其目标函数可表示为 : 旅行商问题限制条件： ① 保证每个节点只有一个入度， ② 保证每个节点只有一个出度 1，2 限制条件保证图中每个城市被通过且只被通过一次。 ③ 保证没有分立的圈存在 查阅资料得， 3 限制条件可以排除经 1，2 条件筛选后不满足题意的几个独立的子环 路的情况，任何一个封闭的子圈必将和 3 限制发生冲突。 ④ 整数规划基本限制条件 基于以上分析，旅行商问题的整数规划数学模型： 3. OPL Development Studio 求解此整数规划问题 旅行商问题的模型文件： /********************************************* * OPL 6.3 Model * Author: Arrow * Creation Date: 2009-12-9 at 下午 03:14:09 *********************************************/ {string} Cities = ...; {string} Uities = ...; int Dis[Cities][Cities] = ...; int Citynumber = ...; dvar int U[Cities]; dvar boolean X[Cities][Cities]; minimize sum( i in Cities , j in Cities ) Dis[i][j]*X[i][j]; subject to{ forall( i in Cities ) sum( j in Cities ) X[i][j] == true; forall( j in Cities ) ctLeaveEachCityOnce: sum( i in Cities ) X[i][j] == true; forall(i in Uities,j in Uities:i!=j) U[i]-U[j]+Citynumber*X[i][j]=Citynumber-1; forall(i in Cities) U[i]=0; } 旅行商问题的数据文件： /********************************************* * OPL 6.3 Data * Author: Arrow * Creation Date: 2009-12-8 at 下午 08:25:04 *******************