This paper presents a likelihood-based algorithm for identifying different phase shift keying (PSK) modulations, i.e., BPSK, QPSK, and 8PSK. This algorithm selects the modulation type that maximizes a loglikelihood function that is based on the known original constellation associated with the constellation of the received signals for the candidate modulation types. However, there are two problems in non-cooperative underwater acoustic Multiple Phase Shift Keying (MPSK) modulation identification based on the likelihood method. One is the original constellation, which as prior information is unknown. The other is the underwater acoustic multipath channel makes the constellation distort seriously. In this paper, we solved these problems by combining sparse bayesian learning (SBL) with expectation maximization (EM). The specific steps are as follows. Firstly, blind channel equalization can be achieved by channel impulse response (CIR), which is estimated by sparse bayesian learning in single input multi output (SIMO) underwater acoustic channel. Subsequently, we used expectation maximization to compensate amplitude attenuation and phase offset, as the original constellation of MPSK is unknown. Finally, modulation can be successfully identified by the Quasi Hybrid Likelihood Ratio Test (QHLRT). The simulation results show that the channel estimation method based on SBL can eliminate the influence of channel effectively, and the EM algorithm can make the received constellation converge to the preset constellation in the case of unknown original transmit constellation, which effectively solves these two problems. We use the proposed SBL-EM-QHLRT method to achieve an identification rate of more than 95% in underwater acoustic multipath channels with Signal to Noise Ratio (SNR) higher than 15 dB, which provides a new idea for modulation identification of non-cooperative underwater acoustic MPSK.