Two different approaches for describing discretetime fractional linear systems are presented. The first is based on the nabla and delta discrete-time derivatives. In this case suitable exponentials are introduced and used to define discrete Laplace transforms. The second approach is based on the bilinear (Tustin) transformations. For both cases, appropriate algorithms for obtaining the impulse, step, and frequency responses are presented. The state-variable representation is also analysed.