Matrices characterizing the orbital angular momentum (OAM) transformations of deterministic and random light beams by commonly used OAM-modulating optical systems are revealed. Such matrices are the counterparts of the Jones and the Mueller matrices in polarization optics defined for the input and output OAM indices and the radial variables. In particular, matrices of systems leading to OAM mode shift, dispersion, and coupling are discussed.