This paper presents a new circular curve fitting approach of articulated manipulators, based on least square. The approach aims at gaining the interpolation of circle from n data points, under the condition that the fitted circle should pass both a starting point and an ending point. First a spherical fitting should be performed, using least squares. Then the circular curve fitting can be resulted from the intersection of the fitted sphere and the plane obtained from 3 points, i. e., a starting point, an ending point and the center of a sphere. The proposed algorithms are shown to be efficient by using MATLAB-based simulation.