Rotational-vibrational spectroscopy of water in solid noble gas matrices has been studied for many decades. Despite that, discrepancies persist in the literature about the assignment of specific bands. We tackle the involved rotational-vibrational spectrum of the water isotopologues (H2O)-O-16, HD16O, and (D2O)-O-16 with an unprecedented combination of experimental high-resolution matrix isolation infrared (MI-IR) spectroscopy and computational anharmonic vibrational spectroscopy by vibrational configuration interaction (VCI) on high-level ab initio potential energy surfaces. With VCI, the average deviation to gas-phase experiments is reduced from >100 to approximate to 1 cm(-1) when compared to harmonic vibrational spectra. Discrepancies between MI-IR and VCI spectra are identified as matrix effects rather than missing anharmonicity in the theoretical approach. Matrix effects are small in Ne (approximate to 1.5 cm(-1)) and a bit larger in Ar (approximate to 10 cm(-1)). Controversial assignments in Ne MI-IR spectra are resolved, for example, concerning the v(3) triad in HDO. We identify new transitions, for example, the V-2 1(01) <- I-10 transition in D2O and H2O or the v(3) 0(00) <- 1(01) transition in D2O, and reassign bands, for example, the band at 3718.9 cm(-1) that is newly assigned as the 1(10) <- 1(11) transition. The identification and solution of discrepancies for a well-studied benchmark system such as water prove the importance of an iterative and one-hand combination of theory and experiment in the field of high-resolution infrared spectroscopy of single molecules. As the computational costs involved in the VCI approach are reasonably low, such combined experimental/theoretical studies can be extended to molecules larger than triatomics.