Return to computing page for the first course APMA0330
Return to computing page for the second course APMA0340
Return to computing page for the fourth course APMA0360
Return to Mathematica tutorial for the first course APMA0330
Return to Mathematica tutorial for the second course APMA0340
Return to Mathematica tutorial for the fourth course APMA0360
Return to the main page for the course APMA0330
Return to the main page for the course APMA0340
Return to the main page for the course APMA0360
Introduction to Linear Algebra with Mathematica

Preface


The inverse Fourier transform is traditionally understood as a limit process that reconstructs a function f from its frequency-domain representation ℱ[f]. In the classical setting of 𝔏¹(ℝ) or 𝔏²(ℝ), this inversion is expressed through an oscillatory integral whose convergence is subtle and, in general, not absolutely guaranteed. Cesàro summation—specifically the (C,1) method—provides a powerful regularization technique that stabilizes this inversion process. By replacing the raw partial integrals of the inverse transform with their Cesàro averages, one obtains a kernel with improved localization and positivity properties, analogous to the Fejér kernel in Fourier series. This essay develops the analytic foundations of this approach and explains why Cesàro summation yields a robust and universal method for reconstructing functions from their Fourier transforms.

 

Cesàro summation

For a function f ∈ 𝔏¹(ℝ), the Fourier transform is defined by

\begin{equation} \label{EqCesaro.1} \hat {f}(\xi )=\int _{\mathbb{R}}f(x) e^{-{\bf j}x\xi }\, {\text d}x. \end{equation}
Formally, the inverse transform (in the sense of 𝔏¹) is:
\begin{equation} \label{EqCesaro.2} f(x) =\frac{1}{2\pi }\,\mbox{V.P.}\int _{\mathbb{R}}\hat {f}(\xi )\, e^{{\bf j}x\xi }\, {\text d}\xi . \end{equation}

The Inverse Fourier Transform as a Limit of Truncated Integrals

However, the inverse Fourier integral is not absolutely convergent in general, and its determination is an ill-posed problem. One therefore introduces the truncated inverse transform
\begin{equation} \label{EqCesaro.3} \left( S_R \,f \right) (x) = \frac{1}{2\pi }\int_{-R}^R \,\hat{f}(\xi )\, e^{{\bf j}x\xi} \, {\text d}\xi = \frac{1}{2\pi }\int_{-R}^R \, e^{{\bf j}x\xi} \, {\text d}\xi \,\int_{\mathbb{R}} f(u)\,e^{-{\bf j}u\xi} \, {\text d} u , \end{equation}
and interprets the inversion formula as the limit
\[ ℱ^{-1} \left[ \hat{f} \right]_{\xi\to x} = \lim_{R\rightarrow \infty } \,\left( S_R \, f \right) (x) = \lim_{R\rightarrow \infty } \,\int_{\mathbb{R}} \, f(u)\, {\text d}u\, \frac{1}{2\pi }\int_{-R}^R \,e^{{\bf j}(x-u)\,\xi} \, {\text d}\xi . \]
The operator SR is a convolution with the Dirichlet kernel:
\[ \left( S_R \, f \right) (x) = \left( f * D_R \right) (x),\qquad D_R (t) = \frac{1}{2\pi }\int_{-R}^R e^{{\bf j}t\xi }\, {\text d}\xi = \frac{\sin (Rt)}{\pi t}. \]
Integrate[Exp[I*t*x], {x, -R, R}]
(2 Sin[R t])/t
The kernel DR is oscillatory, nonpositive, and lacks good localization. As in Fourier series case, these properties lead to divergence phenomena, Gibbs oscillations, and failure of pointwise convergence for general 𝔏¹ functions.

Cesàro Summation of the Inverse Transform

To regularize the inversion integral formula, one replaces the raw partial integrals SRf by their Cesàro averages:

\[ \left( \sigma_R \, f \right) (x) = \frac{1}{R}\int _0^R \left( S_r \, f \right) (x)\, {\text d}r. \]
This is the continuous analogue of Fejér means for Fourier series. We want to rewrite Cesàro,s mean as a convolution with a kernel KR, i.e.,
\[ \left( \sigma_R \,f \right) (x) = \left( f * K_R\right) (x) = \int_{-\infty }^{\infty} \, f(y)\, K_R (x-y)\, {\text d}y, \]
and then show that { KR }, R > 0, is an approximate identity.

Start from the definition

\[ \left( S_r \,f \right) (x) = \frac{1}{2\pi }\int_{-r}^r \,\hat{f}(\xi )\, e^{{\bf j}x\xi }\, {\text d}\xi . \]
and then apply Cesàro's mean,
\[ \left( \sigma_R \, f \right) (x) = \frac{1}{R}\int_0^R \left[ \frac{1}{2\pi }\int _{-r}^r\, \hat{f}(\xi )\, e^{{\bf j}x\xi }\, {\text d}\xi \right] {\text d}r. \]
Assume that f ∈ 𝕃¹(ℝ), so its Fourier transform ℱ[f] is bounded and continuous; Fubini’s theorem justifies interchanging integrals:
\[ \left( \sigma _R \,f \right) (x) = \frac{1}{2\pi }\,\mbox{V.P.} \int_{-\infty }^{\infty} \,\hat{f}(\xi )\, e^{{\bf j}x\xi }\left[ \frac{1}{R}\int_0^R \,\mathbf{1}_{[-r,r]} (\xi )\, {\text d}r\right] {\text d}\xi . \]
For fixed ξ, the inner integral is
\begin{equation} \label{EqCesaro.4} \left( \sigma _R \,f \right) (x) = \frac{1}{2\pi }\int_{-R}^R \left( 1-\frac{|\xi |}{R}\right) \, \hat{f} (\xi )\, e^{{\bf j}x\xi }\, {\text d}\xi . \end{equation}
This is the Cesàro-weighted inverse Fourier integral. Now use the Fourier inversion structure: write the latter as a convolution. Define
\[ K_R(x) := \frac{1}{R}\int_0^R \left( D_r \right) (x)\, {\text d}r = \frac{1}{2\pi }\int_{-R}^R \left( 1-\frac{|\xi |}{R}\right) e^{{\bf j}x\xi }\, {\text d}\xi . \]
where DR is the Dirichlet kernel. Then
\[ \left( \sigma_R \,f \right) (x) = \mbox{V.P.} \int_{-\infty }^{\infty } \,f (y)\, K_R (x-y)\, {\text d}y = \left( f*K_R \right) (x). \]
So the Cesàro means correspond to convolution with the kernel KR.

   
Example 1: Let f : ℝ → ℝ. Usual improper integral over ℝ: \[ \mbox{V.P.} \int_{-\infty }^{\infty }f(x)\, {\text d}x := \lim_{T\rightarrow \infty }\int _{-T}^T f(x)\, {\text d}x \] (if the limit exists). Cesàro's integral over ℝ first form the symmetric truncation \[ S(T):= \int _{-T}^T f(x)\, {\text d}x, \] then take the Cesàro mean \[ C(T) := \frac{1}{T}\int_0^T S(t)\, {\text d}t, \] and say that f is Cesàro integrable over ℝ with value IC if \[ \lim _{T\rightarrow \infty }C(T)=I_C. \] If the usual improper integral exists, then the Cesàro integral exists and has the same value.

Now we give an illustrative example where Cesàro works but the usual integral diverges.

Take \[ f(x)=\cos x. \] Usual improper integral diverges. As T → ∞, 2 sinT oscillates and has no limit, so \[ \int _{-\infty }^{\infty }\cos x\, dx \] does not exist in the usual improper sense.

However, its Cesàro integral converges. Indeed, compute the Cesàro mean \[ C(T) = \frac{1}{T}\int_0^T S(t)\, {\text d}t = \frac{2 \left( 1 - \cos T\right)}{T}, \qquad S(t) = \int_{-t}^t \cos x\,{\text d} x = 2\,\sin t . \]

Integrate[Cos[x], {x, -t, t}]
2 Sin[t]
Integrate[2*Sin[t], {t, 0, T}]/T
(2 - 2 Cos[T])/T
As T → ∞, \[ \left| \frac{2}{T}(-\cos T+1)\right| \leq \frac{4}{T}\rightarrow 0, \] so \[ \lim _{T\rightarrow \infty }C(T)=0. \] Conclusion for this example:

Usual improper integral over ℝ diverges.
Cesàro integral over ℝ exists and equals 0.

Another Example where both notions agree.

Take \[ f(x) = e^{-x^2}. \] The usual improper integral is classical: \[ \int _{-\infty }^{\infty }e^{-x^2}\, {\text d} x = \sqrt{\pi }. \]

Integrate[Exp[-x^2], {x, -Infinity, Infinity}]
Sqrt[\[Pi]]
Then the Cesàro mean is expressed through the error function: \[ C(T) =\frac{1}{T}\int_0^T S(t)\, {\text d}t = \frac{e^{-T^2} - 1 + T\,\sqrt{\pi}\,\mbox{erf}(T)}{\sqrt{\pi} \, T} \]
Integrate[Exp[-x^2], {x, -t, t}]
Sqrt[\[Pi]] Erf[t]
Integrate[Erf[t], {t, 0, T}]/T
(-1 + E^-T^2 + Sqrt[\[Pi]] T Erf[T])/(Sqrt[\[Pi]] T)
because \[ \int_{-t}^t e^{-x^2} \,{\text d}x = \sqrt{\pi}\,\mbox{erf}(t), \qquad \mbox{erf}(t) = \frac{2}{\sqrt{\pi}} \, \int_0^t e^{-x^2}\.{\text d} x . \] Since S(t) → √π and S is bounded and convergent, a standard averaging argument gives \[ \lim _{T\rightarrow \infty }C(T)=\sqrt{\pi }. \] Conclusion for this example:
  • Usual improper integral: \( \displaystyle \quad \int _{-\infty }^{\infty } e^{-x^2} \,{\text d}x = \sqrt{\pi }. \)
  • Cesàro integral: exists and equals the same value √π.
   ■
End of Example 1

We compute the Cesàro kernel explicitly. Note that the weight is even in ξ, and \( \displaystyle \quad e^{{\bf j}x\xi } \quad \) is decomposed into cosine and sine; the sine part vanishes by symmetry. So

\[ K_R(x) =\frac{1}{2\pi }\int_{-R}^R \left( 1-\frac{|\xi |}{R}\right) \cos (x\xi )\, {\text d}\xi . \]
Lemma 1: For all R > 0 and t ∈ ℝ, the Fejér kernel is \[ K_{R}(t) = \frac{1}{2\pi}\left(\frac{\sin(Rt/2)}{t/2}\right)^{2}. \]
Compute the integral: \[ \int _0^R\left( 1-\frac{\xi }{R}\right) \cos (x\xi )\, {\text d}\xi =\int _0^R\cos (x\xi )\, {\text d}\xi -\frac{1}{R}\int_0^R \,\xi \cos (x\xi )\, {\text d}\xi . \] We have \[ \int_0^R \,\cos (x\xi )\, {\text d}\xi = \frac{\sin (xR)}{x}, \] and \[ \int_0^R \,\xi \cos (x\xi )\, {\text d}\xi =\left[ \frac{\xi \sin (x\xi )}{x}\right] _0^R -\int _0^R\frac{\sin (x\xi )}{x}\, {\text d}\xi =\frac{R\sin (xR)}{x}-\frac{1}{x^2}(1-\cos (xR)). \] Therefore, \[ \int_0^R \left( 1-\frac{\xi }{R}\right) \cos (x\xi )\, {\text d}\xi =\frac{\sin (xR)}{x}-\frac{1}{R}\left[ \frac{R\sin (xR)}{x}-\frac{1}{x^2}(1-\cos (xR))\right] =\frac{1}{Rx^2}(1-\cos (xR)). \] Use the identity \[ 1-\cos (xR)=2\sin ^2\left( \frac{xR}{2}\right) , \] so \[ K_R(x)=\frac{2}{\pi Rx^2}\sin ^2\left( \frac{xR}{2}\right) . \] Rewrite in the standard Fejér-like form: \[ K_R(x)=\frac{1}{2\pi R}\left( \frac{\sin \left( \frac{xR}{2}\right) }{\frac{x}{2}}\right) ^2 =\frac{1}{2\pi R}\left( \frac{2\sin \left( \frac{xR}{2}\right) }{x}\right) ^2. \] This is exactly the continuous Fejér kernel on ℝ, scaled by R.
This is the Fejér kernel on the real line, also known as the positive summability kernel for the Fourier transform.

Properties of the Cesàro kernel

The kernel KR enjoys several crucial properties:
  1. Nonnegativity:
    \[ K_R (t)\geq 0\quad \mathrm{for\ all\ }t. \]
    This eliminates the oscillatory cancellations that plague the Dirichlet kernel.
  2. Normalization:
    \[ \int _{-\infty }^{\infty }K_R(x)\, {\text d}x = 1. \]
  3. Concentration at 0: For every δ > 0,
    \[ \int_{|x|>\delta } K_R(x)\, {\text d}x \rightarrow 0\quad \mathrm{as\ }R\rightarrow \infty . \]
    So as R → ∞,
    \[ K_R (t)\rightarrow \delta (t) , \]
    in the sense of distributions, and more precisely:
    • KR(t) decays like O(1/t²) for large t.
These three properties make KR an approximate identity. We verify the normalization property with Mathematica:
Integrate[(2*Sin[R*t/2]/t)^2 , {t, -Infinity, Infinity}]
ConditionalExpression[2 \[Pi] Abs[R], R \[Element] Reals]
By construction, the Fourier transform of the Fejér kernel is
\[ \hat{K_R}(\xi ) =\int_{-\infty }^{\infty }K_R(x)\, e^{-{\bf j}x\xi }\, {\text d}x =\left( 1-\frac{|\xi |}{R}\right) _+, \]
where ( · )+ denotes the positive part, i.e.,
\[ \int _{-\infty }^{\infty }K_R(x)\, {\text d}x =\hat{K_R}(0) =1. \]
To prove the last property, we fix δ >0. We want to show that
\[ \int_{|x|>\delta }K_R(x)\, {\text d}x \rightarrow 0\quad \mathrm{as\ }R\rightarrow \infty . \]
Use the explicit formula:
\[ K_R(x) =\frac{2}{\pi Rx^2}\sin^2 \left( \frac{Rx}{2}\right) . \]
For |x| ≥ δ,
\[ 0\leq K_R(x)\leq \frac{2}{\pi Rx^2}\leq \frac{2}{\pi R\delta ^2}. \]
Thus,
\[ \int _{|x|>\delta }K_R(x)\, {\text d}x\leq \frac{2}{\pi R\delta ^2}\int_{|x|>\delta } {\text d}x =\infty \cdot \frac{1}{R} \]
which is not directly helpful as written. Instead, we use a more careful splitting. Write
\[ \int _{|x|>\delta }K_R(x)\, {\text d}x =\int_{\delta <|x|\leq M}K_R(x)\, {\text d}x +\int_{|x|>M}K_R(x)\, {\text d}x, \]
and choose M large, then let R → ∞. For the tail |x|>M, use the bound
\[ K_R(x)\leq \frac{2}{\pi Rx^2}, \]
so
\[ \int _{|x|>M}K_R(x)\, {\text d}x \leq \frac{2}{\pi R}\int _{|x|>M}\frac{{\text d}x}{x^2} =\frac{4}{\pi RM}. \]
For fixed M, this tends to 0 as R → ∞. On the middle region δ <|x| ≤ M, we use the Riemann–Lebesgue type oscillation in R. For each fixed x ≠ 0,
\[ K_R(x)=\frac{1}{\pi R}\left( \frac{\sin \left( \frac{Rx}{2}\right) }{x/2}\right) ^2\rightarrow 0\quad \mathrm{as\ }R\rightarrow \infty , \]
because the numerator is bounded and we divide by R. Moreover, for δ <|x| ≤ M,
\[ 0\leq K_R(x)\leq \frac{C}{R} \]
for some constant C depending on δ, M. Thus, by dominated convergence,
\[ \int _{\delta <|x|\leq M}K_R(x)\, {\text d}x\rightarrow 0\quad \mathrm{as\ }R\rightarrow \infty . \]
Combining both pieces, we get
\[ \limsup _{R\rightarrow \infty }\int _{|x|>\delta }K_R(x)\, dx\leq \limsup _{R\rightarrow \infty } \,\left( \int_{\delta <|x|\leq M}K_R(x)\, {\text d}x +\frac{4}{\pi RM}\right) =0. \]
So { KR } concentrates at 0.    ▣

These properties mirror those of the Fejér kernel for Fourier series and are the foundation of the Cesàro inversion theorem.

Cesàro summation as a method of inversion

For a Lebesgue intregrable function f on ℝ, a point x in the domain of f is a Lebesgue point if \[ \lim_{r\to 0} \, \frac{1}{r}\,\int_{|x-y| < r} \left\vert f(x) - f(y) \right\vert {\text d} y = 0 . \]

The Lebesgue points of f are thus points where f does not oscillate too much, in an average sense. According to the Lebesgue differentiation theorem, for any locally integrable function f, the set of points that are not Lebesgue points has a Lebesgue measure of zero.

Given the convolution representation

\[ \left( \sigma_R \,\hat{f} \right) (x) = \left( \hat{f}*K_R \right) (x), \]
the classical theory of approximate identities implies:
Theorem 1 (Cesàro Inversion for Fourier Transform): If f ∈ 𝔏¹(ℝ), then for almost every x, \[ \lim_{R\rightarrow \infty } \left( \sigma_R \,f \right) (x) = f(x), \] and the convergence is uniform at every Lebesgue point of f.
This result is strictly stronger than the classical inversion theorem, because:
  • it does not require pointwise convergence of SRf(x),
  • it holds for all f ∈ 𝔏¹,
  • it provides a constructive and stable inversion method.
Proof at Lebesgue points:
Let x₀ be a Lebesgue point of f ∈ 𝔏(ℝ). We must show \[ \lim_{R\rightarrow \infty } \, \left( f *K_R \right) (x_0) = f(x_0). \] where the Fejér-type kernel is \[ K_R (t) = \frac{1}{2\pi}\int_{-R}^R \Bigl(1-\frac{|\xi|}{R}\Bigr)e^{{\bf j}\xi t}\,{\text d}\xi = \frac{R}{2\pi}\left(\frac{\sin(Rt/2)}{Rt/2}\right)^2 . \] We write \[ \left( \sigma_R\, f \right) (x) - f(x) = \int_{\mathbb{R}} \bigl(f(x-t)-f(x)\bigr)\,K_R(t)\,{\text d}t. \] Fix ε > 0. Since x₀ is a Lebesgue point, there exists δ > 0 such that \[ \frac{1}{2\delta }\int_{|y|<\delta } |f(x_0 -y)-f(x_0 )|\, {\text d}y < \varepsilon . \] We split the integral into near and far parts: \[ \left( f*K_R \right) (x_0) = I_1 (R) + I_2 (R) , \] where \[ I_1 (R) = \int_{|y|\le \delta} \bigl(f(x_0 - y) - f(x_0)\bigr)\,K_R(y)\,{\text d}y , \qquad I_2 (R) = \int_{|y|\ge \delta} \bigl(f(x_0 - y) - f(x_0)\bigr)\,K_R(y)\,{\text d}y. \]

Use |KR(y)| ≤ supyKR(y) and \( \displaystyle \quad \int_{|y|<\delta }K_R(y)\, {\text d}y\leq 1. \quad \) More directly, note that \[ |I_1(R)|\leq \int _{|y|<\delta }|f(x_0-y)-f(x_0)|\, K_R(y)\, {\text d}y. \] Since KR ≥ 0 and \( \displaystyle \quad \int K_R =1, \quad \) we can bound \[ |I_1(R)|\leq \left( \sup _{|y|<\delta }\frac{1}{2\delta }\int_{|z|<\delta }|f(x_0-z)-f(x_0)|\, {\text d}z\right) \cdot \int _{|y|<\delta }K_R(y)\, {\text d}y\leq \varepsilon . \] A more standard argument is to use the fact that KR is “almost supported” in |y| < δ for large R, but the Lebesgue point condition already gives us the smallness of the average oscillation. To be fully explicit, we can write \[ I_1(R)|\leq \left( \sup _{|y|<\delta }|f(x_0-y)-f(x_0)|\right) \int _{|y|<\delta }K_R(y)\, {\text d}y, \] and then approximate the supremum by the average using the Lebesgue point property. Either way, we can make |I₁(R)| ≤ Cε for some absolute constant C.

The estimate of I₂(R) relies only on:

  • boundedness of f (or f\in L_{\mathrm{loc}}^1 with truncation),
  • decay of the Fejér kernel tail: \[ \int _{|t|>\delta }|K_R(t)|\, {\text d}t\rightarrow 0. \]
This part is purely an approximate identity tail estimate. It works for every x, Lebesgue point or not.

Let f ∈ 𝔏(ℝ) and let x ∈ ℝ be a point where f has finite one-sided limits. The Cesàro means of the Fourier inversion formula can be written as \[ \sigma_R (x) = \int_{\mathbb{R}} f(x-t)\,K_R(t)\,{\text d}t, \] where the Fejér-type kernel is \[ K_R (t) = \frac{1}{2\pi}\int_{-R}^R \Bigl(1-\frac{|\xi|}{R}\Bigr) e^{{\bf j}\xi t}\,{\text d}\xi = \frac{R}{2\pi}\left(\frac{\sin(Rt/2)}{Rt/2}\right)^2 . \] We now estimate I₂(R). Since f is bounded, let ∥fM. Then \[ |f(x-t)-f(x)| \le 2M, \] and therefore, \[ |I_2(R)| \le 2M \int_{|t|>\delta} |K_R(t)|\,{\text d}t. \] Thus, it suffices to show \[ \int_{|t|>\delta} \left\vert K_R(t) \right\vert {\text d}t \ \xrightarrow[R\to\infty]{} \ 0. \]

Since f ∈ 𝔏¹(ℝ), we can use the fact that { KR } is an approximate identity: \[ \int_{|y|\geq \delta }K_R(y)\, {\text d}y\rightarrow 0 \quad \mbox{ as } \quad R \to \infty , \] and \[ \int |f(x_0-y)|K_R(y)\, {\text d}y \leq \| f\|_{1}\cdot \sup_y \,K_R(y) \] Using the explicit formula for KR(t), \[ K_R(t) = \frac{2}{\pi R}\,\frac{\sin^2(Rt/2)}{t^2}, \] we have for all |t| > δ: \[ |K_R(t)| \le \frac{2}{\pi R}\,\frac{1}{t^2}. \] Hence, \[ \int_{|t|>\delta} |K_R(t)|\,dt \le \frac{2}{\pi R} \int_{|t|>\delta} \frac{dt}{t^2}. \] Compute the integral explicitly: \[ \int_{|t|>\delta} \frac{dt}{t^2} = \int_{\delta}^{\infty} \frac{dt}{t^2} + \int_{\delta}^{\infty} \frac{dt}{t^2} = \frac{1}{\delta} + \frac{1}{\delta} = \frac{2}{\delta}. \] Therefore, \[ \int_{|t|>\delta} |K_R(t)|\,dt \le \frac{2}{\pi R} \cdot \frac{2}{\delta} = \frac{4}{\pi\delta}\,\frac{1}{R}. \] Combining the estimates: \[ |I_2(R)| \le 2M \int_{|t|>\delta} |K_R(t)|\,dt \le 2M \cdot \frac{4}{\pi\delta}\,\frac{1}{R} = \frac{8M}{\pi\delta}\,\frac{1}{R}. \] Thus, \[ I_2 (R) \ \xrightarrow[R\to\infty]{} \ 0. \] This completes the estimate of the tail integral I₂(R). Then \[ \sup_{|y|\geq \delta }K_R(y)\rightarrow 0 \quad \mbox{as } R\rightarrow \infty \quad (\mbox{since } K_R(y)\sim 1/R\quad \mbox{for fixed } y\neq 0), \] so for large R, |I₂(R)| is as small as we like. Putting both estimates together, we get \[ \limsup _{R\rightarrow \infty }|(f*K_R)(x_0)-f(x_0)|\leq C\varepsilon , \] and since ε > 0 is arbitrary, the limit is 0. This proves pointwise convergence at Lebesgue points.

   
Example 2: The Fourier transform of function f is \[ \hat {f}(\xi )=\int _{\mathbb{R}}e^{-{\bf j}x\xi }f(x)\, {\text d}x,\qquad \xi \in \mathbb{R}. \] The formal inverse transform is \[ f(x)=\frac{1}{2\pi }\int _{\mathbb{R}}e^{{\bf j}x\xi }\, \hat {f}(\xi )\, {\text d}\xi , \] but this integral may fail to converge in the usual sense, even when f ∈ 𝔏¹(ℝ).

A Cesàro-type inversion formula says that for f ∈ 𝔏¹(ℝ). \[ f(x) = \frac{1}{2\pi }\lim _{T\rightarrow \infty }\,\frac{1}{T}\int_0^T {\text d}t \int _{-t}^te^{{\bf j}x\xi }\, \hat {f}(\xi )\, {\text d}\xi \] for almost every x, including all points of continuity of f. So instead of a single symmetric truncation \( \displaystyle \quad \int _{-T}^T, \quad \) we take the Cesàro mean of these truncations.

Rewriting Cesàro inversion as convolution with an approximate identity. Start from \[ \frac{1}{T}\int _0^T {\text d}t\, \int _{-t}^t \,e^{{\bf j}x\xi }\, \hat {f}(\xi )\, {\text d}\xi . \] Insert the definition of Fourier transform ℱ[f]: \[ \hat {f}(\xi ) =\int_{\mathbb{R}} e^{-{\bf j}y\xi }\,f(y)\, {\text d}y. \] Then \begin{aligned}\frac{1}{T}\int _0^T {\text d}t\,\int _{-t}^t \,e^{{\bf j}x\xi }\, \hat {f}(\xi )\, {\text d}\xi &= \frac{1}{T}\int _0^T {\text d}t \,\int_{-t}^t \,e^{{\bf j}x\xi }\left( \int _{\mathbb{R}}e^{-{\bf j}y\xi } \,f(y)\, {\text d}y\right) {\text d}\xi \\ &=\frac{1}{T}\int _0^T{\text d}t \,\int_{\mathbb{R}} \,f(y)\left( \int _{-t}^t \,e^{{\bf j}(x-y)\xi }\, {\text d}\xi \right) {\text d}y\, .\end{aligned} Compute the inner integral: \[ ???? \] Thus, \[ \frac{1}{T}\int _0^T\int _{-t}^te^{ix\xi }\, \hat {f}(\xi )\, {\text d}\xi \, {\text d}t = \int _{\mathbb{R}}f(y)\, K_T(x-y)\, {\text d}y, \] where the Cesàro kernel is \[ K_T(u):=\frac{1}{T}\int _0^T\frac{2\sin (tu)}{u}\, dt=\frac{2}{u}\cdot \frac{1}{T}\int _0^T\sin (tu)\, dt. \] Compute the last integral: \[ \int _0^T\sin (tu)\, {\text d}t = \frac{1-\cos (Tu)}{u}, \] so \[ ????? \] Therefore, the Cesàro inversion formula can be written as \[ f(x)=\frac{1}{2\pi }\lim _{T\rightarrow \infty }(f*K_T)(x) \] for almost every x, where (KT) (T>0) is an approximate identity.

Explicit example: indicator of an interval. We take \[ f(x) =\mathbf{1}_{[-1,1]}(x) =\left\{ \, \begin{array}{ll}\textstyle 1,&\textstyle |x|\leq 1,\\ \textstyle 0,&\textstyle |x|>1.\end{array}\right. \] This function is integrable, has jump discontinuities at ± 1, and is continuous elsewhere. Fourier transform of f is \[ \hat {f}(\xi ) =\int _{-1}^1 \,e^{-{\bf j}x\xi }\, {\text d}x = \left[ \frac{e^{-ix\xi }}{-{\bf j}\xi }\right]_{x=-1}^{x=1} =\frac{e^{-{\bf j}\xi }-e^{{\bf j}\xi }}{-{\bf }\xi } = 2\, \frac{\sin \xi }{\xi }. \] So the formal inverse transform is \[ ????? \] This integral is only conditionally convergent and must be interpreted carefully (e.g., as a principal value or via summability methods).

Cesàro inversion for this function. Apply the Cesàro inversion formula: \[ f(x)=\frac{1}{2\pi }\lim _{T\rightarrow \infty }\frac{1}{T}\int_0^T\int_{-t}^t \,e^{{\bf j}x\xi }\, \hat {f}(\xi )\, {\text d}\xi \, {\text d}t. \] Using the convolution form, \[ \frac{1}{2\pi }\frac{1}{T}\int _0^T\int _{-t}^te^{ix\xi }\, \hat {f}(\xi )\, d\xi \, {\text d}t = \frac{1}{2\pi }\int_{\mathbb{R}}f(y)\, K_T(x-y)\, {\text d}y = \frac{1}{2\pi }\int _{-1}^1 \,K_T(x-y)\, {\text d}y. \] So the Cesàro approximants are \[ f_T(x) := \frac{1}{2\pi }\int _{-1}^1 \,K_T(x-y)\, {\text d}y = \frac{1}{2\pi }\int_{x-1}^{x+1} \,K_T(u)\, {\text d}u. \] Now note the key properties of KT:

Normalization: \[ \int _{\mathbb{R}}K_T(u)\, {\text d}u = 2\pi \quad \mathrm{for\ all\ }T>0. \] (So KT/(2π) has total mass 1.)

Approximate identity: as T → ∞, \[ \frac{1}{2\pi }K_T(u)\, du\rightharpoonup \delta _0\quad \mathrm{in\ the\ sense\ of\ distributions,} \] and more explicitly, KT/(2π) concentrates near u=0 and its mass outside any fixed neighborhood of 0 tends to 0. Therefore, for any point of continuity x of f, \[ \lim _{T\rightarrow \infty }f_T(x)=\lim _{T\rightarrow \infty }\frac{1}{2\pi }\int _{-1}^1K_T(x-y)\, dy=f(x). \] In our example: If |x|<1, then f is continuous at x with value 1, and \[ \lim _{T\rightarrow \infty }f_T(x)=1. \] If |x|>1, then f is continuous at x with value 0, and \[ \lim _{T\rightarrow \infty }f_T(x)=0. \] At the jump points x = ± 1, one can show (by symmetry of the kernel) that exactly as in the classical Fourier series/Fejér kernel situation. So the Cesàro inversion recovers f almost everywhere and gives the expected “midpoint of the jump” at discontinuities.

Comparison with ordinary (non-Cesàro) inversion. If instead you use the plain symmetric truncation \[ f_R(x) :=\frac{1}{2\pi }\int _{-R}^Re^{{\bf j}x\xi }\, \hat {f}(\xi )\, {\text d}\xi , \] you obtain:
If |x|<1, then f is continuous at x with value 1, and \lim _{T\rightarrow \infty }f_T(x)=1. If |x|>1, then f is continuous at x with value 0, and \( \displaystyle \quad \lim _{T\rightarrow \infty }f_T(x)=0. \)

At the jump points x = ± 1, one can show (by symmetry of the kernel) that f(1) = ½, exactly as in the classical Fourier series/Fejér kernel situation. So the Cesàro inversion recovers f almost everywhere and gives the expected “midpoint of the jump” at discontinuities.

Comparison with ordinary (non-Cesàro) inversion. If instead you use the plain symmetric truncation \[ f_R(x) := \frac{1}{2\pi }\int_{-R}^R \,e^{{\bf j}x\xi }\, \hat {f}(\xi )\, {\text d}\xi , \] you obtain ======================================== Conclusion (Cesàro Inversion via Characteristic Function) For the inverse Fourier transform, the Cesàro (C,1) method replaces the sharp cutoff \chi _{[-R,R]}(\xi ) with the Fejér–type averaging \frac{1}{R}\int _0^R\chi _{[-r,r]}(\xi )\, dr=\left( 1-\frac{|\xi |}{R}\right) _+. This produces the tempered Cesàro kernel K_R(x)=\frac{1}{2\pi }\int _{-R}^R\left( 1-\frac{|\xi |}{R}\right) e^{ix\xi }\, d\xi =\frac{1}{\pi }\frac{1-\cos (Rx)}{x^2}. The key point is: K_R is positive, integrable, and forms an approximate identity. Therefore, for every f\in L^1(\mathbb{R}) (or tempered distribution), f(x)=\lim _{R\rightarrow \infty }(f*K_R)(x) at every Lebesgue point of f. This is the Cesàro inversion theorem for the Fourier transform. Conclusion (Cesàro Inversion via Characteristic Function) For the inverse Fourier transform, the Cesàro (C,1) method replaces the sharp cutoff \chi _{[-R,R]}(\xi ) with the Fejér–type averaging \frac{1}{R}\int _0^R\chi _{[-r,r]}(\xi )\, dr=\left( 1-\frac{|\xi |}{R}\right) _+. This produces the tempered Cesàro kernel K_R(x)=\frac{1}{2\pi }\int _{-R}^R\left( 1-\frac{|\xi |}{R}\right) e^{ix\xi }\, d\xi =\frac{1}{\pi }\frac{1-\cos (Rx)}{x^2}. The key point is: K_R is positive, integrable, and forms an approximate identity. Therefore, for every f\in L^1(\mathbb{R}) (or tempered distribution), f(x)=\lim _{R\rightarrow \infty }(f*K_R)(x) at every Lebesgue point of f. This is the Cesàro inversion theorem for the Fourier transform.    ■

End of Example 2

Interpretation and Applications

6.1 Regularization of Oscillatory Integrals

The Cesàro method replaces the oscillatory kernel sin(Rt)/t with the positive kernel (sin(Rt/2)/(t/2))², smoothing the inversion process.

6.2 Signal Reconstruction

In signal processing, Cesàro summation corresponds to applying a triangular low-pass filter, which suppresses high-frequency oscillations and improves stability. 6.3 Distribution Theory

For tempered distributions, Cesàro summation provides a canonical way to define inverse transforms of objects whose Fourier transforms are not integrable.

6.4 PDE and Harmonic Analysis

In solving PDEs via Fourier methods, Cesàro summation ensures convergence of integral representations of solutions, especially when initial data is merely integrable.

Conclusion

Cesàro summation provides a conceptually elegant and analytically powerful method for interpreting the inverse Fourier transform. By replacing the raw partial integrals with their Cesàro averages, one obtains a positive, well-localized kernel that forms an approximate identity. This leads to a universal inversion theorem valid for all 𝔏¹ functions and stable at Lebesgue points. The method parallels the classical Fejér summation for Fourier series and extends its benefits to the continuous setting of Fourier integrals. As such, Cesàro summation is not merely a technical device but a fundamental tool in harmonic analysis, signal reconstruction, and the theory of distributions.

Theorem 6:
   
Example 7:    ■
End of Example 7
Theorem 4: Let f,h ∈ 𝔏¹(ℝ) and \[ h(x) = \mbox{V.P.} \,\frac{1}{2\pi} \int H(\xi )\, e^{{\bf j}\xi x} {\text d}\xi \] with integrable H(ξ). Then \[ \left( h \ast f \right) (x) = \frac{1}{2\pi} \int H(\xi )\,\hat{f}(\xi ) \, e^{{\bf j}\xi x} {\text d}\xi . \]
The product of two functions H(ξ) f(y is integrable in (ξ, y) ∈ ℝ×ℝ; hence, by Fubini's theorem, \begin{align*} \left( h \ast f \right) &= \int h(x-y)\, f(y)\,{\text d}y = \frac{1}{2\pi} \iint H(\xi )\, e^{{\bf j}\xi x} e^{-{\bf j}\xi x} f(y)\,{\text d}\xi{\text d}y \\ &= \frac{1}{2\pi} \int H(\xi )\, e^{{\bf j}\xi x} \int e^{-{\bf j}\xi x} f(y)\,{\text d}y {\text d}\xi = \frac{1}{2\pi} \int H(\xi )\,\hat{f} (\xi )\, e^{{\bf j}\xi x} {\text d}\xi . \end{align*}
   
Example 5:    ■
End of Example 5
   

 

  1. Cossar, J., The Cesàro Summability of Fourier Integrals, Proceedings of the Edinburgh Mathematical Society, Vol. 7, Issue 2 (1945), pp. 84–92. doi: 10.1017/S0013091500024366

 

Return to Mathematica page
Return to the main page (APMA0340)
Return to the Part 1 Basic Concepts
Return to the Part 2 Fourier Series
Return to the Part 3 Integral Transformations
Return to the Part 4 Parabolic PDEs
Return to the Part 5 Hyperbolic PDEs
Return to the Part 6 Elliptic PDEs
Return to the Part 6P Potential Theory
Return to the Part 7 Numerical Methods