We consider bounded linear operators from the view of functional reproducing property. For some bounded linear operators associated with orthogonal polynomials we define an inner product space associated with a kernel constructed with orthogonal polynomials, show that it is a functional reproducing kernel Hilbert space (FRKHS) associated with these bounded linear operators and give decay rate for the best FRKHS approximation with a [Formula: see text]-functional associated with the FRKHS. On this basis, we provide a learning rate for kernel regularized regression whose hypothesis space is the defined FRKHS. As applications, we define some concrete FRKHSs associated with polynomial operators such as the Bernstein–Durrmeyer operators, the de la Vallée Poussin operators on both the unit sphere [Formula: see text] and the unit ball [Formula: see text]. We show that these polynomial operators have reproducing property with respect to the corresponding concrete FRKHSs and show the learning rate for the kernel regularized regression. In short, we provide a way of constructing FRKHS operators with Fourier multipliers and show a learning framework from the view of operator approximation.