In this work, we develop a perturbation theory for the approximation of stability spectra (Lyapunov exponents and Sacker-Sell spectrum) for discrete time dynamical systems. The approach is based upon QR based methods that transform linear time varying discrete time systems to an upper triangular...