To improve prediction accuracy and provide sufficient time to control decision-making, a decomposition-based multi-step forecasting model for rabbit house environmental variables is proposed. Traditional forecasting methods for rabbit house environmental parameters perform poorly because the coupling relationship between sequences is ignored. Using the STL algorithm, the proposed model first decomposes the non-stationary time series into trend, seasonal, and residual components and then predicts separately based on the characteristics of each component. LSTM and Informer are used to predict the trend and residual components, respectively. The aforementioned two predicted values are added together with the seasonal component to obtain the final predicted value. The most important environmental variables in a rabbit house are temperature, humidity, and carbon dioxide concentration. The experimental results show that the encoder and decoder input sequence lengths in the Informer model have a significant impact on the model’s performance. The rabbit house environment’s multivariate correlation time series can be effectively predicted in a multi-input and single-output mode. The temperature and humidity prediction improved significantly, but the carbon dioxide concentration did not. Because of the effective extraction of the coupling relationship among the correlated time series, the proposed model can perfectly perform multivariate multi-step prediction of non-stationary time series.