A comprehensive time scale calculation method for hydrogen maser and cesium atomic clock based on the improved exponential filtering and Vondrak_Cepek joint smoothing method is proposed. It is based on the minimum error method, which dynamically estimates the frequency drift parameters of the hydrogen maser and improves the accuracy of predicting the clock bias of the hydrogen atomic clock; With the method of improved quadratic exponential smoothing, the hydrogen masers time scale is generated, and the weighted average method is used to generate the cesium atomic clocks' time scale. At the same time, the Vondrak_Cepek filter is designed to combine the long-term and short-term stability advantages of the two types of time scales to improve the comprehensive time scale performance. Experimental results show that the proposed method produces hydrogen cesium with a comprehensive time scale stability of up to 1.60 X10^-15h^-1, stable up to 3 x10^-15 d^-1, outperforming the time-scale performance generated by three classic methods: ALGOS, AT1, and Kalman filtering.