This paper presents a subspace system identification algorithm for time-invariant structural system. System matrices are estimated from input-output data. The main computational tools are the LQ and the singular value decomposition in subspace system identification. The accuracy and stability of estimated system matrices are influenced by the parameters of Hankel matrix. For stable SI, the optimal parameters of Hankel matrix are determined by considering the computational cost and the accuracy of system matrix.