Majorization 📖 Wikipedia

In mathematics, majorization is a preorder on vectors of real numbers. For two such vectors, ${\displaystyle \mathbf {x} ,\ \mathbf {y} \in \mathbb {R} ^{n}}$, we say that ${\displaystyle \mathbf {x} }$ weakly majorizes (or dominates) ${\displaystyle \mathbf {y} }$ from below, commonly denoted ${\displaystyle \mathbf {x} \succ _{w}\mathbf {y} ,}$ when

${\displaystyle \sum _{i=1}^{k}x_{i}^{\downarrow }\geq \sum _{i=1}^{k}y_{i}^{\downarrow }}$ for all ${\displaystyle k=1,\,\dots ,\,n}$,

where ${\displaystyle x_{i}^{\downarrow }}$ denotes the ${\displaystyle i}$^th largest entry of ${\displaystyle \mathbf {x} }$. If ${\displaystyle \mathbf {x} ,\mathbf {y} }$ further satisfy ${\displaystyle \sum _{i=1}^{n}x_{i}=\sum _{i=1}^{n}y_{i}}$, we say that ${\displaystyle \mathbf {x} }$ majorizes (or dominates) ${\displaystyle \mathbf {y} }$, commonly denoted ${\displaystyle \mathbf {x} \succ \mathbf {y} }$.

Both weak majorization and majorization are partial orders for vectors whose entries are non-decreasing, but only a preorder for general vectors, since majorization is agnostic to the ordering of the entries in vectors, e.g., the statement ${\displaystyle (1,2)\prec (0,3)}$ is simply equivalent to ${\displaystyle (2,1)\prec (3,0)}$.

Specifically, ${\displaystyle \mathbf {x} \succ \mathbf {y} \wedge \mathbf {y} \succ \mathbf {x} }$ if and only if ${\displaystyle \mathbf {x} ,\mathbf {y} }$ are permutations of each other. Similarly for ${\displaystyle \succ _{w}}$.

Majorizing also sometimes refers to entrywise ordering, e.g. the real-valued function f majorizes the real-valued function g when ${\displaystyle f(x)\geq g(x)}$ for all ${\displaystyle x}$ in the domain, or other technical definitions, such as majorizing measures in probability theory.^[1]

For ${\displaystyle \mathbf {x} ,\ \mathbf {y} \in \mathbb {R} ^{n},}$ we have ${\displaystyle \mathbf {x} \prec \mathbf {y} }$ if and only if ${\displaystyle \mathbf {x} }$ is in the convex hull of all vectors obtained by permuting the coordinates of ${\displaystyle \mathbf {y} }$. This is equivalent to saying that ${\displaystyle \mathbf {x} =\mathbf {D} \mathbf {y} }$ for some doubly stochastic matrix ${\displaystyle \mathbf {D} }$.^{[2]: Thm. 2.1} In particular, ${\displaystyle \mathbf {x} }$ can be written as a convex combination of ${\displaystyle n}$ permutations of ${\displaystyle \mathbf {y} }$.^[3] In other words, ${\displaystyle \mathbf {x} }$ is in the permutahedron of ${\displaystyle \mathbf {y} }$.

Figure 1 displays the convex hull in 2D for the vector ${\displaystyle \mathbf {y} =(3,\,1)}$. Notice that the center of the convex hull, which is an interval in this case, is the vector ${\displaystyle \mathbf {x} =(2,\,2)}$. This is the "smallest" vector satisfying ${\displaystyle \mathbf {x} \prec \mathbf {y} }$ for this given vector ${\displaystyle \mathbf {y} }$. Figure 2 shows the convex hull in 3D. The center of the convex hull, which is a 2D polygon in this case, is the "smallest" vector ${\displaystyle \mathbf {x} }$ satisfying ${\displaystyle \mathbf {x} \prec \mathbf {y} }$ for this given vector ${\displaystyle \mathbf {y} }$.

Each of the following statements is true if and only if ${\displaystyle \mathbf {x} \succ \mathbf {y} }$.

• From ${\displaystyle \mathbf {x} }$ we can produce ${\displaystyle \mathbf {y} }$ by a finite sequence of "Robin Hood operations" where we replace two elements ${\displaystyle x_{i}}$ and ${\displaystyle x_{j}<x_{i}}$ with ${\displaystyle x_{i}-\varepsilon }$ and ${\displaystyle x_{j}+\varepsilon }$, respectively, for some ${\displaystyle \varepsilon \in (0,x_{i}-x_{j})}$.^[2]: 11
• For every convex function ${\displaystyle h:\mathbb {R} \to \mathbb {R} }$, ${\displaystyle \sum _{i=1}^{d}h(x_{i})\geq \sum _{i=1}^{d}h(y_{i})}$.^{[2]: Thm. 2.9}
  • In fact, a special case suffices: ${\displaystyle \sum _{i}{x_{i}}=\sum _{i}{y_{i}}}$ and, for every t, ${\displaystyle \sum _{i=1}^{d}\max(0,x_{i}-t)\geq \sum _{i=1}^{d}\max(0,y_{i}-t)}$.^[4]
• For every ${\displaystyle t\in \mathbb {R} }$, ${\displaystyle \sum _{j=1}^{d}|x_{j}-t|\geq \sum _{j=1}^{d}|y_{j}-t|}$.^{[5]: Exercise 12.17}

• 

  Each vector ${\displaystyle \mathbf {x} }$ can be plotted as a concave curve by connecting ${\displaystyle (0,0),(1,x_{1}^{\downarrow }),(2,x_{1}^{\downarrow }+x_{2}^{\downarrow }),\dots ,(n,x_{1}^{\downarrow }+x_{2}^{\downarrow }+\dots +x_{n}^{\downarrow })}$. Then ${\displaystyle \mathbf {x} \succ \mathbf {y} }$ is equivalent to the curve of ${\displaystyle \mathbf {x} }$ being higher than that of ${\displaystyle \mathbf {y} }$.

Among non-negative vectors with three components, ${\displaystyle (1,0,0)}$ and permutations of it majorize all other vectors ${\displaystyle (p_{1},p_{2},p_{3})}$ such that ${\displaystyle p_{1}+p_{2}+p_{3}=1}$. For example, ${\displaystyle (1,0,0)\succ (1/2,0,1/2)}$. Similarly, ${\displaystyle (1/3,1/3,1/3)}$ is majorized by all other such vectors, so ${\displaystyle (1/2,0,1/2)\succ (1/3,1/3,1/3)}$.

This behavior extends to general-length probability vectors: the singleton vector majorizes all other probability vectors, and the uniform distribution is majorized by all probability vectors.

A function ${\displaystyle f:\mathbb {R} ^{n}\to \mathbb {R} }$ is said to be Schur convex when ${\displaystyle \mathbf {x} \succ \mathbf {y} }$ implies ${\displaystyle f(\mathbf {x} )\geq f(\mathbf {y} )}$. Hence, Schur-convex functions translate the ordering of vectors to a standard ordering in ${\displaystyle \mathbb {R} }$. Similarly, ${\displaystyle f(\mathbf {x} )}$ is Schur concave when ${\displaystyle \mathbf {x} \succ \mathbf {y} }$ implies ${\displaystyle f(\mathbf {x} )\leq f(\mathbf {y} ).}$

An example of a Schur-convex function is the max function, ${\displaystyle \max(\mathbf {x} )=x_{1}^{\downarrow }}$. Schur convex functions are necessarily symmetric that the entries of it argument can be switched without modifying the value of the function. Therefore, linear functions, which are convex, are not Schur-convex unless they are symmetric. If a function is symmetric and convex, then it is Schur-convex.

Majorization can be generalized to the Lorenz ordering, a partial order on distribution functions. For example, a wealth distribution is Lorenz-greater than another if its Lorenz curve lies below the other. As such, a Lorenz-greater wealth distribution has a higher Gini coefficient, and has more income disparity.^[6]

The majorization preorder can be naturally extended to density matrices in the context of quantum information.^[5][7] In particular, ${\displaystyle \rho \succ \rho '}$ exactly when ${\displaystyle \mathrm {spec} [\rho ]\succ \mathrm {spec} [\rho ']}$ (where ${\displaystyle \mathrm {spec} }$ denotes the state's spectrum).

Similarly, one can say a Hermitian operator, ${\displaystyle \mathbf {H} }$, majorizes another, ${\displaystyle \mathbf {M} }$, if the set of eigenvalues of ${\displaystyle \mathbf {H} }$ majorizes that of ${\displaystyle \mathbf {M} }$.

• Muirhead's inequality
• Karamata's Inequality
• Schur-convex function
• Schur–Horn theorem relating diagonal entries of a matrix to its eigenvalues.
• For positive integer numbers, weak majorization is called Dominance order.
• Leximin order

• J. Karamata. "Sur une inegalite relative aux fonctions convexes." Publ. Math. Univ. Belgrade 1, 145–158, 1932.
• G. H. Hardy, J. E. Littlewood and G. Pólya, Inequalities, 2nd edition, 1952, Cambridge University Press, London.
• Inequalities: Theory of Majorization and Its Applications Albert W. Marshall, Ingram Olkin, Barry Arnold, Second edition. Springer Series in Statistics. Springer, New York, 2011. ISBN 978-0-387-40087-7
• A tribute to Marshall and Olkin's book "Inequalities: Theory of Majorization and its Applications"
• Matrix Analysis (1996) Rajendra Bhatia, Springer, ISBN 978-0-387-94846-1
• Topics in Matrix Analysis (1994) Roger A. Horn and Charles R. Johnson, Cambridge University Press, ISBN 978-0-521-46713-1
• Majorization and Matrix Monotone Functions in Wireless Communications (2007) Eduard Jorswieck and Holger Boche, Now Publishers, ISBN 978-1-60198-040-3
• The Cauchy Schwarz Master Class (2004) J. Michael Steele, Cambridge University Press, ISBN 978-0-521-54677-5

• Majorization in MathWorld
• Majorization in PlanetMath

• OCTAVE/MATLAB code to check majorization