Multilinear model approach turns out to be an ideal candidate for dealing with nonlinear systems control problem. However, how to identify the optimal active state subspace of each linear subsystem is an open problem due to that the closed-loop performance of nonlinear systems interacts with these subspaces ranges. In this paper, a new systematic method of integrated state space partition and optimal control of multi-model for nonlinear systems based on hybrid systems is initially proposed, which can deal with the state space partition and associated optimal control simultaneously and guarantee an overall performance of nonlinear systems consequently. The proposed method is based on the framework of hybrid systems which synthesizes the multilinear model, produced by nonlinear systems, in a unified criterion and poses a two-level structure. At the upper level, the active state subspace of each linear subsystem is determined under the optimal control index of a hybrid system over infinite horizon, which is executed off-line. At the low level, the optimal control is implemented online via solving the optimal control of hybrid system over finite horizon. The finite horizon optimal control problem is numerically computed by simultaneous method for speeding up computation. Meanwhile, the model mismatch produced by simultaneous method is avoided by using the strategy of receding-horizon. Simulations on CSTR (Continuous Stirred Tank Reactor) confirm that a superior performance can be obtained by using the presented method. (C) 2014 Elsevier Ltd. All rights reserved.