The aromaticity of the ring in variously substituted aniline/anilinium/anilide derivatives in their H-bonded complexes with various Broensted acids and bases was a subject of an analysis based on 332 experimental geometries retrieved from the Cambridge Structural Database and geometries optimized at the B3LYP/6-311+G** and MP2/aug-cc-pVDZ levels of theory. Ab initio modeling was applied to the para-substituted aniline, anilinium cation, and anilide anion derivatives (X = NO, NO2, CN, CHO, H, CH3, OCH3, and OH) and their H-bonded complexes (only for X = NO, NO2, CHO, H, and OH) with B (B = F- and CN-) or HB (HB = HF and HCN). In both cases, the harmonic oscillator model of aromaticity index (HOMA) was used, whereas for computational geometries, additionally, the magnetism-based indices NICS, NICS(1), and NICS(1)zz were also applied (NICS = nucleus-independent chemical shift). There is an equivalent prediction of aromaticity by NICSs and HOMA and approximate monotonic dependences of HOMA and NICS on the C-N bond length. The strongest changes in aromaticity estimated by HOMA and NICSs were found for aniline derivatives with NH2...B and anilide derivatives without and with NH-...HB interactions. The changes observed for two other kinds of interactions, NH2...HB and NH3+...base (for anilinium cations), are much smaller. For all four kinds of interactions, the relationships between ipso-bond angle, mean ipso-ortho bond length, and C-N bond length follow the Bent-Walsh rule.