User:Matteo.focardi/sandbox

2020 Mathematics Subject Classification: Primary: 15Axx [MSN][ZBL]

A formula aimed at expressing the determinant of the product of two matrices $A\in\mathrm{M}_{m,n}(\mathbb{R})$ and $B\in\mathrm{M}_{n,m}(\mathbb{R})$, in terms of the sum of the products of all possible higher order minors of $A$ with corresponding minors of the same order of $B$. More precisely, if $\alpha=(1,\ldots,m)$ and $\beta$ denotes any multi-index $(\beta_1,\ldots,\beta_m)$ with $1\leq \beta_1<\ldots<\beta_m\leq n$ of length $m$, then \[ \det(AB)=\sum_\beta\det A_{\alpha\,\beta}\det B_{\beta\,\alpha}, \] where $A_{\alpha\,\beta}=(a_{\alpha_i\beta_j})$ and $B_{\beta\,\alpha}=(a_{\beta_j\alpha_i})$. In case $m>n$, no such $\beta$ exists and the right-hand side above is set to be $0$ by definition.

Note that if $n=m$ the formula reduces to \[ \det (AB)=\det A\,\det B. \] More generally, if $A\in\mathrm{M}_{m,n}(\mathbb{R})$, $B\in\mathrm{M}_{n,q}(\mathbb{R})$ and $N\leq\min\{m,q\}$, then any minor of order $N$ of the product matrix $AB$ can be expressed as follows by Cauchy-Binet's formula \[ \det C_{\alpha,\,\gamma}=\sum_\beta\det A_{\alpha\,\beta}\det B_{\beta,\,\gamma}, \] where $\alpha=(\alpha_1\ldots,\alpha_N)$, $1\leq\alpha_1<\ldots<\alpha_N\leq m$, $\gamma=(\gamma_1,\ldots,\gamma_N)$, $1\leq\gamma_1<\ldots<\gamma_N\leq q$, and $(\beta_1,\ldots,\beta_m)$ with $1\leq \beta_1<\ldots<\beta_m\leq n$.

A number of interesting consequence of Cauchy-Binet's formula is listed below. First of all, an inequality for the rank of the product matrix follows straightforwardly, i.e., \[ \mathrm{rank}(AB)\leq\min\{\mathrm{rank}A,\mathrm{rank}B\}. \] Moreover, if $m=2$, $\mathbf{a}$, $\mathbf{b}\in\mathbb{R}^n$ are two vectors, by taking $$A=\begin{pmatrix} a_{1}&\dots&a_{n}\\ b_{1}&\dots&b_{n}\\ \end{pmatrix} \quad\text{and}\quad B=\begin{pmatrix} a_{1}&b_{1}\\ \dots&\dots\\ a_{n}&b_{n}\\ \end{pmatrix} $$ by Cauchy-Binet's formula \[ \sum_{1\leq i<j\leq n}\begin{vmatrix} a_{i}&a_{j}\\ b_{i}&b_{j}\\ \end{vmatrix}^2= \begin{vmatrix} \|\mathbf{a}\|^2&\mathbf{a}\cdot \mathbf{b}\\ \mathbf{a}\cdot \mathbf{b}&\|\mathbf{b}\|^2\\ \end{vmatrix}, \] in turn implying Cauchy-Schwartz's inequality, here $\|\cdot\|$ and $\cdot$ are the Euclidean norm and scalar product, respectively.

Let us finally take $B=A^T$, then $\det(A_{\alpha\beta})=\det(A^T_{\beta,\alpha})$, so that by Cauchy-Binet's formula \[\label{p} \det(A^T\,A)=\sum_\beta(\det(A_{\alpha\beta}))^2. \] This is a generalization of the Pythagorean formula, corresponding to $m=1$. Indeed, if $\mathcal{A}:\mathbb{R}^n\to\mathbb{R}^m$ is the linear map associated to $A$ and $Q\subset\mathbb{R}^n$ is the unitary cube, the $n$-th dimensional volume of the parallelepiped $\mathcal{A}(Q)$, recall that $n\leq m$, is $\sqr{\det(A^T\,A)}$ by polar decomposition of $A$.

Then, formula (1) above expresses that the square of the $n$-th dimensional volume of $\mathcal{A}(Q)$ can be obtained by summing up the volumes of the projections on all coordinates $n$ planes.