Posts

Showing posts from February, 2023

Pointwise convergence of Fourier series

Image
Our main question has been: What type of functions can be represented as a Fourier series? We have seen two types of convergence of S[f] to f.  Uniform convergence: If fC2π1 (continuously differentiable), then S[f] converges uniformly to f. L2 convergence: If fR2π2 (square-integrable), then S[f] converges to f in the L2 norm. These are ``global'' convergence properties, so to speak, in the sense that they characterize the function f for the whole period [π,π]. Now, we consider ``local'' convergence properties such as the problem of pointwise convergence.  To begin with, we have the following result: Theorem (du Bois-Reymond) Let f be a periodic continuous function. Then, the Fourier series of f may not converge pointwise. But we don't really need the continuity of f for its Fourier series to converge to f(x0) at x0. Thus, let's define the following. Defi...

L2 convergence of Fourier Series

Image
Let us approximate a function f(x) by a linear combination of {eikx|k=n,,0,,n}. We show that the partial Fourier sum is the ``best'' of such approximators. Then, we prove that the Fourier series of the function f converges to f in the L2 norm. Lemma (``best'' approximator) Let any nN, consider the set of 2n+1 functions of the form φk(x)=eikx   (|k|n). Let vn be an arbitrary linear combination of these functions: vn(x)=|k|nγkφk(x) where γkC are constants. Suppose fC2π0. Then, fSnfvn where Sn=Sn[f] is the partial sum of the Fourier series S[f]. In other words, Sn[f] is the best approximator of f among all the linear combinations of {φk}|k|n with respect to the L2 norm. Proof .  \[\begin{eqnarray*} \|f - v_n\|^2 & = & \|(f ...

Kernels arising from the Fourier series

Image
Here, we present a different way of looking at the Fourier series using  kernels . We see that the partial Fourier sum of a function f(x) can be expressed as the convolution between f and the Dirichlet kernel. We have already seen that the partial Fourier sum converges uniformly to f if f(x) is continuously differentiable. Similarly, we show that the convolution of the function f and the Fejér kernel converges uniformly to f, but this time, if f(x) is continuous. These kernels are used to prove the L2 convergence of the Fourier series . Definition (Kernel function) Let X be a set, and K:X×XR be a two-variable function. K(x,y) is said to be a  kernel function  (or  integral kernel  or  nucleus ) if it satisfies the following conditions: (symmetry) For any x,yX, K(x,y)=K(y,x).  (positive semi-definiteness) For any nN, {x1,x2,,xn}X...

Vector space of functions

Image
We study the collection of functions R2π2 (square-integrable functions with period 2π) as a vector space. We define the L2 norm and L2 inner product on this vector space so that we can investigate the ``geometric'' structure of the space of functions.  Let us show that R2π2, the set of square-integrable functions with period 2π, is a vector space over C. First, we need to define addition and scalar multiplication. Let f,gR2π2. We define f+gR2π2 by (eq:add)(f+g)(x)=f(x)+g(x), xR. Note that the ``+'' on the left-hand side is defined between the two functions f and g, whereas the ``+'' on the right-hand side is the addition between two complex numbers f(x) and g(x). Next, we define scalar multiplication. Let αC and fR2π2. We define \(\alpha ...

Fourier series of piece-wise smooth functions

Image
Our goal in this post is to prove the following: If f is ``piece-wise smooth,'' then the (complex) Fourier series of f converges uniformly to f. Note that the condition in this theorem does not involve the Fourier coefficients of f; hence, it is more direct. That is, we can tell if the Fourier series of a function converges uniformly from the property of the function alone. Let C2πm be the collection of all functions of class Cm such that all derivatives up to the m-th order have period 2π .  We say a function f is piece-wise smooth if fC2π0 (i.e., continuous) and fR2π2 (i.e., the derivative is square-integrable, not necessarily continuous). In the following, we will consider the complex Fourier series: fk=ckeikx where the Fourier coefficients are ck=12πππf(x)eikxdx. Lemma (Bessel's inequality) If \(f\in\mathca...

Complex Fourier series

Image
The Fourier series we have studied so far is a special case of the complex Fourier series. In many cases, the complex Fourier series is more convenient than the (real) Fourier series. The complex Fourier series is defined as a series of the form k=ckeikx where i=1 is the imaginary unit and ckC. Note that the variable x is still real, xR.  In the following, the collections R2π1 and R2π2 also include complex-valued functions. If fR2π1, the complex Fourier coefficients ck of f are defined by (eq:ck)ck=12πππf(x)eikxdx. Here, the right-hand side is the definite integral of a complex-valued function. See also : Calculus of complex-valued functions For m,nZ, if mn, then \[\begin{eqnarray*} \int_{-\pi}^{\pi}e^{imx}e^{-inx}dx &=& \int_{-\pi}^{\pi}e^...

Fourier sine series, Fourier cosine series

Image
When a function is an odd function, its Fourier series contains only the sine functions. Such a Fourier series is called a Fourier sine series . When a function is an even function, its Fourier series contains only the cosine functions (and possibly a constant). Such a Fourier series is called a Fourier cosine series . Example . Consider the function fP2π defined by f(x)=x   (π<x<π). Note that f is an odd function (i.e., f(x)=f(x)). The Fourier coefficients are obtained as (exercise!) an=1π02πf(x)cos(nx)dx=1πππf(x)cos(nx)dx=0bn=1π02πf(x)sin(nx)dx=1πππf(x)sin(nx)dx.=(1)n+12n. Thus, we have x2[sinx1sin(2x)2+sin(3x)3+(1)n+1sin(nx)n+] f...

Collections of periodic functions

Image
To study the Fourier series, it is convenient to define collections of functions with a period of 2π. Definition (Collection of periodic functions) We denote by P2π the collection (set) of all functions on R that have the period of 2π and that have at most finitely many discontinuous points in each period. For each fP2π, its value may not be defined at its discontinuous points. Definition (Collections of integrable functions) We denote by R2π1 the collection of all fP2π such that f absolutely integrable on the period, that is, 02π|f(x)|dx<+. We denote by R2π2 the collection of all fP2π such that f is square-integrable, that is, 02π|f(x)|2dx<+. In these definitions, the integral may be improper. Remark . The R in R2π1 and \(\mathcal{R}_{2\...

Fourier series and the heat equation

Image
Joseph Fourier introduced the Fourier series to solve the heat equation in the 1810s. In this post, we show how the Fourier transform arises naturally in a simplified version of the heat equation. Suppose we have the unit circle S made of a metal wire. Pick an arbitrary point A on the circle. Any point P on the circle is identified by the distance x from A to P along the circle in the counter-clockwise direction (i.e., x is the angle of the section between A and P in radian). Let u(t,x) represent the temperature at position x and time t. The temperature distribution at t=0 is given by u(0,x)=f(x). Assuming no radiation of heat out of the metal wire, u(t,x) for t>0 and 0x2π is determined by the following partial differential equation  (PDE) called the  heat equation : γut=κ2ux2 and the initial condition \[u(0,x) = f(x...

Uniform convergence of Fourier series

Image
We study some of the basic properties of the Fourier series. In particular, we show that if the series of the Fourier coefficients of f converges absolutely, then Fourier series S[f] converges uniformly to f at the continuous points of f. Lemma (Condition for the uniform convergence of Fourier series) Let fR2π1 with its Fourier coefficients satisfying the conditions n=0|an|<+,n=1|bn|<+. Then, the Fourier series S[f] converges uniformly with respect to xR. In particular, S[f] is continuous. Proof . Consider the following positive (non-negative) term series: 12a0+n=1(|an|+|bn|). By assumption, this series converges. Furthermore, this is a dominating series of S[f]. Thus, S[f] converges uniformly. All terms of S[f] are continuous (they are jus...

Fourier series: Introduction

Image
The theory of the  Fourier series is based on a wild assumption: any "well-behaved" periodic function can be represented as a linear combination of sine and cosine functions, and that expression is unique for the given function. This theory (eventually) provided much of the foundations of modern mathematics. But it is also of tremendous practical importance. Trigonometric functions are periodic. But what is a periodic function in general? Definition (Periodic function) A function f(x) is said to be a periodic function if there exists a real number T>0 such that f(x+T)=f(x). In this case, T is called a period of the function. Note that if T is a period of f(x), then 2T, 3T, ,nT (nN) is are also periods of f(x). In fact,  f(x+nT)=f(x+(n1)T+T)=f(x+(n1)T)==f(x). The smallest (non-zero) period is called the fundamental period . Remark . If we say simply a period, we usually mean ...

Arzelà's theorem

Image
Here we list two useful theorems without proofs. If you learn the theory of measure and Lebesgue integral, you will encounter more general versions of these theorems. However, the following theorems do not require the Lebesgue integral but hold for the Riemann integral. Theorem (Arzelà's convergence theorem) Suppose that the sequence of continuous functions fn=fn(x)   (n=1,2,) is uniformly bounded on the closed bounded interval I=[a,b]. That is, there exists an M>0 such that |fn(x)|M   (xI,n=1,2,). In addition, suppose that {fn} converges to a continuous function f0=f0(x) at each xI. That is, limnfn(x)=f0(x)   (xI). Then, we have limnabfn(x)dx=abf0(x)dx. In other words, the limit as n and the integral are ``commutative'' (i.e., the order of limit and integration can be swapped). Perhaps it is clearer if we write \[\...

Special solution of inhomogeneous linear differential equations

Image
The general solution of an inhomogeneous linear differential equation can be obtained as the sum of a special solution and the general solution of the corresponding homogeneous differential equation. We study the method of variation of parameters in particular. Consider the linear differential equation of the form F(D)y=q(x) where F(t) is a polynomial and q(x) is a function. When q(x)=0, this is a homogeneous linear differential equation. Here, we assume q(x)0. Suppose we can factorize F(t) as F(t)=G(t)H(t). If y=y(x) is a solution of F(D)y=q, then z=H(D)y is a solution of G(D)z=q. Thus, the given differential equation F(D)y=q is decomposed into two parts: G(D)z=q (a linear differential equation of z, given q), H(D)y=z (a linear differential equation of y, given z), and we can process one after the other. Thus, by factorizing the polynomial F(D), we only need to consider the case w...

Homogeneous linear differential equations with constant coefficients

Image
Consider the homogeneous linear differential equation (Eq:code)y(n)+an1y(n1)++a1y+a0y=0 where a0,a1,,an1R are constants.  Using the differential operator D=ddx and its polynomial (Eq:epoly)E=Dn+an1Dn1++a1D+a0, (Eq:code) can be expressed as Ey=0. Now, consider the polynomial of variable t F(t)=tn+an1tn1++a1t+a0. Then, we have E=F(D) and (Eq:code) is expressed as F(D)y=0. We can use the properties of polynomials to solve this type of differential equation. We need some results from algebra. Definition (Relatively prime, coprime) Polynomials F1(t) and F2(t) are said to be relatively prime  or coprime if the equations F1(t)=0 and F2(t)=0 do not have common solutions within C. Lemma The polynomials F1(t) and F2(t) in t are relatively prime if and only if th...