Здравствуйте!
Столкнулся с программой расчета определителя матрицы.
Пример принципа расчета на примере матрицы 2*2.

Попытался сопоставить этот метод с другими - не сопоставляется.
Если считать стандартным простым методом, то можно получить: 1*3-2*4=-5, что сходится с ранее полученным ответом.
Этот принцип расчета встретил в следующей программе.
Текст программы:


В ней написано "прямой ход Гаусса". В методе Гаусса написано, что в прямом ходе производятся преобразования строк системы, приводя ее к ступенчатой или треугольной форме, т.е. матрица должна приводится к ступенчатому виду. В этом же случае я не совсем понимаю, как производится это преобразование. Кто-нибудь мог бы мне объяснить?

Метод Гаусса

Приводите матрицу к треугольной, тогда определитель - это произведение чисел на главной диагонали.

Сообщение отредактировал Taradov Alexander - Oct 7 2011, 12:59

Читал статью в Википедии.
Если свести текст программы в виде алгоритма, то получается вот такой:



А как тут реализуется приведение матрицы к треугольному виду?

Для каждого столбца матрицы обнуляются элементы столбца ниже главной диагонали гауссовским исключением.

Если удобнее представить в "школьной" формулировке, то из каждого уравнения системы выражается одна неизвестная и подставляется во все остальные уравнения. Последнее уравнение таким образом содержит только одну неизвестную, после вычисления ее значения происходит обратная подстановка (снизу вверх) этого значения во все остальные уравнения. Последний шаг для поиска определителя не нужен, только для решения СЛАУ.

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



P.S. Taradov Alexander, Вы обогнали с ответом. Я прекрасно помню школьный курс - я не совсем понимаю для чего деление и умножение в алгоритме и как при помощи него реализуется преобразование матрицы к треугольному виду.

Деление не 2-х элементов, а всей первой строки на a11, это и есть выделение переменной при a11 в свободном виде. В дальнейшем происходит "подстановка" во вторую строку.

PS: Запишите СЛАУ с этой матрицей и проследите за своими дейстивиясми при решении, они 1:1 повторят действия алгоритма.

Сообщение отредактировал Taradov Alexander - Oct 7 2011, 15:06

a11 a12 a13
a21 a22 a23
a31 a32 a33


k=-a31/a21

a11 a12 a13
a21 a22 a23
a31+a21*k a32+a22*k a33+a23*k

a31+a21*k=0

Ну и тд

А для чего запись A[i,j]:=A[i,j]-r; в конце цикла подстановки первой строки во все последующие?

Это и есть подстановка. Именно это действие и должно занулить a21 в первом примере.

Сообщение отредактировал Taradov Alexander - Oct 7 2011, 15:28

Вроде как разобрался.
Сначала детерминнант приняли равным 1. Затем начали цикл в котором определитель - произведение чисел на главной диагонали. В первом цикле поизводится деление первой строки на (k+1)-ый элемент, т.е. каждый раз в каждом цике определенная строка будет приниматься первой и для нее будет выделяться дробь (как показал thermit), которая будет подставляться в последующие уравнения во втором цикле.
Спасибо за помощь!

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

Можно в матлабе попробовать. Если совпадет, то там вроде есть ссылка на алгоритм расчета

Если экзекуция Гаусса переставила строки чётное число раз. Если нечётное - то минус произведение. Т.е. ещё нужно умножить на чётность перестановки строк в методе Гаусса.

При численном вычислении определителя (а это сумма эн-факториал произведений - по определению) результат, очень даже возможно, не влезет ни в какую ограниченную плавающую точку, поэтому в хороших подпрограммах вычисления определителей используют пару чисел в плавающей точке (а,b) для представления результата;, так, что |A| = a * 10^b.