An explicit and computable solution for coupled partial difference systems appearing when one considers difference approximations of coupled systems of partial differential equations is given. Our approach is based on the consideration of a discrete matrix separation of variables method. By using discrete Fourier series we avoid the problem of solving algebraic linear systems.