Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

TSP Wirh Matlab

Download as pdf or txt
Download as pdf or txt
You are on page 1of 9

Int. J. Adv. Appl. Math. and Mech.

2(3) (2015) 258 - 266

(ISSN: 2347-2529)

Journal homepage: www.ijaamm.com

International Journal of Advances in Applied Mathematics and Mechanics

Travelling salesman problem with MATLAB programming


Article category

Kirtiwant P Ghadle , Yogesh M Muley


Department of Mathematics, Dr. Babasaheb Ambedkar Marathwada University,Aurangabad - 431004, India.

Received 03 February 2015; accepted (in revised version) 04 March 2015

Abstract:

Travelling Salesman Problem (TSP) is very similar to the Assignment Problem (AP) except there is an additional
restriction i.e. Ci j = , if i = j . In TSP a salesman has to visit n cities. The objective is to select the sequence in
which the cities are visited in such a way that total travelling time is minimized; many times AP does not satisfy
the additional restriction, which becomes difficult to solve. In this article we present modified ones algorithm
with MATLAB programming to solve TSP.

MSC:
Keywords:

90B06 90C05 90C08 90C90


ROA method Assignment Problem TSP MATLAB code

c 2015 IJAAMM all rights reserved.

1.

Introduction

MATLAB is powerful software package. The name MATLAB stands for MATrix LABoratory, it deals with matrix
(array). MATLAB can be used for math computations, modeling and simulations, data analysis and processing,
visualization and graphics, and algorithm development, and has many built in tools for solving problems. It is also
possible to write programs in MATLAB, which are essentially groups of commands that are executed sequentially
[1, 2].
MATLAB is widely used in universities and college introductory and advanced courses in applied mathematics,
science and in engineering. Such as, Keng introduced the boundary element method for solving Laplaces equation
using MATLAB [3].Dunn and Harman presented the Graphics User Interface (GUI) features of MATLAB in designing
programs for calculus [4].Colgan reported the successful use of MATLAB in the teaching of core mathematics [5].
Also, Tonkes et al. managed to overcome several inherent difficulties involved in using MATLAB for numerical
analysis to develop an innovative learning model for computational mathematics [6] .
Travelling Salesman Problem is well known in operation research for minimized travelling cost/ distance. Some
of linear programming concept used with MATLAB, YIN ZANG has described implementation of a primal dual
infeasible - interior point algorithm for large scale linear programming under the MATLAB environment [7].
STURM has shown how to solve optimization problem with linear, quadratic and semidefiniteness constraints
using MATLAB toolbox [8].
In this article, we shall illustrate TSP with MATLAB coding by using revised ones assignment problem.
Corresponding author.
E-mail address: drkp.ghadle@gmail.com

258

Kirtiwant Ghadle, Yogesh Muley / Int. J. Adv. Appl. Math. and Mech. 2(3) (2015) 258 - 266

2.

Mathematical Formulation of Traveling Salesman Problem (TSP)[9]

Let 1, 2, . . . , n be the labels of the n cities and C = Ci , j be an n n cost matrix where Ci , j denotes the cost of traveling from city i to city j. Then, the general formulation of the traveling salesman problem (TSP), as described by
Assignment Problem, is shown below.

If Ci , j = C j ,i , the problem is called symmetric traveling salesman problem (STSP).

3.

Assignment based formulation [9]

Starting from his home, a salesman wishes to visit each of (n 1) other cities and return home at minimal cost. He
must visit each city exactly once and it costs Ci j to travel from city i to city j. We may be tempted to formulate his
problem as the assignment problem:

1; if he goes form city i to city j.
xi j =
0; otherwise
Then the mathematical formulation of the assignment problem is,
n
X

M i n i m i z e Z (X ) =

ci j xi j

(1)

i =1

Subject to constraints,
n
X

xi j = 1, f o r i = 1..n

j =1
n
X

xi j = 1, f o r j = 1..n

i =1

xi j 0, f o r a l l i , j .
and xi j = 0 o r 1, f o r a l l i , j .

4.

Algorithm-Revised ones assignment method (ROA) [9, 10]


Divide each row by

n
X

ai j .

i =1

Divide each column by

n
X

bi j .

j =1

Check optimality (draw lines to covers all one) if it is equal to n then go to last step.
If drawn lines < n , then consider all element less than 1.5 in matrix as 1, and again check optimality.
If still drawn lines < n , then choose smallest element d i from uncovered row or column and divide this d i
from uncovered row or column.
Repeat second last step until optimum solution is attained.

259

Travelling salesman problem with MATLAB programming

5.

MATLAB programming concept [1, 2]

Some basic commands in MATLAB Programming for Assignment Problem and TSP as follows,
1. nnmaxr: Find maximum of each row.
EX.1. >> x= randi(3,3)
x=
3 1 3
3 3 3
2 3 3
>>max=nnmaxr(x)
max =
3
3
3
2. nnminr: Find minimum of each row.
3. Size: [M, N] = size(X) for matrix X, returns the number of rows and columns in X as separate output variables.
>> [M, N] =size(x) (see Ex.1)
M=3
N=3
4. length: length(X) returns the length of vector X. It is equivalent to MAX (SIZE(X)) for non-empty arrays and 0
for empty ones.
>> l=length(x) (see Ex.1)
l=3
5. find: [I,J] = find(X,...) returns the row and column indices instead of linear indices into X. This syntax is especially useful when working with sparse matrices. If X is an N-dimensional array where N > 2, then J is a linear
index over the N-1 trailing dimensions of X.
>> x= randi(3,3)
x=
2 2 1
1 3 2
1 3 2
>> [r c] = find(x==1)
r=
2
3
1
c=
1
1
3
6. sum: S = sum(X) is the sum of the elements of the vector X. If X is a matrix, S is a row vector with the sum over
each column.
>> x= randi(3,3)
x=
1 3 1
2 1 2
1 2 1
>> s=sum(x)
s=4 6 4
>> s1=sum(s)
s1 = 14
>> s= sum(sum(x))
s = 14
260

Kirtiwant Ghadle, Yogesh Muley / Int. J. Adv. Appl. Math. and Mech. 2(3) (2015) 258 - 266

7. for loop:
Repeat statements a specific number of times. The general form of a for statement is:
for variable = expr, statement, ..., statement END
>> x= randi(3,3)
x=
2 2 3
3 1 1
3 1 3
>> [M, N] = size(x)
M=3
N=3
multiply each matrix element by 2
>> for r=1:M
for c=1:N
x(r,c)=x(r,c)*2;
end
end
>> x
x=
4 4 6
6 2 2
6 2 6
8. If Condition:
If conditionally execute statements.
The general form of the if statement is
if expression
statements
ELSEIF expression
statements
ELSE
statements
END
9. Tic/toc: tic is used before program and toc is used after program, it is used for to find out total execution time.

6.

Applications [11]

The Ashtavinayaka Yatra is the most important worship yatra in Maharashtra. Ashtavinayaka means "Eight
Ganeshas" when translated in Sanskrit. Lord Ganesha is the Hindu God of unity, remover of obstacles, learning and
prosperity.
This yatra covers the eight holy temples of Lord Ganesha which are located around Pune . Each temple has its
own individuality and mythology. The deities in each temple are distinct from each other. Some of these deities
are described as Swayambhu figurines. This is to point to that no human created these deities; they were formed
completely by nature.
According to Hindu Shastra, you will have to visit the Moreshwar temple at Morgaon first and then in sequence visit
the temples at Siddhatek, Pali, Mahad, Thevur, Lenyandri, Ozar, and Ranjangaon and in last visit the Moreshwar
temple again. This will complete your Ashtavinayak-Yatra.

261

Travelling salesman problem with MATLAB programming

By using Revised Ones Assignment Method (ROA) and MATLAB coding (see Appendix 1) the final solution obtained
was,

262

Kirtiwant Ghadle, Yogesh Muley / Int. J. Adv. Appl. Math. and Mech. 2(3) (2015) 258 - 266

The solution consist in three cycle which are (1,8),(8,4),(4,9),(9,7),(7,1) ; (2,5),(5,2) ; (3,6),(6,3) and minimum
distance is 675 km. if graph command is ignored for TSP 0.085230 seconds required to execute in MATLAB ( Intel
core i3, 4 GB RAM, 32 bit Windows 7 system).
If we again find next minimum values from column 7th, 5th and 6th , which will cover whole journey of Ashtavinayak
Yatra in 987 km.

7.

Conclusions

In this article, a new approach is used for solving TSP with MATLAB Programming. Due to additional restriction in TSP it is difficult to solve by using AP, even though online Assignment Problem (Hungarian Method http://d3meta.com/hungarian/HungarianMethodSolver.html) solver fail to solve this problem. But using some
tricks and ROA method in MATLAB optimum solution has found within 0.085230 seconds.

263

Travelling salesman problem with MATLAB programming

References
[1] Amous Gilat,MATLAB: An Introduction with Applications, John Wiley and Sons, Inc.2004.
[2] Stormy Attaway, MATLAB: A Practical Introduction to Programming and Problem Solving, Elsevier, Inc 2009.
[3] Keng - Cheng Ang, Introducing the boundary element method with MATLAB, Int. J. Math. Edu. Sci. and Technol.
39(4) (2008) 505-519.
[4] P. Dunn and C. Harman, Calculus demonstrations using Matlab, Int. J. Math. Edu. Sci. and Technol. 33 (2002)
584-596.
[5] L. Colgan, Matlab in first - year engineering mathematics, Int. J. Math. Edu. Sci. and Technol. 31 (2000) 15-25.
[6] E. Tonkes, B. Loch, and A. Stac, An innovative learning model for computation in first year mathematics, Int. J.
Math. Edu. Sci. and Technol. 36 (2005) 751-759.
[7] Yin Zang, Solving large - scale linear programs by interior point methods under the MATLAB environment,
Optimization Method and Software 10 (1998) 1-31.
[8] Jos. F. Sturm, Using SeDuMi 1.02, a MATLAB toolbox for optimization over symmetric cones, Optimization
Method and Software 11 & 12 (1999) 625-653.
[9] Ghadle Kirtiwant P, Muley Yogesh M , Revised Ones Assignment Method for Solving Assignment Problem, Journal of Statistics and Mathematics, Volume 4, Issue 1 (2013) 147-150.
[10] Hadi Basirzadeh, Ones Assignment Method for solving Assignment Problems, Applied Mathematical Sciences
6(47) (2012) 2345-2355.
[11] Ghadle Kirtiwant P, Muley Yogesh M , An Application of Assignment Problem in Traveling Salesman Problem
(TSP), Int. Journal of Engineering Research and Applications 4(1) (2014) 169-172.

Appendix : TSP with MATLAB programming


Main Program:
Program for Travelling Salesman Problem by using Revised Ones Assignment method.
Author :Kirtiwant P Ghadle, Yogesh M Muley
clc;
tic;
x=[ 1000 225 304 236 213 339 187 197 226
225 1000 140 153 15 175 84 160 110
304 140 1000 152 132 41 121 190 108
236 153 152 1000 143 188 70 73 63
213 15 132 143 1000 166 74 145 102
339 175 41 188 166 1000 157 226 144
187 84 121 70 74 157 1000 81 43
197 160 190 73 145 226 81 1000 90
226 110 108 63 102 144 43 90 1000]
xnv=x;
x1=x;
minr= nnminr(x);
[r c ] = s i z e (x );
for i=1:r
for j=1:c
x(i,j)=x(i,j)/minr(i);
end
end
x;
minc=nnminr(x);
for i=1:r
for j=1:c
x(j,i)=x(j,i)/minc(i);
end
end
x;
count=0;
for l=1:r
for i=1:r
for j=1:c
264

Kirtiwant Ghadle, Yogesh Muley / Int. J. Adv. Appl. Math. and Mech. 2(3) (2015) 258 - 266

if x(i,j)==1
count = count+ 1;
end
end
if count <= 1
for I=i
for j=1:c
if x(I,j)==1
x(I,j)= x(I,j)*10;
end
for i=1:r
[y c y r ] = f i nd (x 0 == 10);
l = length(yc);
for m =1:l
x(yr(m),yc(m)) = 10;
for j=1:c
if x(yr(m),j) =10
x(yr(m),j)=0;
end
if x(j,yc(m)) =10
x(j,yc(m))=0;
end
end
continue
end
x;
end
x;
end
end
end
count=0;
end
end
x;
for i= 1:r
for j=1:c
if x(i,j) =0
x(i,j)=1;
end
end
end
x;
count=0;
for i=1:r
for j=1:c
if x(i,j)==1
count = count+ 1;
end
end
if count > 1
for I=i
for j=1:c
x(I,j)= x(I,j)*10;
[y c y r ] = f i nd (x 0 == 10, 1);
l = length(yc);
for m =1:l
x(yr(m),yc(m)) = 10;
for j=1:c
if x(yr(m),j) =10
x(yr(m),j)=0;
end

265

Travelling salesman problem with MATLAB programming

if x(j,yc(m)) =10
x(j,yc(m))=0;
end
end
continue
end
x;
end
x;
end
end
count=0;
end
x;
for i= 1:r
for j=1:c
if x(i,j) =0
x(i,j)=1;
end
end
end
x;
z=xnv.*x
Minvalue=(sum(sum(z)))
plot(z)
toc;

266

You might also like