Refinement equations and refinable functions

[23/100%]
Show
LaTeX

Suppose that φ∈L2(R)\varphi \in L^{2}(\mathbb {R}) is refinable, with refinement coefficients (ck)k∈Z∈ℓ1(Z)\left(c_{k}\right)_{k \in \mathbb {Z}} \in \ell^{1}(\mathbb {Z}). Show that TmφT_{m} \varphi is refinable for every m∈Zm \in \mathbb {Z}, and the refinement coefficients for TmφT_{m} \varphi are (ck+m)k∈Z\left(c_{k+m}\right)_{k \in \mathbb {Z}}.

?

The hat function on [0,2][0,2] is B1(x)=max⁡{1−∣x−1∣,0}B_{1}(x) = \max \left\{ 1-\left|x-1\right|, 0\right\}. Show directly that B1B_{1} is refinable.

?

The Shannon scaling function is the sinc function φ(x)=sin⁡πxπx\varphi (x) = \frac{\sin \pi x}{\pi x}. Show that φ\varphi is refinable, and show that the symbol for the refinement equation is m0=χ[−14,14]m_{0} = \chi_{\left[-\frac{1}{4}, \frac{1}{4}\right]} (extended 1-periodically to R\mathbb {R}).

?

The following refinement equation (based on dilation by 3) has a continuous solution:

φ(x)=12φ(3x)+12φ(3x−1)+φ(3x−2)+12φ(3x−3)+12φ(3x−4). \varphi (x) = \frac{1}{2} \varphi (3 x)+\frac{1}{2} \varphi (3 x-1)+\varphi (3 x-2)+\frac{1}{2} \varphi (3 x-3)+\frac{1}{2} \varphi (3 x-4).

What function satisfies this refinement equation?

?

Justify the factorization performed in equation (12.32) in the proof of Lemma 12.15.

?

Fix φ∈L2(R)\varphi \in L^{2}(\mathbb {R}) and m(ξ)=12∑k∈Zcke−2πikξm(\xi ) = \frac{1}{2} \sum_{k \in \mathbb {Z}} c_{k} e^{-2 \pi i k \xi } with (ck)k∈Z∈ℓ2(Z)\left(c_{k}\right)_{k \in \mathbb {Z}} \in \ell^{2}(\mathbb {Z}), and suppose that

φ^(ξ)=m(ξ/2)φ^(ξ/2) a.e.  \widehat{\varphi }(\xi ) = m(\xi / 2) \widehat{\varphi }(\xi / 2) \text{ a.e. }

Show that if either:

?
(a)

{Tkφ}k∈Z\left\{ T_{k} \varphi \right\}_{k \in \mathbb {Z}} is a Bessel sequence, or

(b)

(ck)k∈Z∈ℓ1(Z)\left(c_{k}\right)_{k \in \mathbb {Z}} \in \ell^{1}(\mathbb {Z}), then φ\varphi is refinable.

Problem 12.20 from Heil C. "A Basis Theory Primer" [Hei11]More problems in Chapter 12 "Wavelet Bases and Frames"

Set B0=χ[0,1]B_{0} = \chi_{[0,1]}, and recursively define the nnth BB-spline BnB_{n} by

Bn=Bn−1∗χ[0,1] B_{n} = B_{n-1} * \chi _{[0,1]}

where * denotes the convolution operation defined in Exercise 12.19.

?
(a)

Show that B1=χ[0,1]∗χ[0,1]B_{1} = \chi_{[0,1]} * \chi_{[0,1]} is the hat function on the interval [0,2][0,2], and find explicit formulas for B2B_{2} and B2′B_{2}^{\prime }.

(b)

Show that BnB_{n} is refinable, integrable, and compactly supported. Show that the symbol for this refinement equation is

m0(ξ)=(1+e−2πiξ2)n+1 m_{0}(\xi ) = \left(\frac{1+e^{-2 \pi i \xi }}{2}\right)^{n+1}
(c)

Find an explicit formula for Bn^\widehat{B_{n}}, and show that Bn^∈L1(R)\widehat{B_{n}} \in L^{1}(\mathbb {R}) for all n>0n > 0.

(d)

Prove that Bn′(x)=Bn−1(x−1)−Bn−1(x)B_{n}^{\prime }(x) = B_{n-1}(x-1)-B_{n-1}(x) for n>1n > 1.

(e)

Show that Bn∈Cn−1(R)B_{n} \in C^{n-1}(\mathbb {R}) for n>0n > 0, and Bn(n−1)B_{n}^{(n-1)} is piecewise linear.

Problem 12.21 from Heil C. "A Basis Theory Primer" [Hei11]More problems in Chapter 12 "Wavelet Bases and Frames"

Suppose that φ∈L2(R)\varphi \in L^{2}(\mathbb {R}) is a compactly supported solution to the finite length refinement equation

φ(x)=∑k=0Nckφ(2x−k) \varphi (x) = \sum _{k = 0}^{N} c_{k} \varphi (2 x-k)

Prove the following statements.

?
(a)

supp⁡(φ)⊆[0,N]\operatorname {supp}\left(\varphi \right) \subseteq [0, N], i.e., φ(x)=0\varphi (x) = 0 for almost every x∉[0,N]x \notin [0, N].

(b)

If φ\varphi is continuous then the vector (φ(1),…,φ(N−1))(\varphi (1), \ldots , \varphi (N-1)) is an eigenvector of the (N−1)×(N−1)(N-1) \times (N-1) matrix

M=[c1c00⋯00c3c2c1⋯00⋮⋮⋮⋱⋮⋮000⋯cNcN−1]=[c2i−j]i,j=1,…,N M = \left[\begin{smallmatrix} c_{1} & c_{0} & 0 & \cdots & 0 & 0 \\ c_{3} & c_{2} & c_{1} & \cdots & 0 & 0 \\ \vdots & \vdots & \vdots & \ddots & \vdots & \vdots \\ 0 & 0 & 0 & \cdots & c_{N} & c_{N-1} \end{smallmatrix}\right] = \left[c_{2 i-j}\right]_{i, j = 1, \ldots , N}

for the eigenvalue 1. Note the double-shift in the rows of MM; for this reason MM is called a two-slanted matrix.

(c)

If φ\varphi is differentiable then φ′\varphi^{\prime } is refinable.

(d)

If φ∈C1(R)\varphi \in C^{1}(\mathbb {R}) then 1 and 1/21 / 2 are both eigenvalues of MM.

(e)

φ\varphi cannot be infinitely differentiable, and in fact can have at most N−2N-2 continuous derivatives.

Suppose φ∈L2(R)\varphi \in L^{2}(\mathbb {R}) is refinable with refinement coefficients (ck)k∈Z∈ℓ2(Z)\left(c_{k}\right)_{k \in \mathbb {Z}} \in \ell^{2}(\mathbb {Z}). Show that if φ\varphi is compactly supported and {Tkφ}k∈Z\left\{ T_{k} \varphi \right\}_{k \in \mathbb {Z}} is orthonormal, then only finitely many of the refinement coefficients are nonzero.

?

Show that equations (12.42) and (12.43) are equivalent when (ck)k∈Z∈ℓ1(Z)\left(c_{k}\right)_{k \in \mathbb {Z}} \in \ell^{1}(\mathbb {Z}).

?

Show that the conclusion of Theorem 12.25 remains valid if we replace the decay hypothesis on the refinement coefficients by ∑k∈Z∣k∣δ∣ck∣<∞\sum_{k \in \mathbb {Z}}\left|k\right|^{\delta }\left|c_{k}\right| < \infty for some δ>0\delta > 0.

?

Suppose that c=(…,0,0,c0,…,cN,0,0,…)c = \left(\ldots , 0,0, c_{0}, \ldots , c_{N}, 0,0, \ldots \right) is a finite sequence with c0,cN≠0c_{0}, c_{N} \neq 0, Show that if the antialiasing requirement (12.42) is satisfied then NN is odd.

?
Problem 12.26 from Heil C. "A Basis Theory Primer" [Hei11]More problems in Chapter 12 "Wavelet Bases and Frames"

Let χ=χ[0,1)\chi = \chi_{[0,1)} be the box function, which satisfies the refinement equation χ(x)=χ(2x)+χ(2x−1)\chi (x) = \chi (2 x)+\chi (2 x-1).

?
(a)

Show that χ^(ξ)=eπiξsin⁡πξπξ\widehat{\chi }(\xi ) = e^{\pi i \xi } \frac{\sin \pi \xi }{\pi \xi }.

(b)

Show that ∣m0(ξ)∣2+∣m0(ξ+12)∣2=1\left|m_{0}(\xi )\right|^{2}+\left|m_{0}\left(\xi +\frac{1}{2}\right)\right|^{2} = 1.

(c)

Give a direct proof of Viète's formula:

∏j=1∞cos⁡(2−jπξ)=sin⁡πξπξ \prod _{j = 1}^{\infty } \cos \left(2^{-j} \pi \xi \right) = \frac{\sin \pi \xi }{\pi \xi }
(d)

Explain how Viète's formula relates to Theorem 12.29.

Remark Viète proved the special case ξ=1/2\xi = 1 / 2, which leads to the formula

2π=222+222+2+22⋯ . \frac{2}{\pi } = \frac{\sqrt{2}}{2} \frac{\sqrt{2+\sqrt{2}}}{2} \frac{\sqrt{2+\sqrt{2+\sqrt{2}}}}{2} \cdots .

Show that Theorem 12.29 continues to hold if instead of the antialiasing condition we only require that ∣m0(ξ)∣2+∣m0(ξ+12)∣2≤1\left|m_{0}(\xi )\right|^{2}+\left|m_{0}\left(\xi +\frac{1}{2}\right)\right|^{2} \leq 1 for ξ∈T\xi \in \mathbf{T}.

?
?
(a)

Show that if ∑k∈Z∣kjck∣<∞\sum_{k \in \mathbb {Z}}\left|k^{j} c_{k}\right| < \infty, then m0m_{0} is jj-times differentiable and its jj th derivative is the continuous function

m0(j)(ξ)=(−2πi)j2∑k∈Zkjcke−2πikx m_{0}^{(j)}(\xi ) = \frac{(-2 \pi i)^{j}}{2} \sum _{k \in \mathbb {Z}} k^{j} c_{k} e^{-2 \pi i k x}

Conclude that

m0(j)(1/2)=0⟺∑k∈Z(−1)kkjck=0 m_{0}^{(j)}(1 / 2) = 0 \quad \Longleftrightarrow \quad \sum _{k \in \mathbb {Z}}(-1)^{k} k^{j} c_{k} = 0

and compare this to the accuracy conditions given in equation (12.46).

(b)

Let BnB_{n} be the nnth BB-spline defined in Exercise 12.20. Show that the accuracy condition m0(j)(1/2)=0m_{0}^{(j)}(1 / 2) = 0 holds for j=0,…,nj = 0, \ldots , n.

Suppose that φ∈L1(R)\varphi \in L^{1}(\mathbb {R}) is refinable, with refinement coefficients (ck)k∈Z∈ℓ1(Z)\left(c_{k}\right)_{k \in \mathbb {Z}} \in \ell^{1}(\mathbb {Z}) that satisfy ∑k∈Zck=2\sum_{k \in \mathbb {Z}} c_{k} = 2. Show that if the "minimal accuracy" condition ∑k(−1)kck=0\sum_{k}(-1)^{k} c_{k} = 0 holds then the periodization of φ\varphi is equal almost everywhere to a constant, specifically

∑j∈Zφ(x+j)=φ^(0) a.e.  \sum _{j \in \mathbb {Z}} \varphi (x+j) = \widehat{\varphi }(0) \quad \text{ a.e. }

(This is easier to do if φ\varphi is assumed to be continuous.)

?

This exercise extends Exercise 12.29, but to avoid issues of convergence we assume that φ\varphi is a compactly supported solution of a refinement equation with finitely many nonzero coefficients ckc_{k}. Suppose that ∑k∈Zck=2\sum_{k \in \mathbb {Z}} c_{k} = 2 and the accuracy conditions ∑k∈Z(−1)kkjck=0\sum_{k \in \mathbb {Z}}(-1)^{k} k^{j} c_{k} = 0 hold for j=0j = 0 and j=1j = 1. Show that

∑j∈Z(j−2a)φ(x+j)=x a.e.  \sum _{j \in \mathbb {Z}}(j-2 a) \varphi (x+j) = x \quad \text{ a.e. }

where

a=∑k∈Z2kc2k=∑k∈Z(2k+1)c2k+1. a = \sum _{k \in \mathbb {Z}} 2 k c_{2 k} = \sum _{k \in \mathbb {Z}}(2 k+1) c_{2 k+1}.
?

The Hilbert transform of f∈L2(R)f \in L^{2}(\mathbb {R}) is the function HfH f whose Fourier transform is (Hf)∧(ξ)=−isign⁡(ξ)f^(ξ)(H f)^{\wedge }(\xi ) = -i \operatorname {sign}\left(\xi \right) \widehat{f}(\xi ), where sign⁡(ξ)=1,0\operatorname {sign}\left(\xi \right) = 1,0, or -1 according to whether ξ>0,ξ=0\xi > 0, \xi = 0, or ξ<0\xi < 0, respectively.

?
(a)

Show that the Hilbert transform is a unitary mapping of L2(R)L^{2}(\mathbb {R}) onto itself, but it does not map L1(R)∩L2(R)L^{1}(\mathbb {R}) \cap L^{2}(\mathbb {R}) onto itself. In particular, if f∈L1(R)∩L2(R)f \in L^{1}(\mathbb {R}) \cap L^{2}(\mathbb {R}) satisfies f^(0)≠0\widehat{f}(0) \neq 0, then Hf∉L1(R)H f \notin L^{1}(\mathbb {R}).

(b)

Suppose that φ∈L2(R)\varphi \in L^{2}(\mathbb {R}) satisfies a refinement equation with coefficients (ck)k∈Z∈ℓ1(Z)\left(c_{k}\right)_{k \in \mathbb {Z}} \in \ell^{1}(\mathbb {Z}). Show that HφH \varphi satisfies the same refinement equation.

Show that if the hypotheses of Theorem 12.25 are satisfied, then the function PP has at most polynomial growth at infinity, i.e., there exists some C>0C > 0 and integer M>0M > 0 such that

∣P(ξ)∣≤C∣ξ∣M,∣ξ∣≥1 \left|P(\xi )\right| \leq C\left|\xi \right|^{M}, \quad \left|\xi \right| \geq 1
?

Assume (ck)k∈Z∈ℓ1(Z)\left(c_{k}\right)_{k \in \mathbb {Z}} \in \ell^{1}(\mathbb {Z}) satisfies ∣m0(0)∣=12∣∑ck∣<1\left|m_{0}(0)\right| = \frac{1}{2}\left|\sum c_{k}\right| < 1. Show that there are no nontrivial integrable solutions to the corresponding refinement equation.

Note: Exercises 12.34-12.36 require some knowledge of distributions, tempered distributions, and the distributional Fourier transform.

?
Problem 12.34 from Heil C. "A Basis Theory Primer" [Hei11]More problems in Chapter 12 "Wavelet Bases and Frames"

The dilation of a distribution μ∈D′(R)\mu \in \mathcal{D}^{\prime }(\mathbb {R}) is defined by duality:

⟨f,D2μ⟩=⟨D1/2f,μ⟩,f∈Cc∞(R) \left\langle f, D_{2} \mu \right\rangle = \left\langle D_{1 / 2} f, \mu \right\rangle , \quad f \in C_{c}^{\infty }(\mathbb {R})

We say that a distribution μ\mu satisfies the refinement equation φ(x)=2φ(2x)\varphi (x) = 2 \varphi (2 x) in a distributional sense if μ=21/2D2μ\mu = 2^{1 / 2} D_{2} \mu.

?
(a)

Show that the function 1/x1 / x satisfies the one-term refinement equation φ(x)=2φ(2x)\varphi (x) = 2 \varphi (2 x). However, 1/x1 / x is not an integrable function and does not determine a tempered distribution.

(b)

The δ\delta distribution is the linear functional on Cc∞(R)C_{c}^{\infty }(\mathbb {R}) defined by ⟨f,δ⟩=f(0)\left\langle f, \delta \right\rangle = f(0) for f∈Cc∞(R)f \in C_{c}^{\infty }(\mathbb {R}). Show that δ\delta satisfies the equation φ(x)=2φ(2x)\varphi (x) = 2 \varphi (2 x) in a distributional sense.

(c)

The jj th distributional derivative δ(j)\delta^{(j)} of δ\delta is defined by ⟨f,δ(j)⟩=(−1)jf(j)(0)\left\langle f, \delta^{(j)} \right\rangle = (-1)^{j} f^{(j)}(0) for f∈Cc∞(R)f \in C_{c}^{\infty }(\mathbb {R}). Show that δ(j)\delta^{(j)} satisfies φ(x)=2j+1φ(2x)\varphi (x) = 2^{j+1} \varphi (2 x) in a distributional sense.

(d)

The principal value of 1/x1 / x is the linear functional pv⁡(1/x)\operatorname {pv}(1 / x) on Cc∞(R)C_{c}^{\infty }(\mathbb {R}) given by

⟨f,pv⁡(1/x)⟩=lim⁡T→∞∫1T<∣t∣<Tf(t)tdt,f∈Cc∞(R) \left\langle f, \operatorname {pv}(1 / x) \right\rangle = \lim _{T \rightarrow \infty } \int _{\frac{1}{T} < \left|t\right| < T} \frac{f(t)}{t} d t, \quad f \in C_{c}^{\infty }(\mathbb {R})

Show that pv⁡(1/x)\operatorname {pv}(1 / x) satisfies φ(x)=2φ(2x)\varphi (x) = 2 \varphi (2 x) in a distributional sense.

Remark δ\delta and pv(1/x)\mathrm{pv}(1 / x) are each tempered distributions, and therefore have distributional Fourier transforms. The Fourier transform of δ\delta is identified with the constant function, i.e., δ^=1\widehat{\delta } = 1. The Fourier transform of pv(1/x)\mathrm{pv}(1 / x) is identified with the function

(pv⁡(1/x))∧(ξ)=−πisign⁡(ξ), (\operatorname {pv}(1 / x))^{\wedge }(\xi ) = -\pi i \operatorname {sign}\left(\xi \right),

where sign⁡(ξ)=1,0\operatorname {sign}\left(\xi \right) = 1,0, or -1 according to whether ξ>0,ξ=0\xi > 0, \xi = 0, or ξ<0\xi < 0, respectively. Up to scale, δ\delta is the unique distributional solution to φ(x)=2φ(2x)\varphi (x) = 2 \varphi (2 x) whose Fourier transform is continuous.

Problem 12.35 from Heil C. "A Basis Theory Primer" [Hei11]More problems in Chapter 12 "Wavelet Bases and Frames"

A tempered distribution μ∈S′(R)\mu \in \mathcal{S}^{\prime }(\mathbb {R}) satisfies the refinement equation in the sense of distributions if we have μ=∑2−1/2ckD2Tkμ\mu = \sum 2^{-1 / 2} c_{k} D_{2} T_{k} \mu, where

⟨f,D2Tkμ⟩=⟨T−kD1/2f,μ⟩,f∈S(R). \left\langle f, D_{2} T_{k} \mu \right\rangle = \left\langle T_{-k} D_{1 / 2} f, \mu \right\rangle , \quad f \in \mathcal{S}(\mathbb {R}).
?
(a)

Show that if the hypotheses of Theorem 12.25 are satisfied, then μ=P∨\mu = \stackrel{\vee }{P} (the inverse distributional Fourier transform of PP) is a tempered distribution that satisfies the refinement equation in the sense of distributions.

(b)

Show that, up to scale, μ\mu is the unique tempered distribution that satisfies the refinement equation and has a Fourier transform that is a continuous function.

(c)

Let θ\theta be any bounded function such that θ(2ξ)=θ(ξ)\theta (2 \xi ) = \theta (\xi ) for all ξ\xi. Show that ν^(ξ)=P(ξ)θ(ξ)\widehat{\nu }(\xi ) = P(\xi ) \theta (\xi ) defines a tempered distribution that satisfies the refinement equation in the sense of distributions.

(d)

Show that a function θ\theta that satisfies θ(2ξ)=θ(ξ)\theta (2 \xi ) = \theta (\xi ) and is continuous at ξ=0\xi = 0 is a constant function. Give examples of nonconstant functions θ\theta such that θ(2ξ)=θ(ξ)\theta (2 \xi ) = \theta (\xi ) for all ξ\xi.

Problem 12.36 from Heil C. "A Basis Theory Primer" [Hei11]More problems in Chapter 12 "Wavelet Bases and Frames"

This exercise considers tempered distributions that are solutions to the refinement equation

φ(x)=∑k=0Nckφ(2x−k), where ∑k=0Nck=2.(12.47) \varphi (x) = \sum _{k = 0}^{N} c_{k} \varphi (2 x-k), \quad \text{ where } \sum _{k = 0}^{N} c_{k} = 2. \tag {12.47}
?
(a)

Let δk\delta_{k} be the point mass at kk distribution, i.e., ⟨f,δk⟩=f(k)\left\langle f, \delta_{k} \right\rangle = f(k) for f∈S(R)f \in \mathcal{S}(\mathbb {R}). Set ν=12∑k=0Nckδk\nu = \frac{1}{2} \sum_{k = 0}^{N} c_{k} \delta_{k}, and show that the distributional Fourier transform of ν\nu is ν^=m0\widehat{\nu } = m_{0}.

(b)

Set

νn=12∑k=0Nckδ2−nk and μn=ν1∗⋯∗νn. \nu _{n} = \frac{1}{2} \sum _{k = 0}^{N} c_{k} \delta _{2^{-n} k} \quad \text{ and } \quad \mu _{n} = \nu _{1} * \cdots * \nu _{n}.

Show that

μn^(ξ)=Pn(ξ)=∏j=1nm0(2−jξ) \widehat{\mu _{n}}(\xi ) = P_{n}(\xi ) = \prod _{j = 1}^{n} m_{0}\left(2^{-j} \xi \right)

and

supp⁡(μn)⊆[0,12N]+⋯+[0,2−nN]⊆[0,N] \operatorname {supp}\left(\mu _{n}\right) \subseteq \left[0, \frac{1}{2} N\right]+\cdots +\left[0,2^{-n} N\right] \subseteq [0, N]
(c)

By Exercise 12.35(a), there is a distributional solution μ\mu to the refinement equation that satisfies μ^=P\widehat{\mu } = P. By Theorem 12.25, μn^=Pn\widehat{\mu_{n}} = P_{n} converges to P=μ^P = \widehat{\mu } uniformly on compact sets. Show that μn→w∗φ\mu_{n} \xrightarrow {\mathrm{w}^{*}} \varphi, i.e.,

lim⁡n→∞⟨f,μn⟩=⟨f,φ⟩,f∈S(R), \lim _{n \rightarrow \infty } \left\langle f, \mu _{n} \right\rangle = \left\langle f, \varphi \right\rangle , \quad f \in \mathcal{S}(\mathbb {R}),

and use this to show that supp⁡(φ)⊆[0,N]\operatorname {supp}\left(\varphi \right) \subseteq [0, N].

(d)

Show that if there exists an integrable solution φ\varphi to the refinement equation (12.47), then φ\varphi is compactly supported and supp⁡(φ)⊆[0,N]\operatorname {supp}\left(\varphi \right) \subseteq [0, N].

(e)

Show that if the antialiasing condition is also satisfied, then there exists a compactly supported solution φ∈L2(R)\varphi \in L^{2}(\mathbb {R}) to the refinement equation (which is therefore integrable by Cauchy-Bunyakovski-Schwarz).

References

[Hei11] Heil C. A Basis Theory Primer. Birkhäuser, 2011.