An Efﬁcient Algorithm for Nontrivial Eigenvectors in Max-Plus Algebra

: The eigenproblem for matrices in max-plus algebra describes the steady state of the system, and therefore it has been intensively studied by many authors. In this paper, we propose an algorithm to compute the eigenvalue and the corresponding eigenvectors of a square matrix in an iterative way. The algorithm is extended to compute the nontrivial eigenvectors for Latin squares in max-plus algebra.


Introduction
Discrete event systems can be used to study processes that are driven by the occurrence of events.The relevant variables represent times at which events are taking place.An important class of systems consists of min-max-plus systems in which the operations-the maximization, the minimization and the addition, can describe discrete event dynamic systems.The union of two sets of equations describes the min-max-plus systems, where the maximization and addition are contained in one set of equations and the minimization and addition are contained in another set of equations.
Max-plus algebra R = R ∪ { } where = −∞, is a set which deals with two main operations-⊕ as maximization and ⊗ as addition.These operations endow R with the structure of a commutative semiring.
There are various applications of max-plus algebra in different areas of mathematics [1].Max-plus algebra is also considered in traffic control, parallel processing systems, telecommunication networks, manufacturing systems, machine scheduling, and control theory [2][3][4].In [5], the author has described the model of all train connections of the whole Dutch railway system with all train types by a max-plus system.Max-plus algebra can be used to model disk events related to synchronization and time delays [6].A new steganography method has been proposed to hide an image into another image using matrix multiplication operations on max-plus algebra [7].In [8] the author has shown how max-plus algebra can be helpful for the dynamic programming of algorithms.The primary reason for the utilization of max-plus algebra in different fields is that many equations which are non-linear in conventional algebra become linear in max-plus algebra [3].
In this paper, we consider the eigenproblems in max-plus algebra.Eigenproblems are the problems of finding an eigenvalue λ and an eigenvector v from an n × n matrix M such that Mv = λv.The algorithm in [9] can solve these problems.
A novel iterative method is proposed for the computation of eigenvalue and the corresponding eigenvectors in max-plus algebra.To evaluate the eigenvalue and the corresponding eigenvectors of discrete event systems, some algorithms have been established in [8][9][10][11].The first algorithm to determine the eigenvalue was given by Karp [12].A new algorithm was proposed to solve systems of max-plus linear inequalities [13].This algorithm uses the technique of cyclic projectors.In [14], there is a method for determining the eigenvectors of monotone and homogenous functions.
The paper is structured as follows.In Section 2, basic results and definitions related to max-plus algebra are described.A new algorithm is proposed to determine the eigenvalue and eigenvectors in an iterative way.In Section 3, we discuss the Latin squares in max-plus algebra and extend the algorithm for Latin squares.Conclusions are made in section 4.

Max-Plus Algebra
In this section, we recall some basic results and definitions related to the max-plus algebra.The symbols (R , ⊕, ⊗) denote the algebraic structure of max-plus algebra where R is the set R ∪ {−∞} and ⊗, ⊕ are the binary operations defined as These operations endow R with a structure of commutative semiring.Both the operations ⊗ and ⊕ satisfy associativity and commutativity.
Also for all m ∈ R So in max-plus algebra, is a null element and 0 is a unit element.For m ∈ R and a non-negative integer n, define for n > 0, 0 for n = 0.
Using the standard multiplication in R, we have In max-plus algebra, R m×n and R m denote the set of all matrices of order m × n with entries in R , and the set of all m × 1 vectors respectively.The set of first n natural numbers is denoted by n, i.e., n = {1, ..., n}.Now let M = (m ij ) and N = (n ij ) be elements of R m×n .The sum of these two matrices Let M ∈ R m×k and N ∈ R k×n be two matrices.The product M ⊗ N = T ∈ R m×n is defined by the formula For M ∈ R n×n , the kth power of M can be defined as The starting point is the following system of equations describing the time evolution of a discrete event system: This can be written as Finally, the previous equation can compactly be written as where For each matrix M ∈ R n×n , we can construct a directed graph (digraph) G(M) = G(N, E), where N is the set of vertices and E represents the set of edges.There exists an arc i → j if m ji = −∞; such an arc will be denoted by (i, j).The weight ω(j, i) for each (i, j)-arc is defined as ω(j, i) = m ji .If m ji = −∞, then there is no (i, j)-arc.A path is a sequence of arcs (a 1 , a 2 ), ..., (a k−1 , a k ).If all vertices a 1 , ..., a k are different then the path is called an elementary path.An elementary path of the form (a 1 , a 2 ), ..., (a k−1 , a 1 ) is known as a circuit.A loop is a circuit consisting of an edge that connects a vertex to itself.The weight of a path p = (a 1 , a 2 ), ..., (a k−1 , a k ), denoted by ω(p), is defined as the sum of weights of each arc.The length of p, denoted by l(p), is equal to the number of arcs in the path p.

The average weight of a path p is defined by ω(p) l(p)
. A critical circuit is a circuit with maximum average weight.If there exists a path from each vertex i to each vertex j then the graph is said to be strongly connected.A matrix M is irreducible if and only if the digraph of M is strongly connected [15].The eigenvalue of a matrix M is equal to the maximal average weight in G(M) [16].The eigenvalue for an irreducible matrix M, if it exists, is unique [16].So, problem of getting the eigenvalue of a matrix M related to the problem of getting the critical circuits in G(M).
In this paper, we find the eigenvalue λ ∈ R of a matrix M ∈ R n×n with the help of the maximal average weight in G(M) as done in [16].Then we define the matrix M λ as Using this relation, an algorithm is derived to calculate the associated eigenvectors of a matrix M. Theorem 1.Let M ∈ R n×n be a matrix and let Proof of Theorem 1.The main task is to prove that M ⊗ v = λ ⊗ v. Indeed, we have The main algorithm of this article stated below (Algorithm 1) gives the eigenvalue and the corresponding eigenvectors.This algorithm consists of the following steps:

Define an arbitrary initial vector u
Each step of the above algorithm is explained with the help of the following example.

Example 1. Consider the following matrix
The corresponding graph G(M) is depicted in Figure 1.Since the average weight of the circuit 1 → 2 → 1 is the maximal circuit mean, the eigenvalue λ = 4. Now we get .
The following sequence is obtained by iterating (3), It follows that s = 2, r = 4.We proceed with the computation of the vector v.
Hence for the eigenvalue λ = 4, the corresponding eigenvector v of the matrix M resulting from the above algorithm is a correct eigenvector.

Latin Square
In this section, the proposed algorithm is used to compute the eigenvalue and the corresponding eigenvectors of Latin squares.A Latin square is a square matrix of size n × n with n different numbers such that each row and each column contains distinct numbers [17].An example of a Latin square is given in the following The purpose behind considering eigenproblems for Latin squares in max-plus algebra is that similar issues have been examined for different matrices like circulant matrices [18,19], Monge matrices [20] and inverse Monge matrices [21].Eigenproblems for these special matrices are more easy to answer.For example, the maximal entry in circulant matrices is its eigenvalue.Now we will discuss two different types of Latin square.
Let L be a Latin square with entries in n.Since each entry in L is finite, the corresponding graph G(L) is strongly connected.Similarly, for a Latin square L with entries in n , exactly one infinite entry appears in each row and column.So there exists a path of length 2 from each vertex r to each vertex s in G(L).Consequently, the graph G(L) is strongly connected.Hence in both cases, a Latin square is an irreducible matrix.In [11] authors show that for a Latin square L, the maximal number in L is equal to the maximal circuit mean in G(L).Therefore the eigenvalue of a Latin square L is equal to the maximal number in L [11].Now using this concept, the proposed algorithm (Algorithm 2) computes the eigenvalue and eigenvectors of Latin square as follows: Algorithm 2 Eigenproblems for Latin Square 1. Compute the eigenvalue λ as the maximal number in L.

Define
Also, recall the power algorithm presented by Subiono [9].This algorithm computes the eigenvalue and corresponding eigenvector for max-plus systems of type (1).
Now we illustrate the Algorithms 2 and 3 by means of an example.In the following example, based on Algorithms 2 and 3, M is a Latin square with entries in n .
After iterating (3), the following sequence is obtained: It follows that s = 2, r = 4.To compute a corresponding eigenvector we proceed with the computation of the vector v.
Hence for the eigenvalue λ = 3, the corresponding eigenvector v of the matrix M resulting from the above algorithm is a correct eigenvector.2. For Algorithm 3, consider the initial vector It follows that q = 2, p = 4.The vector v resulting from Algorithm 3, is given as Remark 1.The main purpose of this article is to present an alternative algorithm for the computation of an eigenvalue and the associated eigenvector in max-plus algebra.This algorithm works in an iterative way.Here we give a computational comparison of the Algorithm 2 with the Algorithm 3 given by Subiono and van dar Woude [9].In the case of a circulant matrix or a Latin square, the maximum value is the eigenvalue.Therefore the Algorithm 2 works quite easily compared to the Algorithm 3 This is because Algorithm 3 ends up in a periodic behavior if there exists a real number c and integers r > s ≥ 0, such that u(r) = c ⊗ u(s), while in the case of Algorithm 2, we do not need any real number c. Also the Algorithm 3 gives an eigenvector v as v = ⊕ r−s j=1 (λ ⊗(r−s−j) ⊗ u(s + j − 1)), however, when using the Algorithm 2 one gets an eigenvector v as v = u * (s) ⊕ ... ⊕ u * (r − 1).Hence the computation of an eigenvector using the Algorithm 2 is quite easy compared to the Algorithm 3.
and a ⊕ b = max(a, b).

Figure 1 .
Figure 1.Graph representation of matrix M.

.
Which is the correct eigenvector.