In the present invention, a correction amount calculation unit 30 acquires operation amount time-series data and multi-point temperature time-series data for a control subject when a target value of an input channel is sequentially changed. An operation amount impact matrix Cmv and a temperature impact matrix Ctemp, in which unit pulse response time-series data, respectively obtained on the basis of the time-series data for the operation amount and the temperature, are obtained. Operation amount time-series data and multi-point temperature time-series data when a disturbance is applied are acquired, and a known operation amount vector Mref and a known temperature vector Tref are obtained. Using a function based on the distribution in respect to an average temperature of predicted multi-point temperatures of the control subject as an evaluation function for minimizing, and setting as a constraint condition the operation amount being kept within a predetermined range, the parameters of the evaluation function and the constraint condition are calculated from the obtained data, and a correction amount vector θ of a target value for minimizing the evaluation function under the constraint condition is calculated. The target value is changed in accordance with the calculated correction amount vector θ to control the control subject.