Определитель квадратной матрицы размеров , заданной над коммутативным кольцом, является элементом кольца . Эта величина определяет многие свойства матрицы , в частности, матрица обратима тогда и только тогда, когда её определитель является обратимым элементом кольца . В случае, когда — поле, определитель матрицы равен нулю тогда и только тогда, когда ранг матрицы меньше , то есть когда системы строк и столбцов матрицы являются линейно зависимыми.
В Европе определители матриц 2 × 2 встречаются у Кардано в XVI веке.
Для старших размерностей определение детерминанта дано Лейбницем в 1693 году.
Первая публикация принадлежит Крамеру. Теория определителей создана Вандермондом, Лапласом, Коши и Якоби. Термин «определитель» в современном его значении ввёл О. Коши (1815), хотя ранее (1801) «детерминантом» К. Гаусс назвал дискриминант квадратичной формы.
Японский математик Сэки Такакадзу ввёл определители независимо в 1683 году[3].
−1, если все индексы различны и подстановка нечётна.
Аксиоматическое построение (определение на основе свойств)[править | править код]
Понятие определителя может быть введено на основе его свойств. А именно, определителем вещественной матрицы называется функция , обладающая следующими тремя свойствами[4]:
— кососимметрическая функция строк (столбцов) матрицы .
Площадь параллелограмма равна модулю определителя матрицы, образованной векторами — сторонами параллелограмма.
Для матрицы определитель вычисляется как:
Эта матрица A может быть рассмотрена как матрица линейного отображения, преобразующего единичный квадрат в параллелограмм с вершинами (0, 0), (a, b), (a + c, b + d), и (c, d).
Абсолютное значение определителя равно площади этого параллелограмма, и, таким образом, отражает коэффициент, на который масштабируются площади при преобразовании A.
Значение определителя со знаком (ориентированная площадь параллелограмма) помимо коэффициента масштабирования также показывает, выполняет ли преобразование A отражение.
Для более удобного вычисления определителя третьего порядка можно воспользоваться правилом Саррюса или правилом треугольника.
Определитель матрицы, составленной из векторов равен их смешанному произведению в правой декартовой системе координат и, аналогично двумерному случаю, представляет собой ориентированный объём параллелепипеда, натянутого на . Геометрически определитель матрицы 3х3 означает объем параллелепипеда, построенного на векторах матрицы.[5]
В общем случае, для матриц более высоких порядков (выше 2-го порядка) определитель можно вычислить, применив следующую рекурсивную формулу:
, где — дополнительный минор к элементу . Эта формула называется разложением по строке.
Легко доказать, что при транспонировании определитель матрицы не изменяется (иными словами, аналогичное разложение по первому столбцу также справедливо, то есть даёт такой же результат, как и разложение по первой строке):
Доказательство
Пусть .
Докажем, что по индукции.
Видно, что для матрицы это верно:
Обобщением вышеуказанных формул является разложение детерминанта по Лапласу (Теорема Лапласа), дающее возможность вычислять определитель по любым строкам (столбцам):
Следующие свойства отражают основные результаты теории определителей, применение которых выходит далеко за пределы этой теории:
(Определитель единичной матрицы равен 1);
(Определитель является однородной функцией степени на пространстве матриц размера );
(Определитель матрицы не меняется при её транспонировании);
(Определитель произведения матриц равен произведению их определителей, и — квадратные матрицы одного и того же порядка);
, причём матрица обратима тогда и только тогда, когда обратим её определитель ;
Существует ненулевое решение уравнения тогда и только тогда, когда (или же должен быть нетривиальным делителем нуля в случае, если — не целостное кольцо).
Определитель как функция строк (столбцов) матрицы[править | править код]
При изучении теории определителей полезно иметь в виду, что в основе этой теории лежит техника манипулирования со строками и столбцами матриц, разработанная К. Ф. Гауссом (преобразования Гаусса). Суть этих преобразований сводится к линейным операциям над строками (столбцами) и их перестановке. Эти преобразования достаточно простым образом отражаются на определителе, и при их изучении удобно «расчленить» исходную матрицу на строки (или столбцы) и считать определитель функцией, определённой над наборами строк (столбцов). Далее буквами обозначаются строки (столбцы) матрицы .
1. Определитель — полилинейная функция строк (столбцов) матрицы. Полилинейность означает линейность функции по каждому аргументу при фиксированных значениях остальных аргументов:
2. Определитель — кососимметрическая функция строк (столбцов) матрицы, то есть при перестановке двух строк (столбцов) матрицы её определитель умножается на −1:
3. Если две строки (столбца) матрицы совпадают, то её определитель равен нулю:
Замечание. Свойства 1-3 являются основными свойствами определителя как функции строк (столбцов), они легко доказываются непосредственно из определения. Свойство 2 (кососимметричность) является логическим следствием свойств 1 и 3. Свойство 3 является логическим следствием свойства 2, если в кольце элемент 2 (то есть 1 + 1) не совпадает с нулём и не является делителем нуля. Из свойств 1 и 3 вытекают также следующие свойства:
4. Общий множитель элементов какой-либо строки (столбца) определителя можно вынести за знак определителя (следствие свойства 1).
5. Если хотя бы одна строка (столбец) матрицы нулевая, то определитель равен нулю (следствие свойства 4).
6. Если две (или несколько) строки (столбца) матрицы линейно зависимы, то её определитель равен нулю (следствие свойств 1 и 3).
7. При добавлении к любой строке (столбцу) линейной комбинации других строк (столбцов) определитель не изменится (следствие свойств 1 и 6).
Фактом, имеющим фундаментальное значение, является универсальность определителя как полилинейной кососимметрической функции полного ранга, аргументами которой являются элементы конечномерного векторного пространства (или -модуля с конечным базисом). Справедлива следующая
Теорема. Пусть — свободный -модуль ранга (-мерное векторное пространство над , если — поле). Пусть — -значная функция на , обладающая свойствами 1-3. Тогда при выборе базиса пространства существует константа такая, что при всех значениях справедливо равенство:
,
где — столбец координат вектора относительно базиса .
Доказательство
Разложим векторы по базису : .
Тогда им будут соответствовать следующие столбцы: .
В силу полилинейности функции
В силу свойства 3, если среди индексов есть совпадающие, то
.
В противном случае, в силу кососимметричности (свойства 2), получается:
Пусть — три вектора в пространстве . Они порождают параллелепипед, вершины которого лежат в точках с радиус-векторами . Этот параллелепипед может быть вырожден, если вектора компланарны (лежат в одной плоскости, линейно зависимы).
Функция ориентированного объёма определяется как объём параллелепипеда, порождённого этими векторами, и взятый со знаком «+», если тройка векторов положительно ориентирована, и со знаком «-», если она ориентирована отрицательно.
Функция полилинейна и кососимметрична. Свойство 3, очевидно, выполнено. Для доказательства полилинейности этой функции, достаточно доказать её линейность по вектору . Если вектора линейно зависимы, значение будет нулевым независимо от вектора , и значит, линейно зависящим от него. Если вектора линейно независимы, обозначим через вектор единичной нормали к плоскости векторов , такой, что . Тогда ориентированный объём параллелепипеда равен произведению площади основания, построенного на векторах и независящего от вектора , и алгебраической величины проекции вектора на нормаль к основанию, которая равна скалярному произведению и является величиной, линейно зависящей от вектора . Линейность по доказана, и аналогично доказывается линейность по остальным аргументам.
Применяя теорему об универсальности определителя как кососимметрической полилинейной функции, получаем, что при выборе ортонормированного базиса пространства
,
где — координаты векторов в выбранном базисе.
Таким образом, определитель матрицы коэффициентов векторов относительно ортонормированного базиса имеет смысл ориентированного объёма параллелепипеда, построенного на этих векторах.
Всё вышесказанное без существенных изменений переносится на пространство произвольной размерности.
Разложение определителя по строке/столбцу и обращение матриц[править | править код]
Формулы разложения определителя по строке/столбцу позволяют сводить вычисление определителей к рекурсивной процедуре, использующей вычисление определителей меньших порядков. Для вывода этих формул сгруппируем и просуммируем в формуле для определителя матрицы , с учётом равенства , все ненулевые члены, содержащие элемент . Эта сумма равна:
,
где — матрица, получаемая из удалением строки с номером и столбца с номером .
Так как произвольный элемент можно перестановками соответствующего столбца вправо и перестановками соответствующей строки вниз переместить в правый нижний угол матрицы, причём дополнительная к нему матрица сохранит свой вид, то сумма всех членов в разложении определителя, содержащих , будет равна
.
Величина называется алгебраическим дополнением элемента матрицы .
Учитывая, что каждый член разложения определителя с ненулевым коэффициентом содержит ровно один элемент из i-ой строки, можно разложить определитель по членам этой строки:
— Формула разложения определителя по i-ой строке
Аналогично, учитывая, что каждый член разложения определителя с ненулевым коэффициентом содержит ровно один элемент из j-го столбца, можно разложить определитель по членам этого столбца:
— Формула разложения определителя по j-ому столбцу
Если элементы k-й строки матрицы скопировать в i-ую строку, её определитель станет равен нулю, а по формуле разложения определителя по i-ой строке получится:
— Формула «фальшивого» разложения определителя по i-ой строке ().
Аналогично для столбцов:
— Формула «фальшивого» разложения определителя по j-ому столбцу ()
Полученные формулы полезно записать в матричном виде. Введём матрицу алгебраических дополнений к элементам матрицы : .
Тогда, согласно с полученными формулами,
.
Следствие 1 (Критерий обратимости матриц). Квадратная матрица является обратимой тогда и только тогда, когда — обратимый элемент кольца , при этом .
Следствие 2. Если произведение матриц равно нулю , и матрица — квадратная, тогда .
Решение систем линейных алгебраических уравнений с помощью определителей[править | править код]
Формула Крамера позволяет выразить решение системы линейных алгебраических уравнений в виде отношения определителей, в знаменателе которого стоит определитель системы, а в числителе — определитель матрицы системы, в которой столбец коэффициентов при соответствующей переменной заменён на столбец из правых частей уравнений.
Формула Крамера. Пусть задана система линейных алгебраических уравнений в матричном виде: , где — матрица коэффициентов системы размера , — столбец из правых частей уравнений системы, и вектор — решение этой системы. Тогда при любых справедливо равенство:
Доказательство
Обозначим через сумму и введём
матрицу и вектор .
Тогда и согласно следствию 2 из предыдущего раздела .
Но так как одна из компонент вектора равна -1, это означает, что . Утверждение доказано, так как
Из этой формулы следует, в частности, что если — не вырожден (не является нулём или делителем нуля), система может обладать не более чем одним решением, а если определитель ещё и обратим, то система обладает единственным решением.
Одной из важнейших теорем в теории определителей является следующая теорема о решениях однородной системы линейных уравнений.
Теорема. Пусть — поле. Однородная система линейных уравнений имеет нетривиальное (ненулевое) решение тогда и только тогда, когда определитель матрицы коэффициентов равен нулю: .
Доказательство
Необходимость условия содержится в следствии 2 предыдущего раздела. Докажем необходимость.
Если матрица — нулевая, любой вектор является решением.
Пусть — максимальный невырожденный минор в матрице имеет размеры . Не ограничивая общности, считаем, что этот минор образован первыми r строками и столбцами (иначе перенумеруем переменные и переставим уравнения в другом порядке.)
Введём вектора и . Тогда первые r уравнений системы в матричном виде записываются так:
Поскольку матрица обратима, любому значению соответствует единственный вектор , удовлетворяющий этим уравнениям. Покажем, что при этом остальные уравнения будут выполнены автоматически. Пусть .
Введём две матрицы:
и .
У матрицы все столбцы являются частями столбцов из матрицы , а последний столбец есть линейная комбинация столбцов матрицы с коэффициентами , поэтому в силу линейности определителя по столбцам есть линейная комбинация определителей миноров матрицы размера . Так как — максимальный по размеру невырожденный минор, все миноры большего размера имеют нулевой определитель, поэтому .
Из соотношения следует, что , где столбец . Поэтому .
Тогда . И так как , то j-ое уравнение системы тоже выполнено. ■
Данная теорема используется, в частности, для нахождения собственных значений и собственных векторов матриц.
Критерий полноты и линейной независимости системы векторов[править | править код]
Тесно связанными с понятием определителя является понятие линейной зависимости и полноты систем векторов в векторном пространстве.
Пусть — поле, — векторное пространство над с конечным базисом . Пусть задан ещё набор из векторов . Их координатами относительно заданного базиса являются коэффициенты разложения . Составим (квадратную) матрицу . Справедлива теорема:
Теорема (Критерий полноты и линейной независимости системы векторов).
(1) Система векторов линейно зависима тогда и только тогда, когда .
(2) Система векторов полна, тогда и только тогда, когда матрица не вырождена ().
Доказательство
(1) Доказательство основано на том, что вектор имеет столбец координат, равный , где .
Если , то . Тогда и если отличен от нуля, то .
Обратно, если , существует ненулевой столбец , такой что . Это означает, что .
(2) Если матрица не вырождена, она обратима. Пусть — произвольный вектор, — столбец его координат, . Тогда . Таким образом, произвольный вектор разложим по системы векторов , что означает её полноту.
Обратно, пусть матрица вырождена. Тогда существует ненулевая строка коэффициентов , такая, что . Это значит, что любой вектор , разложимый по системе векторов , удовлетворяет ограничению . Если какой-то коэффициент отличен от нуля, то базисный вектор не разложим по этой системе векторов, и значит, она не полна.
■
Следствие. В векторном пространстве , обладающем конечным базисом из векторов:
(1) всякая система, состоящая менее, чем из векторов, не полна;
(2) всякая система, состоящая более, чем из векторов, линейно зависима;
(3) всякий базис пространства содержит ровно векторов.
Таким образом, размерность векторного пространства с конечным базисом корректно определена.
Прямые методы вычисления определителя могут быть основаны непосредственно на его определении, как суммы по перестановкам, или на разложении Лапласа по определителям меньшего порядка. Однако такие методы очень неэффективны, так как требуют О(n!) операций для вычисления определителя -го порядка. В то же время они универсальны, применимы в тех случаях, когда элементы матрицы не являются числами (функции, многочлены, дифференциальные формы чётной степени и т. д.), и не требуют выполнения операций деления.
Можно вычислить определитель и выполнив меньшее количество операций умножения; точнее, можно для матрицы обойтись лишь количеством слагаемых, равным количеству разбиений числа вместо , например 5 вместо 6 для [6]:
.
Один из наиболее быстрых численных методов вычисления определителя заключается в простой модификации метода Гаусса. Следуя методу Гаусса, произвольную матрицу можно привести к ступенчатому виду (Верхнетреугольная матрица), используя лишь две следующие операции над матрицей — перестановку двух строк и добавление к одной из строк матрицы другой строки, умноженной на произвольное число. Из свойств определителя следует, что вторая операция не изменяет определителя матрицы, а первая лишь меняет его знак на противоположный. Определитель матрицы, приведённой к ступенчатому виду, равен произведению элементов на её диагонали, так как она является треугольной, поэтому определитель исходной матрицы равен:
где — число перестановок строк, выполненных алгоритмом, а — ступенчатая форма матрицы , полученная в результате работы алгоритма. Сложность этого метода, как и метода Гаусса, составляет , для его реализации необходимо использование операции деления.
Определитель можно вычислить, зная LU-разложение матрицы. Если , где и — треугольные матрицы, то . Определитель треугольной матрицы равен просто произведению её диагональных элементов.
Если доступен алгоритм, выполняющий умножение двух матриц порядка за время, где , для некоторого , то определитель матрицы порядка может быть вычислен за время .[7] В частности это означает, что, используя для умножения матриц алгоритм Копперсмита — Винограда, определитель можно вычислить за время .
↑Houston, Robin; Goucher, Adam P. & Johnston, Nathaniel (2023), A New Formula for the Determinant and Bounds on Its Tensor and Waring Ranks, arΧiv:2301.06586 [math.CO].
↑J. R. Bunch and J.E. Hopcroft. Triangular factorization and inversion by fast matrix multiplication, Mathematics of Computation, 28 (1974) 231—236.