We extend the framework of Neuts' matrix analytic approach to a reflected process generated by a discrete time multidimensional Markov additive process. This Markov additive process has a general background state space and a real vector valued additive component, and generates a multidimensional reflected process. Our major interest is to derive a closed form formula for the stationary distribution of this reflected process. To this end, we introduce a real valued level, and derive new versions of the Wiener-Hopf factorization for the Markov additive process with the multidimensional additive component. In particular, itis represented by moment generating functions, and we consider the domain for it to be valid.</br>Our framework is general...