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

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

Теоретический раздел
Часовой пояс: UTC + 3 часа [ Летнее время ]
новый онлайн-сервис
число, сумма и дата прописью

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


Псевдообратная матрица

Псевдообратная матрица


Обратная матрица в отличие от полуобратной имеет в силу определения очевидные свойства:


\begin{pmatrix}A\cdot A^{-1}\end{pmatrix}^{\ast}=A\cdot A^{-1},\qquad \begin{pmatrix}A^{-1}\cdot A\end{pmatrix}^{\ast}=A^{-1}\cdot A,

так как единичная матрица E=AA^{-1}=A^{-1}A, разумеется, эрмитова. В определении полуобратной матрицы имеется некоторый произвол (см. п.3 замечаний 4.6), которым можно воспользоваться так, чтобы полуобратная матрица обладала аналогичными свойствами.


Пусть A — произвольная матрица размеров m\times n. Полуобратная матрица A^{\sim1} размеров n\times m называется псевдообратной для матрицы A, если матрицы AA^{\sim1} и A^{\sim1}A эрмитовы, т.е. псевдообратная матрица A^{\sim1} определяется четырьмя условиями


AA^{\sim1}A=A;\qquad A^{\sim1}AA^{\sim1}=A^{\sim1};
(4.17)

\begin{pmatrix}AA^{\sim1}\end{pmatrix}^{\ast}=AA^{\ast};\qquad \begin{pmatrix}A^{\sim1}A\end{pmatrix}^{\ast}=A^{\sim1}A.
(4.18)

Покажем, что псевдообратная матрица A^{\sim1} существует для любой матрицы A. Действительно, если A=O — нулевая матрица размеров m\times n, то A^{\sim1}=O^T — нулевая размеров n\times m, что следует из равенств (4.17).


Пусть матрица A — ненулевая. Тогда матрица A^{\sim1}, удовлетворяющая равенствам (4.17), имеет вид (4.14):


A^{\sim1}=T\cdot\! \begin{pmatrix}\dfrac{E_r}{U}\end{pmatrix}\!\cdot\! \begin{pmatrix}E_r\mid V\end{pmatrix}\!\cdot S.
(4.19)

Покажем, что выбором матриц U и V в формуле (4.19) можно получить матрицу, удовлетворяющую условиям (4.18). В самом деле, запишем скелетное разложение (4.10) матрицы A:


A=S^{-1}\Lambda T^{-1}= S^{-1}\! \begin{pmatrix}E_r\!\!&\vline\!\!&O\\\hline O\!\!&\vline\!\!&O\end{pmatrix}\!T^{-1}= S^{-1}\cdot\! \begin{pmatrix}\dfrac{E_r}{O} \end{pmatrix}\!\cdot\! \begin{pmatrix}E_r\mid O\end{pmatrix}\!\cdot T^{-1}.

Найдем произведение


\begin{gathered}AA^{-1}= S^{-1}\cdot\!\begin{pmatrix} \dfrac{E_r}{O} \end{pmatrix}\!\cdot\! \begin{pmatrix}E_r\mid O\end{pmatrix}\!\cdot\underbrace{T^{-1}\cdot T}_{E_n}\cdot\! \begin{pmatrix}\dfrac{E_r}{U} \end{pmatrix}\!\cdot\! \begin{pmatrix}E_r\mid U\end{pmatrix}\!\cdot S=\\[2pt] =S^{-1}\cdot\! \begin{pmatrix}\dfrac{E_r}{O}\end{pmatrix}\!\cdot \underbrace{\begin{pmatrix}E_r\mid O\end{pmatrix}\!\cdot \begin{pmatrix}\dfrac{E_r}{U} \end{pmatrix}}_{E_r}\cdot\! \begin{pmatrix} E_r\mid V\end{pmatrix}\!\cdot S=S^{-1}\cdot\! \begin{pmatrix}\dfrac{E_r}{O} \end{pmatrix}\!\cdot\! \begin{pmatrix}E_r\mid V\end{pmatrix}\!\cdot S= S^{-1}\cdot\! \begin{pmatrix} E_r\!\!&\vline\!\!&V\\ \hline O\!\!&\vline\!\!& O\end{pmatrix}\!\cdot S.\end{gathered}

Подставим его в первое из равенств (4.18):


\left[S^{-1}\cdot\! \begin{pmatrix}E_r\!\!&\vline\!\!&V\\ \hline O\!\!&\vline\!\!& O\end{pmatrix}\!\cdot S\right]^{\ast}= S^{-1}\cdot\! \begin{pmatrix}E_r\!\!&\vline\!\!&V\\ \hline O\!\!&\vline\!\!& O\end{pmatrix}\!\cdot S.

Используя свойства операции сопряжения, а также п. 1 замечаний 4.2, получаем


S^{\ast}\cdot\! \begin{pmatrix} E_r\!\!&\vline\!\!&O\\ \hline V^{\ast}\!\!&\vline\!\!& O \end{pmatrix}\!\cdot\! \begin{pmatrix}S^{-1}\end{pmatrix}^{\ast}= S^{-1}\cdot\! \begin{pmatrix} E_r\!\!&\vline\!\!&V\\ \hline O\!\!&\vline\!\!& O \end{pmatrix}\!\cdot S.

Умножая обе части равенства на матрицу S слева и на матрицу S^{\ast} справа, приходим к равенству


\begin{pmatrix}SS^{\ast}\end{pmatrix}\!\cdot\! \begin{pmatrix} E_r\!\!&\vline\!\!&O\\ \hline V^{\ast}\!\!&\vline\!\!& O \end{pmatrix}= \begin{pmatrix} E_r\!\!&\vline\!\!&V\\ \hline O\!\!&\vline\!\!& O \end{pmatrix}\!\cdot\! \begin{pmatrix}SS^{\ast}\end{pmatrix}\!.

Подставим в это равенство матрицу SS^{\ast}, предварительно разбив ее на блоки SS^{\ast}=\begin{pmatrix} S_1\!\!&\vline\!\!&S_2\\ \hline S_3\!\!&\vline\!\!& S_4 \end{pmatrix} квадратными матрицами S_1 и S_4 порядков r и (m-r) и прямоугольными матрицами S_2 и S_2 размеров r\times(m-r) и (m-r)\times r соответственно. Выполняя умножение блочных матриц, получаем


\begin{gathered} \begin{pmatrix} S_1\!\!&\vline\!\!&S_2\\ \hline S_3\!\!&\vline\!\!& S_4\end{pmatrix}\!\cdot\! \begin{pmatrix} E_r\!\!&\vline\!\!&O\\ \hline V^{\ast}\!\!&\vline\!\!& O\end{pmatrix}= \begin{pmatrix} S_1+S_2V^{\ast}\!\!&\vline\!\!&O\\ \hline S_3+S_4V^{\ast}\!\!&\vline\!\!& O\end{pmatrix}\!,\\[2pt] \begin{pmatrix} E_r\!\!&\vline\!\!&V\\ \hline O\!\!&\vline\!\!& O\end{pmatrix}\!\cdot\! \begin{pmatrix} S_1\!\!&\vline\!\!&S_2\\ \hline S_3\!\!&\vline\!\!& S_4\end{pmatrix}= \begin{pmatrix} S_1+VS_3\!\!&\vline\!\!&S_2+VS_4\\ \hline O\!\!&\vline\!\!& O\end{pmatrix}\!.\end{gathered}

Равенство полученных блочных матриц обеспечивается условием


V=-S_2S_{4}^{-1}, поскольку S_3=S_2^{\ast},~S_1=S_1^{\ast},~ S_4=S_4^{\ast}

в силу эрмитовости матрицы SS^{\ast}, а \begin{pmatrix}S_4^{-1}\end{pmatrix}^{\ast}= \begin{pmatrix}S_4^{\ast}\end{pmatrix}^{-1}.


Аналогичным образом можно показать, что второе из равенств (4.18) выполняется, если положить


U=-T_4^{-1}T_3,

где T_3,T_4 — блоки размеров (n-r)\times r и (n-r)\times(n-r) матрицы T^{\ast}T=\begin{pmatrix} T_1\!\!&\vline\!\!&T_2\\ \hline T_3\!\!&\vline\!\!& T_4\end{pmatrix}.


Таким образом, для любой матрицы существует псевдообратная матрица и притом только одна.




Замечания 4.7


1. Если матрица A обратимая, то обратная матрица A^{-1}, как следует из п. 1 замечаний 4.6, совпадает с псевдообратной, т.е. A^{-1}=A^{\sim1}.


2. Из невырожденности матриц S и T следует, что при любом раз биении эрмитовых матриц


SS^{\ast}=\begin{pmatrix} S_1\!\!&\vline\!\!&S_2\\ \hline S_3\!\!&\vline\!\!& S_4\end{pmatrix}!.,\qquad T^{\ast}T=\begin{pmatrix} T_1\!\!&\vline\!\!&T_2\\ \hline T_3\!\!&\vline\!\!& T_4\end{pmatrix}\!.

на квадратные блоки S_1,\,S_4,\,T_1,\,T_4, существуют обратные матрицы S_1^{-1},\,S_4^{-1},\,T_1^{-1},\,T_4^{-1}.


3. Имеются другие определения псевдообратной матрицы, равносильные приведенному выше. Например,


A^{\sim1}= \lim_{\varepsilon\to0}\Bigl(A^{\ast}A+\varepsilon^2\cdot E\Bigr)^{-1}A^{\ast}= \lim_{\varepsilon\to0}A^{\ast}\Bigl(A^{\ast}A+\varepsilon^2\cdot E\Bigr)^{-1}.

4. В общем случае произведение псевдообратных матриц некоммутативно (AB)^{\sim1}\ne B^{\sim1}\cdot A^{\sim1}.




Свойства псевдообратной матрицы


Операция псевдообращения матриц обладает следующими свойствами:


\begin{aligned} \bold{1.}&~~ (A^{\sim1})^{\sim1}=A;\\[2pt] \bold{2.}&~~ (A^{\ast})^{\sim1}=(A^{\sim1})^{\ast};\\[2pt] \bold{3.}&~~ A^{\sim1}=(A^{\ast}A)^{\sim1}A^{\ast}=A^{\ast}(AA^{\ast})^{\sim1};\\[2pt] \bold{4.}&~~ (AA^{\sim1})^2=AA^{\sim1};\\[2pt] \bold{5.}&~~ (A^{\sim1}A)^2=A^{\sim1}A. \end{aligned}


Эти свойства доказываются по определению (4.17), (4.18). Докажем, например, свойство 3 (первое равенство). По определению псевдообратной матрицы имеем:


\Bigl[(A^{\ast}A)^{\sim1}A^{\ast}\Bigr]\cdot A\cdot\Bigl[(A^{\ast}A)^{\sim1}A^{\ast}\Bigr]= (A^{\ast}A)^{\sim1}\cdot(A^{\ast}A)\cdot(A^{\ast}A)^{\sim1}A^{\ast}= (A^{\ast}A)^{\sim1}A^{\ast}.

Следовательно, A=\bigl[(A^{\ast}A)^{\sim1}A^{\ast}\bigr]^{\sim1}. Тогда по свойству 1: A^{\sim1}=(A^{\ast}A)^{\sim1}A^{\ast}.




Способы нахождения псевдообратной матрицы


Пусть дана ненулевая матрица A размеров m\times n. Требуется найти псевдообратную матрицу A^{\sim1}.


Первый способ. Для нахождения псевдообратной матрицы (4.19) нужно выполнить следующие действия.


1. Составить блочную матрицу \begin{pmatrix}A\!\!&\vline\!\!&E_m\\\hline E_n\!\!&\vline\!\!&{} \end{pmatrix}, приписывая к матрице A слева и снизу единичные матрицы соответствующих размеров. Правый нижний блок этой матрицы может быть произвольным, так как не участвует в дальнейших преобразованиях.


2. Элементарными преобразованиями над первыми m строками и первыми n столбцами привести блочную матрицу к виду \begin{pmatrix}\Lambda\!\!& \vline\!\!&S\\\hline T\!\!&\vline\!\!&{}\end{pmatrix}, где \Lambda — матрица размеров m\times n простейшего вида (4.8), т.е. \Lambda=\begin{pmatrix} E_r\!\!&\vline\!\!&O\\\hline O\!\!&\vline\!\!&O\end{pmatrix}, в которой E_r — единичная матрица r-го порядка (1\leqslant r\leqslant\min\{m;n\}), O — нулевые матрицы соответствующих размеров.


3. Найти произведения SS^{\ast},\,T^{\ast}T и представить их в виде блочных матриц


SS^{\ast}=\begin{pmatrix}S_1\!\!&\vline\!\!&S_2\\\hline S_3\!\!&\vline\!\!&S_4\end{pmatrix}\!,\quad T^{\ast}T=\begin{pmatrix}T_1\!\!&\vline\!\!&T_2\\\hline T_3\!\!&\vline\!\!&T_4\end{pmatrix}\!,
(4.20)

выделяя блоки S_2,\,S_4,\,T_3,\,T_4 размеров r\times(m-r),\,(m-r)\times(m-r),\,(n-r)\times r,\,(n-r)\times(n-r) соответственно.


4. Вычислить матрицы U=-T_4^{-1}T_3,~V=-S_2S_4^{-1}..


5. Получить псевдообратную матрицу


A^{\sim1}=T\cdot\! \begin{pmatrix}\dfrac{E_r}{U}\end{pmatrix}\!\cdot\! \begin{pmatrix}E_r\mid V\end{pmatrix}\!\cdot S.
(4.21)

Замечание 4.8. Если r=m или r=n в (4.20) будут отсутствовать блоки S_2,\,S_4 или T_3,\,T_4, так как в формуле (4.21) будут отсутствовать блоки V или U соответственно (см. п.5 замечаний 4.6). В частных случаях, когда строки матрицы A линейно независимы (r=m) или столбцы матрицы A линейно независимы (r=n), псевдообратную матрицу A^{\sim1} можно найти проще, чем описанным выше способом (см. далее частные случаи нахождения псевдообратной матрицы).


Второй способ. Для нахождения псевдообратной матрицы используем скелетное разложение (4.10).


1, 2. Выполнить два первых пункта, указанных в первом способе. Получим матрицы S,\,T и \Lambda, удовлетворяющие условию


\Lambda=S\cdot A\cdot T,

где S,\,T — элементарные преобразующие матрицы порядков тип соответственно; \Lambda=\begin{pmatrix}E_r\!\!&\vline\!\!&O\\\hline O\!\!&\vline\!\!&O\end{pmatrix} — матрица простейшего вида; E_r — единичная r-го порядка (1\leqslant r\leqslant \min\{m;n\}), O — нулевые матрицы соответствующих размеров.


3. Найти обратные матрицы S^{-1} и T^{-1}.


4. Записать матрицы B=S^{-1}\cdot\! \begin{pmatrix}\dfrac{E_r}{O}\end{pmatrix} и C=\begin{pmatrix}E_r\mid O\end{pmatrix}\!\cdot T^{-1}. Матрица В составлена из первых r столбцов матрицы S^{-1}, а матрица C — из первых r строк матрицы T^{-1}.


5. Получить псевдообратную матрицу по формуле


A^{\sim1}=C^{\ast}\cdot(CC^{\ast})^{-1}\cdot(B^{\ast}B)^{-1}\cdot B^{\ast}.
(4.22)

Докажем, что второй способ дает ту же псевдообратную матрицу, что и первый. В самом деле, имея в виду скелетное разложение матицы A:


A=BC;\quad B=S^{-1}\cdot\! \begin{pmatrix}\dfrac{E_r}{O}\end{pmatrix}\!,\quad C=\begin{pmatrix}E_e\mid O\end{pmatrix}\!\cdot T^{-1},
найдем произведение
\begin{gathered}\begin{pmatrix}B^{\ast}B\end{pmatrix}^{-1}\cdot B^{\ast}\cdot S^{-1}= \left[\begin{pmatrix}E_r&\mid O\end{pmatrix}\!\cdot(S^{-1})^{\ast}\cdot S^{-1}\cdot\! \begin{pmatrix} \dfrac{E_r}{O}\end{pmatrix}\right]^{-1}\cdot\! \begin{pmatrix}E_r\mid O\end{pmatrix} \!\cdot (S^{-1})^{\ast}\cdot S^{-1}=\\[2pt] =\left[\begin{pmatrix}E_r\mid O\end{pmatrix}\!\cdot\! \begin{pmatrix} SS^{\ast}\end{pmatrix}^{-1}\!\cdot\! \begin{pmatrix} \dfrac{E_r}{O}\end{pmatrix} \right]^{-1}\cdot\! \begin{pmatrix}E_r\mid O\end{pmatrix}\!\cdot\! \begin{pmatrix} SS^{\ast}\end{pmatrix}^{-1} \end{gathered}

Обращая блочную матрицу SS^{\ast}=\begin{pmatrix}S_1\!\!&\vline\!\!&S_2\\\hline S_3\!\!&\vline\!\!&S_4\end{pmatrix} по формуле Фробениуса (4.4), имеем


\begin{pmatrix} S_1\!\!&\vline\!\!&S_2\\\hline S_3\!\!&\vline\!\!&S_4\end{pmatrix}^{-1}= \begin{pmatrix} M\!\!&\vline\!\!&-MS_2S_4^{-1}\\\hline -S_4^{-1}S_3M\!\!& \vline\!\!&N \end{pmatrix}\!,

где M=\begin{pmatrix}S_1-S_2S_4^{-1}S_3\end{pmatrix}^{-1},~ N=\begin{pmatrix}S_4-S_3S_1^{-1}S_2\end{pmatrix}^{-1}. Подставляя (SS^{\ast})^{-1} и няя умножение блочных матриц в квадратных скобках, имеем


\begin{gathered}\begin{pmatrix} B^{\ast}B\end{pmatrix}^{-1}\cdot B^{\ast}\cdot S^{-1}=M^{-1}\cdot\! \begin{pmatrix}E_r\mid O\end{pmatrix}\!\cdot\! \begin{pmatrix} M\!\!&\vline\!\!&-MS_2S_4^{-1}\\\hline -S_4^{-1}S_3M\!\!&\vline\!\!&N\end{pmatrix}=\\[2pt] M^{-1}\cdot\! \begin{pmatrix}M\mid -MS_2S_4^{-1}\end{pmatrix}= \begin{pmatrix}E_r\mid -S_2S_4^{-1}\end{pmatrix}= \begin{pmatrix}E_r\mid V\end{pmatrix}\!.\end{gathered}

Таким образом, \begin{pmatrix}E_r\mid V\end{pmatrix}\!S, где V=-S_2S_4^{-1}.


Аналогично показывается, что T\cdot\! \begin{pmatrix}\dfrac{E_r}{U}\end{pmatrix}= C^{\ast}(CC^{\ast})^{-1} для U=-T_4^{-1}T_3. Поэтому формулы (4.22) и (4.21) дают одну и ту же псевдообратную матрицу.




Пример 4.11. Для матрицы A=\begin{pmatrix}1&-1&2&0\\-1&2&-3&1\\ 0&1&-1&1 \end{pmatrix} найти псевдообратную.


Решение. Первый способ. 1, 2. Первые два пункта алгоритма выполнены при решении примера 1.38, где получены матрицы


\Lambda=\begin{pmatrix}1&0\!\!&\vline\!\!&0&0\\ 0&1\!\!&\vline\!\!&0&0\\\hline 0&0\!\!&\vline\!\!&0&0\end{pmatrix}\!,\quad S=\begin{pmatrix}1&0&0\\1&1&0\\-1&-1&1 \end{pmatrix}\!,\quad T=\begin{pmatrix}1&1&-1&-1\\ 0&1&1&-1\\ 0&0&1&0\\ 0&0&0&1 \end{pmatrix}\!,

удовлетворяющие равенству \Lambda=SAT. Следовательно, r=\operatorname{rg}A= \operatorname{rg}\Lambda=2.


3. Найдем произведения


\begin{gathered}SS^{\ast}= \begin{pmatrix}1&0&0\\1&1&0\\-1&-1&1\end{pmatrix} \!\cdot\! \begin{pmatrix}1&1&-1\\0&1&-1\\ 0&0&1 \end{pmatrix}= \begin{pmatrix}1&1\!\!&\vline\!\!&-1\\ 1&2\!\!&\vline\!\!&-2\\\hline -1&-2\!\!&\vline\!\!&3\end{pmatrix}\!;\\[5pt] TT^{\ast}= \begin{pmatrix}1&0&0&0\\1&1&0&0\\-1&1&1&0\\-1&-1&0&1\end{pmatrix}\!\cdot\! \begin{pmatrix}1&1&-1&-1\\ 0&1&1&-1\\ 0&0&1&0\\ 0&0&0&1\end{pmatrix}= \begin{pmatrix} 1&1\!\!&\vline\!\!&-1&-1\\ 1&2\!\!&\vline\!\!&0&-2\\\hline -1&0\!\!&\vline\!\!&3&0\\ -1&-2\!\!&\vline\!\!&0&3 \end{pmatrix}\end{gathered}

Следовательно, S_2=\begin{pmatrix}-1\\-2\end{pmatrix}\!,~S_4=(3),~T_3=\begin{pmatrix}-1&0\\ -1&-2\end{pmatrix}\!,~T_4=\begin{pmatrix}3&0\\0&3\end{pmatrix}

.

4. Находим произведения


\begin{gathered}U=-T_4^{-1}T_3= -\frac{1}{9}\! \begin{pmatrix} 3&0\\0&3\end{pmatrix}\!\cdot\! \begin{pmatrix}-1&0\\ -1&-2\end{pmatrix}= \begin{pmatrix} 1/3&0\\1/3&2/3\end{pmatrix}\!,\\[5pt] V=-S_2S_4^{-1}=-\begin{pmatrix}-1\\-2\end{pmatrix}\!\cdot\! \begin{pmatrix}\dfrac{1}{3}\end{pmatrix}= \begin{pmatrix}1/3\\2/3\end{pmatrix}\!. \end{gathered}

5. По формуле (4.21) получим псевдообратную матрицу


\begin{aligned}A^{-1}&= T\cdot\! \begin{pmatrix}\dfrac{E_r}{U}\end{pmatrix}\!\cdot\! \begin{pmatrix}E_r\mid V\end{pmatrix}\!\cdot S= \begin{pmatrix}1&1&-1&-1\\ 0&1&1&-1\\ 0&0&1&0\\ 0&0&0&1\end{pmatrix}\!\cdot\! \begin{pmatrix}1&0\\0&1\\\hline 1/3&0\\1/3& 2/3 \end{pmatrix}\!\cdot\! \begin{pmatrix}1&0\!\!&\vline\!\!&1/3\\ 0&1\!\!&\vline\!\!&2/3 \end{pmatrix} \!\cdot\! \begin{pmatrix} 1&0&0\\1&1&0\\ -1&-1&1\end{pmatrix}=\\[2pt] &=\begin{pmatrix} 1/3&1/3\\ 0&1/3\\1/3&0\\1/3&2/3\end{pmatrix}\!\cdot\! \begin{pmatrix}2/3&-1/3&1/3\\ 1/3&1/3&2/3 \end{pmatrix}= \begin{pmatrix}1/3&0&1/3\\ 1/9&1/9&2/9\\1/9&1/9&2/9\\2/9&-1/9&1/9\\ 4/9&1/9&5/9 \end{pmatrix}\!.\end{aligned}

Псевдообратная матрица найдена.


Второй способ. 1, 2. Первые два пункта алгоритма выполнены при решении примера 1.38.


3. Найдем обратные матрицы


S^{-1}=\begin{pmatrix}1&0&0\\1&1&0\\-1&-1&0\\-1&-1&1\end{pmatrix}^{-1}= \begin{pmatrix}1&0&0\\-1&1&0\\0&1&1\end{pmatrix}\!;\quad T^{-1}=\begin{pmatrix}1&1&-1&-1\\0&1&1&-1\\0&0&1&0\\0&0&0&1\end{pmatrix}^{-1}= \begin{pmatrix}1&-1&2&0\\ 0&1&-1&1\\ 0&0&1&0\\ 0&0&0&1\end{pmatrix}\!.

4. Запишем матрицы

\begin{gathered}B=S^{-1}\cdot\! \begin{pmatrix}\dfrac{E_r}{O}\end{pmatrix}= \begin{pmatrix}1&0&0\\-1&1&0\\ 0&1&1\end{pmatrix}\!\cdot\! \begin{pmatrix}1&0\\0&1\\\hline 0&0\end{pmatrix}= \begin{pmatrix}1&0\\-1&1\\0&1\end{pmatrix}\!,\\[2pt]  C=\begin{pmatrix}E_r\mid O\end{pmatrix}\!\cdot T^{-1}= \begin{pmatrix}1&0\!\!&\vline\!\!&0&0\\ 0&1\!\!&\vline\!\!&0&0 \end{pmatrix}\!\cdot\! \begin{pmatrix}1&-1&2&0\\0&1&-1&1\\ 0&0&1&0\\ 0&0&0&1\end{pmatrix}= \begin{pmatrix}1&-1&2&0\\0&1&-1&1\end{pmatrix}\!.\end{gathered}

5. Псевдообратную матрицу получим по формуле (4.22). Для этого сначала найдем произведения и обратные матрицы


\begin{aligned}B^{\ast}B&= \begin{pmatrix}1&-1&0\\0&1&1\end{pmatrix}\!\cdot\! \begin{pmatrix}1&0\\-1&1\\0&1\end{pmatrix}= \begin{pmatrix}2&-1\\-1&2\end{pmatrix}\!,\quad \begin{pmatrix}B^{\ast}B\end{pmatrix}^{-1}= \frac{1}{3}\! \begin{pmatrix}2&1\\1&2\end{pmatrix}= \begin{pmatrix}2/3&1/3\\1/3&2/3\end{pmatrix}\!;\\[5pt] CC^{\ast}&= \begin{pmatrix}1&-1&2&0\\0&1&-1&1\end{pmatrix}\!\cdot\! \begin{pmatrix}1&0\\-1&1\\2&-1\\0&1\end{pmatrix}= \begin{pmatrix}6&-3\\-3&3\end{pmatrix}\!,\quad \begin{pmatrix}CC^{\ast}\end{pmatrix}^{-1}= \frac{1}{9}\! \begin{pmatrix} 3&3\\3&6\end{pmatrix}= \begin{pmatrix}1/3&1/3\\1/3&2/3\end{pmatrix}\!, \end{aligned}

а затем псевдообратную матрицу


\begin{gathered}A^{-1}= C^{\ast}\begin{pmatrix}CC^{\ast}\end{pmatrix}^{-1}\!\begin{pmatrix}B^{\ast}B\end{pmatrix}^{-1}B^{\ast}= \begin{pmatrix}1&0\\-1&1\\2&-1\\0&1\end{pmatrix}\!\! \begin{pmatrix}1/3&1/3\\1/3&2/3\end{pmatrix}\!\cdot\! \begin{pmatrix} 2/3&1/3\\1/3&2/3 \end{pmatrix}\!\! \begin{pmatrix}1&-1&0\\0&1&1\end{pmatrix}=\\[2pt] = \begin{pmatrix}1/3&1/3\\0&1/3\\1/3&0\\1/3&2/3\end{pmatrix}\!\cdot\! \begin{pmatrix}2/3&-1/3&1/3\\ 1/3&1/3&2/3\end{pmatrix}= \begin{pmatrix}1/3&0&1/3\\ 1/9&1/9&2/9\\ 2/9&-1/9&1/9\\ 4/9&1/9& 5/9\end{pmatrix}\!,\end{gathered}

что совпадает с полученным ранее результатом.




Третий способ. Для последовательного нахождения строк псевдообратной матрицы используется метод Гревилля.


Пусть a_k — k-й столбец матрицы A\,(1\leqslant k\leqslant n); A_k — матрица, образованная первыми к столбцами матрицы A (A_1=a_1,~A_n=A); A_k^{-1} — соответствующие псевдообратные матрицы (k=1,2,\ldots,n); b_k — последняя строка в матрице A_k^{-1} ( k=1,2,\ldots,n); o — как и ранее, нулевой столбец.


1. Положить k=1. Если a_1\ne o, найти A_1^{-1}=a_1^{-1}=\frac{1}{a_1^{\ast}a_1}\cdot a_1^{\ast}. Если a_1=o, то A_1^{-1}=o^T.


2. Положить k=k+1. Найти:


а) d_k=A_{k-1}^{-1}\cdot a_k;


б) c_k= a_k-A_{k-1}d_k;


в) если c_k\ne o,то b_k=c_k^{-1}=(c_k^{\ast}c_k)^{-1}c_k^{\ast}; если c_k=o,то b_k=(1+d_k^{\ast}d_k)^{-1}d_k^{\ast}A_{k-1}^{-1};


г) B_k=A_{k-1}^{-1}-d_kb_k;


д) A_k^{-1}=\begin{pmatrix}\dfrac{B_k}{b_k}\end{pmatrix}.


3. Если k=n, процесс завершить: A^{-1}=A_k^{-1}. Иначе перейти к п.2.




Пример 4.12. Для матрицы A=\begin{pmatrix}1&-1&0\\-1&2&1\\ 2&-3&-1\\0&-1&1 \end{pmatrix} найти псевдообратную методом Гревилля.


Решение. Так как данная матрица действительная, операция сопряжения совпадает с операцией транспонирования.


1. Положим k=1. Так как первый столбец a_1=\begin{pmatrix}1&-1&2&0 \end{pmatrix}^T ненулевой, то


A_1^{-1}= \frac{1}{a_1^Ta_1}\cdot a_1^T= \frac{1}{6}\cdot \begin{pmatrix}1&-1&2&0 \end{pmatrix}= \begin{pmatrix}\dfrac{1}{6}&-\dfrac{1}{6}&\dfrac{1}{3}&0 \end{pmatrix}\!.

2. Положим k=k+1=2.


a) d_2=A_1^{-1}\cdot a_2= \begin{pmatrix}\dfrac{1}{6}&-\dfrac{1}{6}&\dfrac{1}{3}&0 \end{pmatrix}\!\cdot\! \begin{pmatrix}-1\\2\\-3\\1\end{pmatrix}= -\frac{3}{2}.;


б) c_2=a_2-A_1d_2= \begin{pmatrix}-1\\2\\-3\\1\end{pmatrix}- \begin{pmatrix}1\\-1\\ 2\\0 \end{pmatrix}\!\cdot\! \begin{pmatrix}-\dfrac{3}{2}\end{pmatrix}= \begin{pmatrix}1/2\\1/2\\0\\1 \end{pmatrix};


в) так как c_2\ne o, то


b_2=c_2^{-1}= \begin{pmatrix} c_2^Tc_2\end{pmatrix}^{-1}c_2^{T}= \left[\begin{pmatrix}\dfrac{1}{2}&\dfrac{1}{2}&0&1\end{pmatrix}\!\! \begin{pmatrix}1/2\\1/2\\0\\1 \end{pmatrix}\right]^{-1}\!\begin{pmatrix}\dfrac{1}{2}&\dfrac{1}{2}& 0&1 \end{pmatrix}= \dfrac{2}{3}\! \begin{pmatrix} \dfrac{1}{2}&\dfrac{1}{2}&0&1\end{pmatrix}= \begin{pmatrix} \dfrac{1}{3}&\dfrac{1}{3}&0&\dfrac{2}{3}\end{pmatrix}\!;

г) B_2=A_1^{-1}-d_2b_2= \begin{pmatrix}\dfrac{1}{6}&-\dfrac{1}{6}& \dfrac{1}{3}& 0\end{pmatrix}- \begin{pmatrix}-\dfrac{3}{2}\end{pmatrix}\!\! \begin{pmatrix}\dfrac{1}{3}& \dfrac{1}{3}&0&\dfrac{2}{3}\end{pmatrix}= \begin{pmatrix}\dfrac{2}{3}&\dfrac{1}{3}& \dfrac{1}{3}&1\end{pmatrix}.


д) A_2^{-1}= \begin{pmatrix}\dfrac{B_2}{b_2}\end{pmatrix}= \begin{pmatrix}2/3&1/3&1/3&1\\ 1/3&1/3&0&2/3\end{pmatrix}.


3. Так как k=2\ne n=3, то перейдем к п.2.


2. Положим k=k+1=3.


а) d_3=A_2^{-1}\cdot a_3= \begin{pmatrix}2/3&1/3&1/3&1\\ 1/3&1/3&0&2/3 \end{pmatrix} \!\cdot\! \begin{pmatrix}0\\1\\-1\\1\end{pmatrix}= \begin{pmatrix}1\\1\end{pmatrix};


б) c_3=a_3-A_2d_3= \begin{pmatrix}0\\1\\-1\\1\end{pmatrix}- \begin{pmatrix}1&-1\\ -1&2\\ 2&-3\\ 0&1\end{pmatrix}\!\cdot\! \begin{pmatrix}1\\1\end{pmatrix}= \begin{pmatrix}0\\1\\-1\\1\end{pmatrix}- \begin{pmatrix}0\\1\\-1\\1\end{pmatrix}= \begin{pmatrix}0\\0\\0\\0\end{pmatrix}\!;


в) так как c_3=o, то


\begin{gathered}b_3=\begin{pmatrix}1+d_3^Td_3\end{pmatrix}^{-1}d_3^TA_2^{-1}= \begin{bmatrix}1+ \begin{pmatrix}1&1 \end{pmatrix}\!\cdot\! \begin{pmatrix}1\\1 \end{pmatrix} \end{bmatrix}^{-1}\!\cdot\!\begin{pmatrix}1&1\end{pmatrix}\!\cdot\! \begin{pmatrix} 2/3&1/3&1/3&1\\ 1/3&1/3&0&2/3\end{pmatrix}=\\[2pt] =\begin{pmatrix} \dfrac{1}{3}&\dfrac{1}{3}\end{pmatrix}\!\cdot\! \begin{pmatrix} 2/3&1/3&1/3&1\\ 1/3&1/3&0&2/3\end{pmatrix}= \begin{pmatrix}\dfrac{1}{3}&\dfrac{2}{9}&\dfrac{1}{9}&\dfrac{5}{9} \end{pmatrix}\!; \end{gathered}

г) B_3=A_2^{-1}-d_3b_3= \begin{pmatrix}2/3&1/3&1/3&1\\ 1/3&1/3&0&2/3\end{pmatrix}- \begin{pmatrix} 1\\1\end{pmatrix}\!\cdot\! \begin{pmatrix}\dfrac{1}{3}&\dfrac{2}{9}&\dfrac{1}{9}& \dfrac{5}{9}\end{pmatrix}= \begin{pmatrix}1/3&1/9&2/9&4/9\\ 0&1/9&-1/9&1/9\end{pmatrix}\!;


д) A_3^{\sim1}= \begin{pmatrix}\dfrac{B_3}{b_3}\end{pmatrix}= \begin{pmatrix} 1/3&1/9&2/9&4/9\\ 0&1/9&-1/9&1/9\\ 1/3&2/9&1/9&5/9\end{pmatrix}\!.


3(2). Так как k=3=n, то процесс завершен и A^{\sim1}=A_3^{\sim1}.




Частные случаи нахождения псевдообратной матрицы


1. Если матрица A=(a_{11}) — число, то


A^{\sim1}= \begin{cases}\dfrac{1}{a_{11}},&a_{11}\ne0,\\[4pt] 0,&a_{11}=0. \end{cases}

2. Если матрица A диагональная A=\operatorname{diag}\!\begin{pmatrix}a_{11},a_{22},\ldots,a_{nn}\end{pmatrix}, то псевдообратная матрицы также диагональная


A^{\sim1}=\operatorname{diag}\begin{pmatrix}a_{11}^{-1}, a_{22}^{-1},\ldots, a_{nn}^{-1}\end{pmatrix},~~\text{where}~~a_{ii}^{-1}=\begin{cases}\dfrac{1}{a_{ii}},&a_{ii}\ne0,\\0,&a_{ii}=0,\end{cases}i=1,2,\ldots,n.
(4.23)

3. Если столбцы матрицы A линейно независимы, то


A^{\sim1}=\begin{pmatrix}A^{\ast}\cdot A\end{pmatrix}^{-1}\cdot A^{\ast}.
(4.24)

4. Если строки матрицы A линейно независимы, то


A^{\sim1}=A^{\ast}\cdot\! \begin{pmatrix}AA^{\ast}\end{pmatrix}^{-1}.
(4.25)

Пример 4.13. Для матриц найти псевдообратные


A=\begin{pmatrix}1&0&0\\0&2&0\\0&0&0\end{pmatrix}\!,\quad B=\begin{pmatrix}1&2&3\end{pmatrix}\!,\quad C=\begin{pmatrix}1&0\\1&1\\1&1 \end{pmatrix}\!.

Решение. Матрица A — диагональная. По формуле (4.23) находим


A^{\sim1}=\begin{pmatrix}1&0&0\\0&1/2&0\\0&0&0\end{pmatrix}\!.

Для матрицы-строки B по формуле (4.25) получаем


B^{\sim1}=\begin{pmatrix}1\\2\\3\end{pmatrix}\!\cdot \begin{bmatrix}\begin{pmatrix} 1&2&3\end{pmatrix}\!\cdot\! \begin{pmatrix}1\\2\\3\end{pmatrix}\end{bmatrix}^{-1}= \begin{pmatrix} 1\\2\\3\end{pmatrix}\!\cdot(14)^{-1}= \begin{pmatrix}1/14\\2/14\\ 3/14 \end{pmatrix}\!.

Столбцы матрицы C линейно независимы. По формуле (4.24) имеем


\begin{aligned}C^{\sim1}&= \begin{bmatrix} \begin{pmatrix}1&1&1\\ 0&1&2 \end{pmatrix}\!\cdot\! \begin{pmatrix} 1&0\\1&1\\1&1\end{pmatrix}\end{bmatrix}^{-1}\cdot\! \begin{pmatrix}1&1&1\\0&1&1\end{pmatrix}= \begin{pmatrix}3&2\\2&2\end{pmatrix}^{-1}\cdot\! \begin{pmatrix}1&1&1\\0&1&1\end{pmatrix}=\\[2pt] &=\frac{1}{2}\! \begin{pmatrix}2&-2\\ -2&3 \end{pmatrix}\!\cdot\! \begin{pmatrix}1&1&1\\0&1&1\end{pmatrix}= \begin{pmatrix}1&0&0\\ -1&1/2&1/2\end{pmatrix}\!.\end{aligned}

Найдем матрицу C^{\sim1} первым способом


1. Составляем блочную матрицу

\begin{pmatrix}C\!\!&\vline\!\!&E_3\\ E_2\!\!&\vline\!\!&{}\end{pmatrix}= \begin{pmatrix}1&0\!\!&\vline\!\!&1&0&0\\ 1&1\!\!&\vline\!\!&0&1&0\\ 1&1\!\!&\vline\!\!& 0&0&1\\\hline 1&0\!\!&\vline\!\!&{}&{}&{}\\ 0&1\!\!&\vline\!\!&{}&{}&{} \end{pmatrix}\!.

2. Приводим блок C к простейшему виду


\begin{pmatrix} C\!\!&\vline\!\!&E_3\\ E_2\!\!&\vline\!\!&{}\end{pmatrix}= \begin{pmatrix}1&0\!\!&\vline\!\!&1&0&0\\ 1&1\!\!&\vline\!\!&0&1&0\\ 1&1\!\!&\vline\!\!& 0&0&1\\\hline 1&0\!\!&\vline\!\!&{}&{}&{}\\ 0&1\!\!&\vline\!\!&{}&{}&{}\end{pmatrix}\sim \begin{pmatrix}1&0\!\!&\vline\!\!&1&0&0\\ 0&1\!\!&\vline\!\!&-1&1&0\\ 0&1\!\!&\vline\!\!& -1&0&1\\\hline 1&0\!\!&\vline\!\!&{}&{}&{}\\ 0&1\!\!&\vline\!\!&{}&{}&{}\end{pmatrix}\sim \begin{pmatrix}1&0\!\!&\vline\!\!&1&0&0\\ 0&1\!\!&\vline\!\!&-1&1&0\\ 0&0\!\!&\vline\!\!& 0&-1&1\\\hline 1&0\!\!&\vline\!\!&{}&{}&{}\\ 0&1\!\!&\vline\!\!&{}&{}&{}\end{pmatrix}

Получили матрицы T=E_2,~ S=\begin{pmatrix}1&0&0\\-1&1&0\\ 0&-1&1\end{pmatrix}. Ранг матрицы C равен количеству ненулевых строк, т.е. r=\operatorname{rg}C=2=n, E_r=E_2.


3. Находим произведения


SS^{\ast}= \begin{pmatrix}1&0&0\\-1&1&0\\ 0&-1&1\end{pmatrix}\!\cdot\! \begin{pmatrix}1&-1&0\\0&1&-1\\ 0&0&1\end{pmatrix}= \begin{pmatrix}1&-1\!\!&\vline\!\!&0\\ -1&2\!\!&\vline\!\!&-1\\\hline 0&-1\!\!&\vline\!\!&2\end{pmatrix}=\begin{pmatrix}S_1\!\!&\vline\!\!&S_2\\\hline S_3\!\!&\vline\!\!&S_4\end{pmatrix}\!;\quad T^{\ast}T=\begin{pmatrix}1&0\\0&1 \end{pmatrix}=(T_1).

В матрице T^{\ast}T блоки T_2,T_3,T_4 отсутствуют (см. замечание 4.8).


4. Из-за отсутствия блоков T_3,T_4 вычислять матрицу U не нужно. Находим только матрицу


V=-S_2S_4^{-1}= -\begin{pmatrix}0\\-1\end{pmatrix}\!\cdot(2)^{-1}= \begin{pmatrix} 0\\1/2\end{pmatrix}\!.

5. Находим псевдообратную матрицу по формуле (4.21), учитывая отсутствие матрицы U (см. п.5 замечаний 4.6):


\begin{aligned}C^{\sim1}&= T\cdot E_2\cdot\! \begin{pmatrix}E_2\mid V\end{pmatrix}\!\cdot S= \begin{pmatrix}1&0\\0&1\end{pmatrix}\!\cdot\! \begin{pmatrix}1&0\\0&1 \end{pmatrix}\!\cdot\! \begin{pmatrix}1&0\!\!&\vline\!\!&0\\ 0&1\!\!&\vline\!\! &1/2 \end{pmatrix}\!\cdot\! \begin{pmatrix}1&0&0\\-1&1&0\\ 0&-1&1\end{pmatrix}=\\[2pt] &= \begin{pmatrix}1&0\!\!&\vline\!\!&0\\ 0&1\!\!&\vline\!\! &1/2\end{pmatrix}\!\cdot\! \begin{pmatrix}1&0&0\\-1&1&0\\ 0&-1&1\end{pmatrix}= \begin{pmatrix}1&0&0\\-1&1/2&1/2 \end{pmatrix}\!.\end{aligned}

Результаты обоих способов нахождения матрицы C^{\sim1} совпадают.




Псевдообращение блочных матриц


Пусть дана блочная матрица M=\begin{pmatrix}A\!\!&\vline\!\!&B\\\hline C\!\!&\vline\!\!&D \end{pmatrix}\!, где A — квадратная невырожденная матрица r-го порядка и r=\operatorname{rg}A= \operatorname{rg}M.


Тогда псевдообратная матрица M^{\sim1} — является блочной матрицей и находится по формуле


M^{\sim1}= \begin{pmatrix}\dfrac{A^{\ast}}{B^{\ast}}\end{pmatrix}\!\cdot\! \begin{pmatrix}AA^{\ast}+BB^{\ast}\end{pmatrix}^{-1}\cdot A\cdot\! \begin{pmatrix}A^{\ast}A+ C^{\ast}C\end{pmatrix}^{-1}\cdot\! \begin{pmatrix}A^{\ast}\mid C^{\ast} \end{pmatrix}\!.
(4.26)

Пример 4.14. Для блочной матрицы M=\begin{pmatrix}1&-1\!\!&\vline\!\!&2&0\\ -1&2\!\!&\vline\!\!&-3&1\\\hline 0&1\!\!&\vline\!\!&-1&1\end{pmatrix}= \begin{pmatrix} A\!\!&\vline\!\!&B\\\hline C\!\!&\vline\!\!&D \end{pmatrix} найти псевдообратную.


Решение. Заметим, что \operatorname{rg}A=2=\operatorname{rg}M. Последовательно вычисляем


\begin{aligned}\begin{pmatrix}AA^{\ast}+BB^{\ast}\end{pmatrix}^{-1}&= \begin{bmatrix}\begin{pmatrix}1&-1\\-1&2\end{pmatrix}\!\cdot\! \begin{pmatrix}1&-1\\-1&2 \end{pmatrix}+ \begin{pmatrix}2&0\\-3&1\end{pmatrix}\!\cdot\! \begin{pmatrix}2&-3\\ 0&1\end{pmatrix}\end{bmatrix}^{-1}=\\[2pt] &=\begin{bmatrix}\begin{pmatrix}2&-3\\-3&5\end{pmatrix}+ \begin{pmatrix}4&-6\\-6&10\end{pmatrix} \end{bmatrix}^{-1}= \begin{pmatrix}6&-9\\ -9&15\end{pmatrix}^{-1}= \frac{1}{3}\! \begin{pmatrix}5&3\\3&2\end{pmatrix}\!;\\[5pt] \begin{pmatrix}A^{\ast}A+C^{\ast}C\end{pmatrix}^{-1}&= \begin{bmatrix}\begin{pmatrix}1&-1\\-1&2\end{pmatrix}\!\cdot\! \begin{pmatrix}1&-1\\-1&2 \end{pmatrix}+ \begin{pmatrix}0\\1 \end{pmatrix}\!\cdot\! \begin{pmatrix}0&1\end{pmatrix}\end{bmatrix}^{-1}=\\[2pt] &=\begin{bmatrix}\begin{pmatrix}2&-3\\-3&5\end{pmatrix}+ \begin{pmatrix}0&0\\0&1\end{pmatrix} \end{bmatrix}^{-1}= \begin{pmatrix}2&-3\\ -3&6\end{pmatrix}^{-1}= \frac{1}{3}\! \begin{pmatrix}6&3\\3&2\end{pmatrix}\!;\end{aligned}

\begin{pmatrix} AA^{\ast}+ BB^{\ast}\end{pmatrix}^{-1}\cdot A\cdot\! \begin{pmatrix} A^{\ast}A+ C^{\ast}C\end{pmatrix}^{-1}= \frac{1}{3}\! \begin{pmatrix}5&3\\ 3&2 \end{pmatrix}\!\cdot\! \begin{pmatrix}1&-1\\-1&2\end{pmatrix}\!\cdot \frac{1}{3}\! \begin{pmatrix}6&3\\3&2\end{pmatrix}= \frac{1}{9}\! \begin{pmatrix}2&1\\1&1\end{pmatrix}\!\cdot\! \begin{pmatrix} 6&3\\3&2\end{pmatrix}= \frac{1}{9}\! \begin{pmatrix} 15&8\\9&5 \end{pmatrix}\!.

По формуле (4.26) имеем


\begin{aligned}M^{\sim1}&= \begin{pmatrix}1&-1\\-1&2\\\hline 2&-3\\0&1\end{pmatrix}\!\cdot\frac{1}{9}\! \begin{pmatrix}15&8\\9&5\end{pmatrix}\!\cdot\! \begin{pmatrix}1&-1\!\!&\vline\!\!&0\\ -1&2\!\!&\vline\!\!&1\end{pmatrix}= \frac{1}{9}\! \begin{pmatrix}1&-1\\-1&2\\\hline 2&-3\\0&1\end{pmatrix}\!\cdot\! \begin{pmatrix} 7&1\!\!&\vline\!\!&8\\ 4&1\!\!&\vline\!\!&5\end{pmatrix}=\\[2pt] &=\frac{1}{9}\! \begin{pmatrix}3&0\!\!&\vline\!\!&3\\ 1&1\!\!&\vline\!\!&2\\\hline 2&-1\!\!&\vline\!\!&1\\ 4&1\!\!&\vline\!\!&5\end{pmatrix}= \begin{pmatrix}1/3&0&1/3\\ 1/9&1/9&2/9\\ 2/9&-1/9&1/9\\ 4/9&1/9&5/9\end{pmatrix}\!,\end{aligned}

что совпадает с результатами примера 4.11.
Математический форум (помощь с решением задач, обсуждение вопросов по математике).
Кнопка "Поделиться"
Если заметили ошибку, опечатку или есть предложения, напишите в комментариях.

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


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

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