Based on the Kepler's equation, a new set of relative motion equations with elliptical reference orbits were derived. The equations are expressed by orbit elements, and can be used directly for formation flying orbit design. Several typical relative motions and the precision of the equations were analyzed. Furthermore, a linear formation and a perpendicular circles formation were suggested. The effectiveness was verified by the simulation results.