MAE/MSE 502, Spring 2016, HW3 Solutionshhuang38/mae502_2016s_hw3_sol.pdf · (c) Fourier Cosine...

5

Click here to load reader

Transcript of MAE/MSE 502, Spring 2016, HW3 Solutionshhuang38/mae502_2016s_hw3_sol.pdf · (c) Fourier Cosine...

Page 1: MAE/MSE 502, Spring 2016, HW3 Solutionshhuang38/mae502_2016s_hw3_sol.pdf · (c) Fourier Cosine expansion is better, although it still has Gibbs ripples in the vicinity of x = 0.5

MAE/MSE 502, Spring 2016, HW3 Solutions

Prob 1

u (x ,t ) = 5 t + 3 + t cos(x) +1

√12sinh(√12 t )sin(2 x)

Prob 2

u (x ,t ) = ∑n=−∞

Cn(0) e−n2Dt ein (x−n2E t)

,

where

Cn(0) =1

2π∫0

e[1−cos(x)]2 e−i n x d x .

Plot:

Page 2: MAE/MSE 502, Spring 2016, HW3 Solutionshhuang38/mae502_2016s_hw3_sol.pdf · (c) Fourier Cosine expansion is better, although it still has Gibbs ripples in the vicinity of x = 0.5
Page 3: MAE/MSE 502, Spring 2016, HW3 Solutionshhuang38/mae502_2016s_hw3_sol.pdf · (c) Fourier Cosine expansion is better, although it still has Gibbs ripples in the vicinity of x = 0.5

Prob 3

(a)

(b)

Page 4: MAE/MSE 502, Spring 2016, HW3 Solutionshhuang38/mae502_2016s_hw3_sol.pdf · (c) Fourier Cosine expansion is better, although it still has Gibbs ripples in the vicinity of x = 0.5

(c)

Fourier Cosine expansion is better, although it still has Gibbs ripples in the vicinity of x = 0.5 due to the discontinuity of f(x) there. Fourier Sine expansion has Gibbs ripples in the vicinities of both x = 0.5 and x = 0, because the odd extension of the original f(x) has a discontinuity at x = 0.

Page 5: MAE/MSE 502, Spring 2016, HW3 Solutionshhuang38/mae502_2016s_hw3_sol.pdf · (c) Fourier Cosine expansion is better, although it still has Gibbs ripples in the vicinity of x = 0.5

Prob 4