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

Из Википедии, бесплатной энциклопедии

В математике определитель , — это скалярная величина являющаяся определенной функцией элементов квадратной матрицы . Определитель матрицы A обычно обозначается det( A ) , det A или | А | . Его значение характеризует некоторые свойства матрицы и линейного отображения , представляемого по заданному базису матрицей. В частности, определитель отличен от нуля тогда и только тогда, когда матрица обратима и соответствующее линейное отображение является изоморфизмом . Определителем произведения матриц является произведение их определителей.

Определитель матрицы 2 × 2 равен

а определитель матрицы 3 × 3 равен

Определитель матрицы размера n × n можно определить несколькими эквивалентными способами, наиболее распространенным из которых является формула Лейбница , которая выражает определитель как сумму ( факториал n ) знаковые произведения элементов матрицы. Его можно вычислить с помощью расширения Лапласа , которое выражает определитель как линейную комбинацию определителей подматриц, или с помощью метода исключения Гаусса , который выражает определитель как произведение диагональных элементов диагональной матрицы , полученной последовательностью элементарных вычислений . операции со строками .

Определители также могут определяться некоторыми из их свойств. А именно, определитель — это уникальная функция, определенная на матрицах размера n × n , которая обладает четырьмя следующими свойствами:

  1. Определитель единичной матрицы равен 1 .
  2. Обмен двух строк умножает определитель на −1 .
  3. Умножение строки на число умножает определитель на это число.
  4. Добавление к строке числа, кратного другой строке, не меняет определитель.

Вышеуказанные свойства, относящиеся к строкам (свойства 2–4), могут быть заменены соответствующими утверждениями, относящимися к столбцам.

Определитель квадратной матрицы также может быть определен непосредственно через соответствующий линейный эндоморфизм, представленный матрицей. Поскольку определитель инвариантен относительно преобразования подобия матриц, он на самом деле зависит только от линейного преобразования и в принципе может быть определен полностью «бескоординатным» способом, то есть без использования какого-либо выбора матричного представления.

Определители встречаются во всей математике. Например, для представления коэффициентов в системе линейных уравнений часто используется матрица , а для решения этих уравнений можно использовать определители ( правило Крамера ), хотя другие методы решения в вычислительном отношении гораздо более эффективны. Определители используются для определения характеристического многочлена квадратной матрицы, корнями которого являются собственные значения . В геометрии знаковый n -мерный объем n - мерного параллелепипеда выражается определителем, а определитель линейного эндоморфизма определяет, как ориентация и n преобразуются -мерный объем под действием эндоморфизма. Это используется в исчислении с внешними дифференциальными формами и определителем Якобиана , в частности, для замены переменных в кратных интегралах .

Матрицы два на два [ править ]

Определитель 2 × 2 матрицы обозначается либо « det », либо вертикальными чертами вокруг матрицы и определяется как

Например,

Первые свойства [ править ]

Определитель имеет несколько ключевых свойств, которые можно доказать путем непосредственной оценки определения для -матрицы, и это продолжает сохраняться для определителей более крупных матриц. Они заключаются в следующем: [1] во-первых, определитель единичной матрицы это 1. Во-вторых, определитель равен нулю, если две строки одинаковы:

Это справедливо и в том случае, если два столбца одинаковы. Более того,

Наконец, если какой-либо столбец умножить на некоторое число (т. е. все записи в этом столбце умножаются на это число), определитель также умножается на это число:

Геометрический смысл [ править ]

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

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

Абсолютное значение ad bc — это площадь параллелограмма и, таким образом, представляет собой масштабный коэффициент, с помощью которого площади преобразуются с A. помощью (Параллелограмм, образованный столбцами A , вообще говоря, является другим параллелограммом, но поскольку определитель симметричен относительно строк и столбцов, площадь будет той же.)

Абсолютное значение определителя вместе со знаком становится ориентированной площадью параллелограмма. Ориентированная область такая же, как и обычная площадь , за исключением того, что она отрицательна, когда угол между первым и вторым вектором, определяющим параллелограмм, поворачивается по часовой стрелке (что противоположно направлению, которое можно было бы получить для единичной матрицы ).

Чтобы показать, что ad bc является областью со знаком, можно рассмотреть матрицу, содержащую два вектора u ≡ ( a , b ) и v ≡ ( c , d ) , представляющие стороны параллелограмма. Подписанную область можно выразить как | ты | | в | sin θ для угла θ между векторами, который равен просто основанию, умноженному на высоту, длину одного вектора, умноженную на перпендикулярную составляющую другого. Из-за синуса это уже площадь со знаком, но ее можно выразить более удобно, используя косинус дополнительного угла к перпендикулярному вектору, например u знак равно (− b , a ) , так что | ты | | в | cos θ ' становится рассматриваемой областью со знаком, которую можно определить по образцу скалярного произведения , равного ad - bc , согласно следующим уравнениям:

Объем этого параллелепипеда есть абсолютное значение определителя матрицы, образованной столбцами, построенными из векторов r1, r2 и r3.

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

объект, известный как бивектор С этими идеями связан . В 2D его можно интерпретировать как ориентированный плоский сегмент , образованный путем представления двух векторов, каждый из которых имеет начало координат (0, 0) и координаты ( a , b ) и ( c , d ) . Величина бивектора (обозначается ( a , b ) ∧ ( c , d ) ) — это площадь со знаком , которая также является определителем ad bc . [2]

Если размера n × n действительная матрица A записана через ее векторы-столбцы , затем

Это значит, что отображает единичный n -куб в n -мерный параллелоэдр , определяемый векторами область

Определитель дает знаковый n -мерный объем этого параллелоэдра: и, следовательно, в более общем смысле описывает n -мерный масштабный коэффициент объема линейного преобразования , производимого A . [3] (Знак показывает, сохраняет ли преобразование ориентацию или меняет ее .) В частности, если определитель равен нулю, то этот параллелоэдр имеет нулевой объем и не является полностью n -мерным, что указывает на то, что размерность образа A меньше n . Это означает , что A производит линейное преобразование, которое не является ни однозначным , ни однозначным , и поэтому не является обратимым.

Определение [ править ]

Пусть A квадратная матрица с n строками и n столбцами, так что ее можно записать как

Записи и т. д. во многих случаях являются действительными или комплексными числами. Как обсуждается ниже, определитель также определен для матриц, элементы которых находятся в коммутативном кольце .

Определитель A обозначается det( A ), или его можно обозначить непосредственно через элементы матрицы, написав заключающие столбцы вместо скобок:

Существуют различные эквивалентные способы определения определителя квадратной матрицы A , то есть матрицы с одинаковым количеством строк и столбцов: определитель может быть определен с помощью формулы Лейбница , явной формулы, включающей суммы произведений определенных элементов матрицы. Определитель также можно охарактеризовать как уникальную функцию, зависящую от элементов матрицы, удовлетворяющих определенным свойствам. Этот подход также можно использовать для вычисления определителей путем упрощения рассматриваемых матриц.

Формула Лейбница [ править ]

Матрицы 3 × 3 [ править ]

Формула Лейбница для определителя матрицы 3 × 3 имеет следующий вид:

В этом выражении каждый термин имеет один множитель из каждой строки, все в разных столбцах, расположенных в порядке возрастания строк. Например, bdi имеет b из второго столбца первой строки, d из первого столбца второй строки и i из третьего столбца третьей строки. Знаки определяются тем, сколько транспозиций факторов необходимо, чтобы расположить факторы в порядке возрастания их столбцов (при условии, что члены располагаются слева направо в порядке возрастания порядка строк): положительный при четном числе транспозиций и отрицательный при нечетное число. В примере с bdi однократное транспонирование bd в db дает dbi, три множителя которого взяты из первого, второго и третьего столбцов соответственно; это нечетное количество транспозиций, поэтому термин имеет отрицательный знак.

Правило Сарруса

Правило Сарруса представляет собой мнемонику расширенной формы этого определителя: сумма произведений трех диагональных линий матричных элементов с северо-запада на юго-восток минус сумма произведений трех диагональных линий с юго-запада на северо-восток. восточные строки элементов, когда рядом с ней пишутся копии первых двух столбцов матрицы, как на иллюстрации. Эта схема вычисления определителя матрицы 3×3 не переносится на более высокие измерения.

n × n матрицы [ править ]

Обобщая вышеизложенное на более высокие измерения, определитель Матрица — это выражение, включающее перестановки и их подписи . Перестановка набора является биективной функцией из этого набора в себя, со значениями исчерпать весь набор. Множество всех таких перестановок, называемое симметричной группой , обычно обозначается . Подпись перестановки является если перестановку можно получить четным числом транспозиций (обменами двух записей); в противном случае это

Дана матрица

формула Лейбница для его определителя, используя сигма-обозначение суммы, равна

Используя обозначение числа Пи для продукта, его можно сократить до

.

Символ Леви -Чивита определяется на n - кортежах целых чисел в как 0 , если два целых числа равны, и в противном случае как сигнатура перестановки, определенной кортежем из n целых чисел. С символом Леви-Чивита формула Лейбница принимает вид

где сумма берется по всем n -кортежам целых чисел в [4] [5]

Свойства определителя [ править ]

Характеристика определителя [ править ]

Определитель можно охарактеризовать следующими тремя ключевыми свойствами. Чтобы сформулировать это, удобно рассмотреть -матрица А как составленная из ее столбцы, так называемые

где вектор-столбец (для каждого i ) состоит из элементов матрицы i -го столбца.

  1. , где является единичной матрицей .
  2. Определитель полилинейный : если j -й столбец матрицы записывается как линейная комбинация двух вектор-столбцов v и w и числа r , то определитель A выражается в виде аналогичной линейной комбинации:
  3. Определитель является чередующимся : всякий раз, когда два столбца матрицы идентичны, ее определитель равен 0:

Если определитель определен с использованием формулы Лейбница, как указано выше, эти три свойства можно доказать путем непосредственного изучения этой формулы. Некоторые авторы также подходят к определителю напрямую, используя эти три свойства: можно показать, что существует ровно одна функция, которая присваивает любому -матрица A — число, удовлетворяющее этим трем свойствам. [6] Это также показывает, что этот более абстрактный подход к определителю дает то же определение, что и подход с использованием формулы Лейбница.

Чтобы убедиться в этом, достаточно разложить определитель по многолинейности по столбцам в (огромную) линейную комбинацию определителей матриц, в которой каждый столбец является стандартным базисным вектором. Эти определители равны либо 0 (по свойству 9), либо ±1 (по свойствам 1 и 12 ниже), поэтому линейная комбинация дает приведенное выше выражение в терминах символа Леви-Чивита. Хотя эта характеристика выглядит менее технической, она не может полностью заменить формулу Лейбница при определении определителя, поскольку без нее неясно существование соответствующей функции. [ нужна цитата ]

Непосредственные последствия [ править ]

Эти правила имеют несколько дополнительных последствий:

  • Определитель является однородной функцией , т. е.
    (для матрица ).
  • Перестановка любой пары столбцов матрицы приводит к умножению ее определителя на −1. Это следует из того, что определитель полилинейный и знакопеременный (свойства 2 и 3 выше):
    Эту формулу можно применять итеративно при замене нескольких столбцов. Например
    В более общем смысле любая перестановка столбцов умножает определитель на знак перестановки.
  • Если некоторый столбец может быть выражен как линейная комбинация других столбцов (т. е. столбцы матрицы образуют линейно зависимое множество), определитель равен 0. В особых случаях это включает в себя: если некоторый столбец таков, что все его элементы равны нулю, то определитель этой матрицы равен 0.
  • Добавление скаляра, кратного одному столбцу, к другому столбцу не меняет значение определителя. Это следствие полилинейности и альтернативности: в результате полилинейности определитель изменяется кратно определителю матрицы с двумя равными столбцами, определитель которой равен 0, поскольку определитель чередующийся.
  • Если является треугольной матрицей , т.е. , в любое время или, альтернативно, всякий раз, когда , то его определитель равен произведению диагональных элементов:
    Действительно, такую ​​матрицу можно сократить, соответствующим образом добавив кратные столбцы с меньшим количеством ненулевых элементов к столбцам с большим количеством элементов, до диагональной матрицы (без изменения определителя). Для такой матрицы использование линейности в каждом столбце сводится к единичной матрице, и в этом случае заявленная формула справедлива по самому первому характеризующему свойству определителей. Альтернативно, эту формулу также можно вывести из формулы Лейбница, поскольку единственная перестановка которая дает ненулевой вклад, является тождественной перестановкой.

Пример [ править ]

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

Вычисление определителя матрицы
Матрица

Получено

добавить второй столбец к первому

прибавьте 3 раза третий столбец ко второму

поменять местами первые два столбца

добавлять умножить второй столбец на первый

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

Объединение этих равенств дает

Транспонировать [ править ]

Определитель транспонирования равен определителю A :

.

Это можно доказать, исследуя формулу Лейбница. [7] Это означает, что во всех упомянутых выше свойствах слово «столбец» можно заменить словом «строка». Например, если рассматривать матрицу размера n × n как состоящую из n строк, то определитель представляет собой n -линейную функцию.

Мультипликативность и группы матриц [ править ]

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

Этот ключевой факт можно доказать, заметив, что для фиксированной матрицы , обе части уравнения знакопеременные и полилинейные как функция, зависящая от столбцов . Более того, они оба принимают значение когда является единичной матрицей. Таким образом, вышеупомянутая уникальная характеристика чередующихся полилинейных отображений подтверждает это утверждение. [8]

Матрица с записями в поле обратима именно в том случае , если его определитель не равен нулю. Это следует из мультипликативности определителя и формулы обратного преобразования с участием сопряженной матрицы, упомянутой ниже. В этом случае определитель обратной матрицы определяется выражением

.

В частности, этим свойством по-прежнему обладают произведения и обратные матрицы с ненулевым определителем (соответственно, определителем). Таким образом, набор таких матриц (фиксированного размера над полем ) образует группу, известную как общая линейная группа (соответственно подгруппа , называемая специальной линейной группой . В более общем смысле слово «специальный» обозначает подгруппу другой матричной группы матриц определителя. Примеры включают специальную ортогональную группу (которая, если n равно 2 или 3, состоит из всех матриц вращения ) и специальную унитарную группу .

Поскольку определитель допускает умножение и обратное, на самом деле это гомоморфизм группы из в мультипликативную группу ненулевых элементов . Этот гомоморфизм сюръективен, и его ядро ​​есть (матрицы с определителем один). Следовательно, по первой теореме об изоморфизме это показывает, что является нормальной подгруппой , и что факторгруппа изоморфен .

Формула Коши-Бине является обобщением этой формулы произведения для прямоугольных матриц. Эту формулу также можно преобразовать в мультипликативную формулу для составных матриц , элементы которых являются определителями всех квадратичных подматриц данной матрицы. [9] [10]

Расширение Лапласа [ править ]

Разложение Лапласа выражает определитель матрицы рекурсивно в терминах определителей меньших матриц, известных как их миноры . Несовершеннолетний определяется как определитель -матрица, которая получается из путем удаления -й ряд и -й столбец. Выражение известен как кофактор . Для каждого , имеет место равенство

которое называется разложением Лапласа по i- й строке . Например, разложение Лапласа по первой строке ( ) дает следующую формулу:

Раскрывая детерминанты этих -матрицы возвращают упомянутую выше формулу Лейбница. Аналогично, разложение Лапласа вдоль -й столбец равенство

Разложение Лапласа можно использовать итеративно для вычисления определителей, но этот подход неэффективен для больших матриц. Однако это полезно для вычисления определителей высокосимметричных матриц, таких как матрица Вандермонда.

-членное разложение n Лапласа по строке или столбцу можно обобщить , записав определитель размера n x n в виде суммы члены , каждый из которых является произведением определителя подматрицы x k и k определителя дополнительной подматрицы ( n−k ) x ( n−k ).

Матрица адъюгатов [ править ]

Сопряженная матрица является транспонированием матрицы кофакторов, то есть

Для каждой матрицы имеется [11]

Таким образом, сопряженную матрицу можно использовать для выражения обратной неособой матрицы :

Блочные матрицы [ править ]

Формула определителя а -матрица, приведенная выше, продолжает выполняться при соответствующих дальнейших предположениях для блочной матрицы , т. е. матрицы, состоящей из четырех подматриц. размера , , и , соответственно. Самая простая такая формула, которую можно доказать с помощью формулы Лейбница или факторизации с дополнением Шура , имеет вид

Если обратима что , то из результатов раздела о мультипликативности следует,

что упрощается до когда это -матрица.

Аналогичный результат имеет место, когда обратима, а именно

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

Если блоки представляют собой квадратные матрицы одинакового размера , дальнейшие формулы справедливы. Например, если и ездить на работу (т.е. ), затем [12]

Эта формула была обобщена на матрицы, состоящие из более чем блоков, опять же при соответствующих условиях коммутативности между отдельными блоками. [13]

Для и , справедлива следующая формула (даже если и не ездить на работу) [ нужна цитата ]

Определяющая теорема Сильвестра [ править ]

Детерминантная теорема Сильвестра утверждает, что для A матрица m × n , а для B матрица n × m (так что A и B имеют размеры, позволяющие их умножать в любом порядке, образуя квадратную матрицу):

где I m и I n единичные матрицы m × m и n × n соответственно.

Из этого общего результата вытекает несколько следствий.

  1. Для случая вектора-столбца c и вектора-строки r , каждый из которых имеет m компонентов, формула позволяет быстро вычислить определитель матрицы, которая отличается от единичной матрицы на матрицу ранга 1:
  2. В более общем смысле, [14] для любой обратимой размера m × m матрицы X ,
  3. Для вектора-столбца и строки, как указано выше:
  4. Для квадратных матриц и одинакового размера матрицы и имеют одинаковые характеристические полиномы (следовательно, одинаковые собственные значения).

Сумма [ править ]

Определитель суммы вообще говоря, не выражается через определители A и B. двух квадратных матриц одинакового размера , Однако для положительных полуопределенных матриц , и одинакового размера,

с следствием [15] [16]
И наоборот, если и являются эрмитовыми , положительно определенными и размером , то определитель имеет вогнутую форму й корень; [17] Из этого следует
по однородности.

Тождество суммы для матриц 2×2 [ править ]

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

Удостоверение личности

Это можно показать, расписав каждый термин по компонентам. . Левая сторона

Расширение дает

Члены, являющиеся квадратичными по считаются , и аналогично для , поэтому выражение можно записать

Тогда мы можем записать перекрестные члены как

который можно признать

что завершает доказательство.

У этого есть приложение матричные алгебры. Например, рассмотрим комплексные числа как матричную алгебру. Комплексные числа имеют представление в виде матриц вида

с и настоящий. С , принимая и в приведенном выше тождестве дает

Этот результат следует только из и .

Свойства определителя по отношению к другим понятиям [ править ]

Собственные значения и характеристический полином [ править ]

Определитель тесно связан с двумя другими центральными понятиями линейной алгебры: собственными значениями и характеристическим многочленом матрицы. Позволять быть -матрица со сложными элементами. Тогда по основной теореме алгебры должно иметь ровно n собственных значений . (Здесь подразумевается, что собственное значение с алгебраической кратностью µ встречается в этом списке µ раз.) Тогда оказывается, что определитель A равен произведению этих собственных значений:

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

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

Характеристический полином определяется как [18]

Здесь, является неопределенным многочлена и - единичная матрица того же размера, что и . С помощью этого полинома определители можно использовать для нахождения собственных значений матрицы : они и являются корнями этого многочлена, т. е. теми комплексными числами такой, что

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

быть позитивным для всех между и . [19]

След [ править ]

След tr ( A ) по определению является суммой диагональных элементов A , а также суммой собственных значений. Таким образом, для комплексных матриц A ,

или, для вещественных матриц A ,

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

определитель A определяется выражением

Например, для n = 2 , n = 3 и n = 4 соответственно,

ср. Теорема Кэли-Гамильтона . Такие выражения выводятся из комбинаторных аргументов, тождеств Ньютона или алгоритма Фаддеева-Леверье . То есть для общего det n A = (−1) н c 0 — постоянный член характеристического многочлена со знаком , определяемый рекурсивно из

В общем случае это также можно получить из [20]

где сумма берется по множеству всех целых чисел k l ≥ 0, удовлетворяющих уравнению

Формулу можно выразить через полный экспоненциальный полином Белла от n аргументов s l = −( l – 1)! tr( А л ) как

Эту формулу также можно использовать для нахождения определителя матрицы A я J с многомерными индексами I знак равно ( я 1 , я 2 , ..., я р ) и J знак равно ( j 1 , j 2 , ..., j р ) . Произведение и след таких матриц естественным образом определяются как

Важное тождество произвольной размерности n можно получить из разложения логарифма в ряд Меркатора , когда разложение сходится. Если каждое собственное значение A меньше 1 по абсолютной величине,

где I — единичная матрица. В более общем смысле, если

разлагается в формальный степенной ряд по s, тогда все коэффициенты s м для m > n равны нулю, а оставшийся многочлен равен det( I + sA ) .

Верхняя и нижняя границы [ править ]

Для положительно определенной матрицы A оператор следа дает следующие точные нижние и верхние границы логарифмического определителя:

с равенством тогда и только тогда, A = I. когда Эту взаимосвязь можно вывести с помощью формулы расхождения Кульбака-Лейблера между двумя многомерными нормальными распределениями.

Также,

Эти неравенства можно доказать, выразив следы и определитель через собственные значения. По сути, они отражают хорошо известный факт, что среднее гармоническое меньше среднего геометрического , которое меньше среднего арифметического , которое, в свою очередь, меньше среднеквадратического .

Производная [ править ]

Формула Лейбница показывает, что определитель вещественных (или аналогично комплексных) квадратных матриц является полиномиальной функцией от к . В частности, оно всюду дифференцируемо . Его производную можно выразить по формуле Якоби : [21]

где обозначает сопряжение . В частности, если обратимо, мы имеем

Выражается в виде записей , это

Еще одна эквивалентная формулировка:

,

используя большое обозначение O. Особый случай, когда , единичная матрица, дает

Это тождество используется при описании алгебр Ли, ассоциированных с некоторыми матричными группами Ли . Например, специальная линейная группа определяется уравнением . Приведенная выше формула показывает, что ее алгебра Ли является специальной линейной алгеброй Ли. состоящая из матриц, имеющих нулевой след.

Написание -матрица как где являются векторами-столбцами длины 3, то градиент по одному из трех векторов можно записать как векторное произведение двух других:

История [ править ]

Исторически определители использовались задолго до появления матриц: изначально определитель определялся как свойство системы линейных уравнений . Определитель «определяет», имеет ли система единственное решение (что происходит именно в том случае, если определитель не равен нулю). В этом смысле определители были впервые использованы в китайском учебнике математики «Девять глав математического искусства» (九章算術, китайские учёные, около III века до нашей эры). В Европе решения линейных систем двух уравнений были выражены Кардано в 1545 году с помощью сущности, похожей на определитель. [22]

Собственно определители возникли отдельно от работы Секи Такакадзу в 1683 г. в Японии и параллельно Лейбница в 1693 г. [23] [24] [25] [26] Крамер (1750) без доказательств сформулировал правило Крамера. [27] И Крамер, и Безу (1779) пришли к определителям благодаря вопросу о плоских кривых , проходящих через заданный набор точек. [28]

Вандермонд (1771) впервые признал детерминанты независимыми функциями. [24] Лаплас (1772) дал общий метод разложения определителя по дополнительным минорам : Вандермонд уже привел частный случай. [29] Сразу после этого Лагранж (1773) рассмотрел определители второго и третьего порядка и применил их к вопросам теории исключения ; он доказал множество частных случаев общих тождеств.

Гаусс (1801 г.) сделал следующий шаг. Как и Лагранж, он широко использовал определители в теории чисел . Он ввел слово «определитель» (Лаплас использовал слово «результат»), но не в его нынешнем значении, а скорее применительно дискриминанту кванта к . [30] Гаусс также пришел к понятию обратных (обратных) определителей и очень близко подошел к теореме умножения. [ нужны разъяснения ]

Следующим важным вкладчиком является Бине (1811, 1812), который формально сформулировал теорему о произведении двух матриц из m столбцов и n строк, которая для частного случая m = n сводится к теореме умножения. В тот же день (30 ноября 1812 г.), когда Бине представил свою статью Академии, Коши также представил доклад на эту тему. (См. формулу Коши – Бине .) Здесь он использовал слово «детерминант» в его нынешнем смысле: [31] [32] обобщил и упростил то, что тогда было известно по этому вопросу, улучшил обозначения и дал теорему умножения с доказательством, более удовлетворительным, чем у Бине. [24] [33] С него начинается теория в ее общности.

Якоби (1841) использовал функциональный определитель, который Сильвестр позже назвал якобианом . [34] В своих мемуарах в «Журнале Крелля» за 1841 год он специально рассматривает этот предмет, а также класс знакопеременных функций, которые Сильвестр назвал альтернативами . Примерно во время последних мемуаров Якоби Сильвестр (1839 г.) и Кэли начали свою работу. Кэли в 1841 году ввел современное обозначение определителя с помощью вертикальных черт. [35] [36]

Изучение специальных форм определителей явилось естественным результатом завершения общей теории. Осесимметричные определители изучались Лебегом , Гессе и Сильвестром; персимметричные определители Сильвестра и Ханкеля ; циркулянты Каталана , , Споттисвуда Глейшера и Скотта; косые определители и пфаффианы , в связи с теорией ортогональных преобразований , Кэли; продолжатели Сильвестра; Вронскианы (так называемые Мьюиром ) Кристоффеля и Фробениуса ; составные определители Сильвестра, Рейсса и Пике; Якобианы и гессианцы Сильвестра; и симметричные определители гош Труди . Из учебников по этому предмету первым был учебник Споттисвуда. В Америке трактаты опубликовали Ханус (1886 г.), Уэлд (1893 г.) и Мьюир / Мецлер (1933 г.).

Приложения [ править ]

Правило Крамера [ править ]

Определители можно использовать для описания решений линейной системы уравнений , записанной в матричной форме как . Это уравнение имеет единственное решение если и только если ненулевое значение. В этом случае решение дается по правилу Крамера :

где – матрица, образованная заменой -й столбец по вектор-столбцу . Это немедленно следует из разложения определителя по столбцу, т. е.

где векторы являются столбцами A . Правило также подразумевается тождеством

Правило Крамера может быть реализовано в время, что сравнимо с более распространенными методами решения систем линейных уравнений, такими как LU , QR или разложение по сингулярным значениям . [37]

Линейная независимость [ править ]

Определители можно использовать для характеристики линейно зависимых векторов: равен нулю тогда и только тогда, когда векторы-столбцы (или, что то же самое, векторы-строки) матрицы линейно зависимы. [38] Например, даны два линейно независимых вектора , третий вектор лежит в плоскости , натянутой на первые два вектора, точно, если определитель -матрица, состоящая из трех векторов, равна нулю. Та же идея используется и в теории дифференциальных уравнений : заданные функции (должно быть раз дифференцируемо ), вронскиан определяется как

Оно не равно нулю (для некоторых ) в заданном интервале тогда и только тогда, когда заданные функции и все их производные до порядка линейно независимы. Если можно показать, что вронскиан равен нулю всюду на интервале, то в случае аналитических функций это означает, что данные функции линейно зависимы. См. вронскиан и линейную независимость . Другим примером использования определителя является результирующая , которая дает критерий, когда два многочлена имеют общий корень . [39]

Ориентация основы [ править ]

Определитель можно рассматривать как присвоение номера каждой последовательности из n векторов в R. н , используя квадратную матрицу, столбцы которой являются заданными векторами. Определитель будет отличен от нуля тогда и только тогда, когда последовательность векторов является базой для R н . В этом случае знак определителя определяет, согласуется ли ориентация базиса с ориентацией стандартного базиса или противоположна ему . В случае ортогонального базиса величина определителя равна произведению длин базисных векторов. Например, ортогональная матрица с элементами в R н представляет собой ортонормированный базис в евклидовом пространстве и, следовательно, имеет определитель ±1 (поскольку все векторы имеют длину 1). Определитель равен +1 тогда и только тогда, когда базис имеет ту же ориентацию. Оно равно −1 тогда и только тогда, когда базис имеет противоположную ориентацию.

В более общем смысле, если определитель A положителен, A , сохраняющее ориентацию представляет собой линейное преобразование (если A — ортогональная матрица 2 × 2 или 3 × 3 , это вращение ), а если он отрицательный, A меняет ориентацию. основы.

объема и якобиана Определитель

Как указывалось выше, абсолютное значение определителя вещественных векторов равно объему параллелепипеда, натянутого на эти векторы. Как следствие, если - линейное отображение, заданное умножением на матрицу , и — любое измеримое подмножество , то объём дан кем-то раз больше объёма . [40] В более общем смысле, если линейное отображение представлен -матрица , тогда - габаритный объем дан кем-то:

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

Нелинейная карта превращает небольшой квадрат (слева, красный) в искаженный параллелограмм (справа, красный). Якобиан в точке дает наилучшую линейную аппроксимацию искаженного параллелограмма вблизи этой точки (справа, полупрозрачный белый цвет), а определитель Якобиана дает отношение площади аппроксимирующего параллелограмма к площади исходного квадрата.

Для общей дифференцируемой функции многое из вышесказанного переносится при рассмотрении матрицы Якоби функции f . Для

матрица Якобиана представляет собой матрицу размера n × n , элементы которой задаются частными производными

Его определитель, определитель Якобиана , появляется в многомерной версии интегрирования путем замены : для подходящих функций f и открытого подмножества U в R н (область определения f ), интеграл по f ( U ) от некоторой другой функции φ : R н Р м дан кем-то

Якобиан также встречается в теореме об обратной функции .

Применительно к области картографии определитель можно использовать для измерения скорости расширения карты вблизи полюсов. [41]

Абстрактные аспекты алгебраические

Определитель эндоморфизма [ править ]

Из приведенных выше тождеств, касающихся определителя произведений и обратных матриц, следует, что подобные матрицы имеют один и тот же определитель: две матрицы A и B подобны, если существует обратимая матрица X такая, что A = X −1 БХ . Действительно, повторное применение приведенных выше тождеств дает

Поэтому определитель также называют инвариантом подобия . Определитель линейного преобразования

для некоторого конечномерного векторного пространства V определяется как определитель описывающей его матрицы относительно произвольного выбора базиса в V . В силу инвариантности подобия этот определитель не зависит от выбора базиса V и, следовательно, зависит только от эндоморфизма T .

Квадратные матрицы над коммутативными кольцами [ править ]

Приведенное выше определение определителя с использованием правила Лейбница работает в более общем смысле, когда элементы матрицы являются элементами коммутативного кольца. , например целые числа , в отличие от поля действительных или комплексных чисел. Более того, характеризация определителя как единственного знакопеременного полилинейного отображения, удовлетворяющего условию по-прежнему сохраняется, как и все свойства, вытекающие из этой характеристики. [42]

Матрица обратима (в том смысле, что существует обратная матрица, элементы которой находятся в ) тогда и только тогда, когда его определитель является обратимым элементом в . [43] Для , это означает, что определитель равен +1 или −1. Такая матрица называется унимодулярной .

Определитель, будучи мультипликативным, определяет групповой гомоморфизм

между общей линейной группой (группой обратимых -матрицы с записями в ) и мультипликативная группа единиц в . Поскольку оно соблюдает умножение в обеих группах, это отображение является групповым гомоморфизмом .

Определитель представляет собой естественное преобразование.

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

держит. Другими словами, отображаемая коммутативная диаграмма коммутирует.

Например, определитель комплексно -сопряженной комплексной матрицы (который также является определителем ее сопряженного транспонирования) является комплексно-сопряженным ее определителем, а для целочисленных матриц: приведение по модулю определителя такой матрицы равен определителю матрицы, приведенному по модулю (последний определитель вычисляется с использованием модульной арифметики ). На языке теории категорий определитель — это естественное преобразование между двумя функторами. и . [44] Добавляя еще один уровень абстракции, это можно выразить, сказав, что определитель является морфизмом алгебраических групп , от общей линейной группы до мультипликативной группы .

Внешняя алгебра [ править ]

Определитель линейного преобразования из -мерное векторное пространство или, в более общем смысле, свободный модуль (конечного) ранга над коммутативным кольцом можно сформулировать в бескоординатной форме, рассматривая внешняя сила из . [45] Карта индуцирует линейное отображение

Как одномерна, карта задается путем умножения на некоторый скаляр, т. е. элемент в . Некоторые авторы, такие как ( Bourbaki 1998 ), используют этот факт, чтобы определить определитель как элемент в удовлетворяющее следующему тождеству (для всех ):

Это определение согласуется с более конкретным определением, зависящим от координат. Это можно показать, используя единственность полилинейной знакопеременной формы на -кортежи векторов в . По этой причине самая высокая ненулевая внешняя мощность (в отличие от определителя, связанного с эндоморфизмом) иногда также называют определителем и аналогично для более сложных объектов, таких как векторные расслоения или цепные комплексы векторных пространств. Второстепенные элементы матрицы также можно отобразить в этом случае, рассматривая более низкие чередующиеся формы. с . [46]

Обобщения и родственные понятия [ править ]

Детерминанты, рассмотренные выше, допускают несколько вариантов: перманент матрицы определяется как детерминант, за исключением того, что факторы происходящие в правиле Лейбница, опущены. Имманант обобщает оба , вводя характер симметрической группы. в правление Лейбница.

Определители для конечномерных алгебр [ править ]

Для любой ассоциативной алгебры конечномерное как векторное пространство над полем , существует определительное отображение [47]

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

,

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

Бесконечные матрицы [ править ]

Для матриц с бесконечным числом строк и столбцов приведенные выше определения определителя не переносятся напрямую. Например, в формуле Лейбница необходимо вычислить бесконечную сумму (все члены которой являются бесконечными произведениями). Функциональный анализ предоставляет различные расширения определителя для таких бесконечномерных ситуаций, которые, однако, работают только для определенных типов операторов.

Определитель Фредгольма определяет определитель для операторов, известных как операторы ядерного класса, посредством соответствующего обобщения формулы

Другое бесконечномерное понятие определителя — функциональный определитель .

Операторы в алгебрах Неймана фон

Для операторов с конечным фактором можно определить положительный вещественный определитель, называемый определителем Фугледа-Кадисона, используя канонический след. Фактически, каждому состоянию следа в алгебре фон Неймана соответствует понятие определителя Фугледа–Кадисона.

Связанные понятия для некоммутативных колец [ править ]

Для матриц над некоммутативными кольцами полилинейность и альтернирующие свойства несовместимы при n ≥ 2 , [48] поэтому в этой ситуации нет хорошего определения определителя.

Для квадратных матриц с элементами в некоммутативном кольце возникают различные трудности с определением определителей аналогично тому, как это делается для коммутативных колец. Формуле Лейбница можно придать смысл при условии, что указан порядок произведения, и аналогично для других определений определителя, но некоммутативность тогда приводит к потере многих фундаментальных свойств определителя, таких как мультипликативное свойство или что определитель не изменится при транспонировании матрицы. Над некоммутативными кольцами не существует разумного понятия полилинейной формы (существование ненулевой билинейной формы [ объяснить ] с регулярным элементом R R в качестве значения некоторой пары аргументов следует, что коммутативен ). Тем не менее, были сформулированы различные понятия некоммутативного определителя, которые сохраняют некоторые свойства определителей, особенно квазидетерминантов и определителя Дьедонне . Для некоторых классов матриц с некоммутативными элементами можно определить определитель и доказать теоремы линейной алгебры, очень похожие на их коммутативные аналоги. Примеры включают q -детерминант на квантовых группах, определитель Капелли на матрицах Капелли и березинский на суперматрицах (т. е. матрицах, элементы которых являются элементами - градуированные кольца ). [49] Матрицы Манина образуют класс, наиболее близкий к матрицам с коммутативными элементами.

Расчет [ править ]

Определители в основном используются в качестве теоретического инструмента. Они редко вычисляются явно в числовой линейной алгебре , где для таких приложений, как проверка обратимости и поиск собственных значений, определитель в основном вытесняется другими методами. [50] Однако вычислительная геометрия часто использует вычисления, связанные с определителями. [51]

Хотя определитель можно вычислить напрямую с помощью правила Лейбница, этот подход крайне неэффективен для больших матриц, поскольку эта формула требует вычисления ( факториал ) продукция для -матрица. Таким образом, количество требуемых операций растет очень быстро: оно порядка . Разложение Лапласа также неэффективно. Поэтому для расчета определителей были разработаны более сложные методы.

Методы декомпозиции [ править ]

Некоторые методы вычисления записав матрицу как произведение матриц, определители которых легче вычислить. Такие методы называются методами декомпозиции. Примеры включают LU-разложение , QR-разложение или разложение Холецкого (для положительно определенных матриц ). Эти методы имеют порядок , что является значительным улучшением по сравнению с . [52]

Например, LU-разложение выражает как продукт

матрицы перестановок (который имеет ровно один в каждом столбце, а в противном случае нули), нижняя треугольная матрица и верхняя треугольная матрица . Определители двух треугольных матриц и могут быть быстро вычислены, поскольку они являются произведениями соответствующих диагональных элементов. Определитель это просто знак соответствующей перестановки (которая для четного числа перестановок и для нечетного числа перестановок). Как только такое LU-разложение известно для , его определитель легко вычисляется как

Дальнейшие методы [ править ]

Приказ Достигнутая методами декомпозиции, была улучшена различными методами. Если две матрицы порядка можно умножить во времени , где для некоторых , то существует алгоритм, вычисляющий определитель за время . [53] Это означает, например, что Существует алгоритм вычисления определителя, основанный на алгоритме Копперсмита-Винограда . С 2016 года этот показатель был снижен до 2,373. [54]

Помимо сложности алгоритма, для сравнения алгоритмов можно использовать и другие критерии. Специально для приложений, касающихся матриц над кольцами, существуют алгоритмы, вычисляющие определитель без каких-либо делений. (Напротив, исключение Гаусса требует деления.) Один из таких алгоритмов, имеющий сложность основан на следующей идее: перестановки (как в правиле Лейбница) заменяются так называемыми закрытыми упорядоченными блужданиями , в которых несколько элементов могут повторяться. В полученной сумме больше членов, чем в правиле Лейбница, но в процессе некоторые из этих произведений можно использовать повторно, что делает его более эффективным, чем наивные вычисления по правилу Лейбница. [55] Алгоритмы также можно оценивать по их битовой сложности , т. е. по тому, сколько бит точности необходимо для хранения промежуточных значений, возникающих при вычислении. Например, метод исключения Гаусса (или LU-разложение) имеет порядок , но длина промежуточных значений в битах может стать экспоненциально большой. [56] Для сравнения, алгоритм Барейсса представляет собой метод точного деления (поэтому он использует деление, но только в тех случаях, когда эти деления могут быть выполнены без остатка) имеет тот же порядок, но битовая сложность примерно равна битовому размеру исходные записи в матричных временах . [57]

Если определитель A и обратное значение A уже вычислены, лемма об определителе матрицы позволяет быстро вычислить определитель A + uv. Т , где u и v — векторы-столбцы.

Чарльз Доджсон (т.е. Льюис Кэрролл из «Алисы в стране чудес ») изобрел метод вычисления определителей, названный конденсацией Доджсона . К сожалению, этот интересный метод не всегда работает в исходном виде. [58]

См. также [ править ]

Примечания [ править ]

  1. ^ Ланг 1985 , §VII.1
  2. ^ Вильдбергер, Норман Дж. (2010). Эпизод 4 (видеолекция). ВильдЛинАлг. Сидней, Австралия: Университет Нового Южного Уэльса . Архивировано из оригинала 11 декабря 2021 г. – на YouTube.
  3. ^ «Определители и объемы» . textbooks.math.gatech.edu . Проверено 16 марта 2018 г.
  4. ^ МакКоннелл (1957). Приложения тензорного анализа . Дуврские публикации. стр. 10–17 .
  5. ^ Харрис 2014 , §4.7
  6. ^ Серж Ланг , Линейная алгебра , 2-е издание, Аддисон-Уэсли, 1971, стр. 173, 191.
  7. ^ Ланг 1987 , §VI.7, Теорема 7.5.
  8. ^ Альтернативно, Бурбаки 1998 , §III.8, предложение 1 доказывает этот результат, используя функториальность внешней степени.
  9. ^ Horn & Johnson 2018 , §0.8.7
  10. ^ Кунг, Рота и Ян 2009 , с. 306
  11. ^ Хорн и Джонсон 2018 , §0.8.2.
  12. ^ Сильвестр, младший (2000). «Определители блочных матриц» . Математика. Газ . 84 (501): 460–467. дои : 10.2307/3620776 . JSTOR   3620776 . S2CID   41879675 .
  13. ^ Сотанафан, Нат (январь 2017 г.). «Определители блочных матриц с некоммутирующими блоками». Линейная алгебра и ее приложения . 512 : 202–218. arXiv : 1805.06027 . дои : 10.1016/j.laa.2016.10.004 . S2CID   119272194 .
  14. ^ Доказательства можно найти по адресу http://www.ee.ic.ac.uk/hp/staff/dmb/matrix/proof003.html.
  15. ^ Линь, Минхуа; Сра, Суврит (2014). «Вполне сильная супераддитивность обобщенных матричных функций». arXiv : 1410.1958 [ мат.FA ].
  16. ^ Паксой; Туркменский; Чжан (2014). «Неравенства обобщенных матричных функций через тензорные произведения» . Электронный журнал линейной алгебры . 27 : 332–341. дои : 10.13001/1081-3810.1622 .
  17. ^ Серр, Дени (18 октября 2010 г.). «Вогнутость 1 / н над HPD n " . MathOverflow .
  18. ^ Ланг 1985 , §VIII.2, Horn & Johnson 2018 , Def. 1.2.3
  19. ^ Horn & Johnson 2018 , Наблюдение 7.1.2, Теорема 7.2.5.
  20. ^ Доказательство можно найти в Приложении B к Кондратюк Л.А.; Криворученко, М.И. (1992). «Сверхпроводящая кварковая материя цветовой группы SU (2)». Журнал физики А. 344 (1): 99–115. Бибкод : 1992ZPhyA.344...99K . дои : 10.1007/BF01291027 . S2CID   120467300 .
  21. ^ Horn & Johnson 2018 , § 0.8.10
  22. ^ Граттан-Гиннесс 2003 , §6.6
  23. ^ Каджори, Ф. История математики с. 80
  24. ^ Перейти обратно: а б с Кэмпбелл, Х.: «Линейная алгебра с приложениями», страницы 111–112. Эпплтон Сенчури Крофтс, 1971 год.
  25. ^ Евс 1990 , с. 405
  26. ^ Краткая история линейной алгебры и теории матриц: «Краткая история линейной алгебры и теории матриц» . Архивировано из оригинала 10 сентября 2012 года . Проверено 24 января 2012 г.
  27. ^ Кляйнер 2007 , с. 80
  28. ^ Бурбаки (1994 , стр. 59)
  29. ^ Мьюир, сэр Томас, Теория детерминантов исторического порядка развития [Лондон, Англия: Macmillan and Co., Ltd., 1906]. ЖФМ   37.0181.02
  30. ^ Кляйнер 2007 , §5.2
  31. ^ Первое использование слова «определитель» в современном смысле появилось в: Коши, Огюстен-Луи «Память о функциях, которые могут получить только два равных значения и противоположные знаки в результате транспозиций, выполняемых между переменными, которые они ферментация», которая была впервые прочитана в Институте Франции в Париже 30 ноября 1812 года и которая впоследствии была опубликована в Journal de l'Ecole Polytechnique , Cahier 17, Tome 10, страницы 29–112 (1815).
  32. ^ Происхождение математических терминов: http://jeff560.tripod.com/d.html.
  33. ^ История матриц и определителей: http://www-history.mcs.st-and.ac.uk/history/HistTopics/Matrices_and_determinants.html
  34. ^ Евс 1990 , с. 494
  35. ^ Каджори 1993 , Том. II, с. 92, нет. 462
  36. ^ История матричных обозначений: http://jeff560.tripod.com/matrices.html.
  37. ^ Хабгуд и Арелл, 2012 г.
  38. ^ Ланг 1985 , §VII.3
  39. ^ Ланг 2002 , §IV.8
  40. ^ Ланг 1985 , §VII.6, Теорема 6.10.
  41. ^ Лэй, Дэвид (2021). Линейная алгебра и ее приложения, 6-е издание . Пирсон. п. 172.
  42. ^ Даммит и Фут 2004 , §11.4
  43. ^ Даммит и Фут 2004 , §11.4, Теорема 30
  44. ^ Мак Лейн 1998 , §I.4. См. также Естественное преобразование § Определитель .
  45. ^ Бурбаки 1998 , §III.8
  46. ^ Ломбарди и Квитте 2015 , §5.2, Бурбаки 1998 , §III.5
  47. ^ Гарибальди 2004 г.
  48. ^ В некоммутативной ситуации леволинейность (совместимость с левым умножением на скаляры) следует отличать от праволинейности. Предполагая, что линейность в столбцах считается леволинейной, для некоммутирующих скаляров a , b можно было бы получить :
    противоречие. Не существует полезного понятия полилинейных функций над некоммутативным кольцом.
  49. ^ Варадараджан, В.С. (2004), Суперсимметрия для математиков: введение , Американское математическое общество, ISBN  978-0-8218-3574-6 .
  50. ^ «... мы отмечаем, что детерминант, хотя и является теоретически удобным понятием, редко находит полезную роль в числовых алгоритмах», см. Trefethen & Bau III 1997 , лекция 1.
  51. ^ Фисикопулос и Пеньяранда 2016 , §1.1, §4.3
  52. ^ Камареро, Кристобаль (05 декабря 2018 г.). «Простые, быстрые и практичные алгоритмы разложения Холецкого, LU и QR с использованием быстрого умножения прямоугольных матриц». arXiv : 1812.02056 [ cs.NA ].
  53. ^ Банч и Хопкрофт, 1974 г.
  54. ^ Фисикопулос и Пеньяранда 2016 , §1.1
  55. ^ Роте 2001
  56. ^ Фан, Синь Гуй; Хавас, Джордж (1997). «О сложности целочисленного исключения Гаусса в наихудшем случае» (PDF) . Материалы международного симпозиума 1997 года по символьным и алгебраическим вычислениям . ИССАК '97. Кихеи, Мауи, Гавайи, США: ACM. стр. 28–31. дои : 10.1145/258726.258740 . ISBN  0-89791-875-4 . Архивировано из оригинала (PDF) 7 августа 2011 г. Проверено 22 января 2011 г.
  57. ^ Фисикопулос и Пеньяранда 2016 , §1.1, Bareiss 1968
  58. ^ Абелес, Франсин Ф. (2008). «Конденсация Доджсона: историческое и математическое развитие экспериментального метода» . Линейная алгебра и ее приложения . 429 (2–3): 429–438. дои : 10.1016/j.laa.2007.11.022 .

Ссылки [ править ]

Исторические справки [ править ]

Внешние ссылки [ править ]