Model of a relay mirror system with source composed of coherent beam array is established. Phase distribution of the source is optimized by using the stochastic parallel gradient descent algorithm. The influence factor induced by the intensity non-uniformity of receiving beam is defined and the evaluation function of uplink propagation process is theoretically analyzed. Under H-V5/7 turbulent model, performances of uplink propagation in the relay mirror systems with coherent beam array source and flat beam source are calculated. Results show that power coupling efficiency of uplink propagation and intensity uniformity of received beam are both effectively improved by using the coherent multi-beams array controlled by the stochastic parallel gradient descent algorithm. The evaluation function of 10 km uplink propagation is improved from 0.6730 to 0.8838 and the evaluation function of 30 km uplink propagation is improved from 0.4266 to 0.8560.