You can do this using the SymPY.jl package. As a good chunk of matrix functionality is parametric on element type, you can do matrix operations on SymPY's symbolic variables.

See example 6 in: https://github.com/andreasnoack/andreasnoack.github.io/blob/...