In this article, we discuss a class of multiquadric quasi‐interpolation operator that is primarily on the basis of Wu–Schaback's quasi‐interpolation operator scriptLscriptD and radial basis function interpolation. The proposed operator possesses the advantages of linear polynomial reproducing property, interpolation property, and high accuracy. It can be applied to construct flexible function approximation and scattered data fitting from numerical experiments. Copyright © 2013 John Wiley & Sons, Ltd.