A novel method is developed incorporating a special input sequence for a general class of observable and reachable bilinear system identification. An order reduction scheme is introduced to determine the system matrices. This method makes a significant improvement to its predecessors by removing the observability requirement for the linear part of the bilinear systems. Numerical examples are given to demonstrate the method developed in this paper.