This paper proposes a new representation of discrete bilinear model by developing its coefficients associated to the input, to the output and to the crossed product on three independent Laguerre orthonormal bases. Compared to classical bilinear model, the resulting model entitled bilinear-Laguerre model ensures a significant parameter number reduction as well as simple recursive representation. However, this reduction is still subject to an optimal choice of the Laguerre poles defining the three Laguerre bases. Therefore, we propose an analytical solution to optimise the Laguerre poles which depend on Fourier coefficients defining the bilinear-Laguerre model, and that are identified using the regularised square error. The identification procedures of the Laguerre poles and Fourier coefficients are combined and carried out on a sliding window to provide an online identification algorithm of the bilinear-Laguerre model. The bilinear-Laguerre model as well as the proposed algorithm are illustrated and tested on a numerical simulation and validated on the continuous stirred tank reactor system.