The design of a linear parameter varying (LPV) controller for the F-16 longitudinal axes over the entire ight envelope, using a blending methodology that lets an LPV controller preserve regional optimal solutions over each parameter subset and reduces computationalcosts for synthesizing an LPV controller, is presented. Three blending LPV controller synthesis methodologies are applied to control F-16 longitudinal axes. By the use of a function substitution method, a quasi-LPV model of the F-16 longitudinal axes is constructed from the nonlinear equations of motion over the entire ight envelope, including nontrim regions, to facilitate synthesis of LPV controllers for the F-16 aircraft. The nonlinear simulations of the blended LPV controller show that the desired performance and robustness objectives are achieved across all altitude variations. Nomenclature C m 0 = pitch moment aerodynamic coef cient Cm q = pitch moment aerodynamic stability derivative C X , C Z = X and Z force aerodynamic coef cients C Xq , C Z q = aerodynamic stability derivatives N c, S = chord length, ft, and reference area, ft 2 g = gravity constant, ft/s 2 m, I yy = mass, slug, and inertial moments, slug ¢ ft 2 V , N q = velocity, ft/s, and dynamic pressure, psi X a:c: , X c:g: = aerodynamic center position, ft, and center of gravity position, ft ®, q = angle of attack, rad, and pitch rate, rad/s ± e , T = elevator de ection, rad, and thrust, lb µ ,°= pitch angle, rad, and ight-path angle, rad