Abstract
In this paper, we characterize the discrete Hahn-classical d-orthogonal polynomial sets (and their iterated) by a difference distributional matrix equation satisfied by the associated d-dimensional functional vector, where the coefficient matrices are explicitly given. As a result we describe this set of polynomials via a structure relation. We illustrate the preceding results with some examples.