Метод Гаусса: основы, применение, примеры решения

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

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

Уравнение в матричной форме

Для начала, давайте рассмотрим систему линейных уравнений с n неизвестными и выразим ее в матричной форме. Пусть у нас есть матрица A размером n×n, вектор X, состоящий из n неизвестных, и вектор B, состоящий из n элементов правой части системы. Наша задача состоит в том, чтобы найти вектор X, удовлетворяющий уравнению AX = B.

Процесс решения методом Гаусса

Процесс решения методом Гаусса состоит из трех основных шагов:

1. Прямой ход: Матрица A и вектор B преобразуются с помощью элементарных преобразований таким образом, чтобы система привелась к треугольному виду. В этом шаге мы выполняем следующие действия:
   • Выбираем первое уравнение и делаем его первый элемент равным единице, нормализуя его.
   • Вычитаем из всех последующих уравнений первое уравнение, умноженное на такое число, чтобы обнулить первый элемент каждого уравнения, кроме первого.
   • Повторяем шаги a) и b) для каждого последующего уравнения.
2. Обратный ход: Начиная с последнего уравнения и идя в обратном порядке, находим значения неизвестных методом обратной подстановки. В этом шаге мы выполняем следующие действия:
   • Находим значение последней неизвестной Xn, равное Bn/A[n,n], где n - размерность матрицы.
   • Подставляем значение Xn в предыдущее уравнение и находим значение предыдущей неизвестной Xn-1.
   • Повторяем шаг b) для каждой предыдущей неизвестной, пока все неизвестные не будут найдены.
3. Проверка и вывод результата: Проверяем полученное решение, подставляя его в исходную систему. Если оно удовлетворяет каждому из уравнений, то мы получили корректный ответ и выводим результат. В противном случае, система несовместна или имеет бесконечное количество решений.

Пример кода на языке Python

Приведенный ниже код на языке Python реализует метод Гаусса для решения системы линейных уравнений:

import numpy as np

# Функция для решения системы линейных уравнений методом Гаусса
def solve_gauss(A, B):
    n = len(A)
  
    # Прямой ход
    for i in range(n):
        # Нормализация i-го уравнения
        div = A[i][i]
        A[i] /= div
        B[i] /= div
        
        # Вычитание i-го уравнения из последующих
        for j in range(i+1, n):
            multi = A[j][i]
            A[j] -= multi * A[i]
            B[j] -= multi * B[i]
  
    # Обратный ход
    X = np.zeros(n)
    for i in range(n-1, -1, -1):
        X[i] = B[i]
        for j in range(i+1, n):
            X[i] -= A[i][j] * X[j]
  
    return X

# Пример использования
A = np.array([[2, 1, -1],
              [1, -1, 1],
              [-3, 2, -5]])

B = np.array([3, 2, -1])

X = solve_gauss(A, B)
print("Решение:", X)

В данном примере мы решаем систему линейных уравнений с матрицей A и вектором B, используя функцию solve_gauss. Результат выводится на экран.

Метод Гаусса является эффективным и широко применяемым методом решения систем линейных уравнений. Он позволяет найти единственное решение или определить несовместность или бесконечное количество решений.