Рассматривается проблема получения численных алгоритмов вычисления конечномерной реализации в пространстве состояний стационарных линейных систем дискретного отображения вход-выход, представленного рекуррентной последовательностью матриц. Предложен метод вычисления такой реализации, основанный на LU-разложении ганкелевой матрицы поведения системы.