An industrial model predictive control method realizing dynamic optimization based on linear programming comprises the following steps: 1) carrying out steady state target optimization: computing the expected target value output in a controlled way in the steady state according to the input quantity of the sampling period; and 2) carrying out dynamic optimization: 2.1) obtaining step response models of the industrial process by a least square method, wherein the step response models are established in N-numbered periods; 2.2) predicting the output quantity at each future moment from the k+1 moment, referring to formula (2); 2.3) ensuring the input and the output in the dynamic process to be in own upper and lower bound ranges, referring to formula (3), setting target constraint, referring to formula (4) and creating the following performance indexes, referring to formula (5); forming the dynamic optimization problem by combining the formulas (2), (3), (4) and (5), solving all the control increments deltau(i) by linear programming, wherein i=0, 1,..., N-m-1, and then substituting the control increments into the formula (2) to compute the predictive output values y(k+1), y(k+2),..., y(N). The method simplifies computation, shortens the computation period and is good in practicability.