Wave dark matter, ψDM, has recently been gaining attention as a dark matter candidate. This is due in part to the model’s ability to predict, thus a falsifiable model. Noticeably it predicts a central prominent core, dubbed as the soliton, in every galaxy and large-amplitude density fluctuations in galaxy halos dubbed as halo granules. In this work we conduct the first full ψDM halo simulations including stars for elliptical galaxies. The results show that the solitonic core is enhanced by about one order of magnitude in the peak density due to the star gravity, thereby increasing its contrast to the halo. Furthermore, the solitonic core grows by absorbing mass from its host halo so as to retain the soliton scaling relation. On the other hand, halo granules become “non-isothermal” after stars have thoroughly interacted with the inner halo, a situation unlike isothermal granules in halos of pure ψDM cosmological simulations. We find that the composite (dark matter+stars) mass density can locally follow the isothermal profile near the stellar half-light radius in most cases. Phase separation that segregates a uniform population of stars into two distinct populations of high and very low velocity dispersions stars is also found in some situations. While halo stars have high velocity dispersion, the cold component is bound inside the soliton resembling stars residing in faint dwarf spheroidal galaxies. The most striking result of this study is that the velocity dispersion of halo stars increases rapidly toward the galactic center by a factor of at least 2 or so inside the half-light radius caused by the deepened soliton potential, a result that compares favorably with observations of elliptical galaxies and bulges in spiral galaxies. Finally, we find that the density within the inner halo (a few ten times of its core radius) follows the same self-similar relation as the soliton does within the mass range of dwarf to small galaxies.