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
The Fourier series expansions are applied to the class of periodic functions defined on finite interval [−ℓ, ℓ]. But not every
function is periodic. Of course, we are free to
extend the function outside the given interval by any definition we want
and no harm is done if we continue the function by periodic repetitions.
However, a function may be given in an infinite range. Fourier had the
ingenious idea to consider such a function as a periodic function whose
period grows to infinity. More precisely, we make ℓ increasingly
larger and larger and investigate what happens in the limit when ℓ
becomes infinite.
We begin our analysis with considering a function f(x) that is defined on finite interval [−ℓ, ℓ]. If this function is absolutely integrable on this interval, its Fourier coefficients are well-defined:
Of course, we have no reason to assume that
an arbitrary function would be reducible to strictly periodic functions.
However, we have seen in the study of the Fourier series that an arbitrary
continuous and piecewise differentiable function of a finite range
could certainly be resolved into an infinite superposition of sine and
cosine functions, whose frequencies were not even arbitrary, but
multiples of one fundamental frequency:
where c✶ denotes the complex conjugate of c (in pure math, it is customary also denoted as \( \displaystyle \quad \overline{c} = \left( a + {\bf j}b \right)^{\ast} = \overline{a + {\bf j}b} = a - {\bf j}b \) ) and VP means the Cauchy principal value, so
This means that the values of the function F{ξ) become
a very dense in ℝ, although n still jumps in integers. The Fourier coefficients
of the standard Fourier series (normalized by setting ℓ - π) can be written in the form
\[
c_n = \frac{1}{2\pi}\,F\left( n \right) ,
\]
which brings into evidence the fact that the Fourier components of
the function f(x) are not some arbitrarily capricious numbers, but
they are connected by the common bond that they belong to a continuous
function F(ξ) that is uniquely tied to the given function f{x). This commonℓ
becomes very large, approaching infinity, because in that case we do
not restrict ourselves to the integer values of ξ, but to all values of aξ.
Now we extend the given function f{x) f from the range [−ℓ, ℓ] to a larger interval {−2ℓ, 2ℓ] by setting f{x) &equive; 0 outside the given range [−ℓ, ℓ]. In this zase, we need to evaluate integrals Fourier coefficients
Example 1:
We start with a linear function f(x) = x on interval [−ℓ, ℓ]. Its Fourier series is
\[
x = \frac{2\ell}{\pi}\,\sum_{n\ge 1} (-1)^{n+1} \frac{1}{n}\,\sin\left( n\frac{\pi}{\ell}x \right) .
\]
So
\[
b_n = (-1)^{n+1} \frac{2\ell}{\pi} \cdot \frac{1}{n} , \qquad \hat{f}_n = (-1)^n {\bf j}\,\frac{\ell}{n\pi} .
\]
Their squares are
\[
b_n^2 = \frac{4 \ell^2}{\pi^2 n^2} , \qquad \left\vert \hat{f}_n \right\vert^2 = = \frac{\ell^2}{n^2 \pi^2} = 4\,b_n^2 , \qquad n=1,2,\ldots .
\]
Integrate[x*Sin[n*Pi*x/a], {x, -a, a}]/a
-((2 a (n \[Pi] Cos[n \[Pi]] - Sin[n \[Pi]]))/(n^2 \[Pi]^2))
Using Parseval's identity, we get
\[
\frac{1}{\ell} \int_{-\ell}^{+\ell} x^2 {\text d} x = \frac{2}{3}\,\ell^2 = \sum_{n\ge 1} \frac{4 \ell^2}{\pi^2 n^2} = \frac{4 \ell^2}{\pi^2} \cdot \frac{\pi^2}{6} .
\]
Sum[1/n^2 , {n, 1, Infinity}]
\[Pi]^2/6
Integrate[x^2 , {x, -a, a}]
(2 a^3)/3
■
End of Example 1
It was Fourier's discovery that a large
class of functions of an infinite domain follow the same pattern,
provided that f(x) remains absolutely integrable in the infinite interval. Let ℓ > 0 and consider a complex‑valued function
This is the standard complex Fourier series on [-ℓ,ℓ].
Embedding a nonperiodic function into this framework involves
\[
f\ :\ \mathbb{R}\rightarrow \mathbb{C}
\]
with sufficient decay and regularity—for instance, assume f ∈ 𝔏¹(ℝ) ∩ ℭ¹(ℝ) and that f(x) → 0 as |x| → ∞. We want to represent f on all of points of ℝ.
Rewrite the series explicitly in a Riemann-sum form:
Interchanging limit and integral is justified by dominated convergence or similar theorems when f has sufficient decay and regularity (e.g., f ∈ đť’®(ℝ)).
Identification of the Fourier transform.
Define the Fourier transform (one of the standard normalizations) by
This is the rigorous, academic transition from Fourier series (discrete spectrum on a finite interval) to the Fourier integral / Fourier transform (continuous spectrum on ℝ), with the Riemann-sum limit ℓ → ∞ as the bridge.
Another derivation
We present another derivation of the Fourier integral from Fourier series of different flavor for pedagogical reasons because education is mostly a repetition. We start with an absolutely integrable function f on interval [−π, π]. This condition, f ∈ 𝔏(−π, π), guarantees existence of Fourier coefficients
Recall that we denote by ⅉ or j the imaginary unit on the complex plane ℂ. This allows us to rewrite the Fourier series in terms of exponential function:
where "V.P." or just "VP" symbolizes the Cauchy principal value series regularization. Therefore, we can assign to every absolutely integrable function the corresponding Fourier series:
Let f(x) be absolutely integrable function on interval [−πℓ, πℓ]. We change variable x = uℓ, where u ∈ [−π, π]. In this case, function g(u) = f(uℓ) = ff(x)i>(x) is defined on interval [−π, π]. Its Fourier series becomes
Stein, E.M. and Shakarchi, R., Fourier Analysis: An Introduction, World Book Publishing Company, 2011. Chapter 2, Section “From Fourier Series to the Fourier Transform”.
Titchmarsh, E.C., Introduction to the Theory of Fourier Integrals, Oxford University Press, 1948.
Zygmund, A., Trigonometric Series, Volumes I&II Combined, Cambridge University Press, 2nd edition, 1988.
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