The main purpose of this article is to decompose the response amplitudes of different resonant modes and to further study the wave energy distributions systematically inside long and narrow rectangular harbors with different bottom slopes when harbor oscillations are induced by normal-incident solitary waves. A series of numerical experiments are carried out using the FUNWAVE 2.0 model. The analysis procedure is mainly based on the improved normal mode decomposition method. Results show that when the incident wave height is small, the resonant wave energy inside harbors is dominated by the lowest few modes, and the higher modes only possess a very small proportion of the resonant energy; when the incident wave height increases, the relative energy distribution becomes uniform, and the proportion of energy in the higher modes increases. In addition, for the same incident wave height, the change of the bottom slope inside the harbor has a negligible effect on the relative energy distribution within the ranges of the variation in bottom slopes and the incident wave heights studied in this article.