Системы уравнений

Система уравнений — набор уравнений с несколькими неизвестными, которые должны выполняться одновременно. Решить систему = найти значения переменных, удовлетворяющие всем уравнениям. Основа линейной алгебры и вычислительной математики.

Две дороги, два путника, разные скорости — когда они встретятся? Задача сводится к системе уравнений. Казалось бы, школьная арифметика. Но те же принципы работают в квантовой механике (решение уравнения Шрёдингера), машинном обучении (градиентный спуск), экономике (модели «затраты-выпуск» Леонтьева).

Вавилоняне и китайцы: 4000 лет истории

Глиняные таблички из Месопотамии (2000-1600 до н.э.) содержат задачи вида: «Длина и ширина дают площадь 60, разность 7. Найти длину и ширину». Современная запись: xy = 60, x - y = 7. Вавилоняне решали подстановкой и геометрическими методами.

Китайская книга «Математика в девяти книгах» (II век до н.э., окончательная редакция III век н.э.) описывает метод «фан-чэн» — прообраз исключения Гаусса. Коэффициенты раскладывали на счётной доске бамбуковыми палочками, вычитали строки, получая треугольную форму.

Европа переоткрыла эти методы через 1500 лет. Рене Декарт (1637) ввёл координаты и алгебраические уравнения для геометрии. Готфрид Лейбниц (1693) использовал определители для систем 3×3. Но систематическая теория — только в XIX веке.

Метод Гаусса: треугольник исключения

Карл Фридрих Гаусс в 1809 году вычислял орбиту астероида Паллада. Метод наименьших квадратов требовал решения систем линейных уравнений — десятки переменных вручную. Гаусс формализовал процедуру последовательного исключения.

Идея: преобразовать систему к треугольному виду. Из первого уравнения исключаем x из всех остальных. Из второго — y из оставшихся. Получаем лестницу.

Сложность O(n³) операций для системы n×n. Для 1000 уравнений — миллиард операций, на современном компьютере за секунды. Вручную Гаусс решал системы 10×10 за часы — гениальная техника счёта.

Контринтуитивный факт: метод Гаусса неустойчив к ошибкам округления. Для плохо обусловленных систем (почти параллельные гиперплоскости) погрешности накапливаются. Частичный выбор главного элемента (pivoting) стабилизирует алгоритм — выбираем максимальный коэффициент, меняем строки.