The nonlinear response function associated with the infrared vibrational echo is calculated for a quantum mechanical model of resonantly coupled, anharmonic oscillators at zero temperature. The classical mechanical response function is determined from the quantum response function by setting variant Planck's over 2pi-->0, permitting the comparison of the effects of resonant vibrational coupling among an arbitrary number of anharmonic oscillators on quantum and classical vibrational echoes. The quantum response function displays a time dependence that reflects both anharmonicity and resonant coupling, while the classical response function depends on anharmonicity only through a time-independent amplitude, and shows a time dependence controlled only by the resonant coupling. In addition, the classical response function grows without bound in time, a phenomenon associated with the nonlinearity of classical mechanics, and absent in quantum mechanics. This unbounded growth was previously identified in the response function for a system without resonant vibrational energy transfer, and is observed to persist in the presence of resonant coupling among vibrations. Quantitative agreement between classical and quantum response functions is limited to a time scale of duration inversely proportional to the anharmonicity.