This paper extends the application of discrete orthogonal Laguerre functions (DOLFs) to the analysis and parameter identification of linear time-invariant discrete-time systems described by difference equations with non-zero initial conditions. The extension is achieved through deriving a relationship between the DOLF expansion coefficients of a discrete function y(k + 1) and those of the delayed function y(k).