This paper identifies a linear time-invariant mathematical model of each stand of a five-stand tandem cold mill to design a robust $H_infty$ thickness controller by applying input and output data sets to N4SID (Numerical algorithms for Subspace State Space System Identification) method. The input-output data sets describe interstand interference in the process of tandem cold rolling and are obtained from a nonlinear simulator of the tandem cold mill. In result, it is shown that the identified model well approximates the nonlinear model than a Taylor linearized model. Furthermore, uncertainties including roll eccentricity and incoming strip variation are quantitatively analyzed from the plot of maximum singular values.