The bilinear model has been frequently used in fields like control theory and economics to model seismic data. In this article, time-functional variance (TFV) noises are embedded into a specific bilinear model. We propose a generalized autoregressive conditional heteroskedasticity-type maximum likelihood estimator (GMLE) with a sieve method and then provide various inferential techniques based on this GMLE. It is shown that under the finite fourth moment of errors, the GMLE is consistent and asymptotically normally distributed. A simulation stu...