5.6 Convergence of Fourier series

Figure 5.4 shows periodic function f with period 2π and its partial sum sequence of Fourier series.

PIC

Figure 5.4: f(x)=x2,x∈[−π,π] and its Fourier series

Theorem 5.6.1 (Pointwisely convergence) Assuming f is piecewisely continuous and piecewisely derivable; and ∀ jump discontinuity x0 of f, f′(x0+),f′(x0−) is finite, then ∀x, its Fourier series S is convergent, and S(x)=f(x+)+f(x−)2

Definition 5.6.2 The inner product of 2 functions in R[−π,π] is defined as ⟨f,g⟩=∫−ππf(x)g(x)dx The L2-norm is defined as ∥f∥L2=⟨f,f⟩=∫−ππ|f(x)|2dx

Definition 5.6.3 (Convergence in the mean square) If ∥fn−f∥L2→0 when n→+∞, then fn converges to f in the mean square, marked as fn→L2f.

Theorem 5.6.4 (Least square property) Assuming SN=A02+∑n=1N(Ancos⁡nx+Bnsin⁡nx) is the partial sum of Fourier series of f, mark WN={α02+∑n=1N(αncos⁡nx+βnsin⁡nx)|α0,αnβn∈R} then L2π2[−π,π]={Periodic square-integrable function with period 2π} is a infinite-dimension linear space, WN⊂L2 is the linear subspace and dim⁡WN=2N+1. Assuming T∗ satisfies minT∈WN∥f−T∥=∥f−T∗∥ then T∗=SN.

Proof Notice that SN=T∗⇔(f−SN)⊥SN⇔(f−SN)⊥{1,cos⁡nx,sin⁡nx|1≤n≤N} Since ⟨f−T,1⟩=⟨f,1⟩−⟨T,1⟩=πA0−πα0=0⇔A0=α0⟨f−T,cos⁡nx⟩=⟨f,cos⁡nx⟩−⟨T,cos⁡nx⟩=πAn−παn=0⇔An=αn

Similarly we have Bn=βn. So SN=T∗. ◻

Corollary 5.6.5 (Bessel inequality) π(A022+∑n=1N(An2+Bn2))≤∫−ππ|f(x)|2dx

Proof Notice that ∥f∥L22=∥f−SN∥2+∥SN∥2≥∥SN∥2=‖A02+∑n=1N(Ancos⁡nx+Bnsin⁡nx)‖2=A024⟨1,1⟩+∑n=1N(An2⟨cos⁡nx,cos⁡nx⟩+Bn2⟨sin⁡nx,sin⁡nx⟩)

So Bessel inequality is proved. ◻

Corollary 5.6.6 Series π(A022+∑n=1+∞(An2+Bn2)) is convergent, and π(A022+∑n=1+∞(An2+Bn2))≤∫−ππ|f(x)|2dx This could be used to judge if a trignometric series is a Fourier series.

Corollary 5.6.7 (Parseval’s identity) Notice that SN→L2f⇔∥SN−f∥L2→0⇔∥SN∥L2→∥f∥L2 we have SN→L2f⇔π(A022+∑n=1+∞(An2+Bn2))=∫−ππ|f(x)|2dx

Theorem 5.6.8 For any f∈L2, Parseval’s identity holds.

Corollary 5.6.9 For any f∈L2, SN→L2f.

Corollary 5.6.10 Assuming f(x)=A02+∑n=1+∞(Ancos⁡nx+Bnsin⁡nx) and g(x)=α02+∑n=1+∞(αncos⁡nx+βnsin⁡nx), then ⟨f,g⟩=π[A0α02+∑n=1+∞(Anαn+Bnβn)]

We can use this to compute Riemann-ζ function. Define ζ(z)=∑n=1+∞1nz, we have ζ(2)=π26, ζ(4)=π490, ...

Applications of Fourier series.

Example 5.6.1 Isoperimetric inequality. Given a closed curve γ where its length is L and the area of its interior is A. Assuming its parametric representative is x,y:[0,L]→R where t∈[0,L] is the arc-length parameter and x,y∈C2[0,L] is periodic with period L, then we have 4πA≤L2 where it takes the equality if and only if γ is a circle.

Assuming x(t)=A02+∑n=1+∞(Ancos⁡2nπtL+Bnsin⁡2nπtL)y(t)=α02+∑n=1+∞(αncos⁡2nπtL+βnsin⁡2nπtL)x′(t)=∑n=1+∞(−2nπLAnsin⁡2nπtL+2nπLBncos⁡2nπtL)y′(t)=∑n=1+∞(−2nπLαnsin⁡2nπtL+2nπLβncos⁡2nπtL)

Since t is the arc-length parameter, we have x′(t)2+y′(t)2=1. Then we have L2=L∫0L(x′(t)2+y′(t)2)dt=L(⟨x′,x′⟩+⟨y′,y′⟩)=L[L24π2L2∑n=1+∞n2(An2+Bn2)+L24π2L2∑n=1+∞n2(αn2+βn2)]=2π2∑n=1+∞n2(An2+Bn2+αn2+βn2)A=12∮(−ydx+xdy)=12∫0L(−y(t)x′(t)+x(t)y′(t))dt=12(−⟨y,x′⟩+⟨x,y′⟩)=π∑n=1+∞n(−αnBn+Anβn)

Notice that L2≥2π2∑n=1+∞n(Bn2+αn2+An2+βn2)≥4π2∑n=1+∞n(−αnBn+Anβn)=A Takes the equality.

  • ∀n≥2, An=Bn=αn=βn=0.
  • ∀n∈N∗, An=βn, αn=−Bn.

Therefore, {x(t)=A02+A1cos⁡2πtL+B1sin⁡2πtLy(t)=α02−B1cos⁡2πtL+A1sin⁡2πtL⇒(x−A02)2+(y−α02)2=A12+B12 meaning it’s a circle.

Example 5.6.2 Seek the periodic solution with period 2π to y″+λy=sin⁡x, i.e. {y″+λy=sin⁡xy(0)=y(2π)y′(0)=y′(2π) This is called the boundary problem.

Assuming y(x)=A02+∑n=1+∞(Ancos⁡nx+Bnsin⁡nx), then y″(x)=−∑n=1+∞n2(Ancos⁡nx+Bnsin⁡nx) After substitution, we have

  • Constant term. λA02=0.
  • cos⁡nx term. −n2An+λAn=0.
  • sin⁡nx term. −n2Bn+λBn=δn1.

Discussion on λ.

  • When λ=0, then A0 is arbitrary; An=0, ∀n≥1; Bn=−δn1; so y(x)=A02−sin⁡x are all periodic solutions with period 2π. Infinite solution, dimension=1, affine space.
  • When λ≠n2(n=0,1,...), then An=0, ∀n≥0; Bn=0, ∀n>1; B1=1λ−1; so y(x)=1λ−1sin⁡x is the unique periodic solution with period 2π.
  • When λ=1, then A0=0; A1 is arbitrary; An=Bn=0, ∀n≥1; yet −B1+B1=1 is impossible; so there doesn’t exist a periodic solution with period 2π.
  • When λ=m2(m>1), then An=0, ∀n≠m; Am is arbitrary; B1=1m2−1; Bn=0, ∀n≠m and n≠1; Bm is arbitrary; so y(x)=sin⁡nxm2−1+Amcos⁡mx+Bmsin⁡mx are all periodic solutions with period 2π. Infinite solution, dimension=2, affine space.

Assuming a linear transformation operator L:Ly=−y″ which operates on the linear space V composed of all periodic functions with period 2π, then the original equation could be rewritten as −Ly+λy=sin⁡x Recall in linear algebra, given x,b∈Rn and a square matrix A∈Mn, for the equation −Ax+λx=b⇔(λI−A)x=b Assuming λ is not the eigenvalue of A, then λI−A is invertible, so ∀b∈Rn, (λI−A)x=b has a unique solution.

Assuming λ is the eigenvalue of A, then (λI−A)x=b has a solution ⇔ b∈Im(λI−A). Assuming A is symmetric, then we have Rn=ker⁡(λI−A)⊕⊥Im(λI−A) Here b∈Im(λI−A)⇔b⊥ all eigenvectors of A corresponded to eigenvalue λ.

L is symmetric since ⟨Lf,g⟩=⟨f,Lg⟩. ∀f,g, i.e. ⟨Lf,g⟩=∫02π(−f″)gdx=−∫02πgdf′=−f′g|02π+∫02πf′g′dx=∫02πf′g′dx⟨f,Lg⟩=∫02π(−g″)fdx=−∫02πfdg′=−g′f|02π+∫02πf′g′dx=∫02πf′g′dx

Assuming λ is eigenvalue of L, i.e. (λI−L)y=λy+y″=0.

  • λ<0, only zero solution in V, so λ is not eigenvalue.
  • λ=0 is the eigenvalue, y≡1 is the base of eigen-subspace, dim=1.
  • λ>0,λ≠n2(n=1,2,...), only zero solution in V, so λ is not the eigenvalue.
  • λ=n2 is the eigenvalue, cos⁡nx,sin⁡nx are the bases of eigen-subspace, dim=2.