This paper presents a developed 2-D analytical subdomain model for design and analysis of the double-stator permanent magnet linear synchronous machine (DS-PMLSM) to account for both the primary and secondary end effects. First, considering the time constraint, polar transformed method is adopted to establish the 2-D subdomain analytical model, the DS-PMLSM is deformed into a ring-segment PMLSM, and the analytical model is calculated in polar coordinates instead of Cartesian coordinates. Then, the subdomain method is adopted to analyze the analytical model by solving Laplace’s equation and Poisson’s equation, and 2-D relative permeance function is used to consider the slot effect. The flux density and back electromotive force are calculated based on the developed analytical model. Finally, the analytical results are verified by the finite-element method and experiment. The results show that the proposed analytical model can accurately predict the electromagnetic performance, which is useful for the initial design of the DS-PMLSM.