Computing Popov Forms of Polynomial Matrices