Hadamard product (matrices)

The Hadamard product operates on identically-shaped matrices and produces a third matrix of the same dimensions.

In mathematics, the Hadamard product (also known as the Schur product^[1] or the entrywise product^[2]) is a binary operation that takes two matrices of the same dimensions, and produces another matrix where each element ij is the product of elements ij of the original two matrices. It should not be confused with the more common matrix product. It is attributed to, and named after, either French mathematician Jacques Hadamard, or German mathematician Issai Schur.

The Hadamard product is associative and distributive, and unlike the matrix product it is also commutative.

Definition

For two matrices, $A,B$ , of the same dimension, $m\times n$ , the Hadamard product, $A\circ B$ , is a matrix, of the same dimension as the operands, with elements given by

(A\circ B)_{i,j}=(A)_{i,j}(B)_{i,j}

For matrices of different dimensions ( $m\times n$ and $p\times q$ , where $m\not =p$ or $n\not =q$ or both) the Hadamard product is undefined.

Example

For example, the Hadamard product for a 3×3 matrix A with a 3×3 matrix B is:

\left({\begin{array}{ccc}\mathrm {a} _{11}&\mathrm {a} _{12}&\mathrm {a} _{13}\\\mathrm {a} _{21}&\mathrm {a} _{22}&\mathrm {a} _{23}\\\mathrm {a} _{31}&\mathrm {a} _{32}&\mathrm {a} _{33}\end{array}}\right)\circ \left({\begin{array}{ccc}\mathrm {b} _{11}&\mathrm {b} _{12}&\mathrm {b} _{13}\\\mathrm {b} _{21}&\mathrm {b} _{22}&\mathrm {b} _{23}\\\mathrm {b} _{31}&\mathrm {b} _{32}&\mathrm {b} _{33}\end{array}}\right)=\left({\begin{array}{ccc}\mathrm {a} _{11}\,\mathrm {b} _{11}&\mathrm {a} _{12}\,\mathrm {b} _{12}&\mathrm {a} _{13}\,\mathrm {b} _{13}\\\mathrm {a} _{21}\,\mathrm {b} _{21}&\mathrm {a} _{22}\,\mathrm {b} _{22}&\mathrm {a} _{23}\,\mathrm {b} _{23}\\\mathrm {a} _{31}\,\mathrm {b} _{31}&\mathrm {a} _{32}\,\mathrm {b} _{32}&\mathrm {a} _{33}\,\mathrm {b} _{33}\end{array}}\right)

Properties

The Hadamard product is commutative, associative and distributive over addition. That is,

{\begin{aligned}A\circ B&=B\circ A,\\A\circ (B\circ C)&=(A\circ B)\circ C,\\A\circ (B+C)&=A\circ B+A\circ C.\end{aligned}}

The identity matrix under Hadamard multiplication of two m-by-n matrices is an m-by-n matrix where all elements are equal to 1. This is different from the identity matrix under regular matrix multiplication, where only the elements of the main diagonal are equal to 1. Furthermore, a matrix has an inverse under Hadamard multiplication if and only if none of the elements are equal to zero.^[3]
For vectors $x$ and $y$ , and corresponding diagonal matrices $D_{x}$ and $D_{y}$ with these vectors as their leading diagonals, the following identity holds:^[4]

x^{*}(A\circ B)y=\mathrm {tr} \left(D_{x}^{*}AD_{y}B^{\mathrm {T} }\right),

where

x^{*}

denotes the conjugate transpose of

x

. In particular, using vectors of ones, this shows that the sum of all elements in the Hadamard product is the trace of

AB^{\mathrm {T} }

. A related result for square

A

and

B

, is that the row-sums of their Hadamard product are the diagonal elements of

AB^{\mathrm {T} }

:^[5]

\sum _{i}(A\circ B)_{i,j}=\left(B^{\mathrm {T} }A\right)_{j,j}.

\sum _{j}(A\circ B)_{i,j}=\left(AB^{\mathrm {T} }\right)_{i,i}.

The Hadamard product is a principal submatrix of the Kronecker product.
The Hadamard product satisfies the rank inequality

\operatorname {rank} (A\circ B)\leq \operatorname {rank} (A)\operatorname {rank} (B)

If $A$ and $B$ are positive-definite matrices, then the following inequality involving the Hadamard product is valid:^[6]

\prod _{i=k}^{n}\lambda _{i}(A\circ B)\geq \prod _{i=k}^{n}\lambda _{i}(AB),\quad k=1,\ldots ,n,

where

\lambda _{i}(A)

is the ith largest eigenvalue of

A

Schur product theorem

The Hadamard product of two positive-semidefinite matrices is positive-semidefinite.^[5] This is known as the Schur product theorem,^[3] after German mathematician Issai Schur. For positive-semidefinite matrices A and B, it is also known that

\det(A\circ B)\geq \det(A)\det(B).\,

^[5]

In programming languages

Hadamard multiplication is built into certain programming languages under various names. In MATLAB, GNU Octave, GAUSS and HP Prime, it is known as array multiplication, or in Julia broadcast multiplication, with the symbol .*.^[7] In Fortran, R,^[8] J and Wolfram Language (Mathematica), it is done through simple multiplication operator *, whereas the matrix product is done through the function matmul, %*%, +/ .* and the . operators, respectively. In Python with the numpy numerical library or the sympy symbolic library, multiplication of array objects as a1*a2 produces the Hadamard product, but with otherwise matrix objects m1*m2 will produce a matrix product. The Eigen C++ library provides a cwiseProduct member function for the Matrix class (a.cwiseProduct(b)), while the Armadillo library use the operator % to make compact expressions (a % b; a * b is a matrix product).

Applications

The Hadamard product appears in lossy compression algorithms such as JPEG. The decoding step involves an entry-for-entry product, i.e., Hadamard product.

Analogous operations

Other Hadamard operations are also seen in the mathematical literature,^[9] namely the Hadamard root and Hadamard power (which are in effect the same thing because of fractional indices), defined for a matrix such that: For

{\begin{aligned}B&=A^{\circ 2}\\B_{ij}&=A_{ij}^{2}\end{aligned}}

and for

{\begin{aligned}B&=A^{\circ {\frac {1}{2}}}\\B_{ij}&=A_{ij}^{\frac {1}{2}}\end{aligned}}

The Hadamard inverse reads:^[9]

{\begin{aligned}B&=A^{\circ -1}\\B_{ij}&=A_{ij}^{-1}\end{aligned}}

A Hadamard division is defined as:^[10]^[11]

{\begin{aligned}C&=A{\oslash }B\\C_{ij}&=A_{ij}/B_{ij}\end{aligned}}

References

↑ Davis, Chandler. "The norm of the Schur product operation." Numerische Mathematik 4.1 (1962): 343-344.
↑ (Horn & Johnson 1985, Ch. 5)
1 2 Million, Elizabeth. "The Hadamard Product" (PDF). Retrieved 2 January 2012.
↑ Horn, Roger A., and Charles R. Johnson. Matrix analysis. Cambridge university press, 2012. p. 479
1 2 3 (Styan 1973)
↑ Hiai, Fumio; Lin, Minghua (February 2017). "On an eigenvalue inequality involving the Hadamard product". Linear Algebra and its Applications. 515: 313–320. doi:10.1016/j.laa.2016.11.017.
↑ "Arithmetic Operators + - * / \ ^ ' -". MATLAB documentation. MathWorks. Retrieved 2 January 2012.
↑ "Matrix multiplication". An Introduction to R. The R Project for Statistical Computing. 16 May 2013. Retrieved 24 August 2013.
1 2 Reams, Robert (1999). "Hadamard inverses, square roots and products of almost semidefinite matrices". Linear Algebra and its Applications. 288: 35–43. doi:10.1016/S0024-3795(98)10162-3.
↑ Wetzstein, Gordon; Lanman, Douglas; Hirsch, Matthew; Raskar, Ramesh. "Supplementary Material: Tensor Displays: Compressive Light Field Synthesis using Multilayer Displays with Directional Backlighting" (PDF). MIT Media Lab.
↑ Cyganek, Boguslaw (2013). Object Detection and Recognition in Digital Images: Theory and Practice. John Wiley & Sons. p. 109.

Linear algebra
Basic concepts	Scalar Vector Vector space Scalar multiplication Vector projection Linear span Linear map Linear projection Linear independence Linear combination Basis Column space Row space Dual space Orthogonality Kernel Eigenvalues and eigenvectors Outer product Inner product space Dot product Transpose Gram–Schmidt process Linear equations
Vector algebra	Cross product Triple product Seven-dimensional cross product
Multilinear algebra	Geometric algebra Exterior algebra Bivector Multivector
Matrices	Block Decomposition Invertible Minor Multiplication Rank Transformation Cramer's rule Gaussian elimination
Numerical	Floating point Matrix Laboratory Numerical stability Basic Linear Algebra Subprograms (BLAS) Sparse matrix Comparison of linear algebra libraries Comparison of numerical analysis software

This article is issued from Wikipedia. The text is licensed under Creative Commons - Attribution - Sharealike. Additional terms may apply for the media files.