In this paper, an efficient and effective procedure is successfully developed for parameter identification of linear time-invariant multi-delay systems. The proposed framework is based on a hybrid of block-pulse functions and Taylor's polynomials. Two upper error bounds corresponding to hybrid functions are established. The excellent properties of these functions together with the associated operational matrices of integration and delay are utilised to transform the original problem into a system of linear algebraic equations. The least squares method is then implemented for estimation of the unknown parameters. Several numerical experiments are investigated to demonstrate the usefulness and effectiveness of the proposed procedure. Easy implementation, simple operations and accurate solutions are the main features of the suggested approximation scheme.