Problems of analysis and parameter estimation of bilinear systems via block-pulse functions are discussed, and recursive algorithms convenient for computers are presented. Compared with Walsh function approaches, the results are much simpler in form. In this method, any positive integer can be chosen as the number of sub-intervals. Consequently, much computing time and storage can be saved, and the recursive algorithms developed also make it possible to identify on-line the parameters of bilinear systems.