Matrix norm

Matrix norm

In mathematics, a matrix norm is a natural extension of the notion of a vector norm to matrices.

Contents

Definition

In what follows, K will denote the field of real or complex numbers. Let K^{m \times n} denote the vector space containing all matrices with m rows and n columns with entries in K. Throughout the article A * denotes the conjugate transpose of matrix A.

A matrix norm is a vector norm on K^{m \times n}. That is, if \|A\| denotes the norm of the matrix A, then,

  • \|A\|> 0 if A\ne0 and \|A\|= 0 iff A = 0
  • \|\alpha A\|=|\alpha| \|A\| for all α in K and all matrices A in K^{m \times n}
  • \|A+B\| \le \|A\|+\|B\| for all matrices A and B in K^{m \times n}.

Additionally, in the case of square matrices (thus, m = n), some (but not all) matrix norms satisfy the following condition, which is related to the fact that matrices are more than just vectors:

  • \|AB\| \le \|A\|\|B\| for all matrices A and B in K^{n \times n}.

A matrix norm that satisfies this additional property is called a sub-multiplicative norm ( in some books, the terminology matrix norm is used only for those norms which are sub-multiplicative). The set of all n-by-n matrices, together with such a sub-multiplicative norm, is an example of a Banach algebra.

Induced norm

If vector norms on Km and Kn are given (K is field of real or complex numbers), then one defines the corresponding induced norm or operator norm on the space of m-by-n matrices as the following maxima:

 \begin{align}
\|A\| &= \max\{\|Ax\| : x\in K^n \mbox{ with }\|x\|= 1\} \\
&= \max\left\{\frac{\|Ax\|}{\|x\|} : x\in K^n \mbox{ with }x\ne 0\right\}.
\end{align}

These are different from the entrywise p-norms and the Schatten p-norms for matrices treated below, which are also usually denoted by  \left \| A \right \| _p .

If m = n and one uses the same norm on the domain and the range, then the induced operator norm is a sub-multiplicative matrix norm.

The operator norm corresponding to the p-norm for vectors is:

 \left \| A \right \| _p = \max \limits _{x \ne 0} \frac{\left \| A x\right \| _p}{\left \| x\right \| _p}.

In the case of p = 1 and p=\infty, the norms can be computed as:

 \left \| A \right \| _1 = \max \limits _{1 \leq j \leq n} \sum _{i=1} ^m | a_{ij} |, which is simply the maximum absolute column sum of the matrix.
 \left \| A \right \| _\infty = \max \limits _{1 \leq i \leq m} \sum _{j=1} ^n | a_{ij} |, which is simply the maximum absolute row sum of the matrix

For example, if the matrix A is defined by

 
      A = \begin{bmatrix}
           3 & 5 & 7 \\
           2 & 6 & 4 \\
           0 & 2 & 8 \\
        \end{bmatrix},

then we have ||A||1 = max(5,13,19) = 19. and ||A|| = max(15,12,10) = 15. Consider another example


      A = \begin{bmatrix}
           2 & 4 & 2 & 1 \\
           3 & 1 & 5 & 2 \\
           1 & 2 & 3 & 3 \\
           0 & 6 & 1 & 2 \\
         \end{bmatrix},

where we add all the entries in each column and determine the greatest value, which results in ||A||1 = max (6,13,11,8) = 13.

We can do the same for the rows and get ||A|| = max(9,11,9,9) = 11. Thus 11 is our max.

In the special case of p = 2 (the Euclidean norm) and m = n (square matrices), the induced matrix norm is the spectral norm. The spectral norm of a matrix A is the largest singular value of A i.e. the square root of the largest eigenvalue of the positive-semidefinite matrix A*A:

\left \| A \right \| _2=\sqrt{\lambda_{\text{max}}(A^{^*} A)}=\sigma_{\text{max}}(A)

where A* denotes the conjugate transpose of A.

More generally, one can define the subordinate matrix norm on K^{m\times n} induced by \|\cdot\|_{\alpha} on Kn, and \|\cdot\|_{\beta} on Km as:

 \left \| A \right \| _{\alpha,\beta} = \max \limits _{x \ne 0} \frac{\left \| A x\right \|_{\beta}}{\left \| x\right \|_{\alpha}}.

Subordinate norms are consistent with the norms that induce them, giving


\|Ax\|_{\beta}\leq \|A\|_{\alpha,\beta}\|x\|_{\alpha}.

Any induced norm satisfies the inequality

\left \| A \right \| \ge \rho(A),

where ρ(A) is the spectral radius of A. For a symmetric or hermitian matrix A, we have equality for the 2-norm, since in this case the 2-norm is the spectral radius of A. For an arbitrary matrix, we may not have equality for any . Take


      A = \begin{bmatrix}
           0 & 1  \\
           0 & 0  \\
         \end{bmatrix},

the spectral radius of A is 0, but A is not the zero matrix, and so none of the induced norms are equal to the spectral radius of A.

Furthermore, for square matrices we have the spectral radius formula:

\lim_{r\rarr\infty}\|A^r\|^{1/r}=\rho(A).

"Entrywise" norms

These vector norms treat an  m \times n matrix as a vector of size mn, and use one of the familiar vector norms.

For example, using the p-norm for vectors, we get:

\Vert A \Vert_{p} = \left( \sum_{i=1}^m \sum_{j=1}^n |a_{ij}|^p \right)^{1/p}. \,

This is a different norm from the induced p-norm (see above) and the Schatten p-norm (see below), but the notation is the same.

The special case p = 2 is the Frobenius norm, and p = ∞ yields the maximum norm.

Frobenius norm

For p = 2, this is called the Frobenius norm or the Hilbert–Schmidt norm, though the latter term is often reserved for operators on Hilbert space. This norm can be defined in various ways:

\|A\|_F=\sqrt{\sum_{i=1}^m\sum_{j=1}^n |a_{ij}|^2}=\sqrt{\operatorname{tr}(AA^{{}^*})}=\sqrt{\sum_{i=1}^{\min\{m,\,n\}} \sigma_{i}^2}

where A* denotes the conjugate transpose of A, σi are the singular values of A, and the trace function is used. The Frobenius norm is very similar to the Euclidean norm on Kn and comes from an inner product on the space of all matrices.

The Frobenius norm is sub-multiplicative and is very useful for numerical linear algebra. This norm is often easier to compute than induced norms.

Max norm

The max norm is the elementwise norm with p = ∞:

 \|A\|_{\text{max}} = \max \{|a_{ij}|\}.

This norm is not sub-multiplicative.

Schatten norms

The Schatten p-norms arise when applying the p-norm to the vector of singular values of a matrix. If the singular values are denoted by σi, then the Schatten p-norm is defined by

 \|A\|_p = \left( \sum_{i=1}^{\min\{m,\,n\}} \sigma_i^p \right)^{1/p}. \,

These norms again share the notation with the induced and entrywise p-norms, but they are different.

All Schatten norms are sub-multiplicative. They are also unitarily invariant, which means that ||A|| = ||UAV|| for all matrices A and all unitary matrices U and V.

The most familiar cases are p = 1, 2, ∞. The case p = 2 yields the Frobenius norm, introduced before. The case p = ∞ yields the spectral norm, which is the matrix norm induced by the vector 2-norm (see above). Finally, p = 1 yields the nuclear norm (also known as the trace norm, or the Ky Fan 'n'-norm), defined as

\|A\|_{*} = \operatorname{trace} \left(\sqrt{A^*A}\right) = \sum_{i=1}^{\min\{m,\,n\}} \sigma_i.

(Here \sqrt{A^*A} denotes a matrix B such that BB = A * A. More precisely, since A * A is a positive semidefinite matrix, its square root is well-defined.)

Consistent norms

A matrix norm \| \cdot \|_{ab} on K^{m \times n} is called consistent with a vector norm \| \cdot \|_{a} on Kn and a vector norm \| \cdot \|_{b} on Km if:

\|Ax\|_b \leq \|A\|_{ab} \|x\|_a

for all A \in K^{m \times n}, x \in K^n. All induced norms are consistent by definition.

Compatible norms

A matrix norm \| \cdot \|_{b} on K^{n \times n} is called compatible with a vector norm \| \cdot \|_{a} on Kn if:

\|Ax\|_a \leq \|A\|_b \|x\|_a

for all A \in K^{m \times n}, x \in K^n. Induced norms are compatible by definition.

Equivalence of norms

For any two vector norms ||·||α and ||·||β, we have

r\left\|A\right\|_\alpha\leq\left\|A\right\|_\beta\leq s\left\|A\right\|_\alpha

for some positive numbers r and s, for all matrices A in K^{m \times n}. In other words, they are equivalent norms; they induce the same topology on K^{m \times n}. This is a special case of the equivalence of norms in finite-dimensional Normed vector spaces.

Moreover, for every vector norm \|\cdot\| on \mathbb{R}^{n\times n}, there exists a unique positive real number k such that l\|\cdot\| is a sub-multiplicative matrix norm for every l \ge k.

A matrix norm ||·||α is said to be minimal if there exists no other matrix norm ||·||β satisfying ||·||β ≤ ||·||α.

Examples of norm equivalence

For matrix A\in\mathbb{R}^{m\times n} the following inequalities hold:[1][2]

  • \|A\|_2\le\|A\|_F\le\sqrt{r}\|A\|_2, where r is the rank of A
  • \|A\|_F \le \|A\|_{*} \le \sqrt{r} \|A\|_F, where r is the rank of A
  • \|A\|_{\text{max}} \le \|A\|_2 \le \sqrt{mn}\|A\|_{\text{max}}
  • \frac{1}{\sqrt{n}}\|A\|_\infty\le\|A\|_2\le\sqrt{m}\|A\|_\infty
  • \frac{1}{\sqrt{m}}\|A\|_1\le\|A\|_2\le\sqrt{n}\|A\|_1.

Here, ||·||p refers to the matrix norm induced by the vector p-norm.

Another useful inequality between matrix norms is

\|A\|_2\le\sqrt{\|A\|_1\|A\|_\infty}.

Notes

  1. ^ Golub, Gene; Charles F. Van Loan (1996). Matrix Computations - Third Edition. Baltimore: The Johns Hopkins University Press, 56-57. ISBN 0-8018-5413-X.
  2. ^ Roger Horn and Charles Johnson. Matrix Analysis, Chapter 5, Cambridge University Press, 1985. ISBN 0-521-38632-2.

References

  • James W. Demmel, Applied Numerical Linear Algebra, section 1.7, published by SIAM, 1997.
  • Carl D. Meyer, Matrix Analysis and Applied Linear Algebra, published by SIAM, 2000. [1]
  • John Watrous, Theory of Quantum Information, 2.4 Norms of operators, lecture notes, University of Waterloo, 2008.
  • Kendall Atkinson, An Introduction to Numerical Analysis, published by John Wiley & Sons, Inc 1989

Wikimedia Foundation. 2010.

Игры ⚽ Нужно сделать НИР?

Look at other dictionaries:

  • Norm — or NORM may refer to: Contents 1 In academia 1.1 In mathematics 2 People 3 Miscellaneous 4 …   Wikipedia

  • Norm (mathematics) — This article is about linear algebra and analysis. For field theory, see Field norm. For ideals, see Norm of an ideal. For group theory, see Norm (group). For norms in descriptive set theory, see prewellordering. In linear algebra, functional… …   Wikipedia

  • Matrix function — In mathematics, a matrix function is a function which maps a matrix to another matrix. Contents 1 Extending scalar functions to matrix functions 1.1 Power series 1.2 Jordan decomposition …   Wikipedia

  • Matrix exponential — In mathematics, the matrix exponential is a matrix function on square matrices analogous to the ordinary exponential function. Abstractly, the matrix exponential gives the connection between a matrix Lie algebra and the corresponding Lie group.… …   Wikipedia

  • Matrix Games — Тип Частная компания Год основания 2001 Расположение Статен …   Википедия

  • Matrix (mathematics) — Specific elements of a matrix are often denoted by a variable with two subscripts. For instance, a2,1 represents the element at the second row and first column of a matrix A. In mathematics, a matrix (plural matrices, or less commonly matrixes)… …   Wikipedia

  • Matrix Games — For the Japanese company, see Matrix Software. Matrix Games Industry Strategy and Wargames for Personal Computers. Founded 1999 Headquarters Staten Island, New York Key …   Wikipedia

  • Matrix multiplication — In mathematics, matrix multiplication is a binary operation that takes a pair of matrices, and produces another matrix. If A is an n by m matrix and B is an m by p matrix, the result AB of their multiplication is an n by p matrix defined only if… …   Wikipedia

  • Matrix calculus — Topics in Calculus Fundamental theorem Limits of functions Continuity Mean value theorem Differential calculus  Derivative Change of variables Implicit differentiation Taylor s theorem Related rates …   Wikipedia

  • Matrix ring — In abstract algebra, a matrix ring is any collection of matrices forming a ring under matrix addition and matrix multiplication. The set of n×n matrices with entries from another ring is a matrix ring, as well as some subsets of infinite matrices …   Wikipedia

Share the article and excerpts

Direct link
Do a right-click on the link above
and select “Copy Link”