Absorption infrared spectra have been computed for a variety of polycyclic aromatic hydrocarbon molecules of the pyrene family, taking into account anharmonicity and temperature effects, rovibrational quantization, and couplings. The energy levels are described by a second-order perturbative expansion of the rovibrational Hamiltonian in the vibrational and rotational quantum numbers, as relevant for a symmetric-top molecule, with ingredients obtained from quantum chemistry calculations. Multicanonical Monte Carlo simulations are carried out to compute bidimensional IR intensity histograms as a function of total energy and vibrational frequency, which then provide the absorption spectrum at arbitrary temperatures via a Laplace transformation. The main spectral features analyzed for neutral, anionic, and cationic pyrene indicate a strong dependence on temperature, in agreement with existing laboratory experiments, and a significant contribution of rotational degrees of freedom to the overall broadenings. The spectral shifts and broadenings reveal some sensitivity of anharmonicities to the charge and protonation states and, in the case of protonated pyrene and pyrenyl cation, on possible isomers and between aromatic and aliphatic C-H bands. Implications of the present work to the general issue of interstellar emission features are discussed.