A new approach to the realization problem for fractional discrete-time linear systems is proposed. A procedure for computation of fractional realizations of given transfer matrices is presented and illustrated by numerical examples.