Explicit solution of second-order delayed discrete equations

A system of inhomogeneous second-order difference equations with linear parts given by noncommutative matrix coefficients are considered. Closed form of its solution is derived by means of newly defined delayed matrix sine/cosine using the Z-transform and determining function.