Блог пользователя 2ndSequence

Absorbing markov chains implementation (Help).
Автор 2ndSequence, история, 3 года назад, По-английски

Greeting,

I was trying to implement that algorithm but i don't know how to calculate F = (I — Q) ^ -1 part, I can not calculate the inverse of the matrix (because i also want to keep the proper fractions as it is)..

Any thoughts? (Or even a blog to read).

  • Проголосовать: нравится
  • +1
  • Проголосовать: не нравится

»
3 года назад, # |
  Проголосовать: нравится 0 Проголосовать: не нравится

the matrix times the comatrix is det times identity. so inverse is comatrix divided by det, it can be computed in O(n^4) (the comatrix). Wikipedia speaks also of LU decomposition maybe it might help and be faster I don't know if it stays in fractions.

→ Ответить
  • »
    »
    3 года назад, # ^ |
      Проголосовать: нравится 0 Проголосовать: не нравится

    Gaussian elimination, also known as row reduction,

    enable to compute the inverse in O(n^3) and stay in the fraction.

    You write identity next to your matrix and you apply it the same operation as to your matrix,

    when your matrix becomes identity, the identity becomes the inverse

    → Ответить
  • »
    »
    3 года назад, # ^ |
      Проголосовать: нравится +5 Проголосовать: не нравится

    Thank you so much for your interest..

    I think i knew about Gaussian but i am pretty sure that it doesn't maintain the proper fraction (not sure... to weak in math but it is an assignment..)

    I will try to use the determinant and see what i can do, also i don't care about complexity, Just maintaining the proper fraction.

    → Ответить