This article is devoted to investigating the evolution of the probability density function for power system excited by fractional stochastic noise. First, the single-machine-infinite-bus (SMIB) power system model excited by fractional Gaussian noise (FGN) is established. Second, we derive the Fokker–Planck–Kolmogorov (FPK) equation for the proposed model and solve the FPK equation using the finite difference method. Finally, the numerical results verify that the addition of FGN would influence dynamical stability of the SMIB power system under certain conditions.