A Liouville-Green (WKB) asymptotic approximation theory is developed for some classes of linear second-order difference equations in Banach algebras. The special case of linear matrix difference equations (or, equivalently, of second-order systems) is emphasized. Rigorous and explicitly computable bounds for the error terms are obtained, and this when both, the sequence index and some parameter that may enter the coefficients, go to infinity. A simple application is made to orthogonal matrix polynomials in the Nevai class.