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
Glossary
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
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 transformCesàro Summation of the Inverse Transform
To regularize the inversion integral formula, one replaces the raw partial integrals SRf by their Cesàro averages:
Start from the definition
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 . \]
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 }. \]
- 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 √π.
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
Properties of the Cesàro kernel
The kernel KR enjoys several crucial properties:-
Nonnegativity:
\[ K_R (t)\geq 0\quad \mathrm{for\ all\ }t. \]This eliminates the oscillatory cancellations that plague the Dirichlet kernel.
-
Normalization:
\[ \int _{-\infty }^{\infty }K_R(x)\, {\text d}x = 1. \]
-
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 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
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
- it does not require pointwise convergence of SRf(x),
- it holds for all f ∈ 𝔏¹,
- it provides a constructive and stable inversion method.
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)| ≤ supy KR(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. \]
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 ∥f∥∞ ≤ M. 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.
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. ■
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 ReconstructionIn 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.
- 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