The characteristics of rainfall variation in water resource divisions in Shandong province are analysed using precipitation data of 41 weather stations in recent 40 years. Mean generation function is used as basic function,which is screened through principal component analysis. Mean-step predictive models are set up in water resource divisions in term of mean generation function. The forecasting results are rather good through testing and fitting calculation.