Все формулы РФ

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

2×2 нормальные уравнения без полной инверсии

Малую систему нормальных уравнений 2×2 можно решить вручную через определитель или исключение, не строя полную обратную матрицу.

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

Линейная алгебра

Формула

$$\begin{aligned} \begin{bmatrix}c_{11} & c_{12}\\ c_{12} & c_{22}\end{bmatrix}\begin{bmatrix}x_1\\x_2\end{bmatrix}=\begin{bmatrix}d_1\\d_2\end{bmatrix},\;\Delta=c_{11}c_{22}-c_{12}^2,\\ x_1=\frac{c_{22}d_1-c_{12}d_2}{\Delta},\quad x_2=\frac{-c_{12}d_1+c_{11}d_2}{\Delta}. \end{aligned}$$
table Явные формулы для 2×2

Срочное вычисление решения без общей линейной алгебры.

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

Обозначения

$c_{11}, c_{12}, c_{22}$
элементы симметричной 2×2 матрицы Gram, скаляры
$d_1, d_2$
правая часть 2×1, скаляры
$\Delta$
детерминант Gram-матрицы, скаляр

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

  • Матрица \begin{bmatrix}c_{11}&c_{12}\\ c_{12}&c_{22}\end{bmatrix} симметрична.
  • Δ\neq 0 (иначе система вырождена).
  • Размерность исходной подзадачи: два неизвестных.

Ограничения

  • При очень малом Δ возможны большие потери точности.
  • Формула подходит только для n=2.
  • В коде для больших n использовать неэффективно.

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

В задачах МНК с двумя неизвестными нормальная система имеет симметричный вид [[s00,s01],[s01,s11]] [a,b]^T = [t0,t1]^T. Если определитель s00s11-s01^2 не равен нулю, решение единственно. Его можно получить правилом Крамера, но на практике проще выполнить исключение: занулить один коэффициент, найти вторую переменную и подставить обратно. Этот прием особенно полезен в школьно-вузовском переходе, когда нужно показать связь регрессии, систем линейных уравнений и геометрии проекций без тяжелой матричной нотации. Важно видеть эту формулу в общей цепочке: исходные данные задают матрицу наблюдений A и правую часть b, затем выбирается способ приблизить b в пространстве столбцов A. 2×2 нормальные уравнения без полной инверсии отвечает за прикладная задача наименьших квадратов, поэтому она не существует отдельно от ранга матрицы, ортогональности остатка и устойчивости вычислений. Если столбцы A хорошо различимы и данные имеют умеренный шум, нормальные уравнения могут дать понятный ручной путь. Если столбцы почти зависимы, лучше пользоваться QR или SVD, потому что они меньше усиливают ошибки округления. После вычисления результата полезно проверить три вещи: размерности всех матриц, величину остатка и связь с соседними формулами раздела. Такой подход превращает формулу из механической записи в рабочий инструмент анализа данных, регрессии, инженерных измерений и численной математики.

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

  1. Найдите Δ.
  2. Подставьте c_{ij} и d_i в формулы.
  3. Вычислите x_1 и x_2.
  4. Проверьте результат подстановкой.

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

Ручное решение малых нормальных систем появилось вместе с ранним применением МНК в астрономии и геодезии, когда вычисления выполнялись таблицами и арифмометрами. До массового использования компьютеров умение свести задачу к небольшой системе и аккуратно решить ее было не технической деталью, а важной частью вычислительной культуры. В XX веке эта тема стала частью стандартной численной линейной алгебры: вычислительные машины сделали возможной массовую обработку переопределенных систем, но одновременно показали, что алгебраически эквивалентные формулы могут вести себя по-разному из-за округления. Поэтому учебники начали разделять теоретический вывод МНК, геометрическое объяснение через проекции и практические алгоритмы QR, Холецкого и SVD. Такой исторический сдвиг важен для пользователя: он объясняет, почему на странице рядом стоят не только “красивая формула”, но и условия применимости, ограничения и типичные ошибки.

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

Метод является сочетанием классического исключения Гаусса, правила Крамера для малых систем и нормальных уравнений МНК; он не имеет отдельного единственного автора. Современная запись является результатом развития метода наименьших квадратов, матричной алгебры и численных методов; поэтому атрибуция здесь распределенная: классические идеи связаны с Гауссом и Лежандром, а устойчивые вычислительные формы — с более поздней численной линейной алгеброй.

Пример

Пусть нормальные уравнения имеют вид 3a+3b=5 и 3a+5b=6. Вычитание первого уравнения из второго дает 2b=1, значит b=0.5. Подстановка в первое уравнение дает 3a+1.5=5, поэтому a=7/6. Если решать через определители, главный определитель равен 3·5-3·3=6, а определители для a и b равны 7 и 3. Тогда a=7/6, b=3/6=1/2. Оба пути дают один и тот же результат, но исключение часто быстрее и меньше нагружает запись. Дополнительная проверка: после получения численного ответа всегда подставь найденный вектор обратно в Ax, вычисли остаток r=b-Ax и сравни его норму с нормой остатка для соседнего пробного решения. Если речь идет о МНК, маленькое изменение параметров не должно уменьшать критерий; если оно уменьшает сумму квадратов, значит нормальные уравнения, QR-шаг или ручное исключение выполнены с ошибкой. Такой контроль особенно полезен в учебных задачах, где итоговое число легко получить, но трудно заметить неверный знак или перепутанный порядок умножения.

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

Не стоит выписывать обратную матрицу там, где достаточно двух шагов исключения. Еще одна частая ошибка — округлять промежуточные коэффициенты слишком рано: для маленьких систем лучше сохранять дроби до последнего шага. Отдельно проверяй размерности: произведения A^T A, A^T b, Q^T b и R x допустимы только при согласованных числах строк и столбцов. Ошибка размерности часто маскируется в ручной записи, но сразу ломает смысл формулы.

Практика

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

Ручное решение первого примера

Условие. c_{11}=5,\ c_{12}=2,\ c_{22}=4,\ d_1=11,\ d_2=10.

Решение. Δ=16,\ x_1=(4·11-2·10)/16=3/2,\ x_2=(-2·11+5·10)/16=7/4.

Ответ. (x_1,x_2)=(3/2,7/4).

Второй пример

Условие. c_{11}=3,\ c_{12}=1,\ c_{22}=2,\ d_1=5,\ d_2=3.

Решение. Δ=5,\ x_1=(2·5-1·3)/5=7/5,\ x_2=(-1·5+3·3)/5=4/5.

Ответ. (x_1,x_2)=(7/5,4/5).

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

  • Anton, Elementary Linear Algebra
  • Strang, Linear Algebra and Its Applications

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

Математика

Обратная матрица 2x2

$A^{-1}=\frac{1}{ad-bc}\begin{pmatrix}d&-b\\-c&a\end{pmatrix}$

Обратная матрица 2x2 существует только при ненулевом определителе. Она обращает действие исходной матрицы: A^{-1}A = I, то есть возвращает исходный вектор.

Математика

Холѐцкое разложение нормальных уравнений

$A^\top A = LL^\top,\; L y=A^\top b,\; L^\top \hat x = y.$

Разложение Холецкого применяет положительную определенность A^T A и заменяет решение нормальных уравнений двумя треугольными системами.