This paper presents the problem of robust and nonfragile stabilization of nonlinear systems described by multivariable Hammerstein models. The objective is focused on the design of a nonfragile feedback controller such that the resulting closed-loop system is globally asymptotically stable with robust
H
∞
disturbance attenuation in spite of controller gain variations. First, the parameters of linear and nonlinear blocks characterizing the multivariable Hammerstein model structure are separately estimated by using a subspace identification algorithm. Second, approximate inverse nonlinear functions of polynomial form are proposed to deal with nonbijective invertible nonlinearities. Thereafter, the Takagi–Sugeno model representation is used to decompose the composition of the static nonlinearities and their approximate inverses in series with the linear subspace dynamic submodel into linear fuzzy parts. Besides, sufficient stability conditions for the robust and nonfragile controller synthesis based on quadratic Lyapunov function,
H
∞
criterion, and linear matrix inequality approach are provided. Finally, a numerical example based on twin rotor multi-input multi-output system is considered to demonstrate the effectiveness.