Re: matrix inversion

In article <fuqgr5$6pk$1@xxxxxxxxxxxxxxxxxx>,
Milad Ehtesham <m83ehtesh@xxxxxxxxx> wrote:
I have the following lower triangular matrix, which clearly
has determinant one and is full rank. but matlab gives zero
as its determinant and it cannot correctly inverted by
matlab. why this happens and how can i invert it correctly?

My hypothesis is precision limits and catestrophic cancellation.

Any 16 x 16 (or smaller) submatrix works; at 17 x 17 or larger
you start getting odd answers -- largely 0s, but some very
large and some very small.

rank() starts going wrong at 12 x 12

The matrix is badly conditioned according to rcond(). Even
just a 4 x 4 submatrix has an rcond of about 2E-6
--
"This quitting thing, it's a hard habit to break once you start."
-- Walter Matthau
.

Relevant Pages

  • matrix inversion
    ... has determinant one and is full rank. ... but matlab gives zero ...
    (comp.soft-sys.matlab)
  • Re: Tell whether a matrix is full rank using a calculator
    ... My calculator does not have a built in function that computes the ... Is there another way to know whether a matrix is full rank? ... eigenvalues), or some other method.... ... determinant of your matrix and see whether or not it is equal to 0. ...
    (sci.math)
  • Re: Tell whether a matrix is full rank using a calculator
    ... Is there another way to know whether a matrix is full rank? ... A square matrix is full rank if and only if its determinant is ... determinant of your matrix and see whether or not it is equal to 0. ... If the calculator can do exact determinants, ...
    (sci.math)
  • Re: a big limit of mathematica?
    ... Mathematica's numerical model has been criticized by several people, ... In fact, the determinant of the 10x10, 20x20 etc all come out exactly 0. ... In Matlab the result will ... the 10X10 determinant is exactly zero. ...
    (sci.math.symbolic)
  • Re: Tell whether a matrix is full rank using a calculator
    ... Is there another way to know whether a matrix is full rank? ... A square matrix is full rank if and only if its determinant is ... determinant of your matrix and see whether or not it is equal to 0. ... If the calculator can do exact determinants, ...
    (sci.math)
Loading