1999ApJ...518..442K

Since it has become possible to discover planets orbiting nearby solar-type stars through very precise Doppler-shift measurements, the role of methods used to analyze such observations has grown significantly. The widely employed model-dependent approach based on the least-squares fit of the Keplerian motion to the radial velocity variations can be, as we show, unsatisfactory. Thus, in this paper, we propose a new method that may be easily and successfully applied to the Doppler-shift measurements. This method allows us to analyze the data without assuming any specific model and yet to extract all significant features of the observations. This very simple idea, based on the subsequent subtraction of all harmonic components from the data, can easily be implemented. We show that our method can be used to analyze real 16 Cygni B Doppler-shift observations with a surprising but correct result that is different from that based on the least-squares fit of a Keplerian orbit. Namely, using frequency analysis, we show that with the current accuracy of observations of this star it is not possible to determine precisely the value of the orbital eccentricity, which is claimed to be as high as 0.6.