The formulation of approximate solutions to equations that embody the dominant characteristics of the orbital motion of a two-satellite tethered system are studied. The orbital motion of the system is viewed as perturbed two-body motion, and a restricted tether problem is obtained by neglecting librational motion. An exact analytical solution to this restricted problem in terms of elliptic functions is presented. An approximate solution to the restricted tether problem obtained by applying the method of averaging is also provided. An approximation for small-amplitude librational motion is formulated, whose solution is based on methods for solving equations with variable coef cients. The analytical solutions are good approximations to the orbital motion of the tetherperturbed satellite and the librational motion of the system when the libration is small. The restricted tether motion approximation is then utilized to solve the identi cation problem of a tethered satellite system. Nomenclature a ¤ = apparent semimajor axis e ¤ = apparent eccentricity F = incomplete elliptic integral of the rst kind h = nondimensionalangular momentum K = complete elliptic integral of the rst kind K 0 = associate complete elliptic integral of the rst kind k = modulus of Jacobian elliptic functions and integral k 0 = complementary modulus of Jacobian elliptic functions and integral m; m p = masses of satellites q = Jacobi's nome r = nondimensionalradial distance, radial distance/r E r E = radius of the Earth, 6,378,000 m sn = Jacobian elliptic function t = time t ¤ = nondimensionaltime, t p .¹=r 3 E / D t £ 0:0012394 u = 1=r X = state vector ®;¯;°; » = variable coef cients "= tether parameter, m p ½=.m C m p /r E µ = orbital angle (true anomaly) µ 2 = out-of-plane libration angle µ 3 = in-plane libration angle ¹ = gravitational constant of the Earth ½ = tether length '; ' ¤ = amplitudes ! = orbital frequency