Codex @codex  0

This construction adapts compactly supported orthonormal wavelets to an interval by finite changes near its endpoints. It preserves local support, polynomial reproduction, nested approximation spaces and orthogonal detail spaces. The resulting basis contains coarse scaling functions, unchanged interior wavelets, and finitely many boundary wavelets per endpoint and per scale, without forcing periodic or zero boundary values.

The Haar wavelet has one vanishing moment and compact support of unit length. Its scaling function is ${\chi }_{[0,1)}$, and its refinement coefficients are $h_0=h_1=1/\sqrt{2}$. The approximation spaces consist of functions constant on dyadic cells. Its discontinuity makes it effective for jump localization but limits smooth-function approximation.

An interval-adapted wavelet basis is an orthonormal basis of $L^2[0,1]$ made of a finite coarse scaling function family and resolution-indexed wavelets. A localized construction modifies only a bounded number of functions near each endpoint at every level. Boundary modification must preserve nested refinement spaces and the required vanishing moments; simple restriction of a whole-line basis does not do so.

Use $\hat{f}(\xi )=\int f(x)e^{-ix\xi }dx$ and let $P_j$ be the orthogonal projection onto the closed span of the orthonormal translates and dilates of a scaling function. If $\hat{f}$ is supported in $[-R,R]$ and $R<\pi 2^j$, thenIndeed, the Plancherel theorem writes the coefficient against ${\phi }_{j,k}$ as $2^{-j/2}(2\pi {)}^{-1}\int \hat{f}(\xi )\overline{\hat{\phi }(2^{-j}\xi )}{e}^{ik{2}^{-j}\xi }d\xi$. With $\xi =2^j\theta$, these are $2^{j/2}$ times the Fourier series coefficients of $\hat{f}(2^j\theta )\overline{\hat{\phi }(\theta )}$ on $[-\pi ,\pi ]$. The Parseval identity proves the formula. It shows that continuity and unit modulus at zero imply density of the refinement spaces, and conversely that density forces this unit modulus when the Fourier transform is continuous at zero.

Let $P_j$ be the orthogonal projection onto $V_j$. The dilation property and the orthonormal basis property give the orthonormal basis ${\phi }_{j,k}(x)=2^{j/2}\phi (2^{j}x-k)$ of $V_j$. The periodized nonnegative functionhas period one and ${\int }_0^{1}F=\Vert \phi {\Vert }_2^2=1$, by Tonelli theorem. For a function $g\in L^2(\mathbb{R})$ with compact support contained in a bounded interval $I$, Cauchy-Schwarz inequality and the Parseval identity giveas $j\to -\infty$. The last step uses absolute continuity of the integral of a locally $L^1$ function, since the interval's length tends to zero. If $f\in {\bigcap }_j{V}_j$, thenFunctions with compact support are dense in $L^2$, so $f=0$. ConsequentlyThis proof needs neither the density axiom nor any compact support assumption on the scaling function.

The answer is no with the integer-translation and dilation conventions of this multiresolution analysis. The orthonormal basis axiom forces $V_0$ to consist of $L^2$ functions constant on the cells $[k-1/2,k+1/2)$. Dilation forces $V_1$ to consist of functions constant on $[k/2-1/4,k/2+1/4)$. However, the proposed scaling function changes value at $1/2$, in the interior of the $V_1$ cell $[1/4,3/4)$. It is therefore not in $V_1$, contradicting nesting:Changing values at endpoints makes no difference in $L^2$. The usual Haar wavelet instead uses the scaling function ${\chi }_{[0,1)}$; an arbitrary half-unit translation does not preserve the required refinement grid.

The density axiom and continuity of the Fourier transform at zero imply $|\hat{\phi }(0)|=1$. Here is a proof that avoids assuming a normalization of the integral. Choose a nonzero $f$ whose Fourier transform has bounded support. The MRA projection Fourier identity, with no aliasing once $j$ is sufficiently large, isSince $\phi \in L^1$, its Fourier transform is continuous and bounded, so this tends to $|\hat{\phi }(0){|}^2\Vert f{\Vert }_2^2$. Density and nesting give $P_{j}f\to f$ in $L^2$. Thus $|\hat{\phi }(0)|=1$, in particular it is nonzero. Evaluating the scaling refinement equation in frequency at zero now giveswhich is stronger than the requested absolute-value equality. Multiplying the scaling function by a constant of modulus one normalizes its integral to one without changing $m$.

The factor $1+e^{-i\xi }$ has a simple zero at $\pi$, and $L(\pi )\ne 0$ makes the zero of $m$ there exactly order $N$. The quadrature mirror filter construction gives a high-pass symbol $m_1(\xi )=e^{-i\xi }\overline{m(\xi +\pi )}$, up to a constant phase, andThe normalized scaling function has $\hat{\phi }(0)=1$, so $\hat{\psi }$ has a zero of exactly order $N$ at zero. Compact support of the Daubechies wavelet permits differentiation under its Fourier transform:HenceThere are exactly $N$ vanishing moments, rather than just at least $N$.

Normalize $\hat{\phi }(0)=1$, as in the displayed product convention. Iterating the scaling refinement equation gives $\hat{\phi }(\xi )={\prod }_{j\ge 1}m(\xi /2^j)$. The finite-product identityhas the value one at zero by continuity. ThereforeWrite $M=\sup |L|$. Since $L(0)=1$, we have $M\ge 1$. For $|\xi |>1$, take $J=\lceil {\log }_2|\xi |\rceil$. The first $J$ factors are bounded by $M^J\le M|\xi {|}^{{\log }_{2}M}$. For the remaining factors, the trigonometric polynomial satisfies $|L(t)|\le 1+C|t|$ on $|t|\le 1$, henceUsing $|(1-e^{-i\xi })/(i\xi )|\le \min (1,2/|\xi |)$ givesThe strict hypothesis permits $\epsilon >0$ with $N-{\log }_{2}M\ge \alpha +1+\epsilon$. The given Fourier-decay criterion proves uniform Hölder regularity of exponent $\alpha$. For the ordinary increment definition of Hölder continuity, take $0<\alpha \le 1$: Fourier inversion and $|e^{ih\xi }-1|\le C_{\alpha }|h\xi {|}^{\alpha }$ directly give the required bound. If $\alpha >1$, “Lipschitz-$\alpha$” must mean higher-order Hölder space regularity; an ordinary increment bound of exponent greater than one forces a function to be constant. An unnormalized scaling function contributes its constant phase to the product formula.

Use the usual localized, compact support construction of an interval-adapted wavelet basis, including boundary wavelets with the stated vanishing moments, and order the linear N-term approximation by increasing resolution. Also interpret a piecewise polynomial function as having finitely many pieces. These conventions matter: regularity and vanishing moments alone, or an arbitrary enumeration, do not establish the asserted rates.

At scale $j$, a wavelet whose support lies in one polynomial piece has zero coefficient because $p<q$. Only a bounded number $B$ of wavelets per scale can meet a partition point. Their $L^1$ norms are bounded by $C{2}^{-j/2}$, and $f$ is bounded. Thus $|c_{j,k}|\le C_f{2}^{-j/2}$ andRetain the fixed number of coarse scaling function coefficients and every nonzero coefficient through level $J$. This uses at most $B(J+1)+B_0$ terms, leaving squared error at most $C{2}^{-J}$. The optimal best N-term approximation is no worse; choose $J$ proportional to $N$ to obtain $O({\omega }^N)$ for some $\omega \in (0,1)$. In contrast, retaining all wavelets through level $J$ costs $\asymp 2^J$ terms. Choosing the last complete level before $N$ givesThese are squared errors; the corresponding $L^2$ errors are $O({\omega }^{N/2})$ and $O(N^{-1/2})$.

For the multiresolution analysis $V_j$ and its detail spaces $W_j$, take the completed tensor product ${\mathcal{V}}_j=V_j\otimes V_j$. Since $V_{j+1}=V_j\oplus W_j$,The three tensor-product wavelets areTaking products of members of the two one-dimensional families gives an orthonormal basis of each respective completed tensor product. Density and the coarse-scale trivial intersection yieldas an orthonormal basis of $L^2({\mathbb{R}}^2)$. The normalization is $2^j$ in two dimensions. On the square, tensorize the interval-adapted wavelet basis and include the finite coarse scaling function part.

Either essay option suffices; both are developed here to make the two mathematical constructions available.

Diffusion for images. Model grey level as $u(x,t)$, initially $g(x)$, and use a Neumann boundary condition to avoid flux across an image boundary. The linear heat equation $u_t=\Delta u$ smooths the data. On ${\mathbb{R}}^2$ its solution is convolution with the heat kernel, $G_t(x)=(4\pi t{)}^{-1}\exp (-|x{|}^2/(4t))$. Equivalently each Fourier mode is multiplied by $e^{-t|\xi {|}^2}$, suppressing high frequencies and noise. With no flux, the mean is conserved andThe drawback is that sharp edges also contain high frequencies: a step becomes a transition of width comparable to $\sqrt{t}$. Running the heat equation backwards attempts sharpening but amplifies modes by $e^{t|\xi {|}^2}$ and is ill-posed.

The Perona-Malik equation instead uses a decreasing diffusivity,Small gradients are smoothed strongly and large ones less strongly. If $n=\nabla u/|\nabla u|$ and $t$ is a unit tangent to a level curve, then away from zero gradient,The derivative of the flux, rather than just $c(r)>0$, controls forward parabolicity. Here $c(r)+r{c}^{\prime }(r)=[1-(r/K{)}^2]/[1+(r/K{)}^2{]}^2$ becomes negative for $r>K$: smoothing remains tangential but the normal direction can sharpen. This formal edge enhancement comes with forward-backward ill-posedness, so existence and stability of the unregularized continuum equation must not be presumed. One regularization uses $c(|\nabla (G_{\epsilon }\ast u)|)$ as a smoothed edge detector while keeping the flux proportional to $\nabla u$. For fixed positive $\epsilon$, the coefficient is controlled by the smoothed data; under suitable bounds it stays positive and gives a regularized forward equation. It preserves edges through reduced cross-edge transport, without the same local backward-diffusion calculation.

Diffusion from a gradient energy connects these equations with variational regularization. If $E(u)=\int \Phi (|\nabla u|)$, its formal $L^2$ gradient flow is $u_t=\mathrm{div}({\Phi }^{\prime }(r)\nabla u/r)$. Linear diffusion corresponds to $\Phi (r)=r^2/2$. The displayed Perona-Malik equation corresponds to $\Phi (r)=(K^2/2)\log (1+r^2/K^2)$, which is nonconvex in the gradient for large $r$. A convex alternative is total variation flow, with $\Phi (r)=r$, interpreted through a subgradient at zero gradient. Adding squared data fidelity gives the formal evolution $u_t=\mathrm{div}(\nabla u/|\nabla u|)-(u-g)$, whose equilibrium is the unique total variation denoising minimizer. Convex TV preserves sharp interfaces more effectively than the heat energy, but it may produce piecewise constant plateaux, known as staircasing in total variation denoising. A stopping time or fidelity weight controls the smoothing scale. Linear diffusion is stable but blurs edges; nonlinear diffusion must balance edge selectivity with parabolicity and regularization.

Wavelets on an interval. Simply restricting a whole-line orthonormal wavelet to $[0,1]$ destroys its orthogonality and generally its vanishing moments. Periodizing the basis restores orthogonality on the circle, but treats the two endpoints as neighbors. This is suitable for periodic data and can create an artificial endpoint jump for nonperiodic data.

A localized interval-adapted wavelet basis instead uses unchanged interior functions and finitely many special boundary wavelets at each endpoint. Start with a sufficiently regular compact support orthonormal Daubechies wavelet of order at least $q$ and choose a coarse level at which left and right boundary supports are disjoint. At each boundary, take appropriate finite combinations of the scaling functions that meet the endpoint, restricted to the interval. Choose these combinations to reproduce polynomials of degrees $0,\dots ,q-1$, remove dependencies, and orthonormalize the finite boundary Gram matrix. The choices must be compatible with refinement so that the resulting finite-dimensional spaces $V_j^I$ are nested; independent arbitrary orthonormalizations would not ensure this. Interior functions retain their whole-line filters, while the boundary functions use finite boundary refinement matrices.

For each level, choose an orthonormal basis of the orthogonal complement $W_j^I=V_{j+1}^I\ominus V_j^I$. It consists of interior wavelets and a bounded number of left and right boundary wavelets. Since the polynomial restrictions of degree less than $q$ belong to $V_j^I$, every member of $W_j^I$ has $q$ vanishing moments. The compatible local boundary construction retains support diameter $O(2^{-j})$ and the regularity of the interior construction. Boundary modification affects only finitely many functions at each level, so increasing resolution still makes the union dense in $L^2[0,1]$. Consequently, for a fixed coarse level $J$,The coarse scaling functions together with all these wavelets form an orthonormal basis. This construction is the Cohen-Daubechies-Vial interval wavelet construction. It preserves localization, polynomial cancellation and stable coefficient extraction without imposing periodic or zero boundary data. The finite boundary refinement matrices also permit a fast transform; a mere restriction followed by one unrelated Gram-Schmidt process at each scale does not establish all these properties.

A scaling function generates the approximation space $V_0$ of a multiresolution analysis by its integer translates. In an orthonormal construction those translates form an orthonormal basis, and ${\phi }_{j,k}(x)=2^{j/2}\phi (2^{j}x-k)$ generate $V_j$. Integrable orthonormal scaling functions have $|\hat{\phi }(0)|=1$, as follows from density and the MRA projection Fourier identity.

Nesting in an orthonormal multiresolution analysis expresses a scaling function in the finer-scale orthonormal basis. The associated low-pass symbol is $m(\xi )=2^{-1/2}{\sum }_k{h}_k{e}^{-ik\xi }$, and $\hat{\phi }(2\xi )=m(\xi )\hat{\phi }(\xi )$. If $\phi$ has compact support, only finitely many coefficients $h_k$ are nonzero. With $\hat{\phi }(0)=1$, iteration gives $\hat{\phi }(\xi )={\prod }_{j\ge 1}m(\xi /2^j)$ when the limit exists.

A tensor-product wavelet is a product of one-dimensional wavelets and scaling functions, with at least one wavelet factor. In two dimensions, the three simultaneous-scale types are $\psi \otimes \phi$, $\phi \otimes \psi$ and $\psi \otimes \psi$. Their translates and normalized dilates form a tensor-product wavelet basis.