The deployment of small cells with carrier aggregation (CA) is a significant feature of fifth generation (5G) mobile communication systems which could be characterized by the multi-dimensional heterogeneity on their diversified requirements upon different resources. Taking the heterogeneity into account, we consider here a joint optimization problem wherein multiple kinds of resources are concurrently allocated to optimize the system throughput utility while enhancing the network energy efficiency (EE) and maintaining the system stability. Especially, for the high-dimensional non-deterministic polynomial (NP)-hard allocation problem embedded, we conduct a mathematical programming model involving nonlinear integer constraints to seek the long-term stable utility on throughput and introduce an iterative optimal modulation and coding scheme-based (optimal MCS-based) heuristic algorithm as an effective solver. In addition, as data traffic and channel condition will be time-varying in the real world, an admission control based on the Lyapunov technique that requires no prior knowledge on channel information is proposed to reduce the system overhead. Finally, not only the performance bound is derived in theory, but also the numerical experiments are conduced to reveal its characteristics with respect to the system parameter V and the EE requirement.