An Alternative Proof of Primitivity of Indecomposable Nonnegative Matrices with a Positive Trace

Similar documents
MATH 5720: Unconstrained Optimization Hung Phan, UMass Lowell September 13, 2018

Dedicated to the late Professors David Hawkins and Hukukane Nikaido

Assignment 1: From the Definition of Convexity to Helley Theorem

In particular, if A is a square matrix and λ is one of its eigenvalues, then we can find a non-zero column vector X with

Factor Price Equalization: Geometrical Conditions

Convex Analysis and Economic Theory AY Elementary properties of convex functions

2 Two-Point Boundary Value Problems

Math 341: Convex Geometry. Xi Chen

Key words. Complementarity set, Lyapunov rank, Bishop-Phelps cone, Irreducible cone

Lebesgue Measure on R n

Functional Analysis HW #3

We describe the generalization of Hazan s algorithm for symmetric programming

EC 521 MATHEMATICAL METHODS FOR ECONOMICS. Lecture 1: Preliminaries

Analysis-3 lecture schemes

Solving Linear Systems

Part III. 10 Topological Space Basics. Topological Spaces

AN ELEMENTARY PROOF OF THE SPECTRAL RADIUS FORMULA FOR MATRICES

THE INVERSE FUNCTION THEOREM

RESEARCH ARTICLE. An extension of the polytope of doubly stochastic matrices

6.254 : Game Theory with Engineering Applications Lecture 7: Supermodular Games

Chapter 1. Preliminaries

Convex Analysis and Economic Theory Winter 2018

Lecture 6 - Convex Sets

Topological properties

Convex Feasibility Problems

A Generalization of Theorems on Inverse-Positive Matrices

(convex combination!). Use convexity of f and multiply by the common denominator to get. Interchanging the role of x and y, we obtain that f is ( 2M ε

Finite dimensional topological vector spaces

Notes for Functional Analysis

Climbing an Infinite Ladder

Integral Jensen inequality

Hilbert Spaces. Contents

Real Symmetric Matrices and Semidefinite Programming

Set, functions and Euclidean space. Seungjin Han

Solving the MWT. Recall the ILP for the MWT. We can obtain a solution to the MWT problem by solving the following ILP:

CHAPTER 2: CONVEX SETS AND CONCAVE FUNCTIONS. W. Erwin Diewert January 31, 2008.

FUNCTIONAL ANALYSIS LECTURE NOTES: COMPACT SETS AND FINITE-DIMENSIONAL SPACES. 1. Compact Sets

SCHUR IDEALS AND HOMOMORPHISMS OF THE SEMIDEFINITE CONE

Chapter 2: Preliminaries and elements of convex analysis

A Parametric Simplex Algorithm for Linear Vector Optimization Problems

WEAK CONVERGENCE OF RESOLVENTS OF MAXIMAL MONOTONE OPERATORS AND MOSCO CONVERGENCE

A Distributed Newton Method for Network Utility Maximization, II: Convergence

Lebesgue Measure on R n

Vector Spaces. Addition : R n R n R n Scalar multiplication : R R n R n.

Math 209B Homework 2

Symmetric Matrices and Eigendecomposition

Immerse Metric Space Homework

On self-concordant barriers for generalized power cones

Convex Analysis and Economic Theory Winter 2018

Boundary Behavior of Excess Demand Functions without the Strong Monotonicity Assumption

NATIONAL BOARD FOR HIGHER MATHEMATICS. Research Scholarships Screening Test. Saturday, February 2, Time Allowed: Two Hours Maximum Marks: 40

i=1 β i,i.e. = β 1 x β x β 1 1 xβ d

Extreme points of compact convex sets

Appendix PRELIMINARIES 1. THEOREMS OF ALTERNATIVES FOR SYSTEMS OF LINEAR CONSTRAINTS

A LITTLE REAL ANALYSIS AND TOPOLOGY

DS-GA 1002 Lecture notes 0 Fall Linear Algebra. These notes provide a review of basic concepts in linear algebra.

Locally convex spaces, the hyperplane separation theorem, and the Krein-Milman theorem

An interior-point gradient method for large-scale totally nonnegative least squares problems

OPTIMISATION 3: NOTES ON THE SIMPLEX ALGORITHM

5 Flows and cuts in digraphs

Continuity of convex functions in normed spaces

Primal-Dual Interior-Point Methods for Linear Programming based on Newton s Method

CHAPTER 7. Connectedness

MATH 223A NOTES 2011 LIE ALGEBRAS 35

Local strong convexity and local Lipschitz continuity of the gradient of convex functions

Linear Algebra: Characteristic Value Problem

Separation of convex polyhedral sets with uncertain data

(c) For each α R \ {0}, the mapping x αx is a homeomorphism of X.

CS 6820 Fall 2014 Lectures, October 3-20, 2014

That is, there is an element

On pseudomonotone variational inequalities

Finite dimensional topological vector spaces

Variational Inequalities. Anna Nagurney Isenberg School of Management University of Massachusetts Amherst, MA 01003

Introduction to Kleene Algebras

MATH FINAL EXAM REVIEW HINTS

A NICE PROOF OF FARKAS LEMMA

Solving Linear Systems

1 Topology Definition of a topology Basis (Base) of a topology The subspace topology & the product topology on X Y 3

Topology. Xiaolong Han. Department of Mathematics, California State University, Northridge, CA 91330, USA address:

Math Camp Lecture 4: Linear Algebra. Xiao Yu Wang. Aug 2010 MIT. Xiao Yu Wang (MIT) Math Camp /10 1 / 88

Climbing an Infinite Ladder

Markov Chains and Stochastic Sampling

Convex Geometry. Carsten Schütt

Unbounded Convex Semialgebraic Sets as Spectrahedral Shadows

Notes taken by Graham Taylor. January 22, 2005

Metric Spaces and Topology

Lecture 5. Ch. 5, Norms for vectors and matrices. Norms for vectors and matrices Why?

CO 250 Final Exam Guide

Notions such as convergent sequence and Cauchy sequence make sense for any metric space. Convergent Sequences are Cauchy

Constructive Proof of the Fan-Glicksberg Fixed Point Theorem for Sequentially Locally Non-constant Multi-functions in a Locally Convex Space

Cover Page. The handle holds various files of this Leiden University dissertation

PM functions, their characteristic intervals and iterative roots

Mathematical Induction. Section 5.1

Introduction to Mathematical Programming IE406. Lecture 3. Dr. Ted Ralphs

Cowles Foundation for Research in Economics at Yale University

Laver Tables A Direct Approach

Exercise Solutions to Functional Analysis

Lecture 2 - Unconstrained Optimization Definition[Global Minimum and Maximum]Let f : S R be defined on a set S R n. Then

LP Relaxations of Mixed Integer Programs

Zangwill s Global Convergence Theorem

Transcription:

An Alternative Proof of Primitivity of Indecomposable Nonnegative Matrices with a Positive Trace Takao Fujimoto Abstract. This research memorandum is aimed at presenting an alternative proof to a well known theorem that an indecomposable nonnegative matrix with at least one positive diagonal element is primitive. Our proof is half topological and half combinatorial. Keywords: Indecomposability, Positive Diagonal Elements, Primitive Matrices 1 Introduction It is well known that a nonnegative matrix which is indecomposable and has at least one positive diagonal element is primitive. This is obvious from another well known characterization of indecomposable imprimitive matrices explained, e.g., in Nikaido[4] (and also [5]). This proposition may date back to A. Brauer. A direct and combinatorial proof is given in a forthcoming book by Krause[1], even in a nonlinear setting. However, it seems difficult to extend the method to spaces of an infinite dimension. In this note, a direct and simple proof is presented, which can be used in nonlinear generalizations as well as for infinite dimensional spaces. (See for the importance of primitivity [2] and [3].) Section 2 explains our definitions and notation, and then in Section 3 the theorem is proved in our own way with a help of two lemmata. In Section 4 some remarks are given when generalizing the theorem to nonlinear mappings and to spaces of an infinite dimension. 2 Notation and definitions Let R n be the Euclidean space of dimension n (n 2), R n + be the nonnegative orthant of R n,anda agivenn n nonnegative matrix. Considering future extensions, we denote the mapping Ax, x R n, by f(x) Ax, andinacomponentwiseway, 1

f i (x) (Ax) i nx a ij x j, where a ij is the (i, j) element of the matrix A, and x j is the j-th entry of the vector x. The symbol N stands for the index set j=1 N {1, 2,...,n}. For vector comparison, we use the following inequality signs. x y iff x i y i for all i N, x > y iff x y and x 6= y, x À y iff x i >y i for all i N. Now we define Definition 1. A nonnegative matrix A is decomposable iff there exist two nonempty subsets I and J of the index set N such that I J =, I J = N, and a ij =0 for i I and j J. Definition 2. A nonnegative matrix A is indecomposable iff it is not decomposable. Alternatively we may define the indecomposability as Definition 2a. A nonnegative matrix A is indecomposable iff for any two nonempty subsets I and J of the index set N such that we have I J =, I J = N, f i (x) >f i (y) in at least one i I for arbitrary two vectors x and y such that x i = y i for i I and x i >y i for j J. Definition 3. A nonnegative matrix A is primitive when there exists a positive integer k such that A k À 0. Otherwise it is imprimitive. 2

3 Propositions Let us consider the following 2 2 matrix A = µ 0 1 1 0 This is indecomposable. When we make the second power of A, A 2 = µ 1 0 0 1 which is decomposable. Thus, we need an additional condition to guarantee the indecomposability of the power of an indecomposable matrix. We first state a lemma. Lemma 1. If A is indecomposable and has at least one positive diagonal element, then A k is also indecomposable for any integer k>0. Proof. We prove this lemma using mathematical induction. First for k =1, the proposition is obvious. We assume this lemma holds for (k 1) with k 2. Let us also assume that f 1 (x) >f 1 (y) for any x and y such that x 1 >y 1 and x i = y i when i>1. (That is, a 11 > 0.) Suppose to the contrary that A k be decomposable and we have two nonempty subsets I and J of the index set N such that., I J =, I J = N, and a (k) ij =0 for i I and j J. Here, a (k) ij means the (i, j) element of A k. Among possible bipartitions between I and J, we adopt as I the one with the minimum number of elements. When 1 I, wesplit A k as A A k 1. Since a 11 > 0, a (k 1) 1j =0 for j J. Now if there is a positive entry among a (k 1) ij for i I and j J, say a h` > 0 for some ` J when h H & I, wewill consider the index set I H and J H. (The index subset H is assumed to be maximal in the sense that a (k 1) ij =0 for any pair (i, j) such that i (I H) and j J.) It is clear that a ih =0for any pair (i, h) such that i I and j H: otherwise we would have a (k) ij > 0 for that i I and some j J since A k = A A k 1, thus yielding a contradiction. Then, from another relationship that A k = A k 1 A, wehave a (k) ih =0 for any pair (i, j) such that i (I H) and h H. This implies that we have two nonempty subsets I 0 (I H) and J 0 (J H) of the index set N such that I 0 J 0 =, I 0 J 0 = N, and a (k) ij =0 for i I 0 and j J 0, is shown to be indecom- a contradiction to the minimality of the index set I. Hence,A k posable when 1 I. 3

When 1 J, wefirst split A k as A k 1 A, and in the second stage as A A k 1,looking for positive entries column-wise rather than row-wise. The above proof can proceed mutatis mutandis. Lemma 2. If A is indecomposable, Ax = λx with λ > 0 implies x À 0, i.e., strictly positive. Proof. Thisiswellknownandevidentfromthedefinition. Now we prove our main proposition. Theorem 1. A given nonnegative matrix A is primitive iff its power A k is indecomposable for any positive integer k. Proof. First, let us prove the only if part. Suppose that there is a positive integer p such that A p becomes decomposable. Then we can find a vector x with at least one zero element such that for any positive integer q, A pq x keeps zero elements in the same positions. This is a contradiction to the fact that A k À 0 for some integer k>0, and so A r À 0 for any integer r k. Now we proceed to the if part. Let us adopt the absolute sum norm k k in the n-dimensional Euclidean space R n, and we consider the subset S {x x R n + and kxk =1}. All we have to show is that starting from an arbitrary point x on the boundary of S, bd(s), a sequence of vectors X {x, Ax, A 2 x,...} comes into the interior of S, int(s). Suppose to the contrary that this sequence remains on the boundary for ever starting from a certain vector x. Indecomposability by Definition 2a requires that in the vector series X, there can be no two points on the same sub-simplex facet of the boundary bd(s). Otherwise x and A k x, for some integer k>0, has the same sign pattern, contradicting the indecomposability of A k, assured by Lemma 1. There are, however, only a finite number of facets on bd(s). Hence, the vector sequence X has to enter the interior of S. Corollary 1. If a given indecomposable matrix A has its trace positive, then it is primitive. Proof. It is clear from Theorem 1 and Lemma 1. 4 Remarks When we generalize the above propositions to nonlinear cases, it does not seem to cause much trouble. It is possible to redefine necessary concepts using f i (x) with the help of a basic assumption of monotonicity, i.e., f(x) f(y) if x y. Some difficulties turn up if we try to extend the results to spaces of an infinite dimension. We have in mind the space of piecewise continuous functions on a compact set, or the space L on a certain compact measurable space. To realize a finite number of facets 4

appearing in the proof of Theorem 1, we need to quantize the compact set by use of a neighborhood of a positive measure References [1] Krause, U.: Positive Discrete Dynamical Systems. Theory, Models, and Applications. (to appear in 2004). [2] Fujimoto, T., Krause, U.: Strong ergodicity for strictly increasing nonlinear operators. Linear Algebra and Its Applications 71, 101-112(1984) [3] Fujimoto, T., Krause, U.: Stable inhomogeneous iterations of nonlinear positive operators on Banach spaces. SIAM Journal on Mathematical Analysis 25, 1195-1202(1994) [4] Nikaido, H.: Convex Structures and Economic Theory. New York: Academic Press, 1968 [5] Nikaido, H.: Introduction to Sets and Mappings in Modern Economics. New York: Academic Press, 1970 5