Nonlinear energy sinks (NESs) are critically important for structural vibration suppression. They can absorb vibrational energy across a broad frequency spectrum, possess strong robustness, and have a relatively small mass. This study addresses the vibration suppression in a piecewise linear stiffness NES system under random excitation. Initially, a theoretical model of the piecewise linear stiffness NES system is developed. The piecewise linear stiffness function is approximated using Legendre polynomial approximation. Following this, the steady-state Fokker–Planck–Kolmogorov (FPK) equation of the system is formulated via the Generalized Harmonic Function Method. The FPK equation is solved using the fourth-order central difference method, and the effectiveness of this approach is validated by comparing the FDM results with numerical simulations. Lastly, the influence of varying system parameters on the stability of the piecewise linear stiffness NES system is analyzed.