Математика / Матрицы, определители

Обратная подстановка в методе Гаусса

Обратная подстановка находит неизвестные после прямого хода метода Гаусса. Она идет снизу вверх по ступенчатой системе: сначала последняя ведущая переменная, затем предыдущие.

Опубликовано: Обновлено:

Формула

$$x_i=\frac{b'_i-\sum_{j=i+1}^{n}u_{ij}x_j}{u_{ii}}$$
Лестница снизу вверх Обратная подстановка идет против направления прямого хода

Прямой ход создает ступени сверху вниз, а обратная подстановка считывает значения снизу вверх, используя уже найденные переменные.

Последние строки после прямого хода обычно содержат меньше неизвестных.

Обозначения

$x_i$
неизвестная, которую находят на текущем шаге, значение переменной
$u_{ij}$
элемент верхнетреугольной или ступенчатой матрицы, коэффициент
$b'_i$
правая часть после прямого хода, число
$u_{ii}$
ведущий коэффициент в строке i, коэффициент

Условия применения

  • Система должна быть приведена к ступенчатому или верхнетреугольному виду.
  • Для ведущей переменной в текущей строке коэффициент u_ii должен быть ненулевым.
  • Все переменные правее текущей ведущей переменной должны быть уже найдены или объявлены свободными.

Ограничения

  • Если в системе есть свободные переменные, обратная подстановка дает общее решение через параметры, а не один числовой ответ.
  • Если после прямого хода появилась противоречивая строка, обратную подстановку выполнять бессмысленно: система не имеет решений.
  • В численных задачах ошибки округления из прямого хода переходят в обратную подстановку, поэтому результат стоит проверять подстановкой в исходную систему.

Подробное объяснение

Обратная подстановка работает потому, что прямой ход превращает систему в лестницу. В нижней ненулевой строке обычно меньше всего неизвестных. Если система имеет единственное решение и матрица приведена к верхнетреугольному виду, последняя строка содержит одну неизвестную. Ее легко найти делением правой части на ведущий коэффициент. После этого найденное значение можно использовать в строке выше, где остается на одну неизвестную меньше.

Формула x_i = (b'_i - sum u_ij x_j)/u_ii отражает этот процесс. В строке i есть ведущая переменная x_i и переменные справа от нее. Значения правых переменных уже известны, потому что мы идем снизу вверх. Сначала их вклад переносится в правую часть, то есть вычитается из b'_i. Затем остается u_ii x_i, и деление на u_ii дает x_i.

Если система имеет бесконечно много решений, обратная подстановка не исчезает, но меняет вид. Неведущие переменные объявляют свободными параметрами, например t или s. Затем ведущие переменные выражают через эти параметры. Такой ответ не хуже числового: он описывает все решения системы. Важно не потерять параметры и не превратить бесконечное множество решений в один случайный набор чисел.

После обратной подстановки полезно проверить ответ в исходной системе. Это особенно важно при длинных вычислениях: ошибка в одном знаке на прямом ходе или при подстановке быстро распространяется вверх по строкам.

Как пользоваться формулой

  1. Убедитесь, что прямой ход завершен и противоречивых строк нет.
  2. Начните с нижней ненулевой строки.
  3. Найдите последнюю ведущую переменную или выразите ее через свободные параметры.
  4. Подставьте найденные значения в строку выше.
  5. Продолжайте подниматься вверх, пока не найдены все ведущие переменные.

Историческая справка

Обратная подстановка является естественным завершением методов исключения. Как только вычислитель получает треугольную систему, нижние уравнения становятся простыми, а найденные значения можно возвращать в предыдущие строки. Такой прием использовался в практических расчетах задолго до современной матричной терминологии. В европейской традиции он стал частью метода Гаусса вместе с прямым исключением. В вычислительной математике XX века обратная подстановка стала отдельной стандартной процедурой: она появляется не только после метода Гаусса, но и в LU-разложении, где сначала решают нижнетреугольную систему прямой подстановкой, а затем верхнетреугольную систему обратной подстановкой.

Историческая линия формулы

У обратной подстановки нет отдельного автора: это базовый вычислительный шаг, возникающий после исключения неизвестных. В учебной традиции она рассматривается как вторая часть метода Гаусса, но исторически является общей техникой работы с треугольными системами.

Пример

После прямого хода получена система x + 2y - z = 4, y + 3z = 5, 2z = 6. Начинаем снизу: из последней строки z = 3. Поднимаемся на строку выше: y + 3z = 5, значит y + 9 = 5 и y = -4. Теперь первая строка: x + 2(-4) - 3 = 4, значит x - 11 = 4 и x = 15. Ответ: x = 15, y = -4, z = 3. Формула обратной подстановки делает то же самое в общем виде: из правой части строки вычитаются уже найденные слагаемые, а затем результат делится на ведущий коэффициент. Проверка в исходной ступенчатой системе подтверждает ответ: каждая строка превращается в верное равенство.

Частая ошибка

Частая ошибка - идти сверху вниз, хотя верхняя строка обычно содержит несколько неизвестных. Вторая ошибка - забыть вычесть уже найденные слагаемые из правой части и просто разделить b'_i на u_ii. Третья ошибка - применять формулу при нулевом ведущем коэффициенте, хотя такая строка либо требует другой ведущей переменной, либо показывает свободную переменную. Еще одна ошибка - при наличии параметра выбрать произвольное значение, хотя задача просит общее решение.

Практика

Задачи с решением

Решить треугольную систему

Условие. Решите систему x + y + z = 6, y - z = 1, 2z = 8.

Решение. Из последней строки z = 4. Из второй строки y - 4 = 1, значит y = 5. Из первой строки x + 5 + 4 = 6, значит x = -3.

Ответ. x = -3, y = 5, z = 4

Общее решение с параметром

Условие. После прямого хода получено x + 2y - z = 3, y + z = 4. Выразите решение через свободную переменную z = t.

Решение. Пусть z = t. Из второй строки y + t = 4, значит y = 4 - t. Из первой строки x + 2(4 - t) - t = 3, поэтому x + 8 - 3t = 3 и x = -5 + 3t.

Ответ. x = -5 + 3t, y = 4 - t, z = t

Дополнительные источники

  • MIT OpenCourseWare 18.06SC Linear Algebra, elimination and back substitution
  • Jim Hefferon, Linear Algebra, Chapter One: Linear Systems
  • OpenStax Precalculus 2e, 9.6 Solving Systems with Gaussian Elimination

Связанные формулы

Математика

Прямой ход метода Гаусса

$R_i\leftarrow R_i-\frac{a_{ik}}{a_{kk}}R_k$

Прямой ход метода Гаусса зануляет коэффициенты под ведущими элементами. В результате система приводится к ступенчатому виду, из которого решение находят обратной подстановкой.

Математика

Ступенчатый вид матрицы

$p_1<p_2<\dots<p_r,\quad a_{ij}=0\ \text{ниже ведущих элементов}$

Ступенчатый вид матрицы - это форма, где ведущие элементы ненулевых строк смещаются вправо при движении вниз, а под каждым ведущим элементом стоят нули.

Математика

Метод Гаусса-Жордана

$\left[A\mid b\right]\sim\left[I\mid x\right]$

Метод Гаусса-Жордана продолжает метод Гаусса до приведенного ступенчатого вида. Если система имеет единственное решение, расширенная матрица превращается в [I|x], и ответ читается сразу.