We investigate the solutions of refinement equations of the form

$$\phi (x)\,=\,\sum\limits_{\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 ${{\text{Q}}_{a,M}}$ defined on ${{L}_{p}}({{\mathbb{R}}^{s}})$ by ${{\text{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 ${{(\text{Q}_{a,M}^{n}{{\phi }_{o}})}_{n=1,2,\ldots }}$ 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}}$ , $\text{Q}_{a,M}^{n}{{\phi }_{0}}$ converges in ${{L}_{p}}({{\mathbb{R}}^{s}})$ with an exponential rate.

We consider the shift-invariant space, 𝕊(Φ), generated by a set Φ = {Φ 1,..., Φ r} of compactly supported distributions on R when the vector of distributions ϕ:= {Φ 1,..., Φ r} T satisfies a system of refinement equations expressed in matrix form as

where a is a finitely supported sequence of r x r matrices of complex numbers. Such multiple refinable functions occur naturally in the study of multiple wavelets.

The purpose of the present paper is to characterize the accuracy of Φ, the order of the polynomial space contained in 𝕊(Φ), strictly in terms of the refinement mask a. The accuracy determines the Lp -approximation order of 𝕊(Φ) when the functions in (Φ) belong to Lp(ℝ) (see Jia [10]). The characterization is achieved in terms of the eigenvalues and eigenvectors of the subdivision operator associated with the mask a. In particular, they extend and improve the results of Heil, Strang and Strela [7], and of Plonka [16]. In addition, a counterexample is given to the statement of Strang and Strela [20] that the eigenvalues of the subdivision operator determine the accuracy. The results do not require the linear independence of the shifts of Φ.