This paper investigates the problem of asymptotic stability for a class of linear shift-invariant discrete systems with multiple independent variables. We first establish the equivalence of this problem and the problem of robust stability for a class of ordinary linear shift-varying discrete systems with the matrix uncertainty set defined by the coefficient matrices of the original system. On the basis of this equivalence, by using the variational method and the Lyapunov Second Method, necessary and sufficient conditions for asymptotic stability are obtained in different forms for the class of systems considered in the present paper. The parametric classes of Lyapunov functions which define the necessary and sufficient conditions of asymptotic stability are determined. We use the piecewise linear polyhedral Lyapunov functions of the infinity vector norm type to derive an algebraic criterion for asymptotic stability of the given class of discrete systems in the form of solvability conditions of a set of matrix equations.