Fourier Transform
Fourier Transform
Fourier Transform
The Fourier transform is a fundamental tool, in particular for the study of partial
derivative equations, because it “diagonalizes” (in a sense to be made precise) the dif-
ferential operators.
where the unknown u is a function (say of class C 1 ) from R to E. The favorable case
is the case where u0 is an eigenvector of A, associated to an eigenvalue λ P C. In that
case, It is easy to see that we get a solution by setting, for any t P R,
uptq “ etλ u0 .
Indeed, the function u thus defined is of class C 8 from R to E, we have up0q “ u0 and,
for any t P R,
u1 ptq “ etλ λu0 “ etλ Au0 “ Auptq.
Since the equation u1 “ Au is linear with respect to u, we get that if u0 is a linear
combination of eigenvectors, that is if
k
ÿ
u0 “ αj e j ,
j“1
The interest of results about diagonalization is to give criteria ensuring that any vector
u0 of E can be written as a linear combination of eigenvectors of A. Thus, if A is diag-
onalizable and if we can determine its eigenvalues and eigenvectors, we can easily solve
Julien Royer
Université Toulouse III - Paul Sabatier
February 11, 2025
1
(1) for any initial data u0 .
It turns out that many models from concrete problems can be written in the form
(1), but in a functionnal space E of infinite dmension (typically L2 pRd q), and a linear
map A that is a differential operator. An important example (for example for the heat
equation, but also for Schrödinger equation or the wave equation) is the Laplace operator
d
ÿ B2
A“∆“ .
k“1
Bx2k
To simplify the rest of the discussion, we assume that d “ 1. The simplest differential
operator on R is
Avpxq “ v 1 pxq.
We look for eigenvectors of A, that is functions v such that v 1 “ λv for some λ P C. The
candidates are obviously the exponential functions x ÞÑ eλx . These functions appears in
the Fourier series theory, that have precisely been introduced to solve the heat equation,
but on a bounded interval. According to the Fourier series theory, if we set
2iπnx
en : x ÞÑ e T , for n P Z,
then the family pen qnPZ is a Hilbert basis of the Hilbert space L2T pRq of locally integrable
and T -periodic functions from R to C, endowed with the norm defined by
żT
2 1 2
}f }L2 pRq “ |f pxq|2 dx.
T T ´ T
2
This holds for T -periodic functions, but we can do the same for the functions defined
on an interval of length T .
Here, we are interested in functions that are not periodic and are defined on all R.
Formally, we want to see what happens if we let T go to `8 in the previous expressions.
Observing that for ϕ P C08 pRq we have
ˆ ˙ ż
2π ÿ 2πn
ϕ ÝÝÝÝÝÑ ϕpξq dξ (2)
T nPZ T T Ñ`8 R
2
(this is essentially a Riemann sum), we want to write f pxq as
ż ˆż ˙
1
e ´iyξ
f pyq dy eixξ dξ. (3)
2π R R
If this has a meaning, and if it really gives an expression for f pxq, then we have
written f as an integral (that can be seen as the limit of a series, see (2), or as a “con-
tinuous sum”) of functions of the form cξ pf qeixξ , that are “eigenvectors” (it would be
necessary to specify the space in which one works) associated to the derivative operator.
In addition the corresponding coefficient cξ pf q can be written as a scalar product in
L2 pRq of f with the conjugate function x ÞÑ eixξ “ e´ixξ . Except that the function
x ÞÑ eixξ is not in L2 pRq. . .
Thus, the definition of the Fourier transform on R (or Rd for any d P N˚ ) is motivated
and guided by good properties of the Fourier series on a bounded interval, but to get
similar results in this setting, a specific analysis is necessary. This is the purpose of this
chapter.
3
This proves that f and g are rapidly decreasing functions, then pf ˚ gq is also a rapidly
decreasing function.
By derivation under the integral sign, we see that for f, g P SpRd q the convolution
pf ˚ gq is of class C 8 on Rd and for β P Nd we have B β pf ˚ gq “ pB β f q ˚ g. The previous
result applied to B β f and g ensures that the derivatives of pf ˚ gq are rapidly decreasing,
and so pf ˚ gq P SpRd q.
Since we will also work in different Lebesgue spaces, we give some useful links between
these different spaces.
Proposition 1.5. (i) For any p P r1, `8s we have SpRd q Ă Lp pRd q.
(ii) For p P r1, `8r, SpRd q is dense in Lp pRd q. More precisely, for p, q P r1, `8r and
f P Lp pRd q X Lq pRd q there exists a sequence pfn qnPN of Schwartz functions that
converges to f in both Lp pRd q and Lq pRd q.
The first point is easy, for the second one we use the fact that these properties of
density are valid for C08 pRd q, which is included in SpRd q.
We finish with the property of integration by parts on Rd . The fact that Schwartz
functions ar rapidly decreasing at infinity ensures that there is no difficulty with the
boundary terms to show the following result.
Proposition 1.6. Let u P C 1 pRd q be a slowly increasing function whose first order
partial derivatives are also slowly increasing. Let v P SpRd q. Then for j P J1, dK the
functions pBj uqv and upBj vq are integrable on Rd and we have
ż ż
pBj uqv dx “ ´ upBj vq dx.
Rd Rd
2 Fourier transform in L1
We begin by defining what will play the role of “Fourier coefficient” for a non periodic
function. It is in L1 pRd q that the following definition has a natural meaning.
Definition 2.1. Let f P L1 pRd q. For ξ P Rd we set
ż
fˆpξq “ e´ix¨ξ f pxq dx. (4)
Rd
4
1
Example 2.2. Let a ą 0. We consider on R the function f : x ÞÑ 2a 1r´a,as . Then
f P L1 pRq, fˆp0q “ 1 and for ξ P R˚ we have
1 a ´ixξ 1 e´iaξ ´ eiaξ
ż
ˆ sinpaξq
f pξq “ e dx “ “ .
2a ´a 2a ´iξ aξ
Example 2.3. Let a ą 0. For x P R we set f pxq “ e´a|x| . Then f P L1 pRq and for ξ P R
we have (exercise)
2a
fˆpξq “ .
a2 ` ξ 2
Remark 2.4. Different conventions co-exist for the definition of fˆ. We can also set
ż
ˆ 1
f pξq “ d e´ix¨ξ f pxq dx.
p2πq 2 Rd
The advantage of this definition is that the expression of the Fourier transform and its
inverse transform (see Proposition 3.1) are symetrical, and there is no 2π factor in the
d
Parseval identity (see Proposition 5.1). However a p2πq 2 factor appears for example for
the convolution (see Proposition 2.7). Another convention is to set
ż
ˆ
f pξq “ e´2iπx¨ξ f pxq dx.
Rd
None of these choices is better than the others, if we remove the 2π factor somewhere it
will appear somewhere else in the theory. Thus, when we use the Fourier transform, it
is better to recall which definition is used and to be consistent afterwards (even if the
p2πqd factor is in general harmless for the studied properties).
We now give a number of basic properties for the Fourier transform.
Proposition 2.5. (i) For f P L1 pRd q, fˆ is bounded and }fˆ}8 ď }f }1 .
5
The proof of the following properties are left as exercises for the reader:
(i) If f is an even (respectively odd) function, then fˆ is even (respectively odd). More
generally, if we note P the operator that associates Pf : x ÞÑ f p´xq to f , then for
f P L1 pRd q we have
PFf “ FPf. (5)
A remarkable property (which will play an important role in applications) is that the
Fourier transform of a convolution product is the usual product of Fourier transforms.
˚ g “ fˆĝ.
fz
Proof. Let ξ P Rd . The map px, yq ÞÑ e´ix¨ξ f px ´ yqgpyq is measurable. By the Fubini-
Tonelli theorem and the properties of the convolution product in L1 pRd q we have
ż ˇ ˇ ż ˆż ˙
ˇ ´ix¨ξ
ˇe f px ´ yqgpyqˇ dx dy ď |f px ´ yq| |gpyq| dy dx ď }f }1 }g}1 ă `8.
ˇ
Rd ˆRd Rd Rd
3 Inversion formula
The property which makes the Fourier transform useful in applications is that we can
recover a function if we know its Fourier transform.
Proposition 3.1. Let f P L1 pRd q such that fˆ P L1 pRd q. Then for almost every x P Rd
we have ż
1
f pxq “ eix¨ξ fˆpξq dξ.
p2πqd Rd
6
We observe that the formula giving the inverse Fourier transform is almost the
same as the one giving the Fourier transform itself (only the sign in the exponential
changes, and with our choice for the definition of F there is an extra factor p2πq´d ).
For g P L1 pRd q and x P Rd we set
ż
1 1
Fgpxq “ d
eix¨ξ gpξq dξ “ pFgqp´xq.
p2πq Rd p2πqd
Using the notation P introduced in proposition 2.6 this can be rewritten as
1 1
F“ PF “ FP.
p2πqd p2πqd
Thus we can write
1 1
f “ FFf “ d
PF Ff “ FF Pf “ FFf. (6)
p2πq p2πqd
By continuity under the integral sign, these two functions are continuous on Rd . More-
over, by the dominated convergence theorem we have, for all x P Rd ,
Then
p2εqd f pyq
ż
1
Fε pxq “ dy “ pχε ˚ f qpxq.
p2πqd
śd 2
Rd j“1 pε ` pxj ´ yj q2 q
7
Since pχε qεPs0,1s defines an approximation of unity, we have Fε P L1 pRd q and
In particular, there exists a sequence pεn qnPN that tends to 0 and such that Fεn tends
to f almost everywhere. This proves that F “ f almost everywhere.
Thus we have proved that under the hypotheses of the proposition 3.1, the expres-
sion obtained formally in (3) holds for almost any x P Rd .
We observe that Proposition 3.1 gives injectivity for the Fourier transform in L1 pRd q.
Proposition 4.1. Let f P SpRd q. Let j P J1, dK. Then for ξ P R we have
ˆ
xj f pξq “ iξj f pξq
Bz
This property can in fact be extended to any function f P L1 pRd q X C 1 pRd q such
that the partial derivative Bxj f is also in L1 pRd q.
Since the inverse Fourier transform has an expression analogous to the Fourier trans-
form itself, it is not surprising to find that we actually have an analogous property in
the other direction. That is, the Fourier transform changes the multiplication by the
variable xj into a derivative with respect to ξj :
Proposition 4.2. Let f P SpRd q. Let j P J1, dK. Then fˆ is derivable with respect to ξj
and for any ξ P Rd we have
x
y ˆ
j f pξq “ iBξj f pξq.
8
Proof. For px, ξq P Rd ˆ Rd we set φpx, ξq “ e´ix¨ξ f pxq. Then φ is differentiable with
respect to ξj and for any px, ξq P pRd q2 we have
ˇ ˇ ˇ
ˇ Bφ ˇ ˇ ˇ
ˇ px, ξq ˇ “ ˇ´ixj e´ixξ f pxqˇˇ “ |xj f pxq| .
ˇ Bξj ˇ
Since x ÞÑ |xj f pxq| is integrable on Rd , we get by differentiation under the integral sign
that fˆ is differentiable with respect to ξj in Rd and, for ξ P Rd ,
ż
ˆ
iBξj f pξq “ e´ix¨ξ xj f pxq dx “ x
y j f pξq.
Rd
Notice that this proof holds for any function f P L1 pRd q such that x ÞÑ xj f pxq is
integrable.
The Schwartz space is the space of functions for which one can derive and multiply by
the variable as many times as one wants. Moreover, these two operations play symmetric
roles with respect to the Fourier transform. The two previous propositions suggest that
it is a space in which the Fourier transform is particularly convenient.
Proposition 4.3. The Fourier transform F defines a bijection from SpRd q to itself,
and its inverse is the restiction of F to SpRd q. Moreover, for f P SpRd q, α, β P Nd and
ξ P Rd we have
ξ α Bξβ Ff pξq “ p´iq|α|`|β| FpBxα xβ f qpξq. (7)
Proof. Let f P SpRd q and α, β P Nd . Applying |β| times Proposition 4.2, we get that
the partial derivative Bξβ fˆ exists and
α β |α|`|β| α β ˆ
x x f pξq “ i
B{ ξ Bξ f pξq.
In particular, the function ξ ÞÑ ξ α B β fˆpξq is bounded. This ensures that fˆ P SpRd q and
proves (7). The restriction FS of F to SpRd q is then injective by Corollary 3.2. Finally,
1
according to (6) we have f “ p2πq 2d FF Pf , so FS is surjective. Thus FS is a bijection
Example 4.4. Let f P SpRd q. Then the Fourier transform of ´∆f is ξ ÞÑ |ξ|2 fˆpξq.
Example 4.5. We want to compute on R the Fourier transform of the Gaussian f : x ÞÑ
x2
e´ 2 . Notice that f P SpRq. In particular fˆ P SpRq. Since f 1 pxq “ ´xf pxq for any
x P R, we have that for ξ P R
9
But we know that ż
x2 ?
fˆp0q “ e´ 2 dx “ 2π,
R
So for any ξ P R we have
? ξ2
fˆpξq “ 2πe´ 2 .
x2
More generally, for σ ą 0, if we consider fσ : x ÞÑ e´ 2σ2 then for ξ P R we have
? 2 2
´ σ 2ξ
σ pξq “ σ 2πe
fx .
5 Fourier transform in L2
So far we have discussed the properties of the Fourier transform in L1 pRd q (because it
is the space in which the definition naturally makes sense) or of its restriction to SpRd q
(where it is the most comfortable). But our favorite space in the analysis of partial
differential equations will be L2 pRd q, because it is a Hilbert space. The problem is that
a function in L2 pRd q is not necessarily integrable, and the Fourier transform as defined
in (4) does not even make sense in this case.
We start by going a little further with the analysis of the Fourier transform for
Schwartz functions. Since SpRd q is included in L2 pRd q, it can be endowed with the
corresponding norm. The Fourier transform of SpRd q is then a quasi-isomery for this
norm:
In particular,
d
}fˆ}L2 “ p2πq 2 }f }L2 .
10
Corollary 5.2. For f, ϕ P SpRd q we have
ż ż
fˆϕ dy “ f ϕ̂ dy.
Rd Rd
Now we use the fact that SpRd q is dense in L2 pRd q to extend the Fourier transform
on SpRd q into an ismoetry of L2 pRd q.
Theorem 5.3. There exists a unique isomorphism F of L2 pRd q which coincides with
the Fourier transform on SpRd q and such that
d
@f P L2 pRd q, }Ff }L2 pRd q “ p2πq 2 }f }L2 pRd q . (8)
Proof. ‚ Assume that two isomorphisms of L2 pRd q satisfy the conditions of the theo-
rem. They are in particular continuous applications on L2 pRd q which coincide on the
dense subset SpRd q. Then they coincide everywhere, which gives uniqueness.
‚ Let f P L2 pRd q. There exists a sequence pfn qnPN of Schwartz functions that converge
to f in L2 pRd q. By Proposition 5.1, the sequence pfx 2 d
n qnPN is a Cauchy sequence in L pR q.
2 d 2 d
Since L pR q is complete, fx n has a limit in L pR q, which we denote by Ff . We can
check that this definition of Ff does not depend on the choice of the sequence pfn qnPN
and that for f P SpRd q the limit obtained this way coincides with the Fourier transform
of f . By linearity of the Fourier transform on SpRd q and by Proposition 5.1, the map
F is linear and satisfies (8). By continuity, the equality (6) still holds in L2 pRd q, which
proves that F is surjective.
and ż ż
pFf qg dy “ f pFgq dy.
Rd Rd
Proposition 5.5. The isomorphism of Theorem 5.3 coincides with the Fourier trans-
form already defined on L1 pRd q X L2 pRd q.
Proof. Let f P L1 pRd q X L2 pRd q. For this proof, we denote by fˆ the Fourier transform
of f in the sense (4), and by Ff the function of L2 pRd q given by Theorem 5.3. Let
pfn qnPN be a sequence of Schwartz functions which converges to f both in L1 pRd q and
L2 pRd q. With the same convention as for f , we have fˆn “ Ff for any n P N. Since F
is continuous on L2 pRd q we have
}Ffn ´ Ff }2 ÝÝÝÝÑ 0.
nÑ`8
In particular, possibly after extracting a subsequence, we can assume that Ffn tends to
Ff for almost every x P Rd . On the other hand we have
11
Recall that for a non integrable function, we cannot use the expression of Ff given
by (4). We at least have the following property, where f˜n is a kind of Fourier transform
for locally integrable functions.
Proposition 5.6. Let f P L2 pRd q. For n P N and ξ P Rd we set
ż
f˜n pξq “ e´ix¨ξ f pxq dx,
Bpnq
6 Example of application
For f P SpRd q, we consider the equation
´∆u ` u “ f, (9)
of unknown u P SpRd q. Taking the Fourier transform, we see that u is a solution if and
only if for any ξ P Rd we have
or equivalently
fˆpξq
ûpξq “ .
1 ` |ξ|2
fˆpξq
We consider the function f˜ which maps ξ P Rd to 1`|ξ|2
. This defines a function
f˜ P SpRd q and we have
d
}f˜}L2 ď }fˆ}L2 “ p2πq 2 }f }L2 .
We then set Rf “ F ´1 f˜. This also defines a function of SpRd q and we have
d
}Rf }L2 “ p2πq´ 2 }f˜}L2 ď }f }L2 .
It is not very difficult to check that the map f ÞÑ Rf is linear. Thus R defines a linear
map from SpRd q to itself, continuous for the L2 norm. It is therefore extended by con-
tinuity into a continuous linear map on L2 pRd q.
p´∆ ` 1qRf “ f.
12
We can then check that for any u P SpRd q we have
Rp´∆ ` 1qu “ u.
A final question about this example. Given f P SpRd q, can we give a slightly more
explicit expression for the solution Rf of the problem (9)? We have to compute the
inverse Fourier transform of the product of two functions. This is where the convolution
product appears (see Proposition 2.7). Thus, if G is a function which Fourier transform
is given by
1
Ĝpξq “ ,
1 ` |ξ|2
then we will have
x pξq “ Ĝpξqfˆpξq,
Rf
and hence
Rf pxq “ pG ˚ f qpxq.
It remains to determine G. At least for the dimension d “ 1, we have already done the
computations at example 2.3 and we get
e´|x|
Gpxq “ .
2
This gives
e´|x´y| e´|y|
ż ż
Rf pxq “ f pyq dy “ f px ´ yq dy.
R 2 R 2
We can now check with an explicit computation (left as an exercise) that if f P SpRq
then for any x P R we have
In dimension 1, the result is not very spectacular, because it is only an ordinary differ-
ential equation with constant coefficients, and we could have got the same result using
only a double variation of constants, but the method described on this simple example
will be useful to understand more subtle problems. . .
13