In this article, we construct a novel self-referential fractal multiquadric function which is symmetric about the origin. The scaling factors are suitably restricted to preserve the differentiability and the convexity of the underlying classical multiquadric function. Based on the translates of a fractal multiquadric function defined on a grid, we propose two fractal quasi-interpolants LCαf\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$L_C^{\\alpha }f$$\\end{document} and LDαf\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$L_D^{\\alpha }f$$\\end{document} to approximate smooth and irregular functions. We study the convergence of LCαf\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$L_C^{\\alpha }f$$\\end{document} and LDαf\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$L_D^{\\alpha }f$$\\end{document} to f using uniform error estimates. We investigate the linear polynomial reproducing property, convexity/concavity and monotonicity features of these quasi-interpolation operators. The advantages of fractal quasi-interpolants over the classical quasi-interpolants are demonstrated by various examples.