The authors study the stability of steady state of the compressible Navier-Stokes-Poisson equations with non-flat doping profile. The authors first show that there exists a unique global classical solution to the compressible Navier-Stokes-Poisson equations near the steady state for general doping profile. When the doping profile is small in some norm, the algebraic decay rates are obtained provided the initial perturbation belongs to \(L^p\) with \(p \in [1, 3/2)\).