In the paper, we study entire solutions of the difference equation ψ (z + h) = M (z) ψ (z), z ∈ C{double-struck}, ψ (z) ε C{double-struck}². In this equation, h is a, fixed positive parameter, and M : C{double-struck} →SL (2, C{double-struck}) is a given matrix function. We assume that M(z) is a 2π-periodic trigonometric polynomial. The main aim is to construct the minimal entire solutions, e.i. the solutions with the minimal possible growth simultaneously as for z → -i∞ so for z → +i∞. We show that the monodromy matrices corresponding to the bases made of the minimal solutions are trigonometric polynomials of the same order as the matrix M. This property relates the spectral analysis of the one dimensional difference Schrödinger equations with the potentials being trigonometric polynomials to an analysis of a finite dimensinal dynamical system.