CSR expansions of matrix powers in max algebra

We study the behavior of max-algebraic powers of a reducible nonnegative n by n matrix A. We show that for t>3n^2, the powers A^t can be expanded in max-algebraic powers of the form CS^tR, where C and R are extracted from columns and rows of certain Kleene stars and S is diadonally similar to a Boolean matrix. We study the properties of individual terms and show that all terms, for a given t>3n^2, can be found in O(n^4 log n) operations. We show that the powers have a well-defined ultimate behavior, where certain terms are totally or partially suppressed, thus leading to ultimate CS^tR terms and the corresponding ultimate expansion. We apply this expansion to the question whether {A^ty, t>0} is ultimately linear periodic for each starting vector y, showing that this question can be also answered in O(n^4 log n) time. We give examples illustrating our main results.Comment: 25 pages, minor corrections, added 3 illustration

Cyclic classes and attraction cones in max algebra

In max algebra it is well-known that the sequence A^k, with A an irreducible square matrix, becomes periodic at sufficiently large k. This raises a number of questions on the periodic regime of A^k and A^k x, for a given vector x. Also, this leads to the concept of attraction cones in max algebra, by which we mean sets of vectors whose ultimate orbit period does not exceed a given number. This paper shows that some of these questions can be solved by matrix squaring (A,A^2,A^4, ...), analogously to recent findings concerning the orbit period in max-min algebra. Hence the computational complexity of such problems is of the order O(n^3 log n). The main idea is to apply an appropriate diagonal similarity scaling A -> X^{-1}AX, called visualization scaling, and to study the role of cyclic classes of the critical graph. For powers of a visualized matrix in the periodic regime, we observe remarkable symmetry described by circulants and their rectangular generalizations. We exploit this symmetry to derive a concise system of equations for attraction cpne, and we present an algorithm which computes the coefficients of the system.Comment: 38 page

On the asymptotic period of powers of a fuzzy matrix

AbstractIn our prior study, we have examined in depth the notion of an asymptotic period of the power sequence of an n×n fuzzy matrix with max-Archimedean-t-norms, and established a characterization for the power sequence of an n×n fuzzy matrix with an asymptotic period using analytical-decomposition methods. In this paper, by using graph-theoretical tools, we further give an alternative proof for this characterization. With the notion of an asymptotic period using graph-theoretical tools, we additionally show a new characterization for the limit behaviour, and then derive some results for the power sequence of an n×n fuzzy matrix with an asymptotic period

On the Tightness of Bounds for Transients of Weak CSR Expansions and Periodicity Transients of Critical Rows and Columns of Tropical Matrix Powers

We study the transients of matrices in max-plus algebra. Our approach is based on the weak CSR expansion. Using this expansion, the transient can be expressed by $\max\{T_1,T_2\}$, where $T_1$ is the weak CSR threshold and $T_2$ is the time after which the purely pseudoperiodic CSR terms start to dominate in the expansion. Various bounds have been derived for $T_1$ and $T_2$, naturally leading to the question which matrices, if any, attain these bounds. In the present paper we characterize the matrices attaining two particular bounds on $T_1$, which are generalizations of the bounds of Wielandt and Dulmage-Mendelsohn on the indices of non-weighted digraphs. This also leads to a characterization of tightness for the same bounds on the transients of critical rows and columns. The characterizations themselves are generalizations of those for the non-weighted case.Comment: 42 pages, 9 figure

Cumulants, lattice paths, and orthogonal polynomials

A formula expressing free cumulants in terms of the Jacobi parameters of the corresponding orthogonal polynomials is derived. It combines Flajolet's theory of continued fractions and Lagrange inversion. For the converse we discuss Gessel-Viennot theory to express Hankel determinants in terms of various cumulants.Comment: 11 pages, AMS LaTeX, uses pstricks; revised according to referee's suggestions, in particular cut down last section and corrected some wrong attribution