With the distributed generation technology widely applied, some system problems such as overvoltages and undervoltages are gradually remarkable, which are caused by distributed generations like wind energy system (WES) and photovoltaic system (PVS) because of their probabilistic output power which relied on natural conditions. Since the impacts of WES and PVS are important in the distribution system voltage quality, we study these in this paper using new models with the probability density function of node voltage and the cumulative distribution function of total losses. We apply these models to solve the IEEE33 distribution system to be chosen in IEEE standard database. We compare our method with the Monte Carlo simulation method in three different cases, respectively. In the three cases, these results not only can provide the important reference information for the next stage optimization design, system reliability, and safety analysis but also can reduce amount of calculation.