Матрицы и операции с ними
Элементы матриц часто обозначают одной буквой с двумя индексами $\alpha_{ij}$, первый из которых обозначает номер строки, а второй — номер столбца. Таким образом, матрица \begin{equation} A = \left( \begin{array}{cccc} \alpha_{11} & \alpha_{12} & ... & \alpha_{1n} \\ \alpha_{21} & \alpha_{22} & ... & \alpha_{2n} \\ ... & ... & ... & ... \\ \alpha_{m1} & \alpha_{m2} & ... & \alpha_{mn} \end{array} \right) = \{ \alpha_{ij} \}_{i,j=1}^{m,n}. \end{equation}
Если строки матриц обозначить латинской буквой с индексом внизу $a_{i}$, $i=\overline{1,m}$, а столбцы — с индексом вверху $a^{j}$, $j=\overline{1,n}$, то матрицу $A$ можно записать как $A = [a_{1}, a_2, ..., a_{m} ] = [a^{1}, a^{2}, ..., a^{n}]$.
Если необходимо подчеркнуть размер матрицы, то пишут $A=[m\times n]$.
\begin{equation}
A = \left(\begin{array}{cccc}
\alpha_{11} & 0 & ... & 0 \\
0 & \alpha_{22} & ... & 0 \\
... & ... & ... & ... \\
0 & 0 & ... & \alpha_{nn}
\end{array}\right) = diag\{ \alpha_{11}, \alpha_{22}, ..., \alpha_{nn} \} .
\end{equation}
Таким образом, если ${\bf A} = \{ \alpha_{ij} \}_{i,j=1}^{m,n}$, ${\bf B} = \{ \beta_{ij} \}_{i,j=1}^{m,n}$, то сумма ${\bf A} + {\bf B} = \{ \alpha_{ij} + \beta_{ij} \}_{i,j=1}^{m,n}$.
Свойства суммы матриц и произведения матрицы на число
- ${\bf A} + {\bf B} = {\bf B} + {\bf A}$
- $({\bf A} + {\bf B}) + {\bf C} = {\bf A} + ({\bf B} + {\bf C})$
- Матрица ${\bf 0}$, все элементы которой равны нулю, называется нулевой и играет роль нуля: ${\bf A} + {\bf 0} = {\bf A}$
- Для любой матрицы ${\bf А}$ существует противоположная $-{\bf A} = (-1){\bf A}$, то есть ${\bf A}+(-{\bf A})=0$
- $(\gamma_{1} +\gamma_2 ){\bf A} = \gamma_{1} {\bf A} + \gamma_2 {\bf A}$
- $\gamma ({\bf A} + {\bf B}) = \gamma {\bf A} + \gamma {\bf B}$
- $\gamma_{1} (\gamma_2 {\bf A}) = (\gamma_{1} \gamma_2) {\bf A}$
- $1\cdot {\bf A} = {\bf A}$.
Непосредственно из данного определения следует, что $\gamma_{ik} = a_{i} b^{k} = (\alpha_{i1} \alpha_{i2} ... \alpha_{in} ) \left(\begin{array}{c} \beta_{1k} \\ \beta_{2k} \\ \vdots \\ \beta_{nk} \end{array}\right) = \alpha_{i1} \beta_{1k} + \alpha_{i2} \beta_{2k} + \cdots + \alpha_{in} \beta_{nk}.$
Заметим, что если ${\bf A} = [a_{1}, a_2, ..., a_{m} ] = [m \times n]$, ${\bf B} = [b^{1}, b^{2}, ..., b^{p} ] = [n \times p]$, то произведение
${\bf AB} = {\bf A}[b^{1}, b^{2}, ..., b^{p}] = [{\bf A}b^{1}, {\bf A}b^{2}, ..., {\bf A}b^{p}], {\bf AB} = [a_{1}, a_2, ..., a_{m}]{\bf B} = [a_{1} {\bf B}, a_2 {\bf B}, ..., a_{m} {\bf B}].$
Пусть неизвестные $\xi_{1}, \xi_2, ..., \xi_{m}$ связаны с системой неизвестных $\eta_{1}, \eta_2, ..., \eta_n$ с помощью линейных соотношений \begin{equation} \left\{\begin{array}{c} \xi_{1} = \alpha_{11} \eta_{1} + \alpha_{12} \eta_2 + \cdots \alpha_{1n} \eta_n, \\ \xi_2 = \alpha_{21} \eta_{1} + \alpha_{22} \eta_2 + \cdots \alpha_{2n} \eta_n, \\ \cdots \cdots \cdots \\ \xi_{m} = \alpha_{m1} \eta_{1} + \alpha_{m2} \eta_2 + \cdots \alpha_{mn} \eta_n, \end{array}\right. \end{equation} где $\alpha_{ij} $ некоторые числа, $i=\overline{1,m}$, $j=\overline{1,n}$.
Переход от $\xi_{1}, \xi_2, ..., \xi_{m}$ к $\eta_{1}, \eta_2, ..., \eta_n$ по приведенным формулам называется линейным преобразованием неизвестных.
Введем матрицу ${\bf A} = [a_{1}, a_2, ..., a_{m} ] = [a^{1}, a^{2}, ..., a^{n} ] = \{ \alpha_{ij} \}_{i,j=1}^{m,n} $ и столбцы $x = \left(\begin{array}{c} \xi_{1} \\ \xi_2 \\ \vdots \\ \xi_{m} \end{array}\right)$, $y = \left(\begin{array}{c} \eta_{1} \\ \eta_2 \\ \vdots \\ \eta_n \end{array}\right)$. Тогда в матричной записи линейные преобразования можно представить в виде \begin{equation} x = {\bf A}y \Leftrightarrow \xi_{i} = \sum_{j=1}^{n}\alpha_{ij} \eta_{j}, i=\overline{1,m} \Leftrightarrow \xi_{i} = a_{i} y, i=\overline{1,m} \Leftrightarrow x = \sum_{j=1}^{n}a^{j} \eta_{j}. \end{equation}
С одной стороны, очевидно, что $x={\bf A}({\bf B}z)$. С другой стороны, $\xi_{i} = \sum_{k=1}^{n} {\alpha_{ik} \eta_{k} } = \sum_{k=1}^{n}\alpha_{ik} \left( \sum_{j=1}^{p}\beta_{kj} \zeta_{j} \right) = \sum_{j=1}^{p}\left( \sum_{k=1}^{n}\alpha_{ik} \beta_{kj} \right) \zeta_{j} = \sum_{j=1}^{p}\gamma_{ij} \zeta_{j}, i=\overline{1,m}$, где $\beta_{kj}$ — элемент матрицы ${\bf B} = [n \times p]$, $\zeta_{j}$ — элемент столбца $z$, $j=\overline{1,p}$.
По определению произведения матриц $\gamma_{ij}$ является элементом матрицы ${\bf AB}$, то есть $x=({\bf AB})z$. Таким образом, для любого столбца $z$ справедливо равенство $x={\bf A}({\bf B}z)=({\bf AB})z$.
Свойства операции произведения матриц
- $({\bf A}+{\bf B}){\bf C} = {\bf AC} + {\bf BC}$, ${\bf A}=[m\times n]$, ${\bf B}=[m\times n]$, ${\bf C}=[n\times p]$
- ${\bf A} ({\bf B}+{\bf C} ) = {\bf AB} + {\bf AC}$, ${\bf A}=[m\times n]$, ${\bf B}=[n\times p]$, ${\bf C}=[n\times p]$
- $({\bf AB}){\bf C}={\bf A}({\bf BC})$, ${\bf A}=[m\times n]$, ${\bf B}=[n\times p]$, ${\bf C}=[p\times s]$
- Если ${\bf E}_n =[n\times n]$, ${\bf E}_{m} =[m\times m]$ две единичные матрицы, то для любой ${\bf A}=[m\times n]$ справедливы равенства ${\bf AE}_n = {\bf E}_{m} {\bf A} = {\bf A}$.
Таким образом, если матрица
\begin{equation}
{\bf A} = \left(\begin{array}{cccc}
\alpha_{11} & \alpha_{12} & ... & \alpha_{1n} \\
\alpha_{21} & \alpha_{22} & ... & \alpha_{2n} \\
... & ... & ... & ... \\
\alpha_{m1} & \alpha_{m2} & ... & \alpha_{mn}
\end{array}\right) = [m\times n],
\end{equation}
то транспонированная матрица
\begin{equation}
{\bf A}^{T} = \left(\begin{array}{cccc}
\alpha_{11} & \alpha_{21} & ... & \alpha_{m1} \\
\alpha_{12} & \alpha_{22} & ... & \alpha_{m2} \\
... & ... & ... & ... \\
\alpha_{1n} & \alpha_{2n} & ... & \alpha_{nm}
\end{array}\right) = [n\times m].
\end{equation}
Следовательно, если ${\bf A}^{T} = \{ \delta_{ks} \}_{k,s=1}^{n,m} $, то $\delta_{ks} =\alpha_{sk} $, $s=\overline{1,m}$, $k=\overline{1,n}$.
Непосредственно из определения следует, что $\left({\bf A}^{T} \right)^{T} = {\bf A}$, $({\bf A}+{\bf B})^{T} = {\bf A}^{T} + {\bf B}^{T}$, $(\gamma {\bf A})^{T} =\gamma {\bf A}^{T}$.
Таким образом, если выполнено равенство ${\bf А} = {\bf А}^{Т}$, то ${\bf А}$ — симметрическая матрица и ее элементы симметричны относительно главной диагонали.
Пусть ${\bf A} = \{ \alpha_{ij} \}_{i,j=1}^{m,n} $, ${\bf B} = \{ \beta_{jl} \}_{j,l=1}^{n,p}$. Тогда ${\bf AB} = \left\{\sum_{j=1}^{n}\alpha_{ij} \beta_{jl} \right\}_{i,l=1}^{m,p} =\left\{\gamma_{il} \right\}_{i,l=1}^{m,p}$. Допустим, что ${\bf B}^{T} = \{ \delta_{lj} \}_{l,j=1}^{p,n} $, ${\bf A}^{T} = \{ \theta_{ji} \}_{j,i=1}^{n,m} $, получаем ${\bf B}^{T} {\bf A}^{T} = \left\{\sum_{j=1}^{n}\delta_{lj} \theta_{ji} \right\}_{l,i=1}^{p,m} = \{ \mu_{li} \}_{l,i=1}^{p,m}$. Но $\mu_{li} = \sum_{j=1}^{n}{\delta_{lj} \theta_{ji} } = \sum_{j=1}^{n}{\beta_{jl} \alpha_{ij} } = \sum_{j=1}^{n} {\alpha_{ij} \beta_{jl} } =\gamma_{il},$ то есть окончательно имеем равенство ${\bf B}^{T} {\bf A}^{T} =({\bf AB})^{T} $.