We investigate the solutions of refinement equations of the form $$ \phi(x)=\sum_{\alpha\in\mathbb Z^s}a(\alpha)\:\phi(Mx-\alpha), $$ where the function $\phi$ is in $L_p(\mathbb R^s)$$(1\le p\le\infty)$, $a$ is an infinitely supported sequence on $\mathbb Z^s$ called a refinement mask, and $M$ is an $s\times s$ integer matrix such that $\lim_{n\to\infty}M^{-n}=0$. Associated with the mask $a$ and $M$ is a linear operator $Q_{a,M}$ defined on $L_p(\mathbb R^s)$ by $Q_{a,M} \phi_0:=\sum_{\alpha\in\mathbb Z^s}a(\alpha)\phi_0(M\cdot-\alpha)$. Main results of this paper are related to the convergence rates of $(Q_{a,M}^n \phi_0)_{n=1,2,\dots}$ in $L_p(\mathbb R^s)$ with mask $a$ being infinitely supported. It is proved that under some appropriate conditions on the initial function $\phi_0$, $Q_{a,M}^n \phi_0$ converges in $L_p(\mathbb R^s)$ with an exponential rate.

In the present paper, we introduce a modification of the Meyer-K\"{o}nig and Zeller (MKZ) operators which preserve the test functions $f_{0}(x)=1$ and $f_{2}(x)=x^{2}$, and we show that this modification provides a better estimation than the classical MKZ operators on the interval $[\frac{1}{2},1)$ with respect to the modulus of continuity and the Lipschitz class functionals. Furthermore, we present the $r-$th order generalization of our operators and study their approximation properties.