Строится система осредненных уравнений планетного движения вокруг центральной звезды. Используется астроцентрическая система координат. Рассматривается двупланетная задача, однако все построения легко обобщаются на случай произвольного числа N планет. Движение исследуется в модифицированных (комплексных) канонических элементах Пуанкаре. Осреднение выполняется методом Хори-Депри по быстрым средним долготам до второго порядка относительно планетных масс. Разложение возмущающей функции строится с помощью коэффициентов Лапласа. Даны некоторые члены разложения возмущающей функции и первые члены разложения осредненного гамильтониана. Результаты настоящей работы могут быть использованы для исследования эволюции орбит с умеренными значениями эксцентриситетов и наклонов в различных планетных системах.