FOURIER SERIES PART III: APPLICATIONS
|
|
- Merryl Joseph
- 6 years ago
- Views:
Transcription
1 FOURIER SERIES PART III: APPLICATIONS We extend the construction of Fourier series to functions with arbitrary eriods, then we associate to functions defined on an interval [, L] Fourier sine and Fourier cosine series and then aly these results to solve BVPs. 1. Fourier series with arbitrary eriods Let f : R R be a iecewise continuous function with eriod 2 ( > ). We would like to reresent f by a trigonometric series. We can reeat what we did in the revious note, when we had = π, and reach the sought reresentation. There is however a simle way of obtaining the same series by introducing the function g(s) = f ( s π ) ( f(x) = g ( πx )). Note that since f C(R), then g C(R) and that since f is 2-eriodic, then ( ) (s + 2π) ( s ) ( s ) g(s + 2π) = f = f π π + 2 = f = g(s). π That is, g is 2π-eriodic. We can therefore associate a Fourier series to g: g(s) a 2 + a n cos ns + b n sin ns. In terms of the function f, we have the association f(x) a 2 + a n cos nπx + b n sin nπx. π The coefficients are given by a = 1 π g(s)ds = 1 π ( s ) f ds = 1 f(x)dx π π π π π a n = 1 π g(s) cos nsds = 1 π ( s ) f cos(ns)ds = 1 π π π π π b n = 1 π g(s) sin nsds = 1 π ( s ) f sin(ns)ds = 1 π π π π ( ) nπx f(x) cos dx ( ) nπx f(x) sin dx The fundamental convergence theorem (Fourier theorem) states Theorem. Let f be a 2-eriodic and iecewise smooth function on R. Then f av (x) = f(x+ ) + f(x ) = a a n cos nπx + b n sin nπx, where the coefficients are given by a = 1 f(x)dx, Date: March 2,
2 2 FOURIER SERIES PART III: APPLICATIONS and for n Z +, a n = 1 ( ) nπx f(x) cos dx, b n = 1 ( ) nπx f(x) sin dx. Again if f is continuous at x, then f(x ) is equal to its Fourier series: f(x ) = a 2 + a n cos nπx + b n sin nπx We also have uniform convergence with the additional condition that f is continuous on R. More recisely, Theorem. Let f be a 2-eriodic and iecewise smooth function on R. Suose that f is continuous in R, then f(x) = a 2 + a n cos nπx and the convergence is uniform. + b n sin nπx. x R Analogous results hold for termwise differentiation of Fourier series and termwise integration. Examle. Consider the function f(x) that is 2-eriodic and is given on the interval [, ] by 2Hx + H if x < /2 ; f(x) = 2Hx + H if /2 x < ; if /2 x where H is a ositive constant. H Figure 1. Grah of function of examle The function f is continuous on R and is iecewise smooth. It is also an even function (hence its b n Fourier coefficients are all zero). The Fourier coefficients of f are a = 2 f(x)dx = 2H /2 ( ) 2x + 1 dx = H 2
3 FOURIER SERIES PART III: APPLICATIONS 3 and for n = 1, 2, 3, we have a n = 2 f(x) cos nπx dx = 2H /2 ( ) 2x + 1 = 2H [( ) ] /2 2x nπx + 1 sin + 4H /2 nπ nπ = 4H (1 π 2 n 2 cos nπ ) 2 Since the function f is continuous and iecewise smooth, we have f(x) = H 4 + 4H 1 cos(nπ/2) π 2 n 2 cos nπx x R cos nπx dx sin nπx dx Furthermore, the convergence is uniform. We can aly termwise differentiation to obtain the Fourier series of the derivative f (where f (x) = 2H/ for < x < /2, f (x) = 2H/ for /2 < x < and f (x) = for /2 < x < ). We have f (x) 4H π 1 cos(nπ/2) n 2. Parseval s Identity sin nπx Parseval s identity is a sort of a generalized ythagorean theorem in the sace of functions. Theorem. (Parseval s Identity) Let f be 2-eriodic and iecewise continuous with Fourier series a 2 + a n cos nπx + b n sin nπx. Then a (a 2 n + b 2 n) = 1 2 f(x) 2 dx. Proof. We will rove the identity when f is continuous and iecewise smooth. In this case the Fourier series (equals f),is uniformly convergent, and termwise integration is allowed. We have f(x) 2 = ( a 2 + a n cos nπx = a 2 f(x) + a n f(x) cos nπx ) + b n sin nπx f(x) + b nf(x) sin nπx We integrate from to and divide by 2. The term by term integration gives 1 2 f(x) 2 dx = a 1 f(x)dx a n f(x) cos nπx 2 dx + b 1 n 2 = a (a 2 n + b 2 2 n). f(x) sin nπx dx.
4 4 FOURIER SERIES PART III: APPLICATIONS The Parseval s identity is used to aroximate the average error when relacing a given function f by its N-th Fourier artial sum S N f. Recall that S N f(x) = a N 2 + a n cos nπx + b n sin nπx The mean square error, when relacing f by S N f, is defined as the number E N given by EN 2 = 1 (f(x) S N f(x)) 2 dx. If we use Parseval s identity to the function f S N f, we find that EN 2 = (a 2 n + b 2 n) n=n+1 Examle. Consider the 2π-eriodic function defined over [ π, π] by f(x) = 1 for < x < π and f(x) = 1 for π < x <. The fourier series of f is 4 sin(2j + 1)x. We would like to find N so that the aroximation π (2j + 1) j= f(x) 4 π 2j+1 N sin(2j + 1)x (2j + 1) guarantees that the mean square error is no more than.1. That is E N < 1 2. From the above discussion, we have EN 2 = 16 1 π 2 (2j + 1) 2. j>(n 1)/2 We can estimate the last series by using the integral test (see figure). We have. y=1/(2x+1) 2 M 1 M M+1 j=m Figure 2. Comarison of integral and sum 1 (2j + 1) 2 dx M 1 (2x + 1) 2 = 1 2(2M 1).
5 FOURIER SERIES PART III: APPLICATIONS 5 Hence, it follows from the above calculations that E 2 N 16 π 2 N. Therefore, in order to have E N <.1, it is enough to take N so that N > 16/π 2. That is N = 163. Parseval s identity can also be used to evaluate series. Examle. We have seen that Parseval s identity gives From this we get x = π 2 4 π π π 2 j= j= cos(2j + 1)x (2j + 1) 2 x [ π, π]. 1 (2j + 1) 4 = 1 π x 2 dx = π2 2π π = π Even and Odd Periodic Extensions We would like to reresent a function f given only on a an interval [, L] by a trigonometric series. For this we extend f to the interval [ L, L] as either an even function or as an odd function then extend it to R as a eriodic function with eriod 2L. The Fourier series of this extension gives the sought reresentation of f. The even extension gives the Fourier cosine series of f and the odd extension gives the Fourier sine series of f. More recisely, let f be a iecewise smooth function on the interval [, L]. Let f even and f odd be, resectively, the even and the odd odd extensions of f to the interval [ L, L]. Now we extend f even to R as a 2L-eriodic function F even and grah of f odd grah of f grah of f even L L L L L we extend f odd to R as a 2L-eriodic function F odd. Note f(x) = F even (x) = F odd (x) x [, L]. The Fourier series of F even and F odd are F even (x) a 2 + a n cos nπx L and F odd (x) b n sin nπx L
6 6 FOURIER SERIES PART III: APPLICATIONS Even eriodic extension F even L L Odd eriodic extension F odd L L The Fourier coefficients are a n b n = 2 L = 2 L L L F even (x) cos nπx L dx = 2 L F odd (x) sin nπx L dx = 2 L L L f(x) cos nπx dx n =, 1, 2, L f(x) sin nπx dx n = 1, 2, 3, L This together with Fourier s Theorem give the following reresentations. Theorem Let f be a iecewise smooth function on the interval [, L]. Then f has the following Fourier cosine series reresentation: x (, L) f av (x) = a 2 + a n cos nπx L, where a n = 2 L L In articular at each oint x where f is continuous we have f(x) = a 2 + a n cos nπx L. f(x) cos nπx L dx. Theorem Let f be a iecewise smooth function on the interval [, L]. Then f has the following Fourier sine series reresentation: x (, L) f av (x) = b n sin nπx L, where b n = 2 L L In articular at each oint x where f is continuous we have f(x) = b n sin nπx L. f(x) sin nπx L dx. Examle 1. Let f(x) = x on the interval [, 1]. The 2-eriodic even extension of f is the triangular wave function and the 2-eriodic odd extension of f is the sawtooth function
7 FOURIER SERIES PART III: APPLICATIONS 7 Even eriodic extension of x 1 1 Odd eriodic extension of x 1 1 If we use the even extension we get the Fourier cosine reresentation of x with coefficients and for n 1, 1 a = xdx = 1 a n = 2 x cos nπxdx = 2(( 1)n 1) 1 n 2 π 2 4 Since a 2j = and a 2j+1 = π 2, we get the Fourier cosine reresentation (2j + 1) 2 x over [, 1] as x = π 2 j= cos[(2j + 1)πx] (2j + 1) 2. If we use the odd extension we get the Fourier sine reresentation of x with coefficients b n = x sin nπxdx = 2( 1)n+1 nπ The Fourier sine reresentation of x over the interval [, 1) is x = 2 ( 1) n+1 sin(nπx). π n Examle 2. Find the Fourier cosine series of f(x) = sin x over the interval [, π]. We have a a 1 = 2 π = 2 π π π sin xdx = 4 π ; sin x cos xdx = 1 π To find a n with n 2, we use the identity π sin(2x)dx =. 2 sin x cos(nx) = sin(n + 1)x sin(n 1)x..
8 8 FOURIER SERIES PART III: APPLICATIONS a n π π = 2 sin x cos(nx)dx = 1 n π = 1 [ cos(n 1)x cos(n + 1)x π n 1 n + 1 = 2(( 1)n 1 1) π(n 2 1) We get the Fourier cosine of sin x on [, π] as sin x = 2 π + 2 π n=2 ( 1) n 1 1 n 2 1 (sin(n + 1)x sin(n 1)x) dx ] π cos(nx) = 2 π 4 π j=1 cos(2jx) (2j) Heat Conduction in a Rod Now we are in a osition to solve BVPs with more general nonhomogeneous terms than the ones considered in Note 4. Consider the following BVP for the temerature function u(x, t) in a rod of length L with initial temerature f(x) and with ends ket at temerature. u t = ku xx < x < L, t > u(, t) =, u(l, t) = t > u(x, ) = f(x) < x < L We can aly the method of searation of variables. The homogeneous art of the BVP is u t = ku xx, u(, t) =, u(l, t) = We have seen that solutions u(x, t) = X(x)T (t) (with searated variables) of the homogeneous art leads to the ODE roblems { X (x) + λx(x) = T (t) + kλt (t) =. X() = X(L) = The eigenvalues and eigenfunctions of the X-roblem (Sturm-Liouville roblem) are λ n = ν 2 n, X n (x) = sin(ν n x), where ν n = nπ L, n Z+ For each n Z +, the corresonding T -roblem has a solution T n (t) = e kν2 n t and a solution of the homogeneous art with searated variables is u n (x, t) = T n (t)x n (x). The rincile of suerosition imlies that any linear combination of these solutions is again a solution of the homogeneous art. Thus, u(x, t) = C n T n (t)x n (x) = C n e kν2 n t sin(ν n x) solves formally the homogeneous art of the BVP. At the oints (x, t) where the series converges and term by term differentiation (once in t and twice in x) is allowed, the function u(x, t) defined by the series is a true solution of the homogeneous art. This will be addressed shortly.
9 FOURIER SERIES PART III: APPLICATIONS 9 For now let us find the constants C n so that the formal solution solves also the nonhomogeneous condition u(x, ) = f(x). That is, we would like the constants C n so that u(x, ) = C n e kν2 n sin(ν n x) = f(x). Thus, after relacing ν n by nπ/l, we get f(x) = C n sin nπx L. This is the Fourier sine reresentation of the function f. Therefore the coefficients are given by C n = 2 L f(x) sin nπx L L dx, n Z+. Now we turn our attention to the series and verify that it indeed converges to a twice differentiable function u on < x < L and t > if f is iecewise smooth. For this we will use the Weierstrass M-test to rove uniform convergence. First, let M > be an uer bound of f (i.e. f(x) M for every x [, L]). We have C n 2 L L f(x) sin(ν n x) dx 2M. It follows that for a given t >, we have C n e kν2 n t sin(ν n x) 2Me kν2 n t, t t, x [, L] Since the numerical series n 2Me kν2 n t converges (use ratio or root tests), then it follows from the Weierstrass M-test that the series C n e kν2 n t sin(ν n x) converges uniformly on the set t t, x L. It follows at once that u is a continuous function. We can reeat the argument for the series giving u t and the series giving u xx. That is, the Weierstrass M-test shows that the series u t = ( kνn)c 2 n e kν2 n t sin(ν n x), and u xx = ( νn)c 2 n e kν2 n t sin(ν n x) converge uniformly on t t, x [, L]. We also have u t = ku xx. Consequently the function u(x, t) given by the above series satisfies the comlete BVP. Examle. Consider the BVP We have u t = u xx < x < π, t > u(, t) =, u(π, t) = t > u(x, ) = 1 < x < π π C n = 2 1 sin(nx)dx = 2(1 ( 1)n ) π πn The solution of the BVP is therefore u(x, t) = 2 1 ( 1) n e n2t sin(nx) π n
10 1 FOURIER SERIES PART III: APPLICATIONS 12 Temerature function u(x,t) at various times 1 t=.1 u axis (Temerature) t=1 t=.5 t=.1 2 t= x axis (Rod) or equivalently u(x, t) = 4 π j= ex [ ((2j + 1) 2 t ] sin(2j + 1)x 2j Wave Proagation in a String Consider the BVP for the vibrations of a string with fixed ends. u tt = c 2 u xx < x < L, t > u(, t) =, u(l, t) = t > u(x, ) = f(x) < x < L u t (x, ) = g(x) < x < L Thus u(x, t) reresents the vertical dislacement at time t of the oint x on the string. The initial osition and initial velocities of the string are given by the functions f(x) and g(x). The homogeneous art (HP) of the BVP is u tt = c 2 u xx, u(, t) =, u(l, t) = The solutions u(x, t) = X(x)T (t) (with searated variables) of the homogeneous art leads to the ODE roblems { X (x) + λx(x) = T (t) + c 2 λt (t) =. X() = X(L) = The eigenvalues and eigenfunctions of the X-roblem (SL roblem) are λ n = ν 2 n, X n (x) = sin(ν n x), where ν n = nπ L, n Z+ The corresonding T -roblem has two indeendent solutions T 1 n(t) = cos(cν n t) and T 2 n(t) = sin(cν n t). For each n Z +, we obtain solutions of (HP) with searated variables u 1 n(x, t) = Tn(t)X 1 n (x) = cos(cν n t) sin(ν n x) and u 2 n(x, t) = Tn(t)X 2 n (x) = sin(cν n t) sin(ν n x).
11 FOURIER SERIES PART III: APPLICATIONS 11 The rincile of suerosition imlies that any linear combination of these solutions is again a solution of (HP). Thus, u(x, t) = = A n Tn(t)X 1 n (x) + B n TnX 2 n (x) [A n cos(cν n t) + B n sin(cν n t)] sin(ν n x) is a formal solution of (HP). Now we use the nonhomogeneous conditions to find the coefficients A n and B n. First, we comute the (formal) derivative of u t u t (x, t) = [cν n B n cos(cν n t) cν n A n sin(cν n t)] sin(ν n x). The conditions u(x, ) = f(x) and u t (x, ) = g(x) lead to f(x) = A n sin(ν n x) = A n sin nπx and L cnπ g(x) = cν n B n sin(ν n x) = L B n sin nπx L. These are the Fourier series sine reresentations of f and g on the interval [, L]. Therefore A n = 2 L L f(x) sin nπx L dx and cnπ L B n = 2 L L g(x) sin nπx L dx By using criteria for uniform convergence of Fourier series (Proositions 1 and 2 of Note 7), it can be shown that if f is continuous and iecewise smooth and if g is iecewise smooth, then the series defining u is uniformly convergent and u(x, t) is a continuous function for t and x L. Moreover, we can show that if f, f, f, g, and g are continuous functions on [, L], the function u(x, t) defined by the infinite series is twice differentiable in (x, t) and term by term differentiations in the series are valid. This give u(x, t) as the (unique) solution of BVP. Remark 1. Many concrete roblems involve functions f that are only continuous and iecewise smooth. The series solution u is then only continuous. It is a continuous solution of the BVP. The roblem is understood in a more general sense: in the sense of distributions ( a notion of generalized functions that is beyond the scoe of this course). Remark 2. In concrete alication roblems, to overcome the lack of differentiability of the series solution u, we can to within any degree of accuracy ϵ, relace the functions f and g by their truncated Fourier series S N f and S N g so that f S N f < ϵ, g S N g < ϵ on [, L] The functions S N f and S N g are infinitely differentiable and the corresonding solution u N (the truncated series of u) is infinitely differentiable.
12 12 FOURIER SERIES PART III: APPLICATIONS Remark 3. By using the rincile of suerosition, this BVP could have been slit into two BVPs: BVP1 (lucked string) and BVP2 (struck string) v tt = c 2 v xx < x < L, t > v(, t) =, v(l, t) = t > v(x, ) = f(x) < x < L v t (x, ) = < x < L w tt = c 2 w xx < x < L, t > w(, t) =, w(l, t) = t > w(x, ) = < x < L w t (x, ) = g(x) < x < L The solutions to BVP1 and BVP2 are, resectively, v = w = A n cos(cν n t) sin(ν n x) B n sin(cν n t) sin(ν n x) The solution to the original BVP is u = v + w. Examle 1. (Plucked string) Consider the BVP u tt = 4v xx < x < 1, t > u(, t) =, u(l, 1) = t > u(x, ) = f(x) < x < 1 u t (x, ) = < x < 1 where 1 Initial osition of the string 5 1 f(x) = { x/5 if x 5 (1 x)/5 if5 x 1 For such an initial osition we have (B n = ) and A n = = 1 25 = 8 π 2 n 2 sin nπ 2 f(x) sin nπx 1 dx x sin nπx 1 dx (1 x) sin nπx 1 dx
13 FOURIER SERIES PART III: APPLICATIONS 13 Hence, A 2j = and A 2j+1 = 8( 1)j π 2. The series solution is (2j + 1) 2 u(x, t) = 8 π 2 cos((2j + 1)t/5) ( 1)j sin((2j + 1)πx/1) (2j + 1) 2 j= = 8 ( cos(3t/5) sin(3x/1) π 2 cos(t/5) sin(x/1) + 9 cos(5t/5) sin(5x/1) cos(7t/5) sin(7x/1) ) + The individual comonents u n (x, t) = cos(nπt/5) sin(nπx/1) are called the harmonics or modes of vibrations. The function u n (x, t) is just a sine function in x being scaled by a cosine function in t with frequency n/1. First mode Second mode Third mode Figure 3. The first three modes of vibrations of the lucked string at various times. Examle 2. (Struck string) Consider the BVP u tt = 4v xx < x < 1, t > u(, t) =, u(l, 1) = t > u(x, ) = < x < 1 u t (x, ) = g(x) < x < 1 where g(x) = 1 if 4 < x < 6, and g(x) = elsewhere. This time A n =, and Thus B n = 5 π 2 n 2 The series solution is u(x, t) = 1 π 2 = 1 π nπ 5 B n = 2 1 = 1 πn 1 ( g(x) sin nπx 1 dx cos 2nπ 5 cos 3nπ 5 ( cos 2nπ ) 3nπ cos = π 2 n 2 sin nπ 2 1 n 2 sin nπ 2 ( sin π 1 sin πt 5 nπ nπt nπx sin sin sin ) sin nπ 1 πx sin 1 1 3π 3πt 3πx sin sin sin π 5πt 5πx sin sin sin π 7πt 7πx sin sin sin )
14 14 FOURIER SERIES PART III: APPLICATIONS 6. Problems Dealing with the Lalace Equation Recall that the Dirichlet roblem in a rectangle is to find a harmonic function u inside the rectangle whose values on the boundary are given. That is u(x, y) = < x < L, < y < H u(x, ) = f 1 (x), u(x, H) = f 2 (x) < x < L u(, y) = g 1 (y), u(l, y) = g 2 (2) < y < H To solve this roblem, we use the rincile of suerosition to decomose it into four simler subroblems as in the figure. u 1 = u 1 = u 1 = u 2 =f 2 (x) u = 1 + u = 2 u = 2 + u = 2 u=f 2 (x) u 1 =f 1 (x) u 2 = u=g 1 (y) u= u=g 2 (y) = u=f 1 (x) u 3 = u 4 = u 3 =g 1 (y) u 3 = u 3 = + u 4 = u 4 = u 4 =g 2 (y) u 3 = u 4 = We can find the solution u(x, y) as u(x, y) = u 1 (x, y) + u 2 (x, y) + u 3 (x, y) + u 4 (x, y). Each of the subroblems can be solved by the method of searation of variables. Now we indicate how to find u 1 (x, y). The solutions with searated variables u 1 (x, y) = X(x)Y (y) of the homogeneous art leads to the ODE roblems { { X (x) + λx(x) = Y and (y) λy (y) = X() = X(L) = Y (H) = where λ is the searation constant. The X-roblem is an SL-roblem whose eigenvalues and eigenfunctions are λ n = ν 2 n, X n (x) = sin(ν n x), where ν n = nπ L, n Z+. For each λ n, the corresonding ODE for the Y -roblem has general solution Y n = A cosh(ν n y) + B sinh(ν n y). The boundary condition Y (H) =, imlies that (u to a multilicative constant), the solution of the Y -roblem is Y n (y) = sinh [ν n (H y)].
15 FOURIER SERIES PART III: APPLICATIONS 15 Hence, the solutions with searated variables of the homogeneous art of the u 1 - roblem are generated by u 1,n (x, y) = sinh [ν n (H y)] sin(ν n x) n Z +. A series solution of the u 1 -roblem is therefore u 1 (x, y) = C n sinh [ν n (H y)] sin(ν n x). Such a solution solves the nonhomogeneous condition u(x, ) = f 1 (x) if and only if f 1 (x) = C n sinh(ν n H) sin(ν n x) = C n sinh nπh nπx sin L L. Thus, C n sinh(ν n H) is the n-th Fourier sine coefficient of f 1 over [, L]: C n sinh nπh L = 2 L L f 1 (x) sin nπx L dx. The functions u 2, u 3, and u 4 can be found in a similar way. Examle 1. Consider the following BVP with mixed boundary conditions. u(x, y) = < x < π, < y < 2π u(x, ) = x, u(x, 2π) = < x < π u x (, y) =, u x (π, y) = 1 < y < 2π We decomose the roblem as shown in the figure (see next age). u= v= w= u x = u= u x =1 = v x = v= v x = + w x = w= w x =1 u=x v=x w= The method of searation of variables for the v-roblem leads to the ODE roblems { { X (x) + λx(x) = Y X () = X and (y) λy (y) = (π) = Y (2π) = where λ is the searation constant. The X-roblem is an SL-roblem whose eigenvalues and eigenfunctions are λ =, X (x) = 1 and for n Z + λ n = n 2, X n (x) = cos(nx). For λ =, the general solution of the ODE for the Y -roblem is Y (y) = A + By and in order to get Y (2π) =, we need A = 2Bπ. Thus Y (y) = (2π y) generates the solutions of the Y -roblem For λ n = n 2, the solutions of the Y -roblem are generated by Y n (y) = sinh [n(2π y)].
16 16 FOURIER SERIES PART III: APPLICATIONS The solutions with searated variables of the homogeneous art of the v-roblem are therefore v (x, y) = 2π y and v n (x, y) = sinh [n(2π y)] cos(nx) for n Z + The series solution is v(x, y) = C (2π y) + C n sinh [n(2π y)] cos(nx). In order for such a series to solve the nonhomogeneous condition v(x, ) = x, we need to have x = 2πC + C n sinh(2nπ) cos(nx). This is the Fourier cosine exansion of x over [, π]. Hence, and for n 1 This gives 2πC = 2 π sinh(2nπ)c n = 2 π π π C 2j = and C 2j+1 = The solution of the v-roblem is v(x, y) = 2π y 4 4 π j= xdx = π C = 1 2 x cos(nx)dx = 2 (( 1)n 1) πn 2. 4 π(2j + 1) 2 sinh [2π(2j + 1)]. sinh [(2j + 1)(2π y)] cos(2j + 1)x sinh [2(2j + 1)π] (2j + 1) 2. Now we solve the w-roblem. The searation of variables for the homogeneous art leads to the ODE roblems. { { X (x) λx(x) = Y X and (y) + λy (y) = () = Y () = Y (2π) = This time it is the Y -roblem that is a Sturm-Liouville roblem with eigenvalues and eigenfunctions λ n = n2 2 2, Y n(y) = sin ny 2, n Z+. For each n, a generator of the solutions of the X-roblem is X n (x) = cosh nx 2. The series solution of the w-roblem is therefore w(x, y) = C n cosh nx ny sin 2 2. To find the coefficients C n so that w x (π, y) 1, we need n w x (x, y) = 2 C n sinh nx ny sin 2 2
17 FOURIER SERIES PART III: APPLICATIONS 17 This gives 1 = n 2 C n sinh nπ 2 sin ny 2 (the Fourier sine exansion of 1 over the interval [, 2π]): Equivalently, n 2 C n sinh nπ 2 = 2 2π C 2j =, C 2j+1 = 2π sin ny 2 dy = 2 (1 ( 1)n ). nπ 8 π(2j + 1) 2 sinh [(2j + 1)π/2]. The solution of the w-roblem is therefore w(x, y) = 8 cosh [(2j + 1)x/2] sin [(2j + 1)y/2] π sinh [(2j + 1)π/2] (2j + 1) 2 j= The solution u of the original roblem is u(x, y) = v(x, y) + w(x, y). Examle 2. Consider the Dirichlet roblem in a disk (written in olar coordinates) We take the function f to be given by u(r, θ) = r < 1, θ [, 2π] u(1, θ) = f(θ) θ [, 2π] u(1,θ )=f(θ) u= f(θ) = { sin θ if θ π if π θ 2π Recall that the Lalace oerator in olar coordinates is = 2 r r r r 2 θ 2 A solution with searated variables u(r, θ) = R(r)Θ(θ) leads to the ODEs Θ (θ) + λθ(θ) = and r 2 R (r) + rr (r) λr(r) = where λ is the searation constant. Note that since u(r, θ + 2π) = u(r, θ), then the function Θ and also Θ need to be 2π-eriodic. Thus to the ODE for Θ we need to add Θ() = Θ(2π) and Θ () = Θ (2π). Hence, the Θ-roblem is a eriodic SL-roblem whose eigenvalues and eigenfunctions are λ =, Θ (θ) = 1,
18 18 FOURIER SERIES PART III: APPLICATIONS and for n Z +, λ n = n 2, Θ 1 n(θ) = cos(nθ), Θ 2 n(θ) = sin(nθ). The ODE for the R-function is a Cauchy-Euler equation with characteristic equation m 2 λ =. For λ = λ =, the general solution of the R-equation are generated by R(r) 1 = 1 and R(r) 2 = ln r. For λ = λ n (we have m = ±n), the solutions of the R-equation are generated by R 1 n(r) = r n and R 2 n(r) = r n. The solutions with searated variables of the Lalace equation u = in the disk are therefore 1, ln r, r n cos(nθ), r n sin(nθ), r n cos(nθ), r n sin(nθ). Since we looking for solutions that are bounded in the disk and since ln r and r n cos(nθ) and r n sin(nθ) are not bounded, then we will discard then when forming u The series solution has the form u(r, θ) = A + (A n r n cos(nθ) + B n r n sin(nθ)) The initial condition u(1, θ) = f(θ) leads to f(θ) = A + A n cos(nθ) + B n sin(nθ). This is the Fourier series of f. I leave it as an exercise for you to verify that the Fourier series of f is f(θ) = 1 π sin θ 2 cos(2jθ) π 4j 2 1. The solution of the Dirichlet roblem is u(r, θ) = 1 π + r 2 sin θ 2 π j=1 j=1 7. Exercises r 2j cos(2jθ) 4j 2 1 Exercise 1. (a) Find the Fourier series of the function with eriod 4 that is defined over [ 2, 2] by f(x) = 4 x2. 2 (b) Use Parseval s equality to evaluate the series 1 n 4. (c) Use the integral test to estimate the mean square error E N when relacing f by its truncated Fourier series S N f. (d) Find N so that E N.1 and then find N so that E N.1 Exercise 2. (a) Find the Fourier series of the function with eriod 4 that is defined over [ 2, 2] by { 1 x if x 2 f(x) = 1 + x if 2 x
19 FOURIER SERIES PART III: APPLICATIONS 19 (b) Use Parseval s equality to evaluate the series j= 1 (2j + 1) 4. (c) Use the integral test to estimate the mean square error E N when relacing f by its truncated Fourier series S N f. (d) Find N so that E N.1 and then find N so that E N.1 Exercise 3. Find the Fourier sine series of f(x) = cos x over [, π] (What is the Fourier cosine series of cos x on [, π]?) Exercise 4. Find the Fourier cosine series of f(x) = sin x over [, π] (What is the Fourier sine series of sin x on [, π]?) Exercise 5. Find the Fourier cosine series of f(x) = x 2 over [, 1]. Exercise 6. Find the Fourier sine series of f(x) = x 2 over [, 1]. Exercise 7. Find the Fourier cosine series of f(x) = x sin x over [, π]. Exercise 8. Find the Fourier sine series of f(x) = x sin x over [, π]. Exercise 9. Solve the BVP u t = u xx, < x < 2, t > u(, t) = u(2, t) =, t > u(x, ) = f(x), < x < 2 where f(x) = { 1 if < x < 1 if 1 < x < 2 Exercise 1. Solve the BVP u t = u xx, < x < 2, t > u(, t) = u(2, t) =, t > u(x, ) = cos(πx), < x < 2 Exercise 11. Solve the BVP u t + u = (.1)u xx, < x < π, t > u x (, t) = u x (π, t) =, t > u(x, ) = sin x, < x < 2 Exercise 12. Consider the BVP modeling heat roagation in a rod where the end oints are ket at constant temeratures T 1 and T 2 : u t = ku xx, < x < L, t > u(, t) = T 1, u(l, t) = T 2, t > u(x, ) = f(x), < x < L Since T 1 and T 2 are not necessarily zero, we cannot aly directly the method of eigenfunctions exansion. To solve such a roblem, we can roceed as follows. 1. Find a function α(x) (indeendent on time t) so that α (x) =, α() = T 1 α(l) = T 2.
20 2 FOURIER SERIES PART III: APPLICATIONS 2. Let v(x, t) = u(x, t) α(x). Verify that if u(x, t) solves the given BVP, then v(x, t) solves the following roblem v t = kv xx, < x < L, t > v(, t) =, v(l, t) =, t > v(x, ) = f(x) α(x), < x < L The v-roblem can be solved by the method of searation of variables. The solution u of the original roblem is therefore u(x, t) = v(x, t) + α(x). Exercise 13. Aly the method of described in Exercise 12 to solve the roblem u t = u xx, < x < 2, t > u(, t) = T 1, u(2, t) = T 2, t > u(x, ) = f(x), < x < 2 in the following cases 1. T 1 = 1, T 2 =, f(x) =. 2. T 1 = 1, T 2 = 1, f(x) =. 3. T 1 =, T 2 = 1, f(x) = 5x. In roblems 14 to 16, solve the wave roagation roblem u tt = c 2 u xx, < x < L, t > u(, t) =, u(l, t) =, t > u(x, ) = f(x), u t (x, ) = g(x) < x < L { x if < x < 1 Exercise 14. c = 1, L = 2, f(x) =, g(x) = 2 x if 1 < x < 2 { x if < x < 1 Exercise 15. c = 1/π, L = 2, f(x) = sin x, g(x) = 2 x if 1 < x < 2 Exercise 16. c = 2, L = π, f(x) = x sin x, g(x) = sin(2x). In exercises 17 to 19, solve the wave roagation roblem with daming u tt + 2au t = c 2 u xx, < x < L, t > u(, t) =, u(l, t) =, t > u(x, ) = f(x), u t (x, ) = g(x) < x < L Exercise 17. c = 1, a =.5, L = π, f(x) =, g(x) = x Exercise 18. c = 4, a = π, L = 1, f(x) = x(1 x), g(x) =. Exercise 19. c = 1, a = π/6, L = 2, f(x) = x sin(πx), g(x) = 1. In exercises 2 to 22, solve the Lalace equation u(x, y) = inside the rectangle < x < L, < y < H subject the the given boundary conditions. Exercise 2. L = H = π, u(x, ) = x(π x), u(x, π) =, u(, y) = u(π, y) =. Exercise 21. L = π, H = 2π, u(x, ) =, u(x, 2π) = x, u x (, y) = sin y, u x (π, y) =. Exercise 22. L = H = 1, u(x, ) = u(x, 1) =, u(, y) = 1, u(1, y) = sin y. Exercise 23. Solve the Lalace equation u(r, θ) = inside the semicircle of radius 2 ( < r < 2, < θ < π ) subject to the boundary conditions u(r, ) = u(r, π) = ( < r < 2) and u(2, θ) = θ(π θ) ( < θ < π)
21 FOURIER SERIES PART III: APPLICATIONS 21 Exercise 24. Solve the Lalace equation u(r, θ) = inside the semicircle of radius 2 ( < r < 2, < θ < π ) subject to the boundary conditions u θ (r, ) = u θ (r, π) = ( < r < 2) and u(2, θ) = θ(π θ) ( < θ < π) Exercise 25. Solve the Lalace equation u(r, θ) = inside the quarter of a circle of radius 2 ( < r < 2, < θ < π/2 ) subject to the boundary conditions u(r, ) = u(r, π/2) = ( < r < 2) and u(2, θ) = θ ( < θ < π/2)
THE METHOD OF SEPARATION OF VARIABLES
THE METHOD OF SEPARATION OF VARIABES To solve the BVPs that we have encountered so far, we will use separation of variables on the homogeneous part of the BVP. This separation of variables leads to problems
More informationWave Equation With Homogeneous Boundary Conditions
Wave Equation With Homogeneous Boundary Conditions MATH 467 Partial Differential Equations J. Robert Buchanan Department of Mathematics Fall 018 Objectives In this lesson we will learn: how to solve the
More informationMath 201 Assignment #11
Math 21 Assignment #11 Problem 1 (1.5 2) Find a formal solution to the given initial-boundary value problem. = 2 u x, < x < π, t > 2 u(, t) = u(π, t) =, t > u(x, ) = x 2, < x < π Problem 2 (1.5 5) Find
More informationMA Chapter 10 practice
MA 33 Chapter 1 practice NAME INSTRUCTOR 1. Instructor s names: Chen. Course number: MA33. 3. TEST/QUIZ NUMBER is: 1 if this sheet is yellow if this sheet is blue 3 if this sheet is white 4. Sign the scantron
More informationPartial Differential Equations for Engineering Math 312, Fall 2012
Partial Differential Equations for Engineering Math 312, Fall 2012 Jens Lorenz July 17, 2012 Contents Department of Mathematics and Statistics, UNM, Albuquerque, NM 87131 1 Second Order ODEs with Constant
More informationMath Assignment 14
Math 2280 - Assignment 14 Dylan Zwick Spring 2014 Section 9.5-1, 3, 5, 7, 9 Section 9.6-1, 3, 5, 7, 14 Section 9.7-1, 2, 3, 4 1 Section 9.5 - Heat Conduction and Separation of Variables 9.5.1 - Solve the
More information17 Source Problems for Heat and Wave IB- VPs
17 Source Problems for Heat and Wave IB- VPs We have mostly dealt with homogeneous equations, homogeneous b.c.s in this course so far. Recall that if we have non-homogeneous b.c.s, then we want to first
More informationSAMPLE FINAL EXAM SOLUTIONS
LAST (family) NAME: FIRST (given) NAME: ID # : MATHEMATICS 3FF3 McMaster University Final Examination Day Class Duration of Examination: 3 hours Dr. J.-P. Gabardo THIS EXAMINATION PAPER INCLUDES 22 PAGES
More informationBoundary-value Problems in Rectangular Coordinates
Boundary-value Problems in Rectangular Coordinates 2009 Outline Separation of Variables: Heat Equation on a Slab Separation of Variables: Vibrating String Separation of Variables: Laplace Equation Review
More informationReview For the Final: Problem 1 Find the general solutions of the following DEs. a) x 2 y xy y 2 = 0 solution: = 0 : homogeneous equation.
Review For the Final: Problem 1 Find the general solutions of the following DEs. a) x 2 y xy y 2 = 0 solution: y y x y2 = 0 : homogeneous equation. x2 v = y dy, y = vx, and x v + x dv dx = v + v2. dx =
More informationMath 2930 Worksheet Final Exam Review
Math 293 Worksheet Final Exam Review Week 14 November 3th, 217 Question 1. (* Solve the initial value problem y y = 2xe x, y( = 1 Question 2. (* Consider the differential equation: y = y y 3. (a Find the
More informationMATH-UA 263 Partial Differential Equations Recitation Summary
MATH-UA 263 Partial Differential Equations Recitation Summary Yuanxun (Bill) Bao Office Hour: Wednesday 2-4pm, WWH 1003 Email: yxb201@nyu.edu 1 February 2, 2018 Topics: verifying solution to a PDE, dispersion
More informationM412 Assignment 5 Solutions
M41 Assignment 5 Solutions 1. [1 pts] Haberman.5.1 (a). Solution. Setting u(x, y) =X(x)Y(y), we find that X and Y satisfy X + λx = Y λy =, X () = X () = Y() =. In this case, we find from the X equation
More informationProblem (p.613) Determine all solutions, if any, to the boundary value problem. y + 9y = 0; 0 < x < π, y(0) = 0, y (π) = 6,
Problem 10.2.4 (p.613) Determine all solutions, if any, to the boundary value problem y + 9y = 0; 0 < x < π, y(0) = 0, y (π) = 6, by first finding a general solution to the differential equation. Solution.
More informationSolutions to the Sample Problems for the Final Exam UCLA Math 135, Winter 2015, Professor David Levermore
Solutions to the Sample Problems for the Final Exam UCLA Math 135, Winter 15, Professor David Levermore Every sample problem for the Midterm exam and every problem on the Midterm exam should be considered
More informationName: Math Homework Set # 5. March 12, 2010
Name: Math 4567. Homework Set # 5 March 12, 2010 Chapter 3 (page 79, problems 1,2), (page 82, problems 1,2), (page 86, problems 2,3), Chapter 4 (page 93, problems 2,3), (page 98, problems 1,2), (page 102,
More informationMath 316/202: Solutions to Assignment 7
Math 316/22: Solutions to Assignment 7 1.8.6(a) Using separation of variables, we write u(r, θ) = R(r)Θ(θ), where Θ() = Θ(π) =. The Laplace equation in polar coordinates (equation 19) becomes R Θ + 1 r
More information6 Non-homogeneous Heat Problems
6 Non-homogeneous Heat Problems Up to this point all the problems we have considered for the heat or wave equation we what we call homogeneous problems. This means that for an interval < x < l the problems
More informationMidterm 2: Sample solutions Math 118A, Fall 2013
Midterm 2: Sample solutions Math 118A, Fall 213 1. Find all separated solutions u(r,t = F(rG(t of the radially symmetric heat equation u t = k ( r u. r r r Solve for G(t explicitly. Write down an ODE for
More informationElementary Analysis in Q p
Elementary Analysis in Q Hannah Hutter, May Szedlák, Phili Wirth November 17, 2011 This reort follows very closely the book of Svetlana Katok 1. 1 Sequences and Series In this section we will see some
More informationENGG 2420: Fourier series and partial differential equations
ENGG 2420: Fourier series and artial differential equations Prof. Thierry Blu e-mail: tblu@ee.cuhk.edu.hk web: www.ee.cuhk.edu.hk/~tblu/ Deartment of Electronic Engineering The Chinese University of Hong
More informationMATH 412 Fourier Series and PDE- Spring 2010 SOLUTIONS to HOMEWORK 6
MATH 412 Fourier Series and PDE- Spring 21 SOLUTIONS to HOMEWORK 6 Problem 1. Solve the following problem u t u xx =1+xcos t
More informationAnalysis III Solutions - Serie 12
.. Necessary condition Let us consider the following problem for < x, y < π, u =, for < x, y < π, u y (x, π) = x a, for < x < π, u y (x, ) = a x, for < x < π, u x (, y) = u x (π, y) =, for < y < π. Find
More informationSeparation of Variables
Separation of Variables A typical starting point to study differential equations is to guess solutions of a certain form. Since we will deal with linear PDEs, the superposition principle will allow us
More informationG: Uniform Convergence of Fourier Series
G: Uniform Convergence of Fourier Series From previous work on the prototypical problem (and other problems) u t = Du xx 0 < x < l, t > 0 u(0, t) = 0 = u(l, t) t > 0 u(x, 0) = f(x) 0 < x < l () we developed
More informationCLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE
CLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE 1. Linear Partial Differential Equations A partial differential equation (PDE) is an equation, for an unknown function u, that
More informationMATH 251 Final Examination December 19, 2012 FORM A. Name: Student Number: Section:
MATH 251 Final Examination December 19, 2012 FORM A Name: Student Number: Section: This exam has 17 questions for a total of 150 points. In order to obtain full credit for partial credit problems, all
More informationMethod of Separation of Variables
MODUE 5: HEAT EQUATION 11 ecture 3 Method of Separation of Variables Separation of variables is one of the oldest technique for solving initial-boundary value problems (IBVP) and applies to problems, where
More informationIntroduction to Sturm-Liouville Theory and the Theory of Generalized Fourier Series
CHAPTER 5 Introduction to Sturm-Liouville Theory and the Theory of Generalized Fourier Series We start with some introductory examples. 5.. Cauchy s equation The homogeneous Euler-Cauchy equation (Leonhard
More information21 Laplace s Equation and Harmonic Functions
2 Laplace s Equation and Harmonic Functions 2. Introductory Remarks on the Laplacian operator Given a domain Ω R d, then 2 u = div(grad u) = in Ω () is Laplace s equation defined in Ω. If d = 2, in cartesian
More informationPHYS 301 HOMEWORK #9-- SOLUTIONS
PHYS 0 HOMEWORK #9-- SOLUTIONS. We are asked to use Dirichlet' s theorem to determine the value of f (x) as defined below at x = 0, ± /, ± f(x) = 0, - < x
More informationLecture Notes 4: Fourier Series and PDE s
Lecture Notes 4: Fourier Series and PDE s 1. Periodic Functions A function fx defined on R is caed a periodic function if there exists a number T > such that fx + T = fx, x R. 1.1 The smaest number T for
More information0 3 x < x < 5. By continuing in this fashion, and drawing a graph, it can be seen that T = 2.
04 Section 10. y (π) = c = 0, and thus λ = 0 is an eigenvalue, with y 0 (x) = 1 as the eigenfunction. For λ > 0 we again have y(x) = c 1 sin λ x + c cos λ x, so y (0) = λ c 1 = 0 and y () = -c λ sin λ
More informationFINAL EXAM, MATH 353 SUMMER I 2015
FINAL EXAM, MATH 353 SUMMER I 25 9:am-2:pm, Thursday, June 25 I have neither given nor received any unauthorized help on this exam and I have conducted myself within the guidelines of the Duke Community
More informationStrauss PDEs 2e: Section Exercise 4 Page 1 of 6
Strauss PDEs 2e: Section 5.3 - Exercise 4 Page of 6 Exercise 4 Consider the problem u t = ku xx for < x < l, with the boundary conditions u(, t) = U, u x (l, t) =, and the initial condition u(x, ) =, where
More informationLAPLACE EQUATION. = 2 is call the Laplacian or del-square operator. In two dimensions, the expression of in rectangular and polar coordinates are
LAPLACE EQUATION If a diffusion or wave problem is stationary (time independent), the pde reduces to the Laplace equation u = u =, an archetype of second order elliptic pde. = 2 is call the Laplacian or
More information4.10 Dirichlet problem in the circle and the Poisson kernel
220 CHAPTER 4. FOURIER SERIES AND PDES 4.10 Dirichlet problem in the circle and the Poisson kernel Note: 2 lectures,, 9.7 in [EP], 10.8 in [BD] 4.10.1 Laplace in polar coordinates Perhaps a more natural
More informationMore on Fourier Series
More on Fourier Series R. C. Trinity University Partial Differential Equations Lecture 6.1 New Fourier series from old Recall: Given a function f (x, we can dilate/translate its graph via multiplication/addition,
More informationChapter 10: Partial Differential Equations
1.1: Introduction Chapter 1: Partial Differential Equations Definition: A differential equations whose dependent variable varies with respect to more than one independent variable is called a partial differential
More informationSolutions to Exercises 8.1
Section 8. Partial Differential Equations in Physics and Engineering 67 Solutions to Exercises 8.. u xx +u xy u is a second order, linear, and homogeneous partial differential equation. u x (,y) is linear
More informationLEGENDRE POLYNOMIALS AND APPLICATIONS. We construct Legendre polynomials and apply them to solve Dirichlet problems in spherical coordinates.
LEGENDRE POLYNOMIALS AND APPLICATIONS We construct Legendre polynomials and apply them to solve Dirichlet problems in spherical coordinates.. Legendre equation: series solutions The Legendre equation is
More informationSC/MATH Partial Differential Equations Fall Assignment 3 Solutions
November 16, 211 SC/MATH 3271 3. Partial Differential Equations Fall 211 Assignment 3 Solutions 1. 2.4.6 (a) on page 7 in the text To determine the equilibrium (also called steady-state) heat distribution
More informationMath 251 December 14, 2005 Final Exam. 1 18pt 2 16pt 3 12pt 4 14pt 5 12pt 6 14pt 7 14pt 8 16pt 9 20pt 10 14pt Total 150pt
Math 251 December 14, 2005 Final Exam Name Section There are 10 questions on this exam. Many of them have multiple parts. The point value of each question is indicated either at the beginning of each question
More informationF(p) y + 3y + 2y = δ(t a) y(0) = 0 and y (0) = 0.
Page 5- Chater 5: Lalace Transforms The Lalace Transform is a useful tool that is used to solve many mathematical and alied roblems. In articular, the Lalace transform is a technique that can be used to
More informationFinal: Solutions Math 118A, Fall 2013
Final: Solutions Math 118A, Fall 2013 1. [20 pts] For each of the following PDEs for u(x, y), give their order and say if they are nonlinear or linear. If they are linear, say if they are homogeneous or
More informationSeparation of Variables. A. Three Famous PDE s
Separation of Variables c 14, Philip D. Loewen A. Three Famous PDE s 1. Wave Equation. Displacement u depends on position and time: u = u(x, t. Concavity drives acceleration: u tt = c u xx.. Heat Equation.
More informationMATH 251 Final Examination May 4, 2015 FORM A. Name: Student Number: Section:
MATH 251 Final Examination May 4, 2015 FORM A Name: Student Number: Section: This exam has 16 questions for a total of 150 points. In order to obtain full credit for partial credit problems, all work must
More informationSeparation of variables in two dimensions. Overview of method: Consider linear, homogeneous equation for u(v 1, v 2 )
Separation of variables in two dimensions Overview of method: Consider linear, homogeneous equation for u(v 1, v 2 ) Separation of variables in two dimensions Overview of method: Consider linear, homogeneous
More informationVibrating-string problem
EE-2020, Spring 2009 p. 1/30 Vibrating-string problem Newton s equation of motion, m u tt = applied forces to the segment (x, x, + x), Net force due to the tension of the string, T Sinθ 2 T Sinθ 1 T[u
More informationc2 2 x2. (1) t = c2 2 u, (2) 2 = 2 x x 2, (3)
ecture 13 The wave equation - final comments Sections 4.2-4.6 of text by Haberman u(x,t), In the previous lecture, we studied the so-called wave equation in one-dimension, i.e., for a function It was derived
More informationMATH 241 Practice Second Midterm Exam - Fall 2012
MATH 41 Practice Second Midterm Exam - Fall 1 1. Let f(x = { 1 x for x 1 for 1 x (a Compute the Fourier sine series of f(x. The Fourier sine series is b n sin where b n = f(x sin dx = 1 = (1 x cos = 4
More informationIn this chapter we study elliptical PDEs. That is, PDEs of the form. 2 u = lots,
Chapter 8 Elliptic PDEs In this chapter we study elliptical PDEs. That is, PDEs of the form 2 u = lots, where lots means lower-order terms (u x, u y,..., u, f). Here are some ways to think about the physical
More informationMath 5587 Midterm II Solutions
Math 5587 Midterm II Solutions Prof. Jeff Calder November 3, 2016 Name: Instructions: 1. I recommend looking over the problems first and starting with those you feel most comfortable with. 2. Unless otherwise
More informationPartial differential equation for temperature u(x, t) in a heat conducting insulated rod along the x-axis is given by the Heat equation:
Chapter 7 Heat Equation Partial differential equation for temperature u(x, t) in a heat conducting insulated rod along the x-axis is given by the Heat equation: u t = ku x x, x, t > (7.1) Here k is a constant
More information1 Partial derivatives and Chain rule
Math 62 Partial derivatives and Chain rule Question : Let u be a solution to the PDE t u(x, t) + 2 xu 2 (x, t) ν xx u(x, t) =, x (, + ), t >. (a) Let ψ(x, t) = x tu(ξ, t)dξ + 2 u2 (x, t) ν x u(x, t). Compute
More informationswapneel/207
Partial differential equations Swapneel Mahajan www.math.iitb.ac.in/ swapneel/207 1 1 Power series For a real number x 0 and a sequence (a n ) of real numbers, consider the expression a n (x x 0 ) n =
More informationA NOTE ON RECURRENCE FORMULA FOR VALUES OF THE EULER ZETA FUNCTIONS ζ E (2n) AT POSITIVE INTEGERS. 1. Introduction
Bull. Korean Math. Soc. 5 (4), No. 5,. 45 43 htt://dx.doi.org/.434/bkms.4.5.5.45 A NOTE ON RECURRENCE FORMULA FOR VALUES OF THE EULER ZETA FUNCTIONS ζ E (n) AT POSITIVE INTEGERS Hui Young Lee and Cheon
More informationFourier Series Tutorial
Fourier Series Tutorial INTRODUCTION This document is designed to overview the theory behind the Fourier series and its alications. It introduces the Fourier series and then demonstrates its use with a
More informationIntroduction and preliminaries
Chapter Introduction and preliminaries Partial differential equations What is a partial differential equation? ODEs Ordinary Differential Equations) have one variable x). PDEs Partial Differential Equations)
More informationPartial Differential Equations
Partial Differential Equations Spring Exam 3 Review Solutions Exercise. We utilize the general solution to the Dirichlet problem in rectangle given in the textbook on page 68. In the notation used there
More informationMath 2930 Worksheet Wave Equation
Math 930 Worksheet Wave Equation Week 13 November 16th, 017 Question 1. Consider the wave equation a u xx = u tt in an infinite one-dimensional medium subject to the initial conditions u(x, 0) = 0 u t
More information2.9 Dirac Notation Vectors ji that are orthonormal hk ji = k,j span a vector space and express the identity operator I of the space as (1.
14 Fourier Series 2.9 Dirac Notation Vectors ji that are orthonormal hk ji k,j san a vector sace and exress the identity oerator I of the sace as (1.132) I NX jihj. (2.99) Multilying from the right by
More informationStatics and dynamics: some elementary concepts
1 Statics and dynamics: some elementary concets Dynamics is the study of the movement through time of variables such as heartbeat, temerature, secies oulation, voltage, roduction, emloyment, rices and
More informationLast Update: April 4, 2011
Math ES Winter Last Update: April 4, Introduction to Partial Differential Equations Disclaimer: his lecture note tries to provide an alternative approach to the material in Sections..6 in the textbook.
More informationMA 201: Differentiation and Integration of Fourier Series Applications of Fourier Series Lecture - 10
MA 201: Differentiation and Integration of Fourier Series Applications of Fourier Series ecture - 10 Fourier Series: Orthogonal Sets We begin our treatment with some observations: For m,n = 1,2,3,... cos
More information(The) Three Linear Partial Differential Equations
(The) Three Linear Partial Differential Equations 1 Introduction A partial differential equation (PDE) is an equation of a function of 2 or more variables, involving 2 or more partial derivatives in different
More informationFourier and Partial Differential Equations
Chapter 5 Fourier and Partial Differential Equations 5.1 Fourier MATH 294 SPRING 1982 FINAL # 5 5.1.1 Consider the function 2x, 0 x 1. a) Sketch the odd extension of this function on 1 x 1. b) Expand the
More informationMathematical Modeling using Partial Differential Equations (PDE s)
Mathematical Modeling using Partial Differential Equations (PDE s) 145. Physical Models: heat conduction, vibration. 146. Mathematical Models: why build them. The solution to the mathematical model will
More informationMathematical Methods: Fourier Series. Fourier Series: The Basics
1 Mathematical Methods: Fourier Series Fourier Series: The Basics Fourier series are a method of representing periodic functions. It is a very useful and powerful tool in many situations. It is sufficiently
More informationBoundary value problems for partial differential equations
Boundary value problems for partial differential equations Henrik Schlichtkrull March 11, 213 1 Boundary value problem 2 1 Introduction This note contains a brief introduction to linear partial differential
More informationA proof for the full Fourier series on [ π, π] is given here.
niform convergence of Fourier series A smooth function on an interval [a, b] may be represented by a full, sine, or cosine Fourier series, and pointwise convergence can be achieved, except possibly at
More informationMath 5440 Problem Set 7 Solutions
Math 544 Math 544 Problem Set 7 Solutions Aaron Fogelson Fall, 13 1: (Logan, 3. # 1) Verify that the set of functions {1, cos(x), cos(x),...} form an orthogonal set on the interval [, π]. Next verify that
More information# Points Score Total 100
Name: PennID: Math 241 Make-Up Final Exam January 19, 2016 Instructions: Turn off and put away your cell phone. Please write your Name and PennID on the top of this page. Please sign and date the pledge
More informationDiffusion on the half-line. The Dirichlet problem
Diffusion on the half-line The Dirichlet problem Consider the initial boundary value problem (IBVP) on the half line (, ): v t kv xx = v(x, ) = φ(x) v(, t) =. The solution will be obtained by the reflection
More informationPartial Differential Equations
Partial Differential Equations Xu Chen Assistant Professor United Technologies Engineering Build, Rm. 382 Department of Mechanical Engineering University of Connecticut xchen@engr.uconn.edu Contents 1
More informationThe One-Dimensional Heat Equation
The One-Dimensional Heat Equation R. C. Trinity University Partial Differential Equations February 24, 2015 Introduction The heat equation Goal: Model heat (thermal energy) flow in a one-dimensional object
More informationMath 251 December 14, 2005 Answer Key to Final Exam. 1 18pt 2 16pt 3 12pt 4 14pt 5 12pt 6 14pt 7 14pt 8 16pt 9 20pt 10 14pt Total 150pt
Name Section Math 51 December 14, 5 Answer Key to Final Exam There are 1 questions on this exam. Many of them have multiple parts. The point value of each question is indicated either at the beginning
More informationu tt = a 2 u xx u tt = a 2 (u xx + u yy )
10.7 The wave equation 10.7 The wave equation O. Costin: 10.7 1 This equation describes the propagation of waves through a medium: in one dimension, such as a vibrating string u tt = a 2 u xx 1 This equation
More informationSuggested Solution to Assignment 6
MATH 4 (6-7) partial diferential equations Suggested Solution to Assignment 6 Exercise 6.. Note that in the spherical coordinates (r, θ, φ), Thus, Let u = v/r, we get 3 = r + r r + r sin θ θ sin θ θ +
More informationTHE WAVE EQUATION. F = T (x, t) j + T (x + x, t) j = T (sin(θ(x, t)) + sin(θ(x + x, t)))
THE WAVE EQUATION The aim is to derive a mathematical model that describes small vibrations of a tightly stretched flexible string for the one-dimensional case, or of a tightly stretched membrane for the
More informationReview Sol. of More Long Answer Questions
Review Sol. of More Long Answer Questions 1. Solve the integro-differential equation t y (t) e t v y(v)dv = t; y()=. (1) Solution. The key is to recognize the convolution: t e t v y(v) dv = e t y. () Now
More information1. Solve the boundary-value problems or else show that no solutions exist. y (x) = c 1 e 2x + c 2 e 3x. (3)
Diff. Eqns. Problem Set 6 Solutions. Solve the boundary-value problems or else show that no solutions exist. a y + y 6y, y, y 4 b y + 9y x + e x, y, yπ a Assuming y e rx is a solution, we get the characteristic
More informationMcGill University April 20, Advanced Calculus for Engineers
McGill University April 0, 016 Faculty of Science Final examination Advanced Calculus for Engineers Math 64 April 0, 016 Time: PM-5PM Examiner: Prof. R. Choksi Associate Examiner: Prof. A. Hundemer Student
More informationElementary theory of L p spaces
CHAPTER 3 Elementary theory of L saces 3.1 Convexity. Jensen, Hölder, Minkowski inequality. We begin with two definitions. A set A R d is said to be convex if, for any x 0, x 1 2 A x = x 0 + (x 1 x 0 )
More informationBessel s Equation. MATH 365 Ordinary Differential Equations. J. Robert Buchanan. Fall Department of Mathematics
Bessel s Equation MATH 365 Ordinary Differential Equations J. Robert Buchanan Department of Mathematics Fall 2018 Background Bessel s equation of order ν has the form where ν is a constant. x 2 y + xy
More informationLecture 6. 2 Recurrence/transience, harmonic functions and martingales
Lecture 6 Classification of states We have shown that all states of an irreducible countable state Markov chain must of the same tye. This gives rise to the following classification. Definition. [Classification
More informationThere are five problems. Solve four of the five problems. Each problem is worth 25 points. A sheet of convenient formulae is provided.
Preliminary Examination (Solutions): Partial Differential Equations, 1 AM - 1 PM, Jan. 18, 16, oom Discovery Learning Center (DLC) Bechtel Collaboratory. Student ID: There are five problems. Solve four
More informationPartial Differential Equations Separation of Variables. 1 Partial Differential Equations and Operators
PDE-SEP-HEAT-1 Partial Differential Equations Separation of Variables 1 Partial Differential Equations and Operators et C = C(R 2 ) be the collection of infinitely differentiable functions from the plane
More informationSpecial Instructions:
Be sure that this examination has 20 pages including this cover The University of British Columbia Sessional Examinations - December 2016 Mathematics 257/316 Partial Differential Equations Closed book
More informationMATH 124B: HOMEWORK 2
MATH 24B: HOMEWORK 2 Suggested due date: August 5th, 26 () Consider the geometric series ( ) n x 2n. (a) Does it converge pointwise in the interval < x
More informationPartial Differential Equations Summary
Partial Differential Equations Summary 1. The heat equation Many physical processes are governed by partial differential equations. temperature of a rod. In this chapter, we will examine exactly that.
More informationMATH 2710: NOTES FOR ANALYSIS
MATH 270: NOTES FOR ANALYSIS The main ideas we will learn from analysis center around the idea of a limit. Limits occurs in several settings. We will start with finite limits of sequences, then cover infinite
More information10.2-3: Fourier Series.
10.2-3: Fourier Series. 10.2-3: Fourier Series. O. Costin: Fourier Series, 10.2-3 1 Fourier series are very useful in representing periodic functions. Examples of periodic functions. A function is periodic
More informationProblem set 3: Solutions Math 207B, Winter Suppose that u(x) is a non-zero solution of the eigenvalue problem. (u ) 2 dx, u 2 dx.
Problem set 3: Solutions Math 27B, Winter 216 1. Suppose that u(x) is a non-zero solution of the eigenvalue problem u = λu < x < 1, u() =, u(1) =. Show that λ = (u ) 2 dx u2 dx. Deduce that every eigenvalue
More informationLecture6. Partial Differential Equations
EP219 ecture notes - prepared by- Assoc. Prof. Dr. Eser OĞAR 2012-Spring ecture6. Partial Differential Equations 6.1 Review of Differential Equation We have studied the theoretical aspects of the solution
More information1. Partial differential equations. Chapter 12: Partial Differential Equations. Examples. 2. The one-dimensional wave equation
1. Partial differential equations Definitions Examples A partial differential equation PDE is an equation giving a relation between a function of two or more variables u and its partial derivatives. The
More informationPhysics 250 Green s functions for ordinary differential equations
Physics 25 Green s functions for ordinary differential equations Peter Young November 25, 27 Homogeneous Equations We have already discussed second order linear homogeneous differential equations, which
More informationSums of independent random variables
3 Sums of indeendent random variables This lecture collects a number of estimates for sums of indeendent random variables with values in a Banach sace E. We concentrate on sums of the form N γ nx n, where
More informationMathematics of Physics and Engineering II: Homework problems
Mathematics of Physics and Engineering II: Homework problems Homework. Problem. Consider four points in R 3 : P (,, ), Q(,, 2), R(,, ), S( + a,, 2a), where a is a real number. () Compute the coordinates
More information= =5 (0:4) 4 10 = = = = = 2:005 32:4 2: :
MATH LEC SECOND EXAM THU NOV 0 PROBLEM Part (a) ( 5 oints ) Aroximate 5 :4 using a suitable dierential. Show your work carrying at least 6 decimal digits. A mere calculator answer will receive zero credit.
More information