University of Florida/Egm4313/s12.teamboss/R7

From testwiki
Revision as of 02:22, 17 October 2020 by imported>MaintenanceBot (Category:Pages with broken file links)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to navigation Jump to search

Template:Big3

File:Toc.JPG

Problem R7.1: Determining Orthogonality

Statement

Determine orthogonality of:

ϕi,ϕj=0  for  ij
ϕj,ϕj=L2  for  i=j

(1.0)


Where: ϕi(x)=sin(ωix), ϕj(x)=sin(ωjx), ωi=iπL, and ωj=jπL.

Solution

The definition of a scalar product is:

f¯,g¯:=0Lf¯(x) g¯(x)dx

(1.1)

So, substituting the problem statement values into Eq. (1.0):

ϕi,ϕj:=0Lsin(ωix) sin(ωjx)dx

(1.2)


and:

ϕj,ϕj:=0Lsin(ωjx) sin(ωjx)dx

(1.3)

Solving Eq. 1.2 first:

ϕi,ϕj=0Lsin(ωix) sin(ωjx)dx
=120Lcos((ωiωj)x)cos((ωi+ωj)x)dx
=12(ωiωj)[sin((ωiωj)L)sin(0)]12(ωi+ωj)[sin((ωi+ωj)L)sin(0)]
=12L(iπjπ)sin(iπjπ)12L(iπ+jπ)sin(iπ+jπ)


Because i and j are both integers, the sin function will always be a periodic integer varying by π, meaning that sin(iπjπ)=0 and sin(iπ+jπ)=0. Thus proving:

ϕi,ϕj=0  for  ij

(1.0.a)



Now solving for Eq. 1.3:

ϕj,ϕj=0Lsin(ωjx) sin(ωjx)dx
=120Lcos((ωjωj)x)cos((ωj+ωj)x)dx
=120Lcos((jπjπ)xL)cos((jπ+jπ)xL)dx
=120Lcos(0)cos(2jxπL)dx
=120L1cos(2jxπL)dx
=L2L2jπsin(2jπ)

Again, because j is an integer, the sin function will always be a periodic integer varying by 2jπ, meaning that sin(2jπ)=0. Thus proving:

ϕj,ϕj=L2  for  i=j

(1.0.b)

Author

Solved and Typed By ---Egm4313.s12.team1.durrance (talk) 16:49, 24 April 2012 (UTC)
Reviewed By - Egm4313.s12.team1.armanious (talk) 23:48, 24 April 2012 (UTC)


Problem R7.2: Truncated Series

Statement

Plot the truncated-series (3) p19-12 with n=5, and for:

t=αp1=α2πcω1=α2Lc

(2.0)

α=.5,1,1.5,2

(2.1)

Solution

Let it be stated first that (3) p.19-12 is equal to:

u(x,t)=j=1najcos(cωjt)sin(ωjx)

(2.2)

And in this series representation, c=3, L=2, and aj. is equaled to this:

aj={0forj=2m(even)4/(π3j3)forj=2m+1(odd)

(2.3)

With this in mind, we are prepared to create these truncated-series.
For α=.5

u(x,t)=[a1cos(cω1α2Lc)sin(ω1x)]

(2.4)

+[a2cos(cω2α2Lc)sin(ω2x)]

+[a3cos(cω3α2Lc)sin(ω3x)]

+[a4cos(cω4α2Lc)sin(ω4x)]

+[a5cos(cω5α2Lc)sin(ω5x)]

These simplify below as shown, after plugging in 3 for c, 10, for L, .5 for alpha, and the value of aj as determined from (2.2).

u(x,t)=[4π3cos(ω1(.5)(2*10))sin(ω1x)]

(2.5)

+[4π3*33cos(ω3(.5)(2*10))sin(ω3x)]
+[4π3*53cos(ω5(.5)(2*10))sin(ω5x)]

Making ω equal to 2π, the graph appears as shown below.

Alpha = .5
For α=1

These simplify below as shown, after plugging in 3 for c, 10, for L, 1 for alpha, and the value of aj as determined from (2.2).

u(x,t)=[4π3cos(ω1(1)(2*10))sin(ω1x)]

(2.6)

+[4π3*33cos(ω3(1)(2*10))sin(ω3x)]
+[4π3*53cos(ω5(1)(2*10))sin(ω5x)]

Making ω equal to 2π, the graph appears as shown below.

Alpha = 1

For α=1.5

These simplify below as shown, after plugging in 3 for c, 10, for L, 1.5 for alpha, and the value of aj as determined from (2.2).

u(x,t)=[4π3cos(ω1(1.5)(2*10))sin(ω1x)]

(2.7)

+[4π3*33cos(ω3(1.5)(2*10))sin(ω3x)]
+[4π3*53cos(ω5(1.5)(2*10))sin(ω5x)]

Making ω equal to 2π, the graph appears as shown below.
Alpha = 1.5
For α=2

These simplify below as shown, after plugging in 3 for c, 10, for L, 2 for alpha, and the value of aj as determined from (2.2).

u(x,t)=[4π3cos(ω1(2)(2*10))sin(ω1x)]

(2.8)

+[4π3*33cos(ω3(2)(2*10))sin(ω3x)]
+[4π3*53cos(ω5(2)(2*10))sin(ω5x)]

Making ω equal to 2π, the graph appears as shown below.

Alpha = 2

Author

Solved and Typed By - Egm4313.s12.team1.silvestri (talk) 19:01, 23 April 2012 (UTC)

Reviewed By - Egm4313.s12.team1.armanious (talk) 19:40, 25 April 2012 (UTC)


Problem R7.3 Scalar Product, Magnitude, and Angles of Functions

Statement

Find the scalar product <f,g>, the magnitude of f and g, and the angle between f and g for (see R7.3 Lect. 11-1 pg. 8):

1) f(x)=cosx,g(x)=x,for2x10
2) f(x)=12(3x21),g(x)=12(5x33x),for1x+1

Solution

The scalar product, the function equivalent of the vector dot product, of two functions can be found in the following way:

<f,g>:=abf(x)g(x)dx

(3.0)

The magnitude of a function is defined as the square root of the scalar product between the function and itself:

||f||:=<f,f>1/2=[abf2(x)dx]1/2

(3.1)

The cosine of the angle between two functions can be found in the following way:

cosθ=<f,g>||f||||g||

(3.2)

Note that these equations are very similar to their vector counterparts.
Part 1

f(x)=cosx,g(x)=x,for2x10

(3.3)

The scalar product of f and g can be found using (3.0):

<f,g>=210xcos(x)dx=[xsin(x)+cos(x)]210

(3.4)

Solving this yields the scalar product:

(10sin10+cos10)(2sin(2)+cos(2))=7.68

(3.5)

To find the magnitude of f, the scalar product of f with itself must first be found:

<f,f>=210cos2(x)dx=210(12cos(2x)+12)dx=[14sin(2x)+12x]210

(3.6)

Solving this yields:

[14sin(2x)+12x]210=(14sin(20)+5)(14sin(4)1)=6.04

(3.7)

Therefore the magnitude of f is found to be:

||f||=<f,f>1/2=6.04=2.46

(3.8)

A similar approach can be taken for the magnitude of g:

<g,g>=210x2dx=[x33]210

(3.9)

[x33]210=1000383=10083=336

(3.10)

||g||=<g,g>1/2=336=18.33

(3.11)

Using the above information with (3.2), the cosine of the angle between f and g is found to be

cosθ=<f,g>||f||||g||=7.68(2.46)(18.33)=0.1705

(3.12)

Therefore, the angle between the functions f and g is

θ=cos1(0.1705)=1.74rad=99.8

(3.13)

Part 2

f(x)=12(3x21),g(x)=12(5x33x),for1x+1

(3.14)

The scalar product of f and g can be found using (3.0):

<f,g>=1114(3x21)(5x33x)dx=1411(15x514x3+3x)dx=14[156x6144x4+32x2]11

(3.15)

Solving this yields the scalar product:

14[52x672x4+32x2]11=14(1212)=0

(3.16)

To find the magnitude of f, the scalar product of f with itself must first be found:

<f,f>=11(12(3x21))2dx=1411(9x46x2+1)dx=14[95x52x3+x]11

(3.17)

Solving this yields:

14[(952+1)(95+21)]=14(45+45)=25

(3.18)

Therefore the magnitude of f is found to be:

||f||=<f,f>1/2=0.4=0.632

(3.19)

A similar approach can be taken for the magnitude of g:

<g,g>=11(12(5x33x))2dx=1411(25x615x4+9x2)dx=14[257x73x5+3x3]11

(3.20)

12(2573+3)=2514

(3.21)

||g||=<g,g>1/2=2514=1.336

(3.22)

Using the above information with (3.2), the cosine of the angle between f and g is found to be

cosθ=<f,g>||f||||g||=0(0.632)(1.336)=0

(3.23)

Because the scalar product of the two functions is zero, the two functions are orthogonal. That is:

θ=π2rad=90

(3.24)

Author

Solved and Typed By - Egm4313.s12.team1.armanious (talk) 23:48, 21 April 2012 (UTC)

Reviewed By - --Egm4313.s12.team1.stewart (talk) 02:22, 25 April 2012 (UTC)


Problem R7.4 The Fourier Series and Plotted Functions

Statement

From Lecture 11 Pg. 8 Do K 2011 p.482 pb.6,9,12,13. Problems 6 and 12 ask to graph the function and the Fourier series of the function below (4.1).

f(x)=|x|

(4.1)

Problems 9 and 13 ask to graph the function and Fourier series of the function below (4.2).

f(x)={xifπ<x<0πxif0<x<π}

(4.2)

Solution

For Question number 6, plotted using Wolfram Alpha. The Function 4.1
For Question number 9, plotted using Matlab. The Function 4.2

Question 12 asks for the Fourier series of equation 4.1, f(x)=|x| which is an even function because f(x)=+f(x), and because it is even bn=0. So the Fourier looks like this:

f(x)=a0+n=1ancosnπxL

(4.3)

In which L=π and:

a0=12ππ0xdx+12π0πxdx

(4.4)

an=1πππf(x)cos(nx)dx

(4.5)

bn=1πππf(x)sin(nx)dx=0

(4.6)

These simplify to:

a0=12ππ22+12ππ22=π2

(4.7)

an=1ππ0xcos(nx)dx+1ππ0xcos(nx)dx

(4.8)

an=4n2π

(4.9)

f(x)=π2+n14n2πcosnπxL

(4.10)

The final answer:

f(x)=π24π(cosx+cos3x9+cos5x25+...)

(4.11)

Plotted with Wolfram Aplpha Fourier Series of Function 4.1

Now the Fourier series of equation 4.2, we need to divide the function into two parts. First let's evaluate f(x)=xifπ<x<0, which is an odd function because f(x)=f(x). For an odd function a0=an=0 so:

f(x)=n=1[bnsin(nπxL)]=n=1[bnsin(nx)]

(4.12)

bn=1ππ0xsin(nx)=1π[sin(nx)nxcos(nx)n2]=cos(πn)n

(4.13)

Which simplifies to:

bn=[1,12,13,14,15,16...]

(4.14)

The Fourier Series for just the top part of equation 4.2

f(x)=sinxsin2x2+sin3x3sin4x4+sin5x5...

(4.15)

Now the bottom part of equation 4.2 is neither even nor odd so we must add all the components and at the end combine the solution with the Fourier series for the top part of equation 4.2

f(x)=a0+n=1ancos(nπxL)+n=1bnsin(nπxL)=a0+n=1ancos(nx)+n=1bnsin(nx)

(4.16)

a0=12π0π(πx)dx=12π[πxx22]=π4

(4.17)

an=1π0π(πx)cos(nx)dx=1π[n(πx)sin(nx)cos(nx)n2]=1cos(πn)πn2

(4.18)

an=[2π,0,26π,0,225π,0,249π...]

(4.19)

bn=1π0π(πx)sin(nx)dx=1π[sin(nx)+n(πx)cos(nx)n2]=1n

(4.20)

bn=[1,12,13,14,15,16...]

(4.21)

The Fourier Series of f(x)=πxif0<x<π.

f(x)=a0+n=1ancos(nx)+n=1bnsin(nx)=π4+[2cosxπ+2cos3x6π+2cos5x25π...]+[1+sinx2+sin2x3+sin3x4+sin4x5+sin5x6...]

(4.22)

Combining with the Fourier Series of f(x)=xifπ<x<0:

f(x)=π4+[2cosxπ+2cos3x6π+2cos5x25π...]+[sinx+sin2x2+sin3x3+sin4x4+sin5x5+...]+[sinxsin2x2+sin3x3sin4x4+sin5x5...]

(4.23)

The total series and final answer:

f(x)=π4+[2cosxπ+2cos3x6π+2cos5x25π...]+[2sinx+2sin3x3+2sin5x5...]

(4.24)

Plotted with Wolfram Alpha Fourier Series of Function 4.2

Author

Solved and Typed By - --Egm4313.s12.team1.stewart (talk) 02:22, 25 April 2012 (UTC)

Reviewed By - Egm4313.s12.team1.essenwein (talk) 16:50, 25 April 2012 (UTC)


Problem R7.5 Approximation by Trigonometric Series

Statement

Consider equation (5.0) below:

<ϕ2j1,ϕ2k1>=0Psinjωxsinkωxdx=0

(5.0)

1) Find the exact integration of equation (5.0) with P=2π, j=2, k=3.
2) Confirm the result with Matlab's trapz command for the trapezoidal rule.

Solution

1) Substituting the given values, we get equation (5.1) below:

<ϕ3,ϕ5>=02πsin2ωxsin3ωxdx=0

(5.1)

Using the product and sum trigonometry identities, we can evaluate the definite integral:

02πsin2ωxsin3ωxdx=1202π(cosωxcos5ωx)dx=0

(5.2)

[12ωsinωx15sin5ωx]02π=0

(5.3)

Also, ω=2πP=1, and substituting this into equation (5.3) allows the definite integral to be solved:

[12sinx15sin5x]02π=0

(5.4)

2) The following is the code for evaluating the definite integral in equation (5.1) in Matlab:

EDU>> X = 0:pi/100:2*pi;
EDU>> Y = sin(2*X).*sin(3*X);
EDU>> Z = trapz(X,Y)

Z =

 -2.2633e-17

EDU>> plot(X,Y)
File:Graph7.5.png

The returned value of the evaluated integral in Matlab is essentially equal to zero, and it can be seen in the graph that the areas under the curve cancel each other out.

Author

Solved and Typed By - Egm4313.s12.team1.wyattling

Reviewed By - Egm4313.s12.team1.rosenberg


Contributing Members

Team Contribution Table
Problem Number Lecture Assigned To Solved By Typed By Proofread By
7.1 R7.1 Lect. 19c-1 pg. 10 Jesse Durrance Jesse Durrance Jesse Durrance George Armanious
7.2 R7.2 Lect. 19c-1 pg. 13 Emotion Silvestri Emotion Silvestri Emotion Silvestri George Armanious
7.3 R7.3 Lect. 11-1 pg. 8 George Armanious George Armanious George Armanious Chris Stewart
7.4 R7.4 Lect. 11-1 pg.8 Chris Stewart Chris Stewart Chris Stewart Eric Essenwein
7.5 R7.5 Lect. 12-1 pg. 8 Wyatt Ling Wyatt Ling Wyatt Ling Steven Rosenberg