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

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

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

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


Решение систем уравнений с помощью полуобратных матриц

Решение систем уравнений с помощью полуобратных матриц


Получим решение системы линейных уравнений при помощи полуобратной матрицы.


Требуется решить неоднородную систему линейных уравнений

Ax=b,
(5.17)

где A — произвольная матрица размеров m\times n. Если матрица системы нулевая (A=O), то система либо несовместна (при b\ne o), либо имеет бесконечное множество решений (при b=o любой подходящий по размерам столбец x является решением). Далее рассматриваем случай ненулевой матрицы A.


Пусть A^{\neg1} — матрица, полуобратная к матрице системы. Используя определение (4.12) полуобратной матрицы, запишем (5.17) в виде


A\cdot A^{\neg1}\cdot A\cdot x=b
(5.18)

Если x — решение системы (5.17), то подставляя Ax=b в левую часть (5.18), приходим к равенству


A\cdot A^{\neg1}\cdot b=b.
(5.19)

Преобразуем (5.19) следующим образом:


\begin{pmatrix}E_m-A\cdot A^{\neg1}\end{pmatrix}\!\cdot b=o.
(5.20)

Следовательно, условие (5.20) — необходимое для совместности системы уравнений (5.17). Из равенств (5.17), (5.19) также следует, что столбец


x=A^{\neg1}\cdot b
(5.21)

является решением системы (5.17). Значит, условие (5.20) — это необходимое и достаточное условие совместности системы (5.17).


Поскольку полуобратная матрица определена неоднозначно (см. п.3 замечаний 4.5), то формула (5.21) фактически задает множество решений системы (5.17). Преобразуем (5.21) так, чтобы была видна структура этого множества, в частности, выясним количество независимых параметров, входящих в (5.15).


Возьмем фиксированную полуобратную матрицу (4.16):


A_0^{\neg1}= T\cdot \Lambda^T\cdot S=T\cdot\! \begin{pmatrix}E_r\!\!&\vline\!\!&O\\\hline O\!\!&\vline\!\!&O\end{pmatrix}\!\cdot S,
(5.22)

где T,\,S — элементарные матрицы порядков n и m соответственно, \Lambda -матрица простейшего вида, эквивалентная матрице A (\Lambda=SAT и A=S^{-1}\Lambda T^{-1}), r=\operatorname{rg}A.


Преобразуем матрицу E_m-AA^{\neg1} в левой части (5.20), подставив выражение (5.22) и A=S^{-1}\Lambda T^{-1}\colon


E_m-AA^{\neg1}= E_m-S^{-1}\Lambda\underbrace{T^{-1}\cdot T}_{E_m} \Lambda^TS= S^{-1}\cdot\! \begin{pmatrix}E_m-\begin{pmatrix}E_r\!\!&\vline\!\!&O\\\hline O\!\!&\vline\!\!&O \end{pmatrix}\end{pmatrix}\!\cdot S=S^{-1}\cdot\! \begin{pmatrix} O\!\!&\vline\!\!&O\\\hline O\!\!&\vline\!\!&E_{m-r}\end{pmatrix}\!\cdot S.

Образуем матрицу \Psi размеров (m-r)\times m из последних (m-r) строк матрицы S\colon


\Psi= \begin{pmatrix}O\mid E_{m-r}\end{pmatrix}\!\cdot S.
(5.23)

При помощи этой матрицы равенство (5.20) принимает вид S^{-1}\Psi b=o. Умножая обе его части слева на невырожденную матрицу S, получаем


\Psi b=o,
(5.24)

Условие (5.24) необходимо и достаточно для совместности системы (5.17). Система (5.24) содержит (m-r) числовых равенств. При r=m матрица \Psi нулевая и система (5.17) совместна.


Используя полуобратную матрицу (5.22), по формуле (5.21) получаем частное решение неоднородной системы


x^H=A_{0}^{\neg1}\cdot b=T\cdot\! \begin{pmatrix}E_r\!\!&\vline\!\!&O\\\hline O\!\!&\vline\!\!&O \end{pmatrix}\!\cdot S\cdot b.

Найдем теперь общее решение однородной системы. Составим матрицу \Phi размера n\times(n-r) из последних (n-r) столбцов матрицы T\colon


\Phi=T\cdot\! \begin{pmatrix}\dfrac{O}{E_{n-r}}\end{pmatrix}\!.
(5.25)

Заметим, что эти столбцы линейно независимы, так как матрица T невырожденная (\operatorname{rg}T=n). Кроме того, они являются решениями однородной системы, так как справедливо:


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

Следовательно, матрица \Phi=\begin{pmatrix} \varphi_1& \cdots&\varphi_{n-r}\end{pmatrix} — фундаментальная, ее столбцы образуют фундаментальную систему решений однородной системы. По теореме 5.4 общее решение системы (5.17) можно представить в виде


x=x^H+x^O= A_{0}^{\neg1}b+\Phi\cdot c= T\cdot\! \begin{pmatrix} E_r\!\!&\vline\!\!&O\\\hline O\!\!&\vline\!\!&O \end{pmatrix}\!\cdot S\cdot b+T\cdot\! \begin{pmatrix}\dfrac{O}{E_{n-r}}\end{pmatrix}\!\cdot c,
(5.26)

где c=\begin{pmatrix}C_1&\cdots&C_{n-r}\end{pmatrix}^T — столбец произвольных постоянных.


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


Теорема 5.5 о совместности неоднородной системы и о структуре ее общего решения. Неоднородная система (5.17) совместна тогда и только тогда, когда столбец свободных членов является решением однородной системы (5.24). Если система (5.17) совместна, то ее общее решение имеет вид (5.26).




Алгоритм применения полуобратной матрицы


Нахождение решения системы (5.17) в виде (5.26) сводится к следующим действиям (рассматривается случай, когда матрица системы ненулевая).


1. Привести матрицу A системы (5.17) к простейшему виду: \Lambda=SAT. При этом находятся элементарные преобразующие матрицы S и T, а также ранг r=\operatorname{rg}A \geqslant1.


2. Проверить условие (5.24) совместности системы. При r=m система совместна. Если r<m, то составить матрицу \Psi=\begin{pmatrix}O\mid E_{m-r}\end{pmatrix}\!\cdot S и проверить условие \Psi b=o. Если условие выполняется, то система совместна. В противном случае система несовместна и процесс решения заканчивается.


3. Найти частное решение неоднородной системы по формуле


x^H=A_0^{\neg1}\cdot b=T\cdot\! \begin{pmatrix} E_r\!\!&\vline\!\!&O\\\hline O\!\!&\vline\!\!&O \end{pmatrix}\!\cdot S\cdot b.

4. Составить фундаментальную матрицу \Phi=T\cdot\! \begin{pmatrix}\dfrac{O}{E_{n-r}}\end{pmatrix}.


5. Записать общее решение системы (5.17) в виде (5.26)


x=x^H+\Phi\cdot c, где c=\begin{pmatrix}C_1&\cdots&C_{n-r}\end{pmatrix}^{T} — столбец произвольных постоянных.

Замечания 5.5


1. Формула (5.25) определяет третий способ нахождения фундаментальной матрицы.


2. Положив x^H=o в формуле (5.26), получаем структуру общего решения однородной системы


x=\Phi\cdot c=T\cdot\! \begin{pmatrix}\dfrac{O}{E_{n-r}}\end{pmatrix}\!\cdot c, где c=\begin{pmatrix}C_1&\cdots&C_{n-r}\end{pmatrix}^{T} — столбец произвольных постоянных.

3. Теорема 5.5 равносильна альтернативе Фредгольма.




Пример 5.6. Решить систему уравнений с помощью полуобратной матрицы


\begin{cases}x_1+ x_2+2x_3+x_4=1,\\2x_1+3x_2+x_4=0,\\ 3x_1+4x_2+2x_3+2x_4=1,\end{cases} представив ее решения в виде (5.26).

Решение. 1. К матрице A системы приписываем справа и снизу единичные матрицы порядков m=3 и n=4 соответственно:


\begin{pmatrix}A\!\!&\vline\!\!&E_3\\\hline E_4\!\!&\vline\!\!&{}\end{pmatrix}= \begin{pmatrix}1&1&2&1\!\!&\vline\!\!&1&0&0\\ 2&3&0&1\!\!&\vline\!\!&0&1&0\\ 3&4&2&2\!\!&\vline\!\!&0&0&1\\\hline 1&0&0&0\!\!&\vline\!\!&{}&{}&{}\\ 0&1&0&0\!\!&\vline\!\!&{}&{}&{}\\ 0&0&1&0\!\!&\vline\!\!&{}&{}&{}\\ 0&0&0&1\!\!&\vline\!\!&{}&{}&{} \end{pmatrix}\!.

Элементарными преобразованиями приводим матрицу A к простейшему виду. Взяв ведущий элемент a_{11}=1, делаем равными нулю остальные элементы в первой строке и в первом столбце блока A:


\begin{pmatrix}A\!\!&\vline\!\!&E_3\\\hline E_4\!\!&\vline\!\!&{}\end{pmatrix}\sim \begin{pmatrix}1&1&2&1\!\!&\vline\!\!&1&0&0\\ 0&1&-4&-1\!\!&\vline\!\!&-2&1&0\\ 0&1&-4&-1\!\!&\vline\!\!&-3&0&1\\\hline 1&0&0&0\!\!&\vline\!\!&{}&{}&{}\\ 0&1&0&0\!\!& \vline\!\!&{}& {}&{}\\ 0&0&1&0\!\!&\vline\!\!&{}&{}&{}\\ 0&0&0&1\!\!&\vline\!\!&{}&{}&{} \end{pmatrix}\sim \begin{pmatrix}1&0&0&0\!\!& \vline\!\!&1&0&0\\ 0&1&-4&-1\!\!&\vline\!\!&-2&1&0\\ 0&1&-4&-1\!\!&\vline\!\!&-3&0&1\\\hline 1&-1&-2&-1\!\!&\vline\!\!&{}&{}&{}\\ 0&1&0&0\!\!&\vline\!\!& {}&{}&{}\\ 0&0&1&0\!\!&\vline\!\!&{}&{}&{}\\ 0&0&0&1\!\!&\vline\!\!&{}&{}&{} \end{pmatrix}\!.

Взяв в полученной матрице ведущий элемент a_{22}=1, делаем равными нулю остальные элементы во второй строке и во втором столбце блока A\colon


\begin{pmatrix}A\!\!&\vline\!\!&E_3\\\hline E_4\!\!&\vline\!\!&{}\end{pmatrix}\sim \begin{pmatrix}1&0&0&0\!\!& \vline\!\!&1&0&0\\ 0&1&-4&-1\!\!&\vline\!\!&-2&1&0\\ 0&1&-4&-1\!\!&\vline\!\!&-3&0&1\\\hline 1&-1&-2&-1\!\!&\vline\!\!&{}&{}&{}\\ 0&1&0&0\!\!&\vline\!\!& {}&{}&{}\\ 0&0&1&0\!\!&\vline\!\!&{}&{}&{}\\ 0&0&0&1\!\!&\vline\!\!&{}&{}&{} \end{pmatrix}\sim \begin{pmatrix} 1&0&0&0\!\!&\vline\!\!&1&0&0\\ 0&1&0&0\!\!&\vline\!\!&-2&1&0\\ 0&0&0&0\!\!&\vline\!\!&-1&-1&1\\\hline 1&-1&-6&-2\!\!&\vline\!\!&{}&{}&{}\\ 0&1&4&1\!\!& \vline\!\!&{}& {}&{}\\ 0&0&1&0\!\!&\vline\!\!&{}&{}&{}\\ 0&0&0&1\!\!&\vline\!\!&{}&{}&{} \end{pmatrix}\!.

Простейший вид матрицы A получен, следовательно,


\Lambda=SAT=\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\\-2&1&0\\-1&-1&1\end{pmatrix}\!,\quad T=\begin{pmatrix}1&-1&-6&-2\\ 0&1&4&1\\ 0&0&1&0\\ 0&0&0&1 \end{pmatrix}\!,\quad r=2.

2. Проверяем условие совместности (5.24). Так как r=2<3=m , то составляем матрицу \Psi=\begin{pmatrix}O\mid E_{m-r}\end{pmatrix}\!\cdot S, выделяя последнюю строку матрицы S\colon


\Psi=\begin{pmatrix}0&0\mid 1\end{pmatrix}\!\cdot\! \begin{pmatrix}1&0&0\\-2&1&0\\ -1&-1&1\end{pmatrix}= \begin{pmatrix}-1&-1&1\end{pmatrix}\!.

Записываем условие \Psi b=o\colon\, \begin{pmatrix}-1&-1&1\end{pmatrix}\!\cdot\! \begin{pmatrix}1\\0\\1\end{pmatrix}=0. Условие выполняется, значит, система совместна.


3. Находим частное решение неоднородной системы


\begin{aligned}x^H&=T\cdot\! \begin{pmatrix}E_r\!\!&\vline\!\!&O\\\hline O\!\!&\vline\!\!& O\end{pmatrix}\!\cdot S\cdot b= \begin{pmatrix}1&-1&-6&-2\\0&1&4&1\\ 0&0&1&0\\ 0&0&0&1 \end{pmatrix}\!\cdot\! \begin{pmatrix}1&0\!\!&\vline\!\!&0\\ 0&1\!\!&\vline\!\!&0\\\hline 0&0\!\!& \vline\!\!&0\\ 0&0\!\!&\vline\!\!&0\end{pmatrix}\!\cdot\! \begin{pmatrix}1&0&0\\-2&1&0\\-1&-1&1\end{pmatrix}\!\cdot\! \begin{pmatrix}1\\0\\1\end{pmatrix}=\\[2pt] &=\begin{pmatrix}1&-1&0\\0&1&0\\0&0&0\\0&0&0\end{pmatrix}\!\cdot\! \begin{pmatrix}1&0&0\\-2&1&0\\-1&-1&1\end{pmatrix}\!\cdot\! \begin{pmatrix}1\\0\\1\end{pmatrix}= \begin{pmatrix}1&-1&0\\0&1&0\\0&0&0\\0&0&0\end{pmatrix}\!\cdot\! \begin{pmatrix}1\\-2\\0\end{pmatrix}= \begin{pmatrix}3\\-2\\0\\0\end{pmatrix}\!.\end{aligned}

4. Записываем фундаментальную матрицу \Phi=T\cdot\! \begin{pmatrix}\dfrac{O}{E_{n-r}}\end{pmatrix}, составляя ее из последних n-r=4-2=2 столбцов матрицы T


\Phi=T\cdot\! \begin{pmatrix}\dfrac{O}{E_{n-r}}\end{pmatrix}= \begin{pmatrix}1&-1&-6&-2\\0&1&4&1\\ 0&0&1&0\\ 0&0&0&1\end{pmatrix}\!\cdot\! \begin{pmatrix}0&0\\0&0\\\hline 1&0\\0&1\end{pmatrix}= \begin{pmatrix} -6&-2\\4&1\\1&0\\0&1\end{pmatrix}\!.

Столбцы этой матрицы образуют фундаментальную систему решений однородной системы.


5. Записываем общее решение системы в виде (5.26)


x=x^H+\Phi\cdot c= \begin{pmatrix}3\\-2\\0\\0\end{pmatrix}+ \begin{pmatrix}-6&-2\\4&1\\1&0\\0&1\end{pmatrix}\!\cdot\! \begin{pmatrix}C_1\\C_2\end{pmatrix}= \begin{pmatrix}3\\-2\\0\\0\end{pmatrix}+ C_1\cdot\! \begin{pmatrix}-6\\4\\1\\0\end{pmatrix}+ C_2\cdot\! \begin{pmatrix}-2\\1\\0\\1\end{pmatrix}\!.

Результат совпадает с решением примера 5.5.




Пример 5.7. Проверить совместность системы уравнений с помощью полуобратной матрицы


\begin{cases}x_1+x_2+2x_3+x_4=1,\\2x_1+3x_2+x_4=0,\\3x_1+4x_2+2x_3+2x_4=2,\end{cases} используя условие (5.24).

Решение. Данная система отличается от решаемой в примере 5.6 свободным членом третьего уравнения. При решении примера 5.6 были получены матрица


S=\begin{pmatrix}1&0&0\\-2&1&0\\-1&-1&1\end{pmatrix}

и ранг r=2 матрицы системы. Так как m-r=3-2=1, то в матрицу \Psi следует включить только одну последнюю строку матрицы S, т.е. \Psi=\begin{pmatrix}-1&-1&1\end{pmatrix}. Вычислим левую часть (5.24):


\Psi\cdot b=\begin{pmatrix}-1&-1&1\end{pmatrix}\!\cdot\! \begin{pmatrix}1\\0\\1 \end{pmatrix}=1.

Следовательно, условие \Psi b=o не выполняется, т.е. система несовместна. Тот же вывод получен в примере 5.2 по теореме 5.2 Кронекера-Капелли.

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

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


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

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