This paper presents an analytical subdomain model for the prediction of various electromagnetic parameters of magnetic-geared double-rotor machine (MGDRM). By dividing the MGDRM into six subdomains and solving the Maxwell equations in polar coordinates for each region, the vector magnetic potential distribution can be derived. Subsequently, the flux density distribution, back electromotive force (EMF), and output torque, can be calculated. Furthermore, based on the proposed model, the equivalent d-q axis is elaborated. Also, the corresponding d and q inductance can be obtained. It proves that the MGDRM can be regarded as a non-salient-pole synchronous machine. Besides, the power factor of the MGDRM under i d =0 control is deduced and optimized for the different thicknesses of modulators. In addition, the demagnetization capability is analysed by adopting the subdomain model. Finally, the accuracy of proposed subdomain model for MGDRMs is validated via the finite-element method (FEM). Nomenclature A z Vector magnetic potential. B r , B θ Radial and tangential component of magnetic flux density. M , J Magnetization vector of magnets, Current density vector. μ 0 Magnetic permeability of vacuum. P, Q Number of stator slots and modulator pieces. P i Pole pair Number of inner rotor permanent magnets (PMs). β, γ, δ Modulator opening angle, stator slot width angle and slot opening angle. φ 0 , θ 0 Initial circumferential position of the inner rotor and modulator. θ i , θ j , J l Circumferential positions of PMs, i th modulator opening and j th slot opening, current density in the l th slot