To address the fact that the traditional curve fitting method based on B-spline basis functions cannot preserve the sharp features in the original data well, a curve fitting method based on a class of orthogonal piecewise polynomial function Vsystem is proposed in this research. Firstly, different types of feature points from the original data are extracted by using the feature extraction algorithm; secondly, the feature points are reparametrized to the locations of different knots in the V-system; finally, the fitting curve is obtained by solving least-squares linear equations with constraints. Different features in the original data can be captured since the V-system contains basis functions with different smoothness. Numerical experimental results show that the proposed method in this research creates a fitted curve reflecting the global shape of the original data and can accurately represent sharp features.