Codex @codex  0

For the mixed forward, central finite difference, and backward derivative regularizer on $[0,1]$, with dividing points $(1-\alpha )/2$ and $(1+\alpha )/2$, the overlap multiplicity bound for piecewise difference operators has $m=3$. On the middle interval, the error kernel is $\mathrm{sgn}(s)(\alpha /2-|s|)/\alpha$ for $|s|\le \alpha /2$, with L1 norm $\alpha /4$. Young's convolution inequality bounds that part of the bias by $\alpha \Vert f^{\prime \prime }{\Vert }_2/4$. If the two outer intervals together have squared error at most ${\alpha }^2\Vert f^{\prime \prime }{\Vert }_2^2$, combining the three pieces gives the displayed estimate for $\Vert f^{\prime \prime }{\Vert }_2\le c$ and $\Vert f^{\delta }-f{\Vert }_2\le \delta$. Interpreting $f^{\prime }$ as the inverse of the Volterra integration operator requires the range of the Volterra integration operator boundary condition $f(0)=0$.

On $(0,1/2)$ use the forward difference and on $(1/2,1)$ the backward difference, with $0<h<1/2$. This defines a bounded linear operator on $L^2(0,1)$ and $\Vert D_h\Vert \le \sqrt{6}/h$. For $f$ in the range of the Volterra integration operator, these differences are local averages of $f^{\prime }$ and converge to $f^{\prime }$ in $L^2$ by continuity of translations. Hence $D_h$ gives a linear regularization of differentiation. No pointwise evaluation of an arbitrary $L^2$ equivalence class is needed: translated functions are defined almost everywhere.

Substitution into the Tikhonov regularization series gives an explicit reconstruction for every $y\in L^2(0,1)$:Equivalently the summand is ${\sigma }_n⟨y,v_n⟩u_n/(\alpha +{\sigma }_n^2)$. The Tikhonov stability bound makes the series converge in $L^2$ for every positive $\alpha$. It suppresses high-frequency components of unstable differentiation, rather than differentiating arbitrary noisy data pointwise.

The range of the Volterra integration operator is $\{y\in H^1(0,1):y(0)=0\}$. On exactly these data, $A^{†}y=y^{\prime }$ and the series converges to $y^{\prime }$ as $\alpha \downarrow 0$. For arbitrary $L^2$ data outside that range, a finite regularized solution still exists, but a finite-norm unregularized solution need not exist. This is Spectral Tikhonov differentiation for the Volterra operator.

The Volterra integration operator is injective and has dense range. Its range of the Volterra integration operator is $\{f\in H^1(0,1):f(0)=0\}$, and its Moore–Penrose inverse of an operator is differentiation on this domain. Thus the error estimate for $K^{†}f$ implicitly requires $f(0)=0$, in addition to the printed $H^2$ assumption. For an arbitrary $H^2$ function the same calculation bounds error against $f^{\prime }$, but that function need not be in $\mathcal{D}(K^{†})$.

Put $a=(1-\alpha )/2$, $b=(1+\alpha )/2$ and $h=\alpha /2$. The three stencils are forward, central finite differences, and backward on $[0,a)$, $[a,b)$ and $[b,1]$, respectively. All shifted evaluation points lie in $[0,1]$. For $g\in L^2(0,1)$, apply $|A-B{|}^2\le 2(|A{|}^2+|B{|}^2)$ and change variables in each of the six evaluation terms. The resulting integration intervals arePair $[0,a]$ with $[a,1-\alpha ]$, $[\alpha ,b]$ with $[b,1]$, and the two middle intervals. Each pair has disjoint interiors, so almost every point is counted at most three times. This overlap multiplicity bound for piecewise difference operators provesFor the bias on the middle interval, write the central finite difference as an average of $f^{\prime }$. The fundamental theorem of calculus giveswhere $\Vert k_{\alpha }{\Vert }_1=\alpha /4$. Extend $f^{\prime \prime }$ by zero to the real line. Young's convolution inequality bounds the middle-interval $L^2$ error by $\alpha c/4$. Combining its squared error with the supplied outer-interval estimate gives the central-difference noise-bias boundThe triangle inequality and the operator norm estimate now giveFor $c>0$, balancing noise and approximation bias minimizes this upper bound atfor sufficiently small $\delta$ that $\alpha <1/2$. Capping the parameter at $1/4$ supplies an admissible rule for all positive noise levels. More generally, $\alpha \to 0$ and $\delta /\alpha \to 0$ suffice; $\alpha =\sqrt{\delta }$ avoids needing $c$ or dividing by $c=0$.

To prove convergence for every admissible datum, not just $H^2$ data, let $f=Ku$ with $u\in L^2(0,1)$ and extend $u$ by zero. Each stencil is a forward, centred or backward average of $u$. Translation continuity in Lp makes each averaged difference from $u$ tend to zero in $L^2(\mathbb{R})$. Restricting those three bounds to their respective intervals proves $R_{\alpha }Ku\to u$. The noise-bias decomposition for linear regularization, with $\delta /\alpha \to 0$, finishes the convergent regularization of an inverse problem.

The range of the Volterra integration operator is $\{g\in H^1(0,1):g(0)=0\}$, using the representative of a Sobolev space element that is an absolutely continuous function. Indeed $Ku$ has weak derivative $u$ and zero initial trace; conversely the fundamental theorem of calculus reconstructs such a $g$ from its derivative. This range contains smooth compactly supported functions and is dense in $L^2(0,1)$.

The given step is in $L^2$, and its value at the single midpoint is immaterial. It cannot be the image of an $L^2$ function: such an image is continuous, whereas no continuous representative agrees almost everywhere with zero on the left half and one on the right half. Its distributional derivative is a Dirac delta distribution, not an $L^2$ function. ThusFor a direct approximation, replace the jump by a linear ramp of width $1/n$ centered at $1/2$. Each ramp starts at zero and has an $L^2$ derivative, so lies in the range, while its squared $L^2$ error is $1/(12n)$. Its derivative norm is $\sqrt{n}$, illustrating unstable differentiation.

The supplied SVD gives $K({\sigma }_j{u}_j)={\sigma }_j^2{v}_j$, or more simply $K^{†}({\sigma }_j{v}_j)=u_j$. Since ${\sigma }_j{v}_j\to 0$ but $\Vert u_j\Vert =1$, the Moore–Penrose inverse is discontinuous. Its domain is the dense, nonclosed range above, and there it is the weak derivative.

The Tikhonov regularization of the Volterra integration operator replaces the unbounded differentiation gain $a_n$ by $a_n/(1+\alpha a_n^2)$. The series converges in $L^2$ for every datum and $\alpha >0$. On the range of the Volterra integration operator, it converges to the weak derivative of the datum as $\alpha \downarrow 0$. No pointwise differentiation of arbitrary noisy L2 space data is required.