In the wake of unceasing development in the post-Newtonian celestial mechanics, at present, the research on the post-Newtonian effects has been exhibited gradually due to the fact that the degree of accuracy of astronomical observation improves unceasingly. Hence several authors devoted their research to this subject and the scopes (Brumberg, (1972, 1985)   , Rubincam (1977)  , Soffel (1987, 1989)   , Iorio (2005)  ). These authors research mainly the post- Newtonian effect of the orbital elements of planets and artificial satellites in the solar system. The largest post-Newtonian effects have been exhibited in the orbits of binary systems, especially in the compact binary systems. So, the research in the post-Newtonian orbital effect of the compact binary stars is important and meaningful. Hence some authors research this subject from theory and observation. Such as, Will (1981, 2006)   , Damour & Deruelle (1985, 1986)   , Schāfer & Wex (1993)  , Wex (1995)  . Calura (1997)  , Iorio (2007)  studied this subject from theory. On the other hand, Burgay et al. (2003)  , Konacki (2003)  , Kramer et al. (2005)  , and Weisberg & Taylor (2005)  researched the post-Newtonian effect for the periastron shift of compact binary stars (PSRJ0737-3039, PSR1913+16 PSR1543+12) from observation. These researches are very interested in the theoretical and observing aspects. This paper presents the post-Newtonian effect on all orbital elements of the compact binary stars on the theoretical aspect. The research method of this paper is different from previous studies that this paper uses the method of perturbation theory in the celestial mechanics.
2. R, S and W Components for the General Relativistic Accelerations in the Two-Body Problem
The relative acceleration of two-body with the Post-Newtonian Parameters is given by Will (1981) 
here , , , , , , ,
here f denotes the true anomaly. is a unit vector in the radial direction and are unit vectors in the orbital plane. is directed along the radial direction, and is perpendicular to . In the equation m denotes Gm and the right side should be multiplied by c−2. G is gravitational constant and c is the speed of light.
In this paper we research the general relativistic effect. In the general relativistic case the Post-Newtonian parameters (Will 1981)  . The Equation (1) can be written
Using the relative expressions the Equation (3) may be written
here boldface denotes vector.
We resolve the acceleration into a radial component , a component , normal to and a component W normal to the orbital plane.
i.e, , (the unit vector normal to the orbital plane).
On comparison with the expression (4), we get three scalar accelerative components R, S and W
Substituting the following formulas of the problem of two body into the above formula (Smart, 1953) 
An independent variable dt is transformed to an independent variable df in the Gaussian perturbation equations (Brouwer & Clemence, 1961  , Roy, 1988)  by using the second formula of the Equations (6), we get the perturbation equations with an independent variable true anomaly f
where is the longitude of periastron, E is the eccentric anomaly and is the mean longitude at epoch.
3. Integration for the Perturbation Equations and Its Perturbation Solutions
Substituting R, S and W for expressions (7) into the perturbation Equation (8) by using Kepler’s third law , we obtain
Integrating the above Equation (9), we obtain the perturbation secular and periodic solutions
where denotes the mean longitude of periastron, E denotes the eccentric anomaly. In the last integral expression, we have used the next integral already:
4. The Secular and Periodic Variation of the Orbital Elements
1) The secular variation per cycle (revolution)
By letting , , , the periodic terms are disappeared and one obtains the secular variables per cycle (revolution):
The author Li (2010)  obtained the formulas for the post-Newtonian effect on the time variation of periastron passage of binary stars. In that paper we change the symbol of the first expression of (35) (in the case of relativity) as the symbol of the present paper, which is
2) The secular variable rate:
where the period P is denoted in unit of day
3) The periodic variation of amplitudes:
In the expressions (10) all terms are the periodic variable terms except for all secular terms. Here we list the maximal and minimal amplitudes of the periodic terms for semi-major axis, a and eccentricity, e from the expressions (10).
For the semi-major axis:
5. Numerical Results for Four Compact Binary Systems
We use the formulae (11) - (13) to calculate the secular of the general relativistic secular effect on the orbital elements of four compact binary systems. It is convenient to reduce the formulas (11) - (13) to practical units m1, m2 and a are denoted by the unit in solar mass , and solar radius , and P is denoted by the unit in day = 86400 s, G = 6.67 × 10−8, (c. g. s), c = 3 × 1010 cm/s. The formulae (11) - (13) can be written by taking the secular effect
here (c, g, s).,
This paper chooses four compact binary systems: PSR1913+16, PSR1543+12, PSRJ0737-3039 and a black hole M33 X-7 as an example. For these compact binary stars, their data for P, a, e, M and m are retrieved from Burgay et al. (2003)  , Konacki et al. (2003)  , Kramer et al. (2005)  , Willems et al. (2004)  and Orosz et al. (2007)  . Their data are listed in Table 1.
Substituting these data in Table 1 into formulas (16) and (13) and (14) - (15), we obtain the numerical results for the periodic and secular variation of the orbital elements of four compact binary stars in Tables 2-4.
Table 1. The data of four compact binary stars.
Table 2. Periodic variation of the amplitudes of the orbits of semi-major axis and eccentricity.
It can be seen from Table 2 that the maximum amplitude of semi-major axis is the black hole binary system M33 X-7 and the maximum eccentricity is PSR 1913 + 16.
It can be seen from Table 3 that the maximum secular variable of longitude per period is PSRJO737-3039. The longest time of the periastron passage is M33-7.
It can be seen from Table 4 that the maximum secular variable rates of longitude is PSRJ0737-3039. The longest time of the periastron passage is M33X-7.
6. Discussion and Conclusions
1) The comparison of the theoretical results with the observable results.
The theoretical results in this paper as compared with the observable results given by several authors for three compact binary stars are listed in the Table 5.
It can be seen from the above Table 5 the theoretical results are very close to the observed results.
Table 3. Secular variation of the orbital elements of four compact binary stars per cycle (Revolution).
[Note] The symbol denotes arc-second: Rev denotes Revolution (cycle).
Table 4. Secular variable rates of the orbital elements of four compact binary stars per year.
Table 5. Comparison of the theoretical results with the observable results.
2) The possibility of observing effects.
In the solar system the advance of perihelion of Mercury may be observed by the recent instrument. We can see from Table 4 that maximal value of the advance of PSRJ0737-3039 is 61016" per year which correspond to 145000" time value of advance of perihelion of Mercury in solar system. The value of the advance of periastron of compact binary star is largest than that of Mercury in solar system. Therefore the effects of the advance of compact binary stars can be observed too.
We have four conclusions:
a) The compact binary stars are the best objects for studying the post-Newto- nian effects on the orbits
b) Although there are no secular variation for the semi-major axis and eccentricity, there is maximal amplitude of the periodic terms for semi-major axis, such as .
c) The longitudes of periastron and the mean longitude at epoch exist both secular and periodic variable terms, and the maximal values for and arrive at 16.95/yr and −29.68/yr for PSRJ0737-3039 respectively.
d) The longest time of periastron passage is 19.47 minute per year for black hole M33 X-7. This corresponds to over 3 seconds per a day.
 Rubincam, D.P. (1977) General Relativity and Satellite Orbit: The Motion of a Test Particle in the Schwarzchild Metric. Celestial Mechanics and Dynamical Astronomy, 15, 21-33.
 Burgay, M., et al. (2003) An Increased Estimate of the Merger Rate of Double Neutron Stars from Observation of a Highly Relativistic System. Nature, 426, 531-533.
 Konacki, M., Wolszczan, A. and Stairs, I.H. (2003) Geodetic Precession and Timing of the Relativistic Binary Pulsars B1543+12 and PSR1913+16. The Astrophysical Journal, 589, 495-502.
 Weisberg, J.M. and Taylor, J.H. (2005) The Relativistic Binary Pulsars B 1913+16. Thirty Years of Observations and Analysis. In: Rasio, F.A. and Stairs, I.H., Eds., ASP Conference Series, 328, 25-31.
 Willems, B., Kalogera, V. and Henninger, M. (2004) Plusar Licks and Spin Tilts in the Close Double Newtron Stars PSR J0737-3039. PSR B 1534+12 and PSR 1913+16. The Astrophysical Journal, 616, 414-438.