Interested Article - Определитель

Определи́тель ( детермина́нт ) в линейной алгебре скалярная величина , которая характеризует ориентированное «растяжение» или «сжатие» многомерного евклидова пространства после преобразования матрицей; имеет смысл только для квадратных матриц . Стандартные обозначения определителя матрицы A {\displaystyle A} det ( A ) {\displaystyle \det(A)} , | A | {\displaystyle |A|} , Δ ( A ) {\displaystyle \Delta (A)} .

Определитель квадратной матрицы A {\displaystyle A} размеров n × n {\displaystyle n\times n} , заданной над коммутативным кольцом R {\displaystyle R} , является элементом кольца R {\displaystyle R} . Эта величина определяет многие свойства матрицы A {\displaystyle A} , в частности, матрица обратима тогда и только тогда, когда её определитель является обратимым элементом кольца R {\displaystyle R} . В случае, когда R {\displaystyle R} поле , определитель матрицы A {\displaystyle A} равен нулю тогда и только тогда, когда ранг матрицы A {\displaystyle A} меньше n {\displaystyle n} , то есть когда системы строк и столбцов матрицы A {\displaystyle A} являются линейно зависимыми .

История

Теория определителей возникла в связи с задачей решения систем линейных уравнений .

К понятию определителя близко подошли авторы древнекитайского учебника « Математика в девяти книгах » .

В Европе определители матриц 2 × 2 встречаются у Кардано в XVI веке. Для старших размерностей определение детерминанта дано Лейбницем в 1693 году. Первая публикация принадлежит Крамеру . Теория определителей создана Вандермондом , Лапласом , Коши и Якоби . Термин «определитель» в современном его значении ввёл О. Коши (1815), хотя ранее (1801) «детерминантом» К. Гаусс назвал дискриминант квадратичной формы.

Японский математик Сэки Такакадзу ввёл определители независимо в 1683 году .

Определения

Через перестановки

Для квадратной матрицы A = ( a i j ) {\displaystyle A=(a_{ij})} размера n × n {\displaystyle n\times n} её определитель det A {\displaystyle \det A} вычисляется по формуле:

det A = α 1 , α 2 , , α n ( 1 ) N ( α 1 , α 2 , , α n ) a 1 α 1 a 2 α 2 a n α n {\displaystyle \det A=\sum _{\alpha _{1},\alpha _{2},\ldots ,\alpha _{n}}(-1)^{N(\alpha _{1},\alpha _{2},\ldots ,\alpha _{n})}\cdot a_{1\alpha _{1}}a_{2\alpha _{2}}\dots a_{n\alpha _{n}}} ,

где суммирование проводится по всем перестановкам α 1 , α 2 , , α n {\displaystyle \alpha _{1},\alpha _{2},\ldots ,\alpha _{n}} чисел 1 , 2 , , n {\displaystyle 1,2,\dots ,n} , а N ( α 1 , α 2 , , α n ) {\displaystyle N(\alpha _{1},\alpha _{2},\ldots ,\alpha _{n})} обозначает число инверсий в перестановке α 1 , α 2 , , α n {\displaystyle \alpha _{1},\alpha _{2},\ldots ,\alpha _{n}} .

Таким образом, в определитель входит n ! {\displaystyle n!} слагаемых, которые также называют «членами определителя».

Эквивалентная формула:

det A = i 1 , i 2 , , i n = 1 n ε i 1 i 2 i n a 1 i 1 a 2 i 2 a n i n {\displaystyle \det A=\sum _{i_{1},i_{2},\ldots ,i_{n}=1}^{n}\varepsilon _{i_{1}i_{2}\dots i_{n}}\cdot a_{1i_{1}}a_{2i_{2}}\dots a_{ni_{n}}} ,

где коэффициент ε i 1 i 2 i n {\displaystyle \varepsilon _{i_{1}i_{2}\dots i_{n}}} символ Леви-Чивиты — равен:

0, если не все индексы i 1 , i 2 , , i n {\displaystyle i_{1},i_{2},\ldots ,i_{n}} различны,
1, если все индексы i 1 , i 2 , , i n {\displaystyle i_{1},i_{2},\ldots ,i_{n}} различны и подстановка ( 1 2 n i 1 i 2 i n ) S n {\displaystyle {\begin{pmatrix}1&2&\dots &n\\i_{1}&i_{2}&\dots &i_{n}\end{pmatrix}}\in S_{n}} чётна,
−1, если все индексы i 1 , i 2 , , i n {\displaystyle i_{1},i_{2},\ldots ,i_{n}} различны и подстановка ( 1 2 n i 1 i 2 i n ) S n {\displaystyle {\begin{pmatrix}1&2&\dots &n\\i_{1}&i_{2}&\dots &i_{n}\end{pmatrix}}\in S_{n}} нечётна.

Аксиоматическое построение (определение на основе свойств)

Понятие определителя может быть введено на основе его свойств. А именно, определителем вещественной матрицы называется функция det : R n × n R {\displaystyle \det :\mathbb {R} ^{n\times n}\rightarrow \mathbb {R} } , обладающая следующими тремя свойствами :

  1. det ( A ) {\displaystyle \det(A)} — кососимметрическая функция строк (столбцов) матрицы A {\displaystyle A} .
  2. det ( A ) {\displaystyle \det(A)} — полилинейная функция строк (столбцов) матрицы A {\displaystyle A} .
  3. det ( E ) = 1 {\displaystyle \det(E)=1} , где E {\displaystyle E} — единичная n × n {\displaystyle n\times n} -матрица.

Значение определителя матрицы

Для матрицы первого порядка значение детерминанта равно единственному элементу этой матрицы:

Δ = | a 11 | = a 11 {\displaystyle \Delta ={\begin{vmatrix}a_{11}\end{vmatrix}}=a_{11}}

Матрицы 2 x 2

Схема расчета определителя матрицы 2×2.
Площадь параллелограмма равна модулю определителя матрицы, образованной векторами — сторонами параллелограмма.

Для матрицы 2 × 2 {\displaystyle 2\times 2} определитель вычисляется как:

Δ = | a c b d | = a d b c {\displaystyle \Delta ={\begin{vmatrix}a&c\\b&d\end{vmatrix}}=ad-bc}

Эта матрица A может быть рассмотрена как матрица линейного отображения , преобразующего единичный квадрат в параллелограмм с вершинами (0, 0) , ( a , b ) , ( a + c , b + d ) , и ( c , d ) .

Абсолютное значение определителя | a d b c | {\displaystyle |ad-bc|} равно площади этого параллелограмма, и, таким образом, отражает коэффициент, на который масштабируются площади при преобразовании A .

Значение определителя со знаком ( ориентированная площадь параллелограмма) помимо коэффициента масштабирования также показывает, выполняет ли преобразование A отражение.

Матрицы 3 x 3

Определитель матрицы 3 × 3 {\displaystyle 3\times 3} можно вычислить по формуле:

Δ = | a 11 a 12 a 13 a 21 a 22 a 23 a 31 a 32 a 33 | = a 11 | a 22 a 23 a 32 a 33 | a 12 | a 21 a 23 a 31 a 33 | + a 13 | a 21 a 22 a 31 a 32 | = {\displaystyle \Delta ={\begin{vmatrix}a_{11}&a_{12}&a_{13}\\a_{21}&a_{22}&a_{23}\\a_{31}&a_{32}&a_{33}\end{vmatrix}}=a_{11}{\begin{vmatrix}a_{22}&a_{23}\\a_{32}&a_{33}\end{vmatrix}}-a_{12}{\begin{vmatrix}a_{21}&a_{23}\\a_{31}&a_{33}\end{vmatrix}}+a_{13}{\begin{vmatrix}a_{21}&a_{22}\\a_{31}&a_{32}\end{vmatrix}}=}
= a 11 a 22 a 33 a 11 a 23 a 32 a 12 a 21 a 33 + a 12 a 23 a 31 + a 13 a 21 a 32 a 13 a 22 a 31 {\displaystyle =a_{11}a_{22}a_{33}-a_{11}a_{23}a_{32}-a_{12}a_{21}a_{33}+a_{12}a_{23}a_{31}+a_{13}a_{21}a_{32}-a_{13}a_{22}a_{31}}

Для более удобного вычисления определителя третьего порядка можно воспользоваться правилом Саррюса или правилом треугольника.

Определитель матрицы, составленной из векторов a , b , c {\displaystyle \mathbf {a} ,\mathbf {b} ,\mathbf {c} } равен их смешанному произведению в правой декартовой системе координат и, аналогично двумерному случаю, представляет собой ориентированный объём параллелепипеда, натянутого на a , b , c {\displaystyle \mathbf {a} ,\mathbf {b} ,\mathbf {c} } . Геометрически определитель матрицы 3х3 означает объем параллелепипеда, построенного на векторах матрицы.

Матрицы N × N

В общем случае, для матриц более высоких порядков (выше 2-го порядка) n × n {\displaystyle n\times n} определитель можно вычислить, применив следующую рекурсивную формулу:

Δ = j = 1 n ( 1 ) 1 + j a 1 j M ¯ j 1 {\displaystyle \Delta =\sum _{j=1}^{n}(-1)^{1+j}a_{1j}{\bar {M}}_{j}^{1}} , где M ¯ j 1 {\displaystyle {\bar {M}}_{j}^{1}} дополнительный минор к элементу a 1 j {\displaystyle a_{1j}} . Эта формула называется разложением по строке .

Легко доказать, что при транспонировании определитель матрицы не изменяется (иными словами, аналогичное разложение по первому столбцу также справедливо, то есть даёт такой же результат, как и разложение по первой строке):

Δ = i = 1 n ( 1 ) i + 1 a i 1 M ¯ 1 i {\displaystyle \Delta =\sum _{i=1}^{n}(-1)^{i+1}a_{i1}{\bar {M}}_{1}^{i}}

Также справедливо и аналогичное разложение по любой строке (столбцу):

Δ = j = 1 n ( 1 ) i + j a i j M ¯ j i {\displaystyle \Delta =\sum _{j=1}^{n}(-1)^{i+j}a_{ij}{\bar {M}}_{j}^{i}}

Обобщением вышеуказанных формул является разложение детерминанта по Лапласу ( Теорема Лапласа ), дающее возможность вычислять определитель по любым k {\displaystyle k} строкам (столбцам):

Δ = 1 j 1 < < j k n ( 1 ) i 1 + + i k + j 1 + + j k M j 1 j k i 1 i k M ¯ j 1 j k i 1 i k {\displaystyle \Delta =\sum _{1\leqslant j_{1}<\ldots <j_{k}\leqslant n}(-1)^{i_{1}+\ldots +i_{k}+j_{1}+\ldots +j_{k}}M_{j_{1}\ldots j_{k}}^{i_{1}\ldots i_{k}}{\bar {M}}_{j_{1}\ldots j_{k}}^{i_{1}\ldots i_{k}}}

Альтернативные методы вычисления

Δ = M ¯ 1 1 M ¯ k k M ¯ 1 k M ¯ k 1 M ¯ 1 , k 1 , k {\displaystyle \Delta ={\frac {{\bar {M}}_{1}^{1}{\bar {M}}_{k}^{k}-{\bar {M}}_{1}^{k}{\bar {M}}_{k}^{1}}{{\bar {M}}_{1,k}^{1,k}}}} .

Основные свойства определителей

Следующие свойства отражают основные результаты теории определителей, применение которых выходит далеко за пределы этой теории:

  1. det E = 1 {\displaystyle \det E=1} (Определитель единичной матрицы равен 1);
  2. det c A = c n det A {\displaystyle \det cA=c^{n}\det A} (Определитель является однородной функцией степени n {\displaystyle n} на пространстве матриц размера n × n {\displaystyle n\times n} );
  3. det A T = det A {\displaystyle \det A^{T}=\det A} (Определитель матрицы не меняется при её транспонировании);
  4. det ( A B ) = det A det B {\displaystyle \det(AB)=\det A\cdot \det B} (Определитель произведения матриц равен произведению их определителей, A {\displaystyle A} и B {\displaystyle B} — квадратные матрицы одного и того же порядка);
  5. det A 1 = ( det A ) 1 {\displaystyle \det A^{-1}=(\det A)^{-1}} , причём матрица A {\displaystyle A} обратима тогда и только тогда, когда обратим её определитель det A {\displaystyle \det A} ;
  6. Существует ненулевое решение уравнения A X = 0 {\displaystyle AX=0} тогда и только тогда, когда det A = 0 {\displaystyle \det A=0} (или же det A {\displaystyle \det A} должен быть нетривиальным делителем нуля в случае, если R {\displaystyle R} — не целостное кольцо).

Определитель как функция строк (столбцов) матрицы

При изучении теории определителей полезно иметь в виду, что в основе этой теории лежит техника манипулирования со строками и столбцами матриц, разработанная К.Ф. Гауссом (преобразования Гаусса). Суть этих преобразований сводится к линейным операциям над строками (столбцами) и их перестановке. Эти преобразования достаточно простым образом отражаются на определителе, и при их изучении удобно "расчленить" исходную матрицу на строки (или столбцы) и считать определитель функцией, определённой над наборами строк (столбцов). Далее буквами A ^ i {\displaystyle {\hat {A}}_{i}} обозначаются строки (столбцы) матрицы A {\displaystyle A} .

1. Определитель — строк (столбцов) матрицы. Полилинейность означает линейность функции по каждому аргументу при фиксированных значениях остальных аргументов:
det ( A ^ 1 , , α A ^ i + β A ^ i , , A ^ n ) = α det ( A ^ 1 , , A ^ i , , A ^ n ) + β det ( A ^ 1 , , A ^ i , , A ^ n ) {\displaystyle \det({\hat {A}}_{1},\ldots ,\alpha {\hat {A}}_{i}+\beta {\hat {A'}}_{i},\ldots ,{\hat {A}}_{n})=\alpha \det({\hat {A}}_{1},\ldots ,{\hat {A}}_{i},\ldots ,{\hat {A}}_{n})+\beta \det({\hat {A}}_{1},\ldots ,{\hat {A'}}_{i},\ldots ,{\hat {A}}_{n})}
2. Определитель — кососимметрическая функция строк (столбцов) матрицы, то есть при перестановке двух строк (столбцов) матрицы её определитель умножается на −1:
det ( , A ^ i , , A ^ j , ) = det ( , A ^ j , , A ^ i , ) {\displaystyle \det(\ldots ,{\hat {A}}_{i},\ldots ,{\hat {A}}_{j},\ldots)=-\det(\ldots ,{\hat {A}}_{j},\ldots ,{\hat {A}}_{i},\ldots)}
3. Если две строки (столбца) матрицы совпадают, то её определитель равен нулю:
A ^ i = A ^ j det ( , A ^ i , , A ^ j , ) = 0 {\displaystyle {\hat {A}}_{i}={\hat {A}}_{j}\Longrightarrow \det(\ldots ,{\hat {A}}_{i},\ldots ,{\hat {A}}_{j},\ldots)=0}

Замечание. Свойства 1-3 являются основными свойствами определителя как функции строк (столбцов), они легко доказываются непосредственно из определения. Свойство 2 (кососимметричность) является логическим следствием свойств 1 и 3. Свойство 3 является логическим следствием свойства 2, если в кольце R {\displaystyle R} элемент 2 (т.е. 1 + 1) не совпадает с нулём и не является делителем нуля. Из свойств 1 и 3 вытекают также следующие свойства:

4. Общий множитель элементов какой-либо строки (столбца) определителя можно вынести за знак определителя (следствие свойства 1).
5. Если хотя бы одна строка (столбец) матрицы нулевая, то определитель равен нулю (следствие свойства 4).
6. Если две (или несколько) строки (столбца) матрицы линейно зависимы, то её определитель равен нулю (следствие свойств 1 и 3).
7. При добавлении к любой строке (столбцу) линейной комбинации других строк (столбцов) определитель не изменится (следствие свойств 1 и 6).

Фактом, имеющим фундаментальное значение, является универсальность определителя как полилинейной кососимметрической функции полного ранга, аргументами которой являются элементы конечномерного векторного пространства V {\displaystyle V} (или R {\displaystyle R} -модуля V {\displaystyle V} с конечным базисом). Справедлива следующая

Теорема. Пусть V {\displaystyle V} — свободный R {\displaystyle R} -модуль ранга n {\displaystyle n} ( n {\displaystyle n} -мерное векторное пространство над R {\displaystyle R} , если R {\displaystyle R} — поле). Пусть φ ( X 1 , X 2 , , X n ) {\displaystyle \varphi (X_{1},X_{2},\dots ,X_{n})} R {\displaystyle R} -значная функция на V n {\displaystyle V^{n}} , обладающая свойствами 1-3. Тогда при выборе базиса E 1 , E 2 , , E n {\displaystyle E_{1},E_{2},\dots ,E_{n}} пространства V {\displaystyle V} существует константа c 0 = φ ( E 1 , E 2 , , E n ) {\displaystyle c_{0}=\varphi (E_{1},E_{2},\dots ,E_{n})} такая, что при всех значениях X 1 , X 2 , , X n {\displaystyle X_{1},X_{2},\dots ,X_{n}} справедливо равенство:
φ ( X 1 , X 2 , , X n ) = c 0 det ( X ^ 1 , X ^ 2 , , X ^ n ) {\displaystyle \varphi (X_{1},X_{2},\dots ,X_{n})=c_{0}\det({\hat {X}}_{1},{\hat {X}}_{2},\dots ,{\hat {X}}_{n})} ,

где X ^ i {\displaystyle {\hat {X}}_{i}} — столбец координат вектора X i {\displaystyle X_{i}} относительно базиса E 1 , E 2 , , E n {\displaystyle E_{1},E_{2},\dots ,E_{n}} .

Одним из важнейших следствий универсальности определителя является следующая теорема о мультипликативности определителя.

Теорема. Пусть A {\displaystyle A} — матрица размера n × n {\displaystyle n\times n} . Тогда det ( A X ) = det A det X {\displaystyle \det(AX)=\det A\cdot \det X} для любой матрицы X {\displaystyle X} размера n × n {\displaystyle n\times n} .

Определитель и ориентированный объём

Пусть a , b , c {\displaystyle {\vec {a}},{\vec {b}},{\vec {c}}} — три вектора в пространстве R 3 {\displaystyle \mathbb {R} ^{3}} . Они порождают параллелепипед, вершины которого лежат в точках с радиус-векторами 0 , a , b , c , a + b , c + a , b + c , a + b + c {\displaystyle {\vec {0}},{\vec {a}},{\vec {b}},{\vec {c}},{\vec {a}}+{\vec {b}},{\vec {c}}+{\vec {a}},{\vec {b}}+{\vec {c}},{\vec {a}}+{\vec {b}}+{\vec {c}}} . Этот параллелепипед может быть вырожден, если вектора a , b , c {\displaystyle {\vec {a}},{\vec {b}},{\vec {c}}} компланарны (лежат в одной плоскости, линейно зависимы).

Функция ориентированного объёма V o l ( a , b , c ) {\displaystyle {\rm {Vol}}({\vec {a}},{\vec {b}},{\vec {c}})} определяется как объём параллелепипеда, порождённого этими векторами, и взятый со знаком "+", если тройка векторов ( a , b , c ) {\displaystyle ({\vec {a}},{\vec {b}},{\vec {c}})} положительно ориентирована, и со знаком "-", если она ориентирована отрицательно. Функция V o l ( a , b , c ) {\displaystyle {\rm {Vol}}({\vec {a}},{\vec {b}},{\vec {c}})} полилинейна и кососимметрична. Свойство 3, очевидно, выполнено. Для доказательства полилинейности этой функции, достаточно доказать её линейность по вектору a {\displaystyle {\vec {a}}} . Если вектора b , c {\displaystyle {\vec {b}},{\vec {c}}} линейно зависимы, значение V o l ( a , b , c ) {\displaystyle {\rm {Vol}}({\vec {a}},{\vec {b}},{\vec {c}})} будет нулевым независимо от вектора a {\displaystyle {\vec {a}}} , и значит, линейно зависящим от него. Если вектора b , c {\displaystyle {\vec {b}},{\vec {c}}} линейно независимы, обозначим через n {\displaystyle {\vec {n}}} вектор единичной нормали к плоскости векторов b , c {\displaystyle {\vec {b}},{\vec {c}}} , такой, что V o l ( n , b , c ) > 0 {\displaystyle {\rm {Vol}}({\vec {n}},{\vec {b}},{\vec {c}})>0} . Тогда ориентированный объём параллелепипеда равен произведению площади основания, построенного на векторах b , c {\displaystyle {\vec {b}},{\vec {c}}} и независящего от вектора a {\displaystyle {\vec {a}}} , и алгебраической величины проекции вектора a {\displaystyle {\vec {a}}} на нормаль к основанию, которая равна скалярному произведению a , n {\displaystyle \langle {\vec {a}},{\vec {n}}\rangle } и является величиной, линейно зависящей от вектора a {\displaystyle {\vec {a}}} . Линейность по a {\displaystyle {\vec {a}}} доказана, и аналогично доказывается линейность по остальным аргументам.

Применяя теорему об универсальности определителя как кососимметрической полилинейной функции, получаем, что при выборе ортонормированного базиса ( e 1 , e 2 , e 3 ) {\displaystyle ({\vec {e}}_{1},{\vec {e}}_{2},{\vec {e}}_{3})} пространства R 3 {\displaystyle \mathbb {R} ^{3}}

V o l ( a , b , c ) = V o l ( e 1 , e 2 , e 3 ) det ( a 1 b 1 c 1 a 2 b 2 c 2 a 3 b 3 c 3 ) = det ( a 1 b 1 c 1 a 2 b 2 c 2 a 3 b 3 c 3 ) {\displaystyle {\rm {Vol}}({\vec {a}},{\vec {b}},{\vec {c}})={\rm {Vol}}({\vec {e}}_{1},{\vec {e}}_{2},{\vec {e}}_{3})\det {\begin{pmatrix}a_{1}&b_{1}&c_{1}\\a_{2}&b_{2}&c_{2}\\a_{3}&b_{3}&c_{3}\end{pmatrix}}=\det {\begin{pmatrix}a_{1}&b_{1}&c_{1}\\a_{2}&b_{2}&c_{2}\\a_{3}&b_{3}&c_{3}\end{pmatrix}}} ,

где a i , b i , c i {\displaystyle a_{i},b_{i},c_{i}} — координаты векторов ( a , b , c ) {\displaystyle ({\vec {a}},{\vec {b}},{\vec {c}})} в выбранном базисе.

Таким образом, определитель матрицы коэффициентов векторов относительно ортонормированного базиса имеет смысл ориентированного объёма параллелепипеда, построенного на этих векторах.

Всё вышесказанное без существенных изменений переносится на пространство R n {\displaystyle \mathbb {R} ^{n}} произвольной размерности.

Разложение определителя по строке/столбцу и обращение матриц

Формулы разложения определителя по строке/столбцу позволяют сводить вычисление определителей к рекурсивной процедуре, использующей вычисление определителей меньших порядков. Для вывода этих формул сгруппируем и просуммируем в формуле для определителя матрицы C = ( c i j ) {\displaystyle C=(c_{ij})} , с учётом равенства ε i 1 i n 1 n = ε i 1 i n 1 {\displaystyle \varepsilon _{i_{1}\dots i_{n-1}n}=\varepsilon _{i_{1}\dots i_{n-1}}} , все ненулевые члены, содержащие элемент c n n {\displaystyle c_{nn}} . Эта сумма равна:

c n n i 1 , i n 1 = 1 n 1 ε i 1 i n 1 n c 1 i 1 c n 1 , i n 1 = c n n i 1 , i n 1 = 1 n 1 ε i 1 i n 1 c 1 i 1 c n 1 , i n 1 = c n n det C n n {\displaystyle c_{nn}\sum _{i_{1},\dots i_{n-1}=1}^{n-1}\varepsilon _{i_{1}\dots i_{n-1}n}c_{1i_{1}}\dots c_{n-1,i_{n-1}}=c_{nn}\sum _{i_{1},\dots i_{n-1}=1}^{n-1}\varepsilon _{i_{1}\dots i_{n-1}}c_{1i_{1}}\dots c_{n-1,i_{n-1}}=c_{nn}\det C_{n}^{n}} ,

где C i j {\displaystyle C_{i}^{j}} — матрица, получаемая из C {\displaystyle C} удалением строки с номером i {\displaystyle i} и столбца с номером j {\displaystyle j} .

Так как произвольный элемент c i j {\displaystyle c_{ij}} можно ( n j ) {\displaystyle (n-j)} перестановками соответствующего столбца вправо и ( n i ) {\displaystyle (n-i)} перестановками соответствующей строки вниз переместить в правый нижний угол матрицы, причём дополнительная к нему матрица сохранит свой вид, то сумма всех членов в разложении определителя, содержащих c i j {\displaystyle c_{ij}} , будет равна

( 1 ) ( n i ) ( 1 ) ( n j ) c i j det C i j = ( 1 ) i + j c i j det C i j = c i j A i j {\displaystyle (-1)^{(n-i)}(-1)^{(n-j)}c_{ij}\det C_{i}^{j}=(-1)^{i+j}c_{ij}\det C_{i}^{j}=c_{ij}A_{i}^{j}} .

Величина A i j = ( 1 ) i + j det C i j {\displaystyle A_{i}^{j}=(-1)^{i+j}\det C_{i}^{j}} называется алгебраическим дополнением элемента матрицы c i j {\displaystyle c_{ij}} .

Учитывая, что каждый член разложения определителя с ненулевым коэффициентом содержит ровно один элемент из i-ой строки, можно разложить определитель по членам этой строки:

det C = j = 1 n c i j A i j {\displaystyle \det C=\sum _{j=1}^{n}c_{ij}A_{i}^{j}} — Формула разложения определителя по i-ой строке

Аналогично, учитывая, что каждый член разложения определителя с ненулевым коэффициентом содержит ровно один элемент из j-го столбца, можно разложить определитель по членам этого столбца:

det C = i = 1 n c i j A i j {\displaystyle \det C=\sum _{i=1}^{n}c_{ij}A_{i}^{j}} — Формула разложения определителя по j-ому столбцу

Если элементы k-й строки матрицы C {\displaystyle C} скопировать в i-ую строку, её определитель станет равен нулю, а по формуле разложения определителя по i-ой строке получится:

0 = j = 1 n c k j A i j {\displaystyle 0=\sum _{j=1}^{n}c_{kj}A_{i}^{j}} — Формула "фальшивого" разложения определителя по i-ой строке ( i k {\displaystyle i\neq k} ).

Аналогично для столбцов:

0 = j = 1 n c i k A i j {\displaystyle 0=\sum _{j=1}^{n}c_{ik}A_{i}^{j}} — Формула "фальшивого" разложения определителя по j-ому столбцу ( j k {\displaystyle j\neq k} )

Полученные формулы полезно записать в матричном виде. Введём матрицу алгебраических дополнений к элементам матрицы C {\displaystyle C} : A = ( A i j ) {\displaystyle A=(A_{i}^{j})} . Тогда, согласно с полученными формулами,

C A T = A T C = ( det C ) E {\displaystyle CA^{T}=A^{T}C=(\det C)\cdot E} .

Следствие 1 (Критерий обратимости матриц). Квадратная матрица C {\displaystyle C} является обратимой тогда и только тогда, когда det C {\displaystyle \det C} — обратимый элемент кольца R {\displaystyle R} , при этом C 1 = ( det C ) 1 A T {\displaystyle C^{-1}=(\det C)^{-1}\cdot A^{T}} .

Следствие 2. Если произведение матриц равно нулю C X = 0 {\displaystyle CX=0} , и матрица C {\displaystyle C} — квадратная, тогда ( det C ) X = 0 {\displaystyle (\det C)\cdot X=0} .

Решение систем линейных алгебраических уравнений с помощью определителей

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

Формула Крамера . Пусть задана система линейных алгебраических уравнений в матричном виде: C X = B {\displaystyle CX=B} , где C = ( c i j ) {\displaystyle C=(c_{ij})} — матрица коэффициентов системы размера n × n {\displaystyle n\times n} , B = ( b 1 , b 2 , , b n ) T {\displaystyle B=(b_{1},b_{2},\dots ,b_{n})^{T}} — столбец из правых частей уравнений системы, и вектор X = ( x 1 , x 2 , , x n ) T {\displaystyle X=(x_{1},x_{2},\dots ,x_{n})^{T}} — решение этой системы. Тогда при любых λ 1 , λ 2 , , λ n R {\displaystyle \lambda _{1},\lambda _{2},\dots ,\lambda _{n}\in R} справедливо равенство:

( λ 1 x 1 + λ 2 x 2 + + λ n x n ) ( det C ) = det ( c 11 c 12 c 1 n b 1 c 21 c 22 c 2 n b 2 c n 1 c n 2 c n n b n λ 1 λ 2 λ n 0 ) {\displaystyle (\lambda _{1}x_{1}+\lambda _{2}x_{2}+\dots +\lambda _{n}x_{n})\cdot (\det C)=-\det {\begin{pmatrix}c_{11}&c_{12}&\dots &c_{1n}&b_{1}\\c_{21}&c_{22}&\dots &c_{2n}&b_{2}\\\dots &\dots &\dots &\dots &\dots \\c_{n1}&c_{n2}&\dots &c_{nn}&b_{n}\\\lambda _{1}&\lambda _{2}&\dots &\lambda _{n}&0\end{pmatrix}}}

Из этой формулы следует, в частности, что если det C {\displaystyle \det C} — не вырожден (не является нулём или делителем нуля), система может обладать не более чем одним решением, а если определитель ещё и обратим, то система обладает единственным решением.

Одной из важнейших теорем в теории определителей является следующая теорема о решениях однородной системы линейных уравнений.

Теорема. Пусть R {\displaystyle R} — поле. Однородная система линейных уравнений C X = 0 {\displaystyle CX=0} имеет нетривиальное (ненулевое) решение тогда и только тогда, когда определитель матрицы коэффициентов равен нулю: det C = 0 {\displaystyle \det C=0} .

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

Критерий полноты и линейной независимости системы векторов

Тесно связанными с понятием определителя является понятие линейной зависимости и полноты систем векторов в векторном пространстве.

Пусть R {\displaystyle R} — поле, V {\displaystyle V} — векторное пространство над R {\displaystyle R} с конечным базисом e 1 , e 2 , , e n {\displaystyle {\vec {e_{1}}},{\vec {e_{2}}},\dots ,{\vec {e_{n}}}} . Пусть задан ещё набор из n {\displaystyle n} векторов v 1 , v 2 , , v n {\displaystyle {\vec {v_{1}}},{\vec {v_{2}}},\dots ,{\vec {v_{n}}}} . Их координатами v i j {\displaystyle v_{ij}} относительно заданного базиса являются коэффициенты разложения v i = v i j e j {\displaystyle {\vec {v_{i}}}=\sum v_{ij}{\vec {e_{j}}}} . Составим (квадратную) матрицу A = ( v i j ) {\displaystyle A=(v_{ij})} . Справедлива теорема:

Теорема (Критерий полноты и линейной независимости системы векторов).

(1) Система векторов v 1 , v 2 , , v n {\displaystyle {\vec {v_{1}}},{\vec {v_{2}}},\dots ,{\vec {v_{n}}}} линейно зависима тогда и только тогда, когда det A = 0 {\displaystyle \det A=0} .
(2) Система векторов v 1 , v 2 , , v n {\displaystyle {\vec {v_{1}}},{\vec {v_{2}}},\dots ,{\vec {v_{n}}}} полна, тогда и только тогда, когда матрица A {\displaystyle A} не вырождена ( det A 0 {\displaystyle \det A\neq 0} ).

Следствие. В векторном пространстве V {\displaystyle V} , обладающем конечным базисом из n {\displaystyle n} векторов:

(1) всякая система, состоящая менее, чем из n {\displaystyle n} векторов, не полна;
(2) всякая система, состоящая более, чем из n {\displaystyle n} векторов, линейно зависима;
(3) всякий базис пространства V {\displaystyle V} содержит ровно n {\displaystyle n} векторов.

Таким образом, размерность векторного пространства V {\displaystyle V} с конечным базисом корректно определена.

Некоторые специальные свойства определителей

Алгоритмическая реализация

  • Прямые методы вычисления определителя могут быть основаны непосредственно на его определении, как суммы по перестановкам, или на разложении Лапласа по определителям меньшего порядка. Однако такие методы очень неэффективны, так как требуют О ( n ! ) операций для вычисления определителя n {\displaystyle n} -го порядка. В то же время они универсальны, применимы в тех случаях, когда элементы матрицы не являются числами (функции, многочлены, дифференциальные формы чётной степени и т. д.), и не требуют выполнения операций деления.
  • Можно вычислить определитель и выполнив меньшее количество операций умножения; точнее, можно для матрицы n × n {\displaystyle n\times n} обойтись лишь количеством слагаемых, равным количеству разбиений числа n {\displaystyle n} вместо n ! {\displaystyle n!} , например 5 вместо 6 для 3 × 3 {\displaystyle 3\times 3} :
    | a 11 a 12 a 13 a 21 a 22 a 23 a 31 a 32 a 33 | = ( a 11 + a 12 ) ( a 21 + a 23 ) ( a 31 + a 32 ) a 12 a 21 ( a 31 + a 32 + a 33 ) a 13 ( a 21 + a 22 + a 23 ) a 31 ( a 11 + a 12 + a 13 ) a 23 a 32 + a 11 a 22 a 33 {\displaystyle {\begin{vmatrix}a_{11}&a_{12}&a_{13}\\a_{21}&a_{22}&a_{23}\\a_{31}&a_{32}&a_{33}\end{vmatrix}}=(a_{11}+a_{12})(a_{21}+a_{23})(a_{31}+a_{32})-a_{12}a_{21}(a_{31}+a_{32}+a_{33})-a_{13}(a_{21}+a_{22}+a_{23})a_{31}-(a_{11}+a_{12}+a_{13})a_{23}a_{32}+a_{11}a_{22}a_{33}} .
  • Один из наиболее быстрых численных методов вычисления определителя заключается в простой модификации метода Гаусса . Следуя методу Гаусса, произвольную матрицу A {\displaystyle A} можно привести к ступенчатому виду (Верхнетреугольная матрица), используя лишь две следующие операции над матрицей — перестановку двух строк и добавление к одной из строк матрицы другой строки, умноженной на произвольное число. Из свойств определителя следует, что вторая операция не изменяет определителя матрицы, а первая лишь меняет его знак на противоположный. Определитель матрицы, приведённой к ступенчатому виду, равен произведению элементов на её диагонали, так как она является треугольной , поэтому определитель исходной матрицы равен:
    det A = ( 1 ) s det A ref , {\displaystyle \det A=(-1)^{s}\cdot \det A_{\mbox{ref}},}
где s {\displaystyle s} — число перестановок строк, выполненных алгоритмом, а A ref {\displaystyle A_{\mbox{ref}}} — ступенчатая форма матрицы A {\displaystyle A} , полученная в результате работы алгоритма. Сложность этого метода, как и метода Гаусса, составляет O ( n 3 ) {\displaystyle O(n^{3})} , для его реализации необходимо использование операции деления.
  • Определитель можно вычислить, зная LU-разложение матрицы. Если A = L U {\displaystyle A=LU} , где L {\displaystyle L} и U {\displaystyle U} — треугольные матрицы, то det A = ( det L ) ( det U ) {\displaystyle \det A=(\det L)(\det U)} . Определитель треугольной матрицы равен просто произведению её диагональных элементов.
  • Если доступен алгоритм, выполняющий умножение двух матриц порядка n {\displaystyle n} за время M ( n ) {\displaystyle M(n)} , где M ( n ) n a {\displaystyle M(n)\geqslant n^{a}} , для некоторого a > 2 {\displaystyle a>2} , то определитель матрицы порядка n {\displaystyle n} может быть вычислен за время O ( M ( n ) ) {\displaystyle O(M(n))} . В частности это означает, что, используя для умножения матриц алгоритм Копперсмита — Винограда , определитель можно вычислить за время O ( n 2.376 ) {\displaystyle O(n^{2.376})} .

Специальные виды определителей

См. также

Примечания

  1. Бронштейн И. Н., Семендяев К. А. Справочник по математике для инженеров и учащихся втузов. — 13-е изд., исправленное. — М. : Наука, 1986.
  2. Э. И. Березкина. Математика древнего Китая. — М. : Наука, 1980.
  3. H. W. Eves. . — Saunders College Publishing, 1990.
  4. Элементы алгебры. — М.: Наука, 1986. — С. 16-23. — Тираж 21 000 экз.
  5. (рус.) . Математика для всех . Дата обращения: 15 декабря 2022. 15 декабря 2022 года.
  6. Houston, Robin; Goucher, Adam P.; Johnston, Nathaniel (2023). "A New Formula for the Determinant and Bounds on Its Tensor and Waring Ranks". arXiv : [ ].
  7. J. R. Bunch and J.E. Hopcroft. Triangular factorization and inversion by fast matrix multiplication, Mathematics of Computation , 28 (1974) 231—236.

Литература

Same as Определитель