My old professor has a book that starts out going through this in detail.
Look at pages 20 - 22 or so to start. Here's the draft of his book
My Full Answer:
If you ask "My question is: how to translate these rules into matrix equations?" I have to say that they already are matrix equations. Let me explain a bit.
In an early course on Special relativity you may actually write out the Lorentz Transformations as a matrix-column vector operation:
$$\tilde{x} = \Lambda x$$
This gets tedious quickly as you consider more general situations, like boosts with rotations, or more general geometry like when there is curvature. So we don't often write out what the components of $\Lambda$ are except in the easy cases.
When multiplying vectors with a matrix you get a vector. But it is often helpful to think about the components one at a time. This is what the indices do. They do more also, because you've chosen some basis to write them down in. When you write down
$$\tilde{x}^\alpha = \Lambda^{\alpha}_{\phantom{\alpha}\beta} x^\beta$$
It sometimes help to think of this just like
$$ y = Ax \Rightarrow y_i = \sum_j A_{i j} x_j$$
It is essentially the same thing. Seeing this equation should not make you think that it is not a matrix equation. It is! $A_{ i j}$ are the components of the matrix (one at a time!). It is an algebraic expression with indices that represents a real life matrix multiplying a real life vector.
Now one big difference with tensors is that with higher rank you need more multiplies of the $\Lambda$ matrix. This gets confusing because with things like
$$\tilde{R}^{\mu \nu \epsilon \phi} = \Lambda^{\mu}_{\phantom{\mu} \alpha} \Lambda^{\nu}_{\phantom{\nu} \beta} \Lambda^{\epsilon}_{\phantom{\epsilon} \gamma} \Lambda^{\phi}_{\phantom{\phi} \delta } R^{\alpha \beta \gamma \delta}$$
You're sort of still talking about matrix multiply but not with a vector. With a rank 4 tensor. So the algebra works out that you have to do 4 sums. There is little in high school or undergraduate math where this is the case out side of SR or GR. In fact the only thing I've ever seen that even comes close is where you talk about transforming a matrix by right and left multiplication of some other matrix: $$\tilde{M} = U.M.V $$ $$ \tilde{M}_{i l} = \sum_j \sum_k U_{i j} V_{k l} M_{j k}$$
Leaving out the sums really only simplifies things when you already know what you're talking about.
If you are trying to write this down in a familiar way using matrix notation with entries and all of that: GOOD LUCK. This is problematic for many reason but not least of all because usually your higher rank tensors are not really 2D things! So writing them out on paper is hard because you have to somehow represent them and the traditional box shaped matrix won't work (the components are multi-dimensional arrays).
Sorry if my explanation was rambling. I can elaborate but I feel that I may have just made things worse. Let me know what was most confusing about my answer and I'll try to fix it.
Best of luck in a confusing subject.