A method for calculating the pulsation response of reciprocating compressor piping systems using convolution integral is proposed. The method considers the interaction between the valve motion of compressors and the pulsation of piping systems. Compressors are simulated by integrating the deferential equations, and pulsations in piping systems are simultaneously calculated by convolution integrals. Impulse responses of piping systems are evaluated as Fourier transforms of the impedances that are calculated by admittance matrix method including attenuations by mean flow. By this approach, Complex piping systems connected to multiple compressors are easily analyzed. Experiments were carried out with a model piping system that had a reciprocating compressor and a vibrating piston. Impedances of the piping system with mean flows were measured using an accelerometer on the vibrating piston and pressure sensors at the pipes. Predicted impedances and simulated pressures for the piping system were compared with the experimental results and good agreement was obtained.