k and v = v 1 j + u 3 i + v 2

Similar documents
DEFINITION The inner product of two functions f 1 and f 2 on an interval [a, b] is the number. ( f 1, f 2 ) b DEFINITION 11.1.

(9) P (x)u + Q(x)u + R(x)u =0

Abstract inner product spaces

MATH34032: Green s Functions, Integral Equations and the Calculus of Variations 1

ENGI 9420 Lecture Notes 7 - Fourier Series Page 7.01

Physics 116C Solution of inhomogeneous ordinary differential equations using Green s functions

Green function and Eigenfunctions

McGill University Math 354: Honors Analysis 3 Fall 2012 Solutions to selected Exercises. g(x) 2 dx 1 2 a

Best Approximation in the 2-norm

Chapter 28. Fourier Series An Eigenvalue Problem.

(4.1) D r v(t) ω(t, v(t))

LECTURE 10: JACOBI SYMBOL

Wave Equation on a Two Dimensional Rectangle

LINEAR ALGEBRA APPLIED

Chapter 6 Techniques of Integration

Math 473: Practice Problems for the Material after Test 2, Fall 2011

AMATH 731: Applied Functional Analysis Fall Additional notes on Fréchet derivatives

Duke Math Meet

Quadratic Residues. Chapter Quadratic residues

NOTES ON HILBERT SPACE

Inner-product spaces

Discrete Least-squares Approximations

Math Theory of Partial Differential Equations Lecture 2-9: Sturm-Liouville eigenvalue problems (continued).

Chapter 6. Riemann Integral

g i fφdx dx = x i i=1 is a Hilbert space. We shall, henceforth, abuse notation and write g i f(x) = f

7.3 Problem 7.3. ~B(~x) = ~ k ~ E(~x)=! but we also have a reected wave. ~E(~x) = ~ E 2 e i~ k 2 ~x i!t. ~B R (~x) = ~ k R ~ E R (~x)=!

Families of Solutions to Bernoulli ODEs

Best Approximation. Chapter The General Case

Theoretical foundations of Gaussian quadrature

c n φ n (x), 0 < x < L, (1) n=1

Orthogonal Polynomials and Least-Squares Approximations to Functions

(PDE) u t k(u xx + u yy ) = 0 (x, y) in Ω, t > 0, (BC) u(x, y, t) = 0 (x, y) on Γ, t > 0, (IC) u(x, y, 0) = f(x, y) (x, y) in Ω.

Sturm-Liouville Theory

The Dirichlet Problem in a Two Dimensional Rectangle. Section 13.5

63. Representation of functions as power series Consider a power series. ( 1) n x 2n for all 1 < x < 1

Advanced Calculus: MATH 410 Uniform Convergence of Functions Professor David Levermore 11 December 2015

Polynomials and Division Theory

8 Laplace s Method and Local Limit Theorems

Supplement 4 Permutations, Legendre symbol and quadratic reciprocity

Analytical Methods Exam: Preparatory Exercises

Fundamental Theorem of Calculus

Mapping the delta function and other Radon measures

Space Curves. Recall the parametric equations of a curve in xy-plane and compare them with parametric equations of a curve in space.

Quadratic reciprocity

Orthogonal Polynomials

QUADRATIC RESIDUES MATH 372. FALL INSTRUCTOR: PROFESSOR AITKEN

USA Mathematical Talent Search Round 1 Solutions Year 25 Academic Year

Review of Gaussian Quadrature method

PRIMES AND QUADRATIC RECIPROCITY

PHYSICS 116C Homework 4 Solutions

Chapter 5. , r = r 1 r 2 (1) µ = m 1 m 2. r, r 2 = R µ m 2. R(m 1 + m 2 ) + m 2 r = r 1. m 2. r = r 1. R + µ m 1

Matrix Algebra. Matrix Addition, Scalar Multiplication and Transposition. Linear Algebra I 24

The Dirac distribution

a 2 +x 2 x a 2 -x 2 Figure 1: Triangles for Trigonometric Substitution

ODE: Existence and Uniqueness of a Solution

Chapter 8: Methods of Integration

Lecture 3. In this lecture, we will discuss algorithms for solving systems of linear equations.

Sections 5.2: The Definite Integral

The Riemann Integral

Taylor Polynomial Inequalities

Partial Derivatives. Limits. For a single variable function f (x), the limit lim

Definite Integrals. The area under a curve can be approximated by adding up the areas of rectangles = 1 1 +

10 Vector Integral Calculus

10 Elliptic equations

Fourier series. Preliminary material on inner products. Suppose V is vector space over C and (, )

LECTURE 3. Orthogonal Functions. n X. It should be noted, however, that the vectors f i need not be orthogonal nor need they have unit length for

1.2. Linear Variable Coefficient Equations. y + b "! = a y + b " Remark: The case b = 0 and a non-constant can be solved with the same idea as above.

UNIFORM CONVERGENCE MA 403: REAL ANALYSIS, INSTRUCTOR: B. V. LIMAYE

Summary: Method of Separation of Variables

Name of the Student:

7. Indefinite Integrals

Integration Techniques

Chapter 6. Infinite series

Integrals - Motivation

MATH 174A: PROBLEM SET 5. Suggested Solution

Review of Riemann Integral

STURM-LIOUVILLE THEORY, VARIATIONAL APPROACH

Chapter 3 Single Random Variables and Probability Distributions (Part 2)

The Regulated and Riemann Integrals

Lecture 1. Functional series. Pointwise and uniform convergence.

The area under the graph of f and above the x-axis between a and b is denoted by. f(x) dx. π O

1 The Lagrange interpolation formula

MA Handout 2: Notation and Background Concepts from Analysis

21.6 Green Functions for First Order Equations

STURM-LIOUVILLE PROBLEMS

Consequently, the temperature must be the same at each point in the cross section at x. Let:

Some integral inequalities on time scales

Chapter 3 The Schrödinger Equation and a Particle in a Box

Hilbert Spaces. Chapter Inner product spaces

Mrgolius 2 In the rticulr cse of Ploue's constnt, we tke = 2= 5+i= 5, nd = ;1, then ; C = tn;1 1 2 = ln(2= 5+i= 5) ln(;1) More generlly, we would hve

The Fundamental Theorem of Calculus

Math 554 Integration

df dt f () b f () a dt

STURM-LIOUVILLE BOUNDARY VALUE PROBLEMS

Bases for Vector Spaces

A BRIEF INTRODUCTION TO UNIFORM CONVERGENCE. In the study of Fourier series, several questions arise naturally, such as: c n e int

Numerical Methods I Orthogonal Polynomials

Math Advanced Calculus II

Improper Integrals, and Differential Equations

Lecture 2e Orthogonal Complement (pages )

Transcription:

ORTHOGONAL FUNCTIONS AND FOURIER SERIES Orthogonl functions A function cn e considered to e generliztion of vector. Thus the vector concets like the inner roduct nd orthogonlity of vectors cn e extended to functions. Inner roduct Consider the vectors u = u 1 i + u 2 j + u 3 k nd v = v 1 i + v 2 j + v 3 k in R 3, then the inner roduct or dot roduct of u nd v is rel numer, sclr, defined s (u,v) = u 1 v 1 + u 2 v 2 + u 2 v 2 = 3 u k v k k=1

The inner roduct ossesses the following roerties (u,v) = u 1 v 1 + u 2 v 2 + u 2 v 2 = 3 u k v k k=1 (i) (u,v) = (v, u) (ii) (ku, v) = k(u, v) (iii) (u, u) = 0 if u = 0 nd (u, u) > 0 if u 0 (iv) (u + v, w) = (u, w) + (v, w)

Suose tht f 1 nd f 2 re iecewise continuous functions defined on n intervl [, ]. The definite integrl on the intervl of the roduct f 1 (x) f 2 (x) ossesses roerties (i) - (iv) ove. Definition: Inner roduct of functions The inner roduct of two functions f 1 nd f 2 on n intervl [, ] is the numer ( f 1, f 2 ) = f 1 (x) f 2 (x) dx

Definition: Orthogonl functions Two functions f 1 nd f 2 re sid to e orthogonl on n intervl [, ] if ( f 1, f 2 ) = f 1 (x) f 2 (x) dx = 0 Exmle: f 1 (x) = x 2 nd f 2 (x) = x 3 re orthogonl on the intervl [ 1, 1] since ( f 1, f 2 ) = 1 1 x 2. x 3 dx = 1 1 6 x6 1 = 0

Orthogonl sets We re rimrily interested in n infinite sets of orthogonl functions. Definition: Orthogonl set A set of rel-vlued functions {φ 0 (x), φ 1 (x), φ 2 (x),...} is sid to e orthogonl on n intervl [, ] if (φ m, φ n ) = φ m (x) φ n (x) dx = 0, m n

Orthonorml sets The norm u of vector u cn e exressed using the inner roduct: (u, u) = u 2 u = (u, u) Similrly the squre norm of function φ n is φ n 2 = (φ n, φ n ), nd so the norm is φ n = (φ n, φ n ). In other words, the squre norm nd the norm of function φ n in n orthogonl set {φ n (x)} re, resectively, φ n 2 = φ 2 n(x) dx nd φ n = φ 2 n(x) dx

Exmle 1: Orthogonl set of functions Show tht the set {1, cos x, cos 2x,...} is orthogonl on the intervl [, ]: φ 0 = 1, φ n = cos nx (φ 0, φ n ) = = 1 sin nx n φ 0 (x) φ n (x) dx = nd for m n, using the trigonometric identity (φ m, φ n ) = φ m (x) φ n (x) dx = cos nx dx = 1 [sin n sin( n)] = 0 n cos mxcos nx dx = 1 [cos(m + n)x + cos(m n)x] dx 2 = 1 sin(m + n)x sin(m n)x + = 0 2 m + n m n

Exmle 2: Norms Find the norms of the functions given in the Exmle 1 ove. φ 0 2 = dx = 2 φ 0 = 2 φ n 2 = cos 2 nx dx = 1 [1 + cos 2nx] dx = 2 φ n = Any orthogonl set of nonzero functions {φ n (x)}, n = 0, 1, 2... cn e normlized, i.e. mde into n orthonorml set. Exmle: An orthonorml set on the intervl [, ]: 1, cos x cos 2x,,... 2

Vector nlogy Suose v 1, v 2, nd v 3 re three mutully orthogonl nonzero vectors in R 3. Such n orthogonl set cn e used s sis for R 3, tht is, ny three-dimensionl vector cn e written s liner comintion u = c 1 v 1 + c 2 v 2 + c 3 v 3 where c i, i = 1, 2, 3 re sclrs clled the comonents of the vector. Ech comonent cn e exressed in terms of u nd the corresonding vector v i : (u,v 1 ) = c 1 (v 1,v 1 ) + c 2 (v 2,v 1 ) + c 3 (v 3,v 1 ) = c 1 v 1 2 + c 2.0 + c 3.0 (u,v 2 ) = c 2 v 2 2 (u,v 3 ) = c 3 v 3 2

Hence nd c 1 = (u,v 1) v 1 2 c 2 = (u,v 2) v 2 2 c 3 = (u,v 3) v 3 2 u = (u,v 1) v 1 2 v 1 + (u,v 2) v 2 2 v 2 + (u,v 3) v 3 2 v 3 = 3 n=1 (u,v n ) v n 2 v n

Orthogonl series exnsion Suose {φ n (x)} is n infinite orthogonl set of functions on n intervl [, ]. If y = f (x) is function defined on the intervl [, ], we cn determine set of coefficients c n, n = 0, 1, 2,... for which f (x) = c 0 φ 0 (x) + c 1 φ 1 (x) + c 2 φ 2 (x) +... + c n φ n (x) +... (1) using the inner roduct. Multilying the exression ove y φ m (x) nd integrting over the intervl [, ] gives f (x) φ m (x) dx = = c 0 φ 0 (x) φ m (x) dx + c 1 = c 0 (φ 0, φ m ) + c 1 (φ 1, φ m ) +... + c n (φ n, φ m ) +... φ 1 (x) φ m (x) dx +... + c n φ n (x) φ m (x) dx +...

By orthogonlity, ech term on r.h.s. is zero excet when m = n, in which cse we hve f (x) φ n (x) dx = c n φ2 n(x) dx The required coefficients re then In other words f (x) = n=0 c n = c n φ n (x) = n=0 f (x) φ n (x) dx, n = 0, 1, 2... φ2 n(x) dx f (x) φ n (x) dx φ n (x) 2 φ n (x) = n=0 ( f, φ n ) φ n (x) 2 φ n(x)

Definition: Orthogonl set / weight function A set of rel-vlued functions {φ 0 (x), φ 1 (x), φ 2 (x),...} is sid to e orthogonl with resect to weight function w(x) on n intervl [, ] if w(x) φ m (x) φ n (x) dx = 0, m n The usul ssumtion is tht w(x) > 0 on the intervl of orthogonlity [, ]. For exmle, the set {1, cos x, cos 2x,...} is orthogonl w.r.t. the weight function w(x) = 1 on the intervl [, ].

If {φ n (x)} is orthogonl w.r.t. weight function w(x) on the intervl [, ], them multilying the exnsion (1), f (x) = c 0 φ 0 (x) + c 1 φ 1 (x)..., yw(x) nd integrting y rts yields c n = f (x) w(x) φ n (x) dx φ n (x) 2 where φ n (x) 2 = w(x) φ 2 n(x) dx

The series f (x) = n=0 c n φ n (x) (2) with the coefficients given either y c n = f (x) φ n (x) dx φ n (x) 2 or c n = f (x) w(x) φ n (x) dx φ n (x) 2 (3) is sid to e n orthogonl series exnsion of f or generlized Fourier series. Comlete sets We shll ssume tht n orthogonl set {φ n (x)} is comlete. Under this ssumtion f cn not e orthogonl to ech φ n of the orthogonl set.

Fourier series Trigonometric series The set of functions 1, cos x, cos 2 x, cos 3 x,..., sin x, sin 2 x, sin 3 x,... is orthogonl on the intervl [, ]. We cn exnd function f defined on [, ] into the trigonometric series f (x) = 0 2 + n cos n x + n sin n x n=1 (4)

Determining the ocefficients 0, 1, 2,..., 1, 2,... : We multily y 1 (the first function in our orthogonl set) nd integrte oth sides of the exnsion (4) from to f (x) dx = 0 dx + n cos n 2 xdx+ n sin n xdx n=1 Since cos(nx/) nd sin(nx/), n 1, re orthogonl to 1 on the intervl, the r.h.s. reduces s follows Solving for 0 yields f (x) dx = 0 2 0 = 1 dx = 0 2 x = 0 f (x) dx (5)

Now, we multily (4) y cos(mx/) nd integrte f (x) cos m xdx = 0 2 + n=1 By orthogonlity, we hve cos m xdx n cos m x cos n xdx+ n cos m xdx= 0, m > 0 cos m x sin n xdx= 0 cos m x cos n xdx= δ mn where the Kronecker delt δ mn = 0 if m n, nd δ mn = 1 if m = n. cos m x sin n xdx

Thus the eqution (6) ove reduces to f (x) cos n xdx= n nd so n = 1 f (x) cos n xdx (6)

Finlly, multilying (4) y sin(mx/), integrting nd using the orthogonlity reltions we find tht n = 1 sin m xdx= 0, m > 0 sin m x cos n xdx= 0 sin m x sin n xdx= δ mn f (x) sin n xdx (7) The trigonometric series (4) with coefficients 0, n, nd n defined y (5), (6) nd (7), resectively re sid to e the Fourier series of the function f. The coefficients otined from (5), (6) nd (7) re referred s Fourier coefficients of f.

Definition: Fourier series The Fourier series of function f defined on the intervl (, ) is given y f (x) = 0 2 + n cos n x + n sin n x where n=1 0 = 1 n = 1 n = 1 (8) f (x) dx (9) f (x) cos n xdx (10) f (x) sin n xdx (11)

Exmle 1: Exnsion in Fourier series f (x) = With = we hve from (9) nd (10) tht 0 = 1 f (x) dx = 1 0 0 dx + 0 n = 1 f (x) cos nx dx = 1 = 1 sin nx ( x) + 1 n 0 n 0 = 1 cos nx n n 0 cos n + 1 = n 2 = 1 ( 1)n n 2 0, < x < 0 x, 0 x < ( x) dx = 1 x x2 0 2 0 dx + ( x) cos nx dx 0 sin nx dx 0 = 2

Similrly, we find from (11) n = 1 n 0 ( x) sin nx dx = 1 n The function f (x) is thus exnded s f (x) = 1 ( 1) n 4 + n=1 n 2 cos nx + 1 n sin nx (12) We lso note tht 1 ( 1) n = 0, n even 2, n odd.

Convergence of Fourier series Theorem: Conditions for convergence Let f nd f e iecewise continuous on the intervl (, ); tht is, let f nd f e continuous excet t finite numer of oints in the intervl nd hve only finite discontinuities t these oints. Then the Fourier series of f on the intervl converges to f (x) t oint of continuity. At oint of discontinuity, the Fourier series converges to the verge f (x+) + f (x ), 2 where f (x+) nd f (x ) denote the limit of f t x from the right nd from the left, resectively.

Exmle 2: Convergence of oint of discontinuity The exnsion (12) of the function (Exmle 1) 0, < x < 0 f (x) = x, 0 x < will converge to f (x) for every x from the intervl (, ) excet t x = 0 where it will converge to f (0+) + f (0 ) 2 = + 0 2 = 2.

Periodic extension Oserve tht ech of the functions in the sis set 1, cos x, cos 2 x, cos 3 x,..., sin x, sin 2 x, sin 3 x,... hs different fundmentl eriod 2/n, n 1, ut since ositive integer multile of eriod is lso eriod, we see tht ll the functions hve in common the eriod 2. Thus the r.h.s. of f (x) = 0 2 + n cos n x + n sin n x n=1 is 2-eriodic; indeed 2 is the fundmentl eriod of the sum. We conclude tht Fourier series not only reresents the function on the intervl (, ) ut lso gives the eriodic extension of f outside this intervl.

We cn now ly the Theorem on conditions for convergence to the eriodic extension or simly ssume the function is eriodic, f (x + T) = f (x), with eriod T = 2 from the outset. When f is iecewise continuous nd the right- nd left-hnd derivtives exist t x = nd x =, resectively, then the Fourier series converges to the verge [ f ( )+ f (+)]/2 t these oints nd lso to this vlue extended eriodiclly to ±3, ±5, ±7, nd so on. Exmle: The Fourier series of the function f (x) in the Exmle 1 converges to the eriodic extension of the function on the entire x-xis. At 0, ±2, ±4,..., nd t ±, ±3, ±5,..., the series converges to the vlues f (0+) + f (0 ) 2 = 2 nd f (+) + f ( ) 2 = 0

Sequence of rtil sums It is interesting to see how the sequence of rtil sums {S N (x)} of Fourier series roximtes function. For exmle S 1 (x) = 4, S 2(x) = 4 + 2 cos x + sin x, S 3(x) = 4 + 2 cos x + sin x + 1 sin 2x 2