The paper deals with mathematical modelling and computer simulation of a gearbox driving system with a double stage gearbox. Mathematical modelling and computer simulations are used for supporting diagnostic inference. Vibration is thought of as a signal of gear condition. It is stressed that vibration generated by gears is influenced by many factors. These factors are divided into four groups: design, production technology, operational, condition change. The condition change of a gearbox is given by gear faults that are divided into single faults such as a tooth crack or breakage or distributed faults as pitting, scuffing, and erosion. The faults are modelled in the case of a crack as a change of tooth stiffness in the case of distributed faults they are given multi-parameter functions. Simulated signals undergo signal analysis by spectrum, cepstrum, time-frequency spectrogram. It has been shown by computer simulation that single and distributed faults are identified by cepstrum. For explicit fault identification time-frequency spectrogram has to be additionally used. The computer simulation results are confirmed by analysis of measured vibration signals received from a gearbox wall/housing. The aim of mathematical modelling and computer simulation, besides finding the relationship between gear condition and vibration signal is in the future to give vibration signals for neural network training.