@PengOne: тем временем вы редактировали ответ. Вы сказали: «поскольку det (A) = 1, собственные значения равны + -1». Кроме того, вычисление жордановой формы происходит намного медленнее, чем вычисление обратной матрицы. Я не знаю о стабильности жордановой формы, но подозреваю, что она сильно зависит от числа условий A, которое, вероятно, будет очень высоким, поскольку матрица велика. Но если все собственные значения равны единице, то A = I + N с N нильпотентами и трюк с сериями работает (хотя в наихудшем случае 5000 слагаемых, что дает ту же сложность O (N ^ 3), что и простая инверсия).

я есть много больших (около 5000 x 5000) матриц, которые мне нужно инвертировать в Matlab. На самом деле мне нужно обратное, поэтому я не могу использовать вместо этого mldivide, который намного быстрее для решения Ax = b только для одного ...

Three.js: Показать мировые координатные оси в углу сцены

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

Улучшение плохо обусловленной матрицы

У меня плохо обусловленная матрица, чьяrcond() [http://www.mathworks.com/help/matlab/ref/rcond.html]близка к нулю, и, следовательно, обратная сторона этой матрицы не оказывается правильной. Я пытался использоватьpinv() ...

Почему выходные данные inv () и pinv () не равны в Matlab и Octave?

Я заметил, что если A матрица NxN и имеет обратную матрицу. Но то, что выводят функции inv () и pinv (), отличается. - Моя среда - Win7x64 SP1, Matlab R2012a, Cygwin Octave 3.6.4, FreeMat 4.2 Посмотрите на примеры из Octave: A = rand(3,3) A = ...

Обратная матрица в R

Мне было интересно, каков ваш рекомендуемый способ вычисления обратной матрицы? Пути, которые я нашел, кажутся неудовлетворительными. Например,

Большая матричная инверсия

Я смотрю на обратную матрицу большого размера, обычного размера 1000 x 1000, но иногда превышающую 100000 x 100000 (которая в настоящее время дает сбой из-за...

Левый обратный в numpy или scipy?

в питоне

инвертирование матрицы 4х4

Я ищу пример реализации кода о том, как инвертировать матрицу 4x4. я знаю, что есть гауссовское уничтожение, декомпозиция LU и т. д., но вместо того, чтобы с...

Простой обратный код матрицы 3х3 (C ++)

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