Computing Matrix Canonical Forms of Ore Polynomials