This article studies model reduction of continuous-time stable positive linear systems under the Hankel norm, H-infinity norm and H-2 norm performance. The reduced-order systems preserve the stability as well as the positivity of the original systems. This is achieved by developing new necessary and sufficient conditions of the model reduction performances in which the Lyapunov matrices are decoupled with the system matrices. In this way, the positivity constraints in the reduced-order model can be imposed in a natural way. As the model reduction performances are expressed in linear matrix inequalities with equality constraints, the desired reduced-order positive models can be obtained by using the cone complementarity linearisation iterative algorithm. A numerical example is presented to illustrate the effectiveness of the given methods.