A solution procedure is proposed to approximate the probability density function (PDF) solution of high-dimensional non-linear systems under Poisson impulses. The PDF solution yields the generalized Fokker-Planck-Kolmogorov (FPK) equation. First a state-space-split method is proposed to reduce the highdimensional generalized FPK equation to a low dimensional equation. After that, the exponential-polynomial closure method is further adopted to solve the reduced FPK equation for the PDF solution. In order to show the effectiveness of the proposed solution procedure, a two-degree-of-freedom coupled pitch-roll ship motion system and a 10-degree-of-freedom mass-spring-damper system are investigated, respectively. Compared to the simulated results, the proposed solution procedure is effective to obtain the PDF solution, especially in the tail region which is very important for reliability analysis.