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

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

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

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


Метод Гаусса решения системы линейных уравнений

Метод Гаусса решения системы линейных уравнений


Пусть дана система (5.3) m линейных уравнений с n неизвестными. Для ее решения нужно выполнить следующие действия.


1. Составить расширенную матрицу (5.7) системы:


\begin{pmatrix}A\mid b\end{pmatrix}= \begin{pmatrix} a_{11}&\cdots&a_{1n}\!\!&\vline\!\!&b_1\\ \vdots&\ddots&\vdots\!\!&\vline\!\!& \vdots\\ a_{m1}&\cdots&a_{mn}\!\!&\vline\!\!&b_m\end{pmatrix}\!.

2. Используя элементарные преобразования над строками матрицы \begin{pmatrix}A\mid b\end{pmatrix}, привести ее к ступенчатому виду. Если базисный минор матрицы A расположен в первых r строках и r столбцах, получится следующий вид:


\begin{pmatrix} \widetilde{A}\mid \widetilde{b}\end{pmatrix}= \begin{pmatrix} 1&\widetilde{a}_{12}&\cdots&\widetilde{a}_{1r}&\cdots&\widetilde{a}_{1n}\!\!&\vline\!\!&\widetilde{b}_1\\ 0&1&\cdots&\widetilde{a}_{2r}&\cdots&\widetilde{a}_{2n}\!\!&\vline\!\!&\widetilde{b}_2\\ \vdots&\vdots&\ddots&\vdots&\ddots&\vdots\!\!&\vline\!\!&\vdots\\ 0&0&\cdots&1&\cdots& \widetilde{a}_{m}\!\!&\vline\!\!&\widetilde{b}_r\\ 0&0&\cdots&0&\cdots&0\!\!&\vline\!\!& \widetilde{b}_{r+1}\\ \vdots&\vdots&\ddots&\vdots&\ddots&\vdots\!\!&\vline\!\!&\vdots\\ 0&0&\cdots&0&\cdots&0\!\!&\vline\!\!&0 \end{pmatrix}\!.
(5.8)

3. Выяснить, совместна система или нет. Для этого определить ранги матриц A и (A\mid b):


\operatorname{rg}A=\operatorname{rg}\widetilde{A}=r — число ненулевых строк в матрице \widetilde{A};


\operatorname{rg}(A\mid b)= \operatorname{rg}(\widetilde{A}\mid \widetilde{b})= \begin{cases}r+1,&\text{if}\quad \widetilde{b}_{r+1}\ne0,\\ \phantom{r}r,&\text{if}\quad \widetilde{b}_{r+1}=0.\end{cases}


Если \operatorname{rg}A\ne \operatorname{rg}(A\mid b) (при \widetilde{b}_{r+1}\ne0), то система не имеет решений. Процесс решения завершен. Если \operatorname{rg}A=\operatorname{rg}(A\mid b) (при \widetilde{b}_{r+1}=0), то система совместна. Процесс решения продолжается.


4. Для совместной системы \operatorname{rg}A=\operatorname{rg}(A\mid b)=r привести матрицу (\widetilde{A}\mid \widetilde{b}) к упрощенному виду. Для этого при помощи элементарных преобразований над строками добиваемся того, чтобы в каждом столбце, входящем в базисный минор, все элементы были равны нулю, за исключением одного, равного единице. Если базисный минор матрицы A расположен в первых r строках и первых r столбцах, то матрица приводится к упрощенному виду:


\begin{pmatrix}A'\mid b'\end{pmatrix}= \begin{pmatrix} 1&0&\cdots&0&a'_{1\,r+1}&\cdots&a'_{1n}\!\!&\vline\!\!&b'_1\\ 0&1&\cdots&0&a'_{2\,r+1}&\cdots&a'_{2n}\!\!&\vline\!\!&b'_2\\ \vdots&\vdots&\ddots&\vdots&\vdots&\ddots&\vdots\!\!&\vline\!\!&\vdots\\ 0&0&\cdots&1&a'_{r\,r+1}&\cdots&a'_{m}\!\!&\vline\!\!&b'_r\\ 0&0&\cdots&0&0&\cdots&0\!\!&\vline\!\!&0\\ \vdots&\vdots&\ddots&\vdots&\vdots&\ddots&\vdots\!\!&\vline\!\!&\vdots\\ 0&0&\cdots&0&0&\cdots&0\!\!&\vline\!\!&0 \end{pmatrix}\!.
(5.9)

Первые четыре пункта составляют прямой ход метода Гаусса. В результате прямого хода исходная система существенно упрощается (имеет вид A'x=b'):


\begin{cases}x_1+a'_{1\,r+1}x_{r+1}+\ldots+a'_{1n}x_n=b'_1,\\ \cdots\cdots\cdots\cdots\cdots\\ x_r+a'_{r\,r+1}x_{r+1}+\ldots+a'_{rn}x_n=b'_r. \end{cases}
(5.10)

5. По упрощенному виду (5.9) разделяем все неизвестные x_1,x_2,\ldots,x_n на базисные и свободные. Неизвестные, которым соответствуют столбцы, входящие в базисный минор, называются базисными переменными, остальные неизвестные — свободными переменными. Для системы (5.10) базисными переменными являются x_1,x_2,\ldots,x_r, свободными переменнымиx_{r+1},x_{r+2},\ldots,x_n . Выражаем в (5.10) базисные переменные через свободные:


\begin{cases}x_1=b'_1-a'_{1\,r+1}x_{r+1}-\ldots-a'_{1n}x_n,\\ \cdots\cdots\cdots\cdots\cdots\\ x_r=b'_r-a'_{r\,r+1}x_{r+1}-\ldots-a'_{rn}x_n.\end{cases}
(5.11)

Если ранг r матрицы системы равен числу n неизвестных (r=\operatorname{rg}A=n), то левый блок матрицы (5.9) будет представлен единичной матрицей E_n:


\begin{pmatrix}A'\mid b'\end{pmatrix}= \begin{pmatrix}1&0&\cdots&0\!\!&\vline\!\!&b'_1\\ 0&1&\cdots&0\!\!&\vline\!\!&b'_2\\ \vdots&\vdots&\ddots&\vdots\!\!&\vline\!\!&\vline\\ 0&0&\cdots&1\!\!&\vline\!\!&b'_n\end{pmatrix}\!.

Все неизвестные x_1,x_2,\ldots,x_n будут базисными и формула (5.11) будет определять единственное решение системы


\begin{cases}x_1=b'_1,\\x_2=b'_2,\\\vdots\\x_n=b'_n.\end{cases}
(5.12)

Если ранг матрицы системы меньше числа неизвестных (\operatorname{rg}A<n), то система имеет бесконечно много решений, задаваемых формулой (5.11), которая обладает следующими свойствами:


– при любых значениях свободных переменных x_{r+1},x_{r+2},\ldots,x_n по формуле (5.11) получаются такие значения базисных переменных, что упорядоченный набор чисел x_1,x_2,\ldots,x_n является решением системы (5.3);


– любое решение x_1,x_2,\ldots,x_n системы (5.3) удовлетворяет равенствам (5.11).


Равенства (5.11), выражающие базисные переменные через свободные, называются общим решением системы (5.3). Решение системы, получающееся по формуле (5.11) при задании конкретных значений свободных переменных, называется частным решением системы (5.3).


Процесс решения совместной системы (5-3) заканчивается получением формулы (5.11) общего решения (в частности, определением единственного решения (5.12)).


Содержание п.5 алгоритма составляет обратный ход метода Гаусса.




Пример 5.3. Решить систему линейных уравнений


\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}

Решение. 1. Составим расширенную матрицу системы


\begin{pmatrix}A\mid B\end{pmatrix}= \begin{pmatrix}1&1&2&1\!\!&\vline&\!\!1\\ 2&3&0&1\!\!&\vline&\!\!0\\ 3&4&2&2\!\!&\vline&\!\!1\end{pmatrix}\!.

2. Используя элементарные преобразования над строками матрицы \begin{pmatrix}A\mid b\end{pmatrix}, приводим ее к ступенчатому виду. Берем в качестве ведущего элемента a_{11}=1\ne0. Ко второй строке прибавим первую, умноженную на (-2), а к третьей — первую, умноженную на (-3):


\begin{pmatrix}A\mid B\end{pmatrix}= \begin{pmatrix}1&1&2&1\!\!&\vline&\!\!1\\ 2&3&0&1\!\!&\vline&\!\!0\\ 3&4&2&2\!\!&\vline&\!\!1\end{pmatrix}\sim \begin{pmatrix}1&1&2&1\!\!&\vline&\!\!1\\ 0&1&-4&-1\!\!&\vline&\!\!-2\\ 0&1&-4&-1\!\!&\vline&\!\!-2\end{pmatrix}\!.

Ведущий элемент a_{22}=1\ne0. К третьей строке прибавим вторую, умноженную на (-1):


\begin{pmatrix}A\mid B\end{pmatrix}\sim \begin{pmatrix}1&1&2&1\!\!&\vline&\!\!1\\ 0&1&-4&-1\!\!&\vline&\!\!-2\\ 0&1&-4&-1\!\!&\vline&\!\!-2\end{pmatrix}\sim \begin{pmatrix}1&1&2&1\!\!&\vline&\!\!1\\ 0&1&-4&-1\!\!&\vline&\!\!-2\\ 0&0&0&0\!\!&\vline&\!\!0\end{pmatrix}= \begin{pmatrix}\widetilde{A}\mid \widetilde{B}\end{pmatrix}\!.

Расширенная матрица системы приведена к ступенчатому виду.


3. Определяем ранги матриц: \operatorname{rg}A= \operatorname{rg}(A\mid b)=2. Согласно теореме 5.2, система совместна.


4. Приводим матрицу к упрощенному виду. Выбираем в качестве базисного минора M_{1\,2}^{1\,2}. К первой строке прибавляем вторую, умноженную на (-1):


\begin{pmatrix}\widetilde{A}\mid \widetilde{B}\end{pmatrix}= \begin{pmatrix} 1&1&2&1\!\!&\vline&\!\!1\\ 0&1&-4&-1\!\!&\vline&\!\!-2\\ 0&0&0&0\!\!&\vline&\!\!0\end{pmatrix}\sim \begin{pmatrix}1&0&6&2\!\!&\vline&\!\!3\\ 0&1&-4&-1\!\!&\vline&\!\!-2\\ 0&0&0&0\!\!&\vline&\!\!0\end{pmatrix}= \begin{pmatrix}A'\mid b'\end{pmatrix}\!.

5. Переменные x_1,\,x_2 — базисные, а x_2,\,x_4 — свободные. Записываем общее решение системы, используя (5.11):


\begin{cases}x_1=3-6x_3-2x_4,\\ x_2=-2+4x_3+x_4.\end{cases}

Система имеет бесконечно много решений. Найдем частное решение. Например, для x_3=x_4=0 получаем x_1=3,~x_2=-2. Следовательно, столбец x=\begin{pmatrix}3&-2&0&0\end{pmatrix}^T — частное решение системы.

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

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


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

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