Two-point correlation functions of spin operators in the minimal models M p,p′ perturbed by the field Φ 13 are studied in the framework of conformal perturbation theory [Nucl. Phys. B 348 (1991) 619]. The first-order corrections for the structure functions are derived analytically in terms of gamma functions. Together with the exact vacuum expectation values of local operators, this gives the short-distance expansion of the correlation functions. The long-distance behaviors of these correlation functions in the case M 2,2n+1 have been worked out using a form-factor bootstrap approach. The results of numerical calculations demonstrate that the short- and long-distance expansions match at the intermediate distances. Including the descendent operators in the OPE drastically improves the convergency region. The combination of the two methods thus describes the correlation functions at all length scales with good precision.