Revision history [back]
 | 1 | initial version |
There is a theorem that says that every doubly stochastic matrix is a linear combination of permutation matrices.
If you want random stochastic matrices, generate random linear combinations of permutations. See where that gets you.
It is called the Birkhoff Von Neumann theorem, by the way.
And there is this, too:
http://www.sagemath.org/doc/reference/matrices/sage/matrix/matrix2.html#sage.matrix.matrix2.Matrix.as_sum_of_permutations
Nathann
