The rolling element bearing is a key part in many mechanical facilities and the diagnosis of its faults is very important in the field of predictive maintenance. Till date, the resonant demodulation technique (envelope analysis) has been widely exploited in practice. However, much practical diagnostic equipment for carrying out the analysis gives little flexibility to change the analysis parameters for different working conditions, such as variation in rotating speed and different fault types. Because the signals from a flawed bearing have features of non-stationarity, wide frequency range and weak strength, it can be very difficult to choose the best analysis parameters for diagnosis. However, the kurtosis of the vibration signals of a bearing is different from normal to bad condition, and is robust in varying conditions. The fast kurtogram gives rough analysis parameters very efficiently, but filter centre frequency and bandwidth cannot be chosen entirely independently. Genetic algorithms have a strong ability for optimization, but are slow unless initial parameters are close to optimal. Therefore, the authors present a model and algorithm to design the parameters for optimal resonance demodulation using the combination of fast kurtogram for initial estimates, and a genetic algorithm for final optimization. The feasibility and the effectiveness of the proposed method are demonstrated by experiment and give better results than the classical method of arbitrarily choosing a resonance to demodulate. The method gives more flexibility in choosing optimal parameters than the fast kurtogram alone.