Docsity
Docsity

Prepare for your exams
Prepare for your exams

Study with the several resources on Docsity


Earn points to download
Earn points to download

Earn points by helping other students or get them with a premium plan


Guidelines and tips
Guidelines and tips

Fourier Series - Seminar in Engineering Analysis - Lecture Slides, Slides of Engineering Analysis

Course title is Seminar in Engineering Analysis. Analytic and numerical methods applied to the solution of engineering problems at an advanced level. Solution methods are demonstrated on a wide range of engineering topics, including structures, fluids, thermal, thermal energy transport, and mechanical systems.

Typology: Slides

2013/2014

Uploaded on 02/01/2014

sashie
sashie 🇮🇳

4.2

(40)

185 documents

1 / 5

Toggle sidebar

This page cannot be seen from the preview

Don't miss anything!

bg1
Fourier Series January 26, 2009
ME 501B – Engineering Analysis 1
Fourier Series
Fourier Series
Larry Caretto
Mechanical Engineering 501B
Seminar in Engineering Analysis
January 26, 2009
2
Outline
Review last class
Fourier series as expansions in periodic
functions
Comparison to eigenfunction expansions
Odd and even functions
Periodic extensions of non-periodic
functions
Complex Fourier series
3
Review Last Lecture
Discussed Sturm-Liouville Problem
Solutions are a set of orthogonal
eigenfunctions, ym(x) that can be used to
express other functions, f(x)
=
=
0
)()(
mmm xyaxf
()
== b
a
mm
b
a
m
mm
m
m
dxxyxyxp
dxxfxyxp
yy
fy
a
)()()(
)()()(
,
),(
p(x) is
weight
function
Have to
compute
am
4
Review Orthogonal Functions
Defined in terms of inner product
Norm of two like eigenfunctions ||yi||
()
iji
b
a
jiji ydxxpxyxyyy
δ
2
*)()()(, ==
()
ij
b
a
jiji dxxpxfxfff
δ
== )()()(, *
Convert orthogonal eigenfunctions to
orthonormal eigenfunctions
Orthonormal eigenfunctions
i
i
iy
y
f=
5
Review Sturm-Liouville
General equation whose solutions provide
orthogonal eigenfunctions
Defined for a x b with p(x), q(x) and r(x)
continuous and p(x) > 0
(Homogenous) differential equation and
boundary conditions shown below
[]
0)()(
)(
=+
+
yxpxq
dx
dy
xr
dx
d
λ
0)(
0)(
21
21
=+
=+
=
=
bx
ax
dx
dy
by
dx
dy
kayk
ll
6
Review Sturm-Liouville Results
Eigenvalues are real
Eigenfunctions defined over a region a
x b form an orthogonal set over that
region.
Eigenfunctions form a complete set over
an infinite-dimensional vector space
We can expand any function over the
region in which the Sturm-Liouville
problem is defined in terms of the
eigenfunctions for that problem
pf3
pf4
pf5

Partial preview of the text

Download Fourier Series - Seminar in Engineering Analysis - Lecture Slides and more Slides Engineering Analysis in PDF only on Docsity!

Fourier SeriesFourier Series

Larry Caretto

Mechanical Engineering 501B

Seminar in Engineering Analysis

January 26, 2009

2

Outline

  • Review last class
  • Fourier series as expansions in periodic functions - Comparison to eigenfunction expansions
  • Odd and even functions
  • Periodic extensions of non-periodic functions
  • Complex Fourier series

3

Review Last Lecture

  • Discussed Sturm-Liouville Problem
  • Solutions are a set of orthogonal eigenfunctions, ym(x) that can be used to express other functions, f(x)

=

= 0

( ) () m

fx am ymx

∫ = = b

a

m m

b

a

m

m m

m m pxy xy xdx

pxy xfxdx

y y

y f a

() () ( )

() () ()

,

( , )

  • p(x) is

weight function

  • Have to

compute a (^) m 4

Review Orthogonal Functions

  • Defined in terms of inner product
  • Norm of two like eigenfunctions ||yi ||

( ) (^) i ij

b

a

yi yj yi xyjxpxdx y δ

  • 2

( ) (^) ij

b

a

fi , fj = ∫ fi ( x ) fj ( x ) p ( x ) dx = δ

  • Convert orthogonal eigenfunctions to orthonormal eigenfunctions
  • Orthonormal eigenfunctions

i

i i (^) y

y f =

5

Review Sturm-Liouville

General equation whose solutions provide orthogonal eigenfunctions

  • Defined for a ≤ x ≤ b with p(x), q(x) and r(x) continuous and p(x) > 0
  • (Homogenous) differential equation and boundary conditions shown below

[ ( ) ( )] 0

qx px y

dx

dy

rx

dx

d

λ ( )^0

1 2

1 2

=

=

x b

xa

dx

dy

yb

dx

dy

kya k

l l

6

Review Sturm-Liouville Results

  • Eigenvalues are real
  • Eigenfunctions defined over a region a ≤ x ≤ b form an orthogonal set over that region.
  • Eigenfunctions form a complete set over an infinite-dimensional vector space
  • We can expand any function over the region in which the Sturm-Liouville problem is defined in terms of the eigenfunctions for that problem

7

Review Eigenfunction Expansions

  • Eigenfunction expansion formula

=

0

m

f x am ym x

∫ = = b

a

m m

b

a

m

m m

m m pxy xy xdx

pxy x f xdx

y y

y f a

( ) ( ) ( )

  • Equation for a (^) m coefficients in eigenfunction expansion of f(x)

8

Review Expansion of f(x) = x

  • Start with general equation for a (^) m
  • Use y (^) m = sin(mπx) over 0 ≤ x ≤ 1 which is a Sturm-Liouville solution
  • Weight function p(x) = 1

π

= π

π

= =

m

x mxdx

mxdx

x mxdx

y y

y f a

m

m m

m m

1

1

0 1

0

2

1

0 2 (^1 )

2

1

sin( )

sin( )

sin( )

,

( , )

⎥ ⎦

⎤ ⎢ ⎣

⎡ = = − + − +L 4

sin( 4 ) 3

sin( 3 ) 2

sin( 2 ) 1

2 sin( ) ( )

x x x x fx x

π π π π π

9

Review Partial Sums – Small

0

0.

0.

0.

0.

0.

0.

0.

0.

0.

1

1.

0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 x

Series sum

Exact 1 term 2 terms 3 terms 5 terms 10 terms

10

Review Partial Sums – Large

0

0.

0.

0.

0.

0.

0.

0.

0.

0.

1

1.

1.

0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 x

Series sums

Exact 10 terms 25 terms 50 terms 100 terms

11

Fourier Series

  • Have same basic idea as eigenfunction

expansions

  • Represent other functions, f(x), as a series of sines and cosines
  • Compute coefficients in a similar way to eigenfunction expansions
  • Fourier series based on periodicity of trigonometric functions

12

Fourier Series

  • Equations for series and coefficients

defined for –L < x < L

=

1

( ) 0 cos sin

n

n n

L

n x

b

L

n x

f x a a

∫ −

L

L

f xdx

L

a ( )

0

∫ −

L

L

n dx

L

n x

f x

L

b

( ) sin

∫ −

L

L

n dx

L

n x

f x

L

a

( ) cos

19

Even Function Cosine Series

  • f(x) is an even function f(–x) = f(x)
  • Cosine only series defined for –L ≤ x ≤ L
  • a (^) n evaluated by integral from 0 to L

=

1

( ) 0 cos

n

n

L

n x

f x a a

π

∫ ∫

L L

L

f xdx

L

f xdx

L

a

0

∫ ∫ ⎟ ⎠

L L

L

n dx L

nx f x L

dx L

n x f x L

a 0

( )cos

( )cos

20

Half-interval series

  • Based on equations for sine series for odd functions and cosine series for even functions
  • Either the sine or cosine series can apply to any function - Can use sines to expand even functions - Can use cosines to expand odd functions
  • Each series defined for 0 ≤ x ≤ L
  • Behavior outside region 0 ≤ x ≤ L depends on the function

21

Half-interval Series Example

  • Last class computed series for f(x) = x in terms of sin(nπx/L) with L =

1 for 0 ≤ x ≤ L

  • Would get same coefficients from

equations for Fourier sine series

  • Get correct result for –L ≤ x ≤ L with periodic extensions

=

⎟ ⎠

⎞ ⎜ ⎝

1

( ) sin n

n L

nx f x b

π

π

− ⎟ = ⎠

⎞ ⎜ ⎝

⎛ π ⎟ = ⎠

⎞ ⎜ ⎝

⎛ π

n

L dx L

nx x L

dx L

nx f x L

b

L L n n

1

0 0

2 1 sin

2 ()sin

2

Fourier Sine Series for f(x) = x

-1.

-1.

-0.

-0.

-0.

-0.

-0.

-0.

-0.

-0.

-0.

0.

0.

0.

0.

0.

0.

0.

0.

0.

0.

1.

1.

-1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2. x

f(x)

Exact 1 term 2 terms 5 terms 10 terms

23

Cosine Series for f(x) = x

  • Set f(x) = x in equations for a 0 and an

0

2

0 0

0

x L

L

xdx L

f xdx L

a

L L L

⎥^ = ⎦

∫ ∫

L

L L

n

L

nx

n

x

L

n x

n

L

L

dx

L

n x

x

L

dx

L

n x

f x

L

a

0

2

0 0

cos sin

cos

( )cos

∫ ∫

24

Cosine Series for f(x) = x

  • Result for a (^) n is zero for even n values

[ ( ) ] [ ( )

]

2 cos 1 4 ( 0 )sin( 0 )

cos cos( 0 ) sin

cos sin

2 2

2

0

2

odd n n

L

n

L n

L n n

L

n n

L
L
L

n x

n

x

L

n x

n

L
L

a

L

n

25

Cosine Series for f(x) = x

  • Final series, shown below, uses n = 2m+1 to get odd values of n only

∑ ∑

=

=

⎟ ⎠

⎞ ⎜ ⎝

⎛ ⎟= − ⎠

⎞ ⎜ ⎝

⎛ = + 1 , 3 , 5 ,K

2 2 1

0 cos

4 1

2

( ) cos n n

n L

nx

n

L L

L

nx fx a a

π

π

π

( )

=

0

2 2

cos 2 1

m L

m x

m

L L

f x

= − L
L

x

L

x

L

L L π x π π

cos 25

cos 9

cos

2

Cosine series for f(x) = x

0.

0.

0.

0.

0.

0.

0.

0.

0.

0.

1.

-1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2. x

f(x)

Exact 1 term 2 terms 5 terms 10 terms

27

Complex Fourier Series

  • Euler relationship e ix = cos(x) + i sin(x)
  • Setting x = -x gives e -ix = cos(-x) + i sin(-x) = cos(x) - i sin(x)
  • Get result for f(x) and coefficient cn

=−∞

n

inx f ( x ) cn e

L

L

inx

n f xe

L

c ()

28

Summary

  • Fourier series is alternative approach to Sturm-Liouville for developing series expansions in sines and cosines
  • Based on periodic functions
  • Can use half-interval expansions to

show application to any function

  • Basic approach similar to eigenfunction

expansions: find the coefficients