Дискуссионный математический форумМатематический форум

Математический форум Math Help Planet

Обсуждение и решение задач по математике, физике, химии, экономике

Теоретический раздел
Часовой пояс: UTC + 4 часа [ Летнее время ]
MathHelpPlanet.com RSS-лента Математического форума

Часовой пояс: UTC + 4 часа [ Летнее время ]


Блочные матрицы и кронекеровские произведение и сумма матриц
ОглавлениеЛинейная алгебра

Блочные матрицы и кронекеровские произведение и сумма матриц


Блочные (клеточные) матрицы и операции над ними


Числовая матрица [math]A[/math] размеров [math]m\times n[/math], разделенная горизонтальными и вертикальными линиями на блоки (клетки), которые представляют собой матрицы, называется блочной (клеточной) матрицей. Элементами блочной матрицы [math]A[/math] являются матрицы [math]A_{ij}[/math] размеров [math]m_i\times n_j[/math], [math]i=1,2,\ldots,p,[/math] [math]j=1,2,\ldots,q[/math], причем [math]m_1+m_2+\ldots+m_p=m[/math] и [math]n_1+n_2+\ldots+n_q=n[/math].


Операции с блочными матрицами выполняются по тем же правилам, что и с числовыми матрицами. Если числовые матрицы [math]A[/math] и [math]B[/math] равных размеров одинаково разбиты на блоки [math]A=\begin{pmatrix}A_{ij}\end{pmatrix}[/math] и [math]B=\begin{pmatrix}B_{ij}\end{pmatrix}[/math], то их сумму [math]C=A+B[/math] можно аналогичным образом разбить на блоки [math]C=\begin{pmatrix}C_{ij}\end{pmatrix}[/math], причем для каждого блока [math]C_{ij}=A_{ij}+B_{ij}[/math]. Если блочную матрицу [math]A=\begin{pmatrix}A_{ij}\end{pmatrix}[/math] умножить на число [math]\lambda[/math], то получим матрицу [math]\lambda A=A\lambda=\begin{pmatrix}\lambda A_{ij}\end{pmatrix}[/math]. При транспонировании блочной матрицы транспонированию подлежит блочная структура и все ее блоки, например


[math]A^T= \begin{pmatrix}A_{11}&\!\!\vline&\!\!A_{12}\\\hline A_{21}&\!\!\vline&\!\!A_{22}\end{pmatrix}^T= \begin{pmatrix} A_{11}^T&\!\!\vline&\!\!A_{21}^T\\\hline A_{12}^T&\!\!\vline&\!\!A_{22}^T\end{pmatrix}\!.[/math]



Пример 1.25. Даны блочные матрицы


[math]A=\begin{pmatrix}2&3&\!\!\vline&\!\!4\\\hline 3&4&\!\!\vline&\!\!5\\ 4&5&\!\!\vline&\!\!6\end{pmatrix}= \begin{pmatrix}A_{11}&\!\!\vline&\!\!A_{12}\\\hline A_{21}&\!\!\vline&\!\!A_{22}\end{pmatrix}\!,\quad B=\begin{pmatrix}1&1&\!\!\vline&\!\!0\\\hline 2&1&\!\!\vline&\!\!2\\ 3&0&\!\!\vline&\!\!1\end{pmatrix}= \begin{pmatrix}B_{11}&\!\!\vline&\!\!B_{12}\\\hline B_{21}&\!\!\vline&\!\!B_{22}\end{pmatrix}\!.[/math]

Найти матрицы [math]C=A+B,~D=5B,~B^T[/math].

Решение. Матрицы [math]A[/math] и [math]B[/math] имеют блоки одинаковых размеров: блоки [math]A_{11}[/math] и [math]B_{11}[/math] имеют размеры [math]m_1\times n_1=1\times2[/math]; блоки [math]A_{12}[/math] и [math]B_{12}[/math][math]m_1\times n_2=1\times1[/math]; блоки [math]A_{21}[/math] и [math]B_{21}[/math][math]m_2\times n_1=2\times2[/math]; блоки [math]A_{22}[/math] и [math]B_{22}[/math][math]m_2\times n_2=2\times1[/math]. Матрица [math]C=A+B[/math] будет иметь такие же по размерам блоки [math]C=\begin{pmatrix}C_{11}&\!\!\vline&\!\!C_{12}\\\hline C_{21}&\!\!\vline&\!\!C_{22}\end{pmatrix}[/math]. Для каждого блока находим


[math]\begin{gathered} C_{11}=A_{11}+B_{11}= \begin{pmatrix}2&3\end{pmatrix}+ \begin{pmatrix}1&1\end{pmatrix}= \begin{pmatrix}3&4\end{pmatrix};\quad C_{12}=A_{12}+B_{12}= \begin{pmatrix}4\end{pmatrix}+ \begin{pmatrix}0\end{pmatrix}= \begin{pmatrix}4\end{pmatrix}\!;\\[3pt] C_{21}=A_{21}+B_{21}= \begin{pmatrix}3&4\\4&5\end{pmatrix}+ \begin{pmatrix}2&1\\3&0\end{pmatrix}= \begin{pmatrix}5&5\\7&5\end{pmatrix}\!;~~ C_{22}=A_{22}+B_{22}= \begin{pmatrix}5\\6\end{pmatrix}+ \begin{pmatrix}2\\1\end{pmatrix}= \begin{pmatrix}7\\7\end{pmatrix}\!.\end{gathered}[/math]

Следовательно, матрица [math]C[/math] будет следующая


[math]C=\begin{pmatrix}3&4&\!\!\vline&\!\!4\\\hline 5&5&\!\!\vline&\!\!7\\ 7&5&\!\!\vline&\!\!7\end{pmatrix}= \begin{pmatrix}C_{11}&\!\!\vline&\!\!C_{12}\\\hline C_{21}&\!\!\vline&\!\!C_{22}\end{pmatrix}\!.[/math]

Матрица [math]D=5B[/math] будет иметь блоки тех же размеров, что и [math]B[/math]:


[math]\begin{gathered}D_{11}=5B_{11}=5\cdot \begin{pmatrix}1&1\end{pmatrix}= \begin{pmatrix}5&5\end{pmatrix}\!;\quad D_{12}=5B_{12}=5\cdot \begin{pmatrix}0\end{pmatrix}= \begin{pmatrix}0\end{pmatrix}\!;\\[3pt] D_{21}=5B_{21}=5\cdot \begin{pmatrix}2&1\\3&0\end{pmatrix}= \begin{pmatrix}10&5\\15&0\end{pmatrix}\!;\quad D_{22}=5B_{22}=5\cdot \begin{pmatrix}2\\1\end{pmatrix}= \begin{pmatrix}10\\5\end{pmatrix}\!. \end{gathered}[/math]

Поэтому матрица [math]D[/math] будет иметь вид

[math]D=\begin{pmatrix}5&5&\!\!\vline&\!\!0\\\hline 10&5&\!\!\vline&\!\!10\\ 15&0&\!\!\vline&\!\!5\end{pmatrix}= \begin{pmatrix}D_{11}&\!\!\vline&\!\!D_{12}\\\hline D_{21}&\!\!\vline&\!\!D_{22}\end{pmatrix}\!.[/math]

Используя правило транспонирования блочных матриц, получаем


[math]B^T= \begin{pmatrix}B_{11}^T&\!\!\vline&\!\!B_{21}^T\\\hline B_{12}^T&\!\!\vline&\!\!B_{22}^T\end{pmatrix}= \begin{pmatrix}1&\!\!\vline&\!\!2&3\\ 1&\!\!\vline&\!\!1&0\\\hline 0&\!\!\vline&\!\!2&1\end{pmatrix}\!.[/math]



Умножение блочных матриц


Рассмотрим теперь операцию умножения блочных матриц [math]A[/math] и [math]B[/math]. Блочные матрицы [math]A[/math] и [math]B[/math] называются согласованными, если разбиение матрицы [math]A=\begin{pmatrix}A_{ik}\end{pmatrix}[/math] на блоки по столбцам совпадает с разбиением матрицы [math]B=\begin{pmatrix}B_{kj}\end{pmatrix}[/math] по строкам, т.е. блоки [math]A_{ik}[/math] имеют размеры [math]m_i\times p_k[/math], а блоки [math]B_{kj}[/math][math]p_k\times n_j[/math] [math](k=1,2,\ldots,s)[/math]. У согласованных блочных матриц блоки [math]A_{ik}[/math] и [math]B_{kj}[/math] являются согласованными матрицами.


Произведением [math]C=A\cdot B[/math] согласованных блочных матриц [math]A[/math] и [math]B[/math] называется блочная матрица [math]C=\begin{pmatrix}C_{ij}\end{pmatrix}[/math], блоки которой вычисляются по следующей формуле


[math]C_{ij}= A_{i1}\cdot B_{1j}+ A_{i2}\cdot B_{2j}+\ldots+ A_{is}\cdot B_{sj}.[/math]

Это означает, что блочные матрицы, разделенные на блоки надлежащим образом, можно перемножать обычным способом. Чтобы получить блок [math]C_{ij}[/math] произведения, надо выделить i-ю строку блоков матрицы [math]A[/math] и j-й столбец блоков матрицы [math]B[/math]. Затем найти сумму попарных произведений соответствующих блоков: первый блок i-й строки блоков умножается на первый блок j-го столбца блоков, второй блок i-й строки блоков умножается на второй блок j-го столбца и т.д., а результаты умножений складываются.




Пример 1.26. Даны блочные матрицы


[math]A= \begin{pmatrix}2&3\!\!&\vline\!\!&4\\\hline 3&4\!\!&\vline\!\!&5\\ 4&5\!\!&\vline\!\!&6\end{pmatrix}= \begin{pmatrix}A_{11}&\!\!\vline&\!\!A_{12}\\\hline A_{21}&\!\!\vline&\!\!A_{22}\end{pmatrix}\!,\quad B=\begin{pmatrix}1&1&\!\!\vline&\!\!0\\ 2&1&\!\!\vline&\!\!2\\\hline 3&0&\!\!\vline&\!\!1\end{pmatrix}= \begin{pmatrix}B_{11}&\!\!\vline&\!\!B_{12}\\\hline B_{21}&\!\!\vline&\!\!B_{22}\end{pmatrix}\!.[/math]

Найти произведение блочных матриц
[math]C=AB[/math]
.

Решение. Матрица [math]A[/math] разбита на блоки: [math]A_{11}[/math] размеров [math]m_1\times p_1=1\times2;[/math] [math]A_{12}[/math][math]m_1\times p_2=1\times1;[/math] [math]A_{21}[/math][math]m_2\times p_1=2\times2;[/math] [math]A_{22}[/math][math]m_2\times p_2=2\times1[/math]. Матрица [math]B[/math] разбита на блоки: [math]B_{11}[/math] размеров [math]p_1\times n_1=2\times2;[/math] [math]B_{12}[/math][math]p_1\times n_2=2\times1;[/math] [math]B_{21}[/math][math]p_2\times n_1=1\times2;[/math] [math]B_{22}[/math][math]p_2\times n_2=1\times1[/math]. Блочные матрицы [math]A[/math] и [math]B[/math] согласованы. Матрица [math]A[/math] разбита по столбцам на два и один (считая слева), матрица [math]B[/math] разбита по строкам на две и одну (считая сверху). Поэтому произведение [math]AB[/math] определено. Матрица [math]C=AB[/math] будет иметь блоки [math]C=\begin{pmatrix}C_{11}&\!\!\vline&\!\!C_{12}\\\hline C_{21}&\!\!\vline&\!\!C_{22}\end{pmatrix}[/math]. Для каждого блока находим


[math]\begin{gathered}C_{11}=A_{11}B_{11}+A_{12}B_{21}= \begin{pmatrix}2&3\end{pmatrix}\!\cdot\! \begin{pmatrix}1&1\\2&1\end{pmatrix}+ \begin{pmatrix}4\end{pmatrix}\!\cdot\! \begin{pmatrix}3&0\end{pmatrix}= \begin{pmatrix}8&5\end{pmatrix}+ \begin{pmatrix}12&0\end{pmatrix}= \begin{pmatrix}20&5\end{pmatrix}\!;\\[3pt] C_{12}=A_{11}B_{12}+A_{12}B_{22}= \begin{pmatrix}2&3\end{pmatrix}\!\cdot\! \begin{pmatrix}0\\2\end{pmatrix}+ \begin{pmatrix}4\end{pmatrix}\!\cdot\! \begin{pmatrix}1\end{pmatrix}= \begin{pmatrix}6\end{pmatrix}+ \begin{pmatrix}4\end{pmatrix}= \begin{pmatrix}10\end{pmatrix}\!;\\[3pt] C_{21}=A_{21}B_{11}+A_{22}B_{21}= \begin{pmatrix}3&4\\4&5\end{pmatrix}\!\! \begin{pmatrix}1&1\\2&1\end{pmatrix}+ \begin{pmatrix}5\\6\end{pmatrix}\!\! \begin{pmatrix}3&0\end{pmatrix}= \begin{pmatrix}11&7\\14&9\end{pmatrix}+ \begin{pmatrix}15&0\\18&0\end{pmatrix}= \begin{pmatrix}26&7\\32&9\end{pmatrix}\!;\\[3pt] C_{22}=A_{21}B_{12}+A_{22}B_{22}= \begin{pmatrix}3&4\\4&5\end{pmatrix}\!\cdot\! \begin{pmatrix}0\\2\end{pmatrix}+ \begin{pmatrix}5\\6\end{pmatrix}\!\cdot\! \begin{pmatrix}1\end{pmatrix}= \begin{pmatrix}8\\10\end{pmatrix}+ \begin{pmatrix}5\\6\end{pmatrix}= \begin{pmatrix}13\\16\end{pmatrix}\!.\end{gathered}[/math]

Следовательно, матрица [math]C[/math] будет иметь вид


[math]C= \begin{pmatrix}20&5\!\!&\vline\!\!&10\\\hline 26&7\!\!&\vline\!\!&13\\ 32&9\!\!&\vline\!\!&16\end{pmatrix}= \begin{pmatrix}C_{11}&\!\!\vline&\!\!C_{12}\\\hline C_{21}&\!\!\vline&\!\!C_{22}\end{pmatrix}\!.[/math]



Замечания 1.7


1. Операции сложения, умножения на число и произведения блочных матриц выполняются по тем же правилам, что и для обычных матриц, только вместо элементов в формулах используются блоки.


2. Выполняя операции над блочными матрицами, всегда можно их рассматривать как числовые, и проводить указанные операции по обычным правилам для числовых матриц. При этом результат операций (числовая матрица) будет один и тот же. Действия с блочными матрицами предпочтительнее, чем с числовыми, в том случае, когда в результате вычислений требуется искать не всю матрицу, а только ее часть — блок.


3. Матрица, у которой большинство элементов отличны от нуля, называется плотной матрицей. Матрица, большинство элементов которой -нули, называется разреженной матрицей. Для разреженных матриц, подавляющее количество элементов которых равно нулю, полезно выделять нулевые блоки с целью уменьшения вычислительных операций.




Пример 1.27. Найти произведение матриц четвёртого порядка


[math]A=\begin{pmatrix}1&0&0&0\\ 0&1&0&0\\ 1&2&1&0\\ 3&4&0&1\end{pmatrix}\!,\quad B=\begin{pmatrix}1&0&1&0\\ 0&1&0&1\\ 0&0&5&6\\ 0&0&7&8\end{pmatrix}\!.[/math]

Решение. Разобьем данные матрицы на блоки размеров [math]2\times2:[/math]


[math]A=\begin{pmatrix}1&0&\!\!\vline&\!\!0&0\\ 0&1&\!\!\vline&\!\!0&0\\\hline 1&2&\!\!\vline&\!\!1&0\\ 3&4&\!\!\vline&\!\!0&1\end{pmatrix}= \begin{pmatrix}E&\!\!\vline&\!\!O\\\hline C&\!\!\vline&\!\!E\end{pmatrix}\!,\quad B=\begin{pmatrix}1&0&\!\!\vline&\!\!1&0\\ 0&1&\!\!\vline&\!\!0&1\\\hline 0&0&\!\!\vline&\!\!5&6\\ 0&0&\!\!\vline&\!\!7&8\end{pmatrix}= \begin{pmatrix}E&\!\!\vline&\!\!E\\\hline O&\!\!\vline&\!\!D\end{pmatrix}\!.[/math]

где [math]C=\begin{pmatrix}1&2\\3&4\end{pmatrix}\!,~ D=\begin{pmatrix}5&6\\7&8\end{pmatrix}\!,~E=\begin{pmatrix}1&0\\0&1\end{pmatrix}[/math] — единичная, [math]O=\begin{pmatrix} 0&0\\0&0\end{pmatrix}[/math] — нулевая. Запишем сначала произведение блочных матриц


[math]AB= \begin{pmatrix}E&\!\!\vline&\!\!O\\\hline C&\!\!\vline&\!\!E\end{pmatrix}\!\cdot\! \begin{pmatrix}E&\!\!\vline&\!\!E\\\hline O&\!\!\vline&\!\!D\end{pmatrix}= \begin{pmatrix}E\cdot E+O\cdot O&\!\!\vline&\!\!E\cdot E+O\cdot D\\\hline C\cdot E+E\cdot O&\!\!\vline&\!\!C\cdot E+E\cdot D\end{pmatrix}= \begin{pmatrix}E&\!\!\vline&\!\!E\\\hline C&\!\!\vline&\!\!C+D\end{pmatrix}\!.[/math]

Следовательно, вместо умножения матриц [math]A[/math] и [math]B[/math] достаточно определить только один блок, сложив матрицы [math]C[/math] и [math]D:[/math]


[math]C+D= \begin{pmatrix}1&2\\3&4\end{pmatrix}+ \begin{pmatrix}5&6\\7&8\end{pmatrix}= \begin{pmatrix}6&8\\10&12\end{pmatrix}\!.[/math]

Осталось записать результат


[math]AB=\begin{pmatrix}1&0&\!\!\vline&\!\!1&0\\ 0&1&\!\!\vline&\!\!0&1\\\hline 1&2&\!\!\vline&\!\!6&8\\ 3&4&\!\!\vline&\!\!10&12\end{pmatrix}= \begin{pmatrix}1&0&1&0\\ 0&1&0&1\\ 1&2&6&8\\ 3&4&10&12\end{pmatrix}\!.[/math]



Кронекеровские произведение и сумма матриц


Пусть даны матрицы [math]A=\begin{pmatrix}a_{ij}\end{pmatrix}[/math] и [math]B[/math] размеров [math]m\times n[/math] и [math]p\times q[/math] соответственно. Числовая матрица размеров [math]mp\times nq[/math], составленная из блоков [math]a_{ij}B:[/math]


[math]A\otimes B= \begin{pmatrix}a_{11}\cdot B&\!\!\vline&\!\!\cdots&\!\!\vline&\!\!a_{1n}\cdot B\\\hline \vdots&\!\!\vline&\!\!\ddots&\!\!\vline&\!\!\vdots\\\hline a_{m1}\cdot B&\!\!\vline&\!\!\cdots&\!\!\vline&\!\!a_{mn}\cdot B \end{pmatrix}[/math]

называется правым кронекеровским произведением матриц [math]A[/math] и [math]B[/math] (или правым прямым произведением матриц).

Пусть [math]A[/math] и [math]B[/math] квадратные матрицы n-го и m-го порядков соответственно. Кронекеровскои суммой матриц [math]A[/math] и [math]B[/math] называется квадратная матрица mn-го порядка


[math]A\oplus B=\begin{pmatrix}E_{m}\otimes A\end{pmatrix}+\begin{pmatrix}B\otimes E_n\end{pmatrix}\!,[/math]

где [math]E_m,\,E_n[/math] — единичные матрицы соответствующих порядков.




Пример 1.28. Даны матрицы


[math]A=\begin{pmatrix}1&2\\3&4\end{pmatrix}\!,\quad B=\begin{pmatrix}1&0&2\\ 0&2&3 \end{pmatrix}\!,\quad C=\begin{pmatrix}5&6\\7&8\end{pmatrix}\!.[/math]

Найти кронекеровское произведение [math]A\otimes B[/math] и кронекеровскую сумму [math]A\oplus B[/math].

Решение. По определению находим


[math]A\otimes B=\begin{pmatrix}1\cdot \begin{pmatrix}1&0&2\\0&2&3\end{pmatrix}\!\!&\vline\!\!& 2\cdot \begin{pmatrix}1&0&2\\0&2&3\end{pmatrix}\\\hline 3\cdot \begin{pmatrix}1&0&2\\0&2&3\end{pmatrix}\!\!&\vline\!\!& 4\cdot \begin{pmatrix}1&0&2\\0&2&3\end{pmatrix}\end{pmatrix}= \begin{pmatrix}1&0&2&2&0&4\\ 0&2&3&0&4&6\\ 3&0&6&4&0&8\\ 0&6&9&0&8&12 \end{pmatrix}\!;[/math]

[math]\begin{gathered}A\oplus C=(E\otimes A)+(C\otimes E)= \begin{pmatrix} 1\cdot \begin{pmatrix}1&2\\3&4\end{pmatrix}\!\!&\vline\!\!& 0\cdot \begin{pmatrix} 1&2\\ 3&4\end{pmatrix}\\\hline 0\cdot \begin{pmatrix}1&2\\3&4 \end{pmatrix}\!\!&\vline\!\!& 1\cdot \begin{pmatrix}1&2\\3&4\end{pmatrix}\end{pmatrix}+ \begin{pmatrix}5\cdot \begin{pmatrix}1&2\\ 3&4\end{pmatrix}\!\!&\vline\!\!& 6\cdot \begin{pmatrix} 1&2\\3&4\end{pmatrix}\\ \hline 7\cdot \begin{pmatrix}1&2\\3&4 \end{pmatrix}\!\!&\vline\!\!& 8\cdot \begin{pmatrix}1&2\\ 3&4 \end{pmatrix}\end{pmatrix}=\\[2pt] =\begin{pmatrix}1&2&0&0\\3&4&0&0\\ 0&0&1&2\\ 0&0&3&4\end{pmatrix}+ \begin{pmatrix}5&0&6&0\\0&5&0&6\\ 7&0&8&0\\ 0&7&0&8\end{pmatrix}= \begin{pmatrix}6&2&6&0\\ 3&9&0&6\\ 7&0&9&2\\ 0&7&3&12 \end{pmatrix}\!.\end{gathered}[/math]

Часовой пояс: UTC + 4 часа [ Летнее время ]


Яндекс.Метрика

Copyright © 2010-2016 MathHelpPlanet.com. All rights reserved