旅行商问题数学建模.docx
《旅行商问题数学建模.docx》由会员分享，可在线阅读，更多相关《旅行商问题数学建模.docx（16页珍藏版）》请在冰豆网上搜索。
旅行商问题数学建模
黄石理工学院
数学建模 大型作业
2011—２0１2学年第1学期
一.摘要
二．旅行问题
1.问题描述
2.符号说明
3.模型设计
4.建模求解
5.模型分析
6.
三．建模过程及心得体会
四．参考文件
一.摘要
本文是一个围绕旅行商问题和背包问题这两个经典问题的论文。
问题一,是一个依赖与每个城市去一次且仅去一次的路线确定问题，问题二类似于问题一。
问题三是一个依赖于可背重量限制的背包问题。
关键词：
HAMＩLＴＯN回路LＩＮGＯ最优旅行路线0－１模型
二.旅行问题
问题描述
某人要在假期内从城市A出发,乘火车或飞机到城市B，C，D,E,F旅游购物。
他计划走遍这些城市各一次且仅一次,最后返回城市A。
已知城市间的路费数据见附表1，请你设计一条旅行路线使得他的总路费最少。
如果临行他因故只能去4个城市，该怎样修订旅行路线?
在城市间旅游时他计划购买照相机,衣服,书籍，摄像机，渔具，白酒，食品,而受航空行李重量的限制以及个人体力所限,所买物品的总重量不能超过15kｇ，各种物品的价格见附表2．请你为他决策买哪些物品,使所买物品价值最大。
附表1:
A
Ｂ
C
D
E
F
A
1０
１5０
B
12
00
C
25
8０
D
3３
１85
E
21
４20
Ｆ
５
42０
０
附表２
照相机
衣服
书籍
摄像机
渔具
白酒
食品
香烟
重量（kg）
1
2
4
3
4
２
2
1
价格（元）
435
60０
模型设计
首先给出一个定义:
设v1，v2,......,vn是图G中的n个顶点,若有一条从某一顶点v1出发，经过各节点一次且仅一次，最后返回出发点ｖ1的回路，则称此回路为ＨＡMILTON回路。
问题1.
分析:
这个优化问题的目标是使旅行的总费用最少,要做的决策是如何设定旅行路线,决策受的约束条件:
每个城市都必须去，但仅能去一次。
按题目所给，将决定变量,目标函数和约束条件用数学符号及式子表示出来，就可得一下模型。
模型建立:
对于6个城市的旅行问题设A，B,C,D,E,F六个城市分别对应v1,v２,v3,ｖ4,v5，v6。
假设
表示从城市i到城市j的费用。
定义0-1整数型变量
=1表示从城市i旅行到城市ｊ，否则
=0。
则旅行问题的数学模型可表示为一个整数规划问题。
min z=
（i
j）
s.t．
=１ （i
j;j＝1，2,……,６）
=１ （i
j;ｉ＝1,２,……,6）
（i
ｊ;i=2,3,……,6;j=2，３,……6）
其中辅助变量
（ｉ=2,3,……，6）可以是连续变化的,虽然这些变量在最优解中取普通的整数值（从而在约束条件中，可以限定这些变量为整数）。
事实上,在最优解中,
=访问城市的顺序数。
模型求解
运用ＬＩNGO,输入程序：
MODＥL:
!
Travelinｇ SalesPｒｏｂleｍ fｏr ｔｈecitiesoｆsixcity;
SETS:
CITY／1..６/:
U;!
U（I）＝sｅquｅncｅ nｏ.ofcｉty;
ＬINK（CITY,CＩTY）:
COST,!
Ｔhecost matrix;
X; !
Ｘ（Ｉ,J）＝1ifｗeｕselｉnkI， J;
EＮＤSＥTS
DATA:
!
Ｃｏst mａtriｘ，it neeｄnｏｔ ｂe symmｅｔｒｉc;
COST=0 12０ 2５0 33０ 210 150
１20 09８ 3５0 ２2５３0０
２５098 0 520 ４30 280
３３0 ３50 5200 270 18５
210 2２5 ４3０ 27０ ０ ４20
1５０300２80 1８5 ４2０ 0；
ＥＮDDATＡ
N=＠SＩZE（CITY）;
ＭIN=＠SUM（LIＮＫ:
COST＊X）;
@FOR（ＣITＹ（K）：
!
Ｉｔmusｔbe ｅｎtereｄ;
@SUM（ＣＩTＹ（I）| I＃NＥ＃Ｋ:
X（I,Ｋ））=1;
！
It mｕsｔbedｅparｔｅd;
@SＵM（ＣITY（J）|J＃NE#K：
X（K, J））＝1;
!
Weakｆormoｆ thｅsubtouｒbreakinｇcoｎstrainｓ;
！
Theｓeａrenotvｅrypoweｒfulｆoｒlargｅ ｐｒｏblems;
＠FＯR（ＣIＴY（Ｊ）|J#ＧT＃ 1 #AＮD#J#NE#K:
U（Ｊ）>=U（K）+X（K,J）-（N-2）*（1－Ｘ（K,J））+（N-３）*X（J,K））;）；
@FOR（LINK：
@ＢIＮ（X））;!
Maｋｅｔｈe X's0／1;
!
Foｒthｅｆirｓtandlaststoｐ ｗｅknow...；
@FＯR（CＩTY（K）｜ K#GT＃1:
U（K） <= N- １-（N-2） *X（1，K）;
Ｕ（K）>=１＋（ N-2）＊ X（ K, 1））；
EＮD
得到结果：
总费用为１163
路线：
A－B－C-Ｆ-Ｄ-E-A
问题２.
临时因故只能去4个城市，那么应该怎样安排旅行路线。
分析：
在B，C,D,E,F五个城市中选4个城市，显然有５种可能,按照问题一中的模型,将费用矩阵稍作修改,运用LINGO分别解除５种可能的费用,进行比较，即得结果。
（1）选取Ｂ,Ｄ,E,F，计算旅行费用:
MODEL：
!
Tｒavｅlｉng SalesPｒoblemfortheｃiｔieｓofsｉｘcity;
SETS:
CＩTY/１..5/:
U;!
Ｕ（I）= seｑｕence no.ｏfcｉty;
LIＮK（ＣIＴY, CIＴY）:
COST, !
Tｈeｃosｔmaｔｒｉｘ;
X;！
X（ I, Ｊ）＝1if ｗｅuｓelinｋI,Ｊ;
ENDSETS
ＤAＴA：
！
Coｓtmatriｘ,it nｅｅdｎoｔｂesyｍmeｔric;
CＯST=0 １20 3３02１0 １5０
１200 ３50２25 30０
3３0 350 0 27０ 18５
210 225２700 ４2０
150 300 185 ４2０0;
ENDDAＴA
Ｎ=@SＩZE（ CITY）;
MIＮ= ＠SUM（LＩNK:
ＣOST*X）;
@FOＲ（ＣＩTY（K）：
!
Ｉｔmｕst be ｅnｔered;
@ＳUＭ（ CITY（I）｜I ＃NE# Ｋ：
X（I,K））=1；
!
Ｉtmｕｓtbe departed；
@SＵM（CIＴY（ Ｊ）｜J#NE#K：
X（K, J））=1;
!
Weaｋfoｒｍof ｔhesuｂtour bｒeａkｉｎgｃonｓtrainｓ；
！
Ｔｈeｓearenoｔverypｏwerfｕlforlarｇe proｂleｍs；
@ＦOR（ CＩTY（J）| Ｊ#GT#1#ANＤ#J＃ＮＥ# K:
U（Ｊ）>=U（K）+Ｘ（K,J）-（N－2）＊（１－X（Ｋ,Ｊ））+（N-３）*X（J,Ｋ））;）;
@FOR（ＬINK:
@BＩN（X））;!
MａketｈeX'ｓ0/１；
!
Ｆoｒthefｉrstandlａstｓｔopwｅ ｋnoｗ...;
@FOＲ（ＣITY（K）| K #GT#1:
U（K） <=N - 1-（N -2）*X（ 1,Ｋ）;
U（K） >=1＋（ N - ２）*X（K,1））;
ENＤ
得到结果：
总费用：
950
路线：
A-B－Ｅ-D-F-A
（2）选取B,C，E,F,计算旅行费用：
MODEL：
!
TravelingＳalｅｓPrｏblemforｔhecitiｅsｏｆsix city；
SETS:
CIＴY／ １..5/:
Ｕ;!
U（Ｉ） =ｓequencｅ no. of ciｔy;
LINK（CＩTY, CITY）:
COSＴ,!
Ｔｈecostmatrix;
X;!
X（Ｉ,J）＝ １ifwｅ ｕseliｎk I,Ｊ;
ＥNDSETS
DＡTA:
！
Coｓtmａｔrｉx,iｔ nｅedｎｏtbｅｓymｍetrｉc；
ＣOST=0 120 ２５0 ２10 15０
１20 0 9８ 2２5 ３00
25098 0 ４30 28０
2102254300 ４20
15０３00280 ４２0 0；
ENDDATA
N=@SIZE（CIＴY）;
MIN=@SUM（LＩNK:
ＣOST*X）；
@FOR（CＩTY（ K）:
!
Itｍusｔbe entｅｒｅd;
＠ＳUM（ ＣITＹ（I）|I #NE＃Ｋ:
X（ I, Ｋ））=1;
!
Itｍuｓtbedeｐaｒted;
＠ＳUM（CIＴY（J）|Ｊ#NE＃Ｋ:
X（Ｋ， J）） =1;
!
Ｗeakfｏrm oftｈeｓubtouｒbreakingcｏnstrainｓ;
！
Ｔheｓearenｏｔveryｐowerfulforlａrgeｐroblems;
＠ＦOR（CITY（J）|J #GT#1#AND＃J#NE#K：
Ｕ（Ｊ）＞＝U（K）+Ｘ（Ｋ,J）-（N-2）*（1-X（Ｋ，Ｊ））＋（N-３）*Ｘ（Ｊ，K））；）；
@FＯR（LＩNK:
@BIN（X））; !
MaketheX's0/1;
!
For thefiｒstanｄ laststop weknｏw．..;
@FＯR（CITY（K）| K #GT# 1:
Ｕ（ K）＜=N-1- （Ｎ－2）*X（１，K）；
Ｕ（K）>=１ ＋（N-2） *X（Ｋ,1））;
END
得到结果：
总费用：
９６３
路线:
A-E－B－Ｃ-Ｆ-A
（3）选取B,C,D,F，计算旅行费用:
MODＥＬ：
!
ＴraｖelingSalｅsＰroｂｌeｍforｔhｅcitiesofｓixcity;
ＳEＴS:
CＩTＹ/１..5/:
U; ！
U（Ｉ） =sｅｑuencｅｎｏ．of citｙ;
ＬINK（CIＴＹ, CITＹ）：
ＣOＳT,!
The costｍatｒix;
X;!
X（I,J） =1 ifweusｅｌinkI,Ｊ;
EＮＤＳETＳ
ＤＡTＡ：
!
Costmａｔｒix, itneｅdｎotｂesｙmmetrｉc;
ＣOST＝0 １2０250 ３30 150
12０ 0 98 350 30０
2５09８ ０ ５2０ 280
330 35０ 520 0 185
15０ ３00 280 １85 0 ;
ENＤＤATA
Ｎ= @SＩＺＥ（ CITY）;
ＭIN= @SＵM（ LINＫ:
ＣＯST＊X）;
@ＦOＲ（CITY（ K）:
!
It must bｅenteｒed;
＠SUM（ ＣITY（ Ｉ）|I＃NE＃K:
X（I, Ｋ））＝1;
!
Iｔｍustｂedepａrted;
＠ＳUＭ（ CIＴY（J）｜J#NＥ#K:
X（K,J）） ＝１;
！
Weakformｏｆtｈesubtourbreakingconｓtrains;
！
Theｓe arenｏtverｙpoｗerｆulfoｒlａrgeprobｌems;
@ＦＯR（CITY（J）|J#GT＃1 #ANＤ＃Ｊ#NE# K:
U（J）>=U（K）+X（K,J）-（N-2）*（1-X（K,J））＋（Ｎ-3）＊Ｘ（J,K））;）；
＠FＯＲ（ＬＩNK:
@BIN（X））;！
Makeｔhe X＇s0/1；
!
Ｆoｒ thefirsｔａndlast stopweknow...;
@ＦOR（ CITＹ（Ｋ）| K ＃GT＃ 1：
U（ K）＜= N-1-（N －2） ＊X（1， K）；
U（Ｋ）>=1+（ N －２）* X（ K, 1））;
ENＤ
得到结果:
总费用:
1013
路线：
A－B-C－Ｆ-D-A
（4）选取路线:
Ｂ,C,Ｄ,E,计算旅行费用：
MODEL：
!
ＴraveｌinｇＳales Pｒoblem forthe cｉｔiｅｓofｓiｘ city；
SETS：
CITY／ １.．5／：
U;!
U（ I）= ｓequｅnｃeno.ofciｔy;
LINＫ（ＣITY, ＣIＴＹ）：
ＣOＳＴ， !
The costmａｔriｘ;
X；！
X（I,J）=１if wｅuselinｋ I,J;
EＮＤSETS
DＡTA:
!
Cosｔmatrix,ｉtnｅｅdｎotｂesyｍｍｅtric;
COSＴ=０ 120２５0 330 21０
120０ ９8 350 2２5
250 98 ０ 520 4３０
３30 3５0 520 0 27０
210 ２2５43０270 0 ;
ＥＮDDAＴＡ
N=@ＳIZE（CITY）;
MＩN=＠ＳUM（ＬINK:
ＣＯST *Ｘ）;
＠FＯR（CITＹ（K）：
!
Ｉtmusｔ beｅntered;
@ＳUＭ（CIＴY（I）|I ＃NＥ#K:
X（I, K）） =１；
!
Iｔ mｕstbeｄｅｐarted;
＠ＳＵＭ（ＣＩＴY（J）|J #NE#K:
X（K，J）） =１；
!
Weak fｏrmｏf the ｓubtｏuｒbreakiｎg conｓｔrainｓ；
!
Thesｅarenoｔverｙｐowｅrfulｆorlargｅproｂｌeｍs；
@ＦOＲ（CITY（J）|J #ＧT#１#AND# J#NE＃K：
Ｕ（J）>=U（K）+X（K,Ｊ）-（Ｎ-2）*（1-Ｘ（Ｋ，J））+（N－3）*X（J,Ｋ））;）；
@FＯＲ（ LINK:
@BＩN（X））；!
Maｋetｈｅ X's0/１;
!
Fｏrthｅｆirsｔ ａnｄ lａstsｔoｐweknｏw...;
＠FＯＲ（ CITY（K）｜K#GT#1:
Ｕ（K）<=N-1－（N-2）* Ｘ（1，K）;
U（K） >＝1 +（Ｎ -２）＊X（ K,1））;
END
得到结果：
总费用:
1１７3
路线:
A-C-B-E-Ｄ-Ａ
（5）选取C,D,E,F,计算旅行费用：
ＭODＥＬ:
！
Traｖｅlｉｎg ＳalesPｒoｂlemfoｒ tｈe ｃitiesofsix city;
SETS:
CIＴY／1..5/:
Ｕ;!
U（I） =sequｅnｃenｏ. of cｉty;
LＩNK（ CITY,CITY）：
COST, ！
Ｔhｅcost matｒix;
X;!
Ｘ（I,Ｊ）＝1 ifｗe ｕselｉnｋI,J；
ENDSETS
DAＴA：
!
Costmaｔrｉｘ,itneedｎｏt besymmetric;
CＯST=0 2５0 330 2１01５０
25０0 ５２0 43０ 2８0
330 52０ 0 270 １85
210 430270 0 4２０
150２80 185 ４20 0 ；
EＮＤＤATA
N＝@SIZE（CITY）；
MIN= @SUM（ＬIＮK：
COSＴ* X）;
@FOR（ CＩTY（ K）:
!
Itｍuｓtbe enteｒed;
＠SUM（ＣITY（ I）｜I＃NE#Ｋ:
X（Ｉ,K）） ＝1；
！
Iｔ ｍustbeｄeparted;
@SUM（CＩTY（ Ｊ）｜J#NE# Ｋ:
X（K,Ｊ）） =１;
！
Ｗeakformｏｆ thｅｓubｔｏｕr brｅakinｇconｓtｒaｉnｓ;
!
Thesearｅnoｔverypowｅrfuｌ ｆｏrlａｒgeｐroｂlems;
@FOR（CITY（Ｊ）| Ｊ#ＧＴ#1 #ANＤ#J#NE＃K：
Ｕ（J）>＝U（K）+X（Ｋ,J）－（N-２）*（1-Ｘ（Ｋ,J））+（N－３）*Ｘ（J,K））;）；
＠ＦOR（ ＬINK:
@BIＮ（ X））;!
MakeｔheＸ's 0/１;
！
Ｆorthefirstaｎｄlａｓtｓｔoｐｗeknow...;
＠FＯR（ CITY（K）|K＃GT＃１：
U（K） ＜=N-1－（Ｎ －2）＊X（ 1,K）;
U（ K）＞=１+（N -２）*X（Ｋ,1））;
END
得到结果：
总费用：
１195
路线：
A-Ｃ－F－Ｄ-E-A
因此，总结以上
（1），
（2）,（3），（４）,（5）五种情况可得:
应该选用路线：
A-Ｂ-Ｅ-Ｄ-F-A。
总费用花费最少,为９50.
问题3.
在城市间旅游时他计划购买照相机、衣服、书籍、摄像机、渔具、白酒、食品、香烟等物品，而受航空行李重量的限制以及个人体力所限，所买物品的总重量不能超过１5kg，各种物品的价格见附表2。
请你为他决策买哪些物品，使所买物品总价值最大?
分析:
解读题意,实际上此题涉及背包问题,题目转化为一个限重15公斤的背包，要求在８件可带可不带的物品中做出合理的方法。
模型建立：
将照相机,衣服，书籍,摄像机,渔具，酒,食品和香烟编号为1，2,３，4,５，６，7，８。
设总容量为b,
为每个物品的重量,
为每个物品各自的价格。
定义一个变量
，当
=1是,表示装第i件物品,当
=0时，则不装（ｉ=1,2,……,8）．
则约束条件为:
maxz＝
=0或1,（i=1，2，……,８）
模型求解：
利用LINGO软件求解,在LINGO中输入以下程序：
moｄｅｌ：
sｅｔs:
M/1..8/:
Ｗ；
price/1.．8/:
p;
number/1．.8/:
x;
endsｅts
dａta:
W=1，２,4,3,4,2，2,1;
P=13０0,2７5０，3２0,43５0,1400,300,１20,600;
enddaｔａ
MAX=@SUM（Ｍ：
X*P）;
@SUM（M:
X*W）<=１5;
@FOR（M:
@ＢＩＮ（X））;
Eｎd
得到结果:
选择照相机，衣服，摄像机,渔具，酒,食品和香烟,最大价值为１0820。
三.建模过程及心得体会
数学建模是一个经历观察、思考、归类、抽象与总结的过程，也是一个信息捕捉、筛选、整理的过程,更是一个思想与方法的产生与选择的过程。
它给学生再现了一种“微型科研”的过程。
数学建模教学有利于激发学生学习数学的兴趣,丰富学生数学探索的情感体验;有利于学生自觉检验、巩固所学的数学知识，促进知识的深化、发展；有利于学生体会和感悟数学思想方法。
同时教师自身具备数学模型的构建意识与能力，才能指导和要求学生通过主动思维,自主构建有效的数学模型，从而使数学课堂彰显科学的魅力。
为了使描述更具科学性，逻辑性,客观性和可重复性，人们采用一种普遍认为比较严格的语言来描述各种现象,这种语言就是数学。
使用数学语言描述的事物就称为数学模型。
有时候我们需要做一些实验,但这些实验往往用抽象出来了的数学模型作为实际物体的代替而进行相应的实验,实验本身也是实际操作的一种理论替代。
数学建模要有团队精神。
数学建模不是一个人的事,是团队的一项活动。
三个人要相互支持，相互鼓励。
不能只管自己的那一部分。
特别是模型的建立，一个人根本不可能想的那么全面，只有大家一起讨论才能把问题搞清楚。
必须合理的安排时间。
做任何事情，合理的时间安排非常重要，建模也是一样,事先要做好一个规划，建模一共分十个板块。
你每天昨晚哪几个板块事先要确定好，这样做才会使自己游刃有余，保证在规定的时间内完成论文,以避免在时间上的不妥,以至于最后无法完成论文。
有些同学经常会借鉴一些别人的论文里的思想，然后变成自己的东西，不过那也是一种能力,不能小瞧。
所以也要多看些别人的论文来调高自己。
四．参考文件。
【１】.运筹学,科学出版社,徐玖平，胡知能,２0０3年11月第一版
【２】.运筹学——数据,模型，决策，科学出版社,徐玖平，胡知能,20０6年3月第一版
升级会员 