The inclusion complexes of β-cyclodextrin (β-CD) with prednisolone 1, ethinyloestradiol 2 and estriol 3 in aqueous solution were investigated using ¹H NMR and molecular modelling. The NMR spectra of the steroids studied in the presence of β-CD are fully assigned and interpreted by means of 2D GCOSY and NOESY spectra. The parallel interpretation of β-CD chemical shift changes and dipolar contacts allows the mode of binding to be established. On the basis of ROESY data, the “low resolution” β-cyclodextrin complexes of 1, 2 and 3 have been determined by multistep restrained molecular dynamics calculations. Calculated structures of β-cyclodextrin complexes with 1, 2 and 3 fully agree with experimental data. Combined approaches allow the distinction of weak nonspecific binding for 1 as compared to stronger, “through cavity”, inclusion established for 2 and 3.