SummaryBlast furnace ironmaking is a complex nonlinear dynamic process, which has strict requirements on the quality of final molten iron. To improve the control effect of molten iron quality, this paper proposes an optimal tracking control method for molten iron quality based on the bilinear subspace identification and Krotov's method. First, a state space model between key process variables and blast furnace molten iron quality, using the bilinear subspace identification method, is established based on actual industrial data. Then, an augmented system is established combining the state vector and the reference trajectory, and the improving feedback controller is designed by solving a set of Riccati equations, which are established based on Krotov's method. In combination with the designed controller and the established system, the specific implementation steps are proposed. Finally, the effectiveness and advancement of the proposed method are verified based on actual industrial data of a large blast furnace.