This study aims to make a numerical simulation program to calculate the vibration level. The level, which is used to evaluate pollution vibration, must be measured by a vibration level meter that satisfies JIS (JIS C 1510) requirements. The level meters are generally made of two circuits to simulate human body responses for frequency and duration. The response of the vertical direction for frequency is different from its horizontal direction. We used the z-transformation method for numerical simulation of those characteristics, and obtained the linear difference equations from the pulse transfer function. On the other hand, the response for duration was simulated by a non-linear circuit which was a CR circuit with a diode. The result was similarly expressed by the linear difference equations. Using the difference equations. we developed a program which could simulate those response characteristics.