It is important to predict the Mach number in transonic wind tunnel system. According to the aerodynamic mechanism, the Mach number is indirectly calculated from the Total Pressure (TP) in the stilling chamber and the Static Pressure (SP) in the test section. The high-dimensional input features and large-scale data are the main difficulties to build the TP and the SP models. The fixed-size LS-SVM is a popular method to build a nonlinear model for a large-scale problem. However, it is difficult to further improve the sparsity in the high-dimensional input space. Based on the multivariate fuzzy Taylor theorem, the Feature Subsets Ensemble (FSE) method is proposed to deal with the high-dimensional problem in this paper. The set of direct, exhaustive, independent feature-space subdivisions forms the basis to develop FSEs. In the EFS sub-models are learned using substantially low-dimensional data sets and characterized by low complexity. The TP and the SP are estimated with the FSE based ensemble fixed-size LS-SVMs. Experiments show that the FSEs speed up both training and testing time that would otherwise be infeasible for individual, Bagging and Random Subspace. FSE models meet the requirements of the forecasting speed, the accuracy and generalization of the Mach number prediction.