In this paper a matrix method for establishing discrete system dynamic equations is presented. Using this method to establish discrete system dynamic equations, netting but the derivative of function, and matrix multiplication are required. This method suits to the holonomic system as well as the linear nonholomic system.