Ordinary Differential Equations Higher-order ODEs Problem Set 5 (Version 1.0) Wayne Hacker Copyright ©Wayne Hacker 2007, 2014. All rights reserved. higher-order odes, problem set 5 Copyright ©Wayne Hacker 2007. All rights reserved. 1 Contents 1 Higher-order ODEs 1.1 2 Definitions, theorems, and terminology for higher-order ODEs . . . . . . 2 1.1.1 Determining when functions are dependent or independent . . . . 2 1.1.2 Verifying solutions to higher-order equations . . . . . . . . . . . . 4 1.2 Differential operators . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5 1.3 Second-order linear homogeneous ODEs . . . . . . . . . . . . . . . . . . 7 1.4 Reduction of order for second-order ODEs . . . . . . . . . . . . . . . . . 8 1.5 Higher-order linear homogeneous ODEs . . . . . . . . . . . . . . . . . . . 9 1.6 Cauchy-Euler equations . . . . . . . . . . . . . . . . . . . . . . . . . . . 10 1.7 The method of undetermined coefficients . . . . . . . . . . . . . . . . . . 11 1.8 Variation of parameters . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12 1.9 Applications: harmonic oscillating systems . . . . . . . . . . . . . . . . . 13 1.9.1 Harmonic motion w/o friction and external forcing . . . . . . . . 13 1.9.2 Damped harmonic motion w/o external forcing . . . . . . . . . . 14 1.9.3 Forced harmonic motion w/o damping . . . . . . . . . . . . . . . 16 higher-order odes, problem set 5 1 1.1 1.1.1 Copyright ©Wayne Hacker 2007. All rights reserved. 2 Higher-order ODEs Definitions, theorems, and terminology for higher-order ODEs Determining when functions are dependent or independent For the following problems determine whether the functions given in the form: {f1 (x), f2 (x)} or {f1 (x), f2 (x), f3 (x)} are independent or dependent on the entire real line (−∞, ∞). Problem 1. Let f1 (x) = x and f2 (x) = |x|. Problem 2. Let f1 (x) = 1, f2 (x) = x, and f3 (x) = x2 . Problem 3. Let f1 (x) = 1, f2 (x) = cos2 x, and f3 (x) = sin2 x. Problem 4. Let f1 (x) = x, f2 (x) = 2x, and f3 (x) = x + 3. Problem 5. Let f1 (x) = 1, f2 (x) = cos(2x), and f3 (x) = cos2 x. Problem 6. Let f1 (x) = ex , f2 (x) = e−x , and f3 (x) = cosh x. higher-order odes, problem set 5 Copyright ©Wayne Hacker 2007. All rights reserved. 3 For the following problems use the Wronskian to determine whether the functions given in the form: {f1 (x), f2 (x)} or {f1 (x), f2 (x), f3 (x)} are independent or dependent on the given interval. Problem 7. Let f1 (x) = 1, f2 (x) = x, and f3 (x) = x2 . Problem 8. Let f1 (x) = x1/2 and f2 (x) = x on (0, ∞). Problem 9. Let f1 (x) = cos x and f2 (x) = sin x on (−∞, ∞). Problem 10. Let f1 (x) = eαx and f2 (x) = xeαx on (−∞, ∞). Problem 11. (i) Show graphically that f1 (x) = x2 and f2 (x) = x|x| are linearly independent on (−∞, ∞). (ii) Even though |x| is not differentiable, you can still compute the Wronskian on the intervals (0, ∞) and (−∞, 0). What result does this give? Does it contradict the result in part (i)? Problem 12. [When are a function and its derivative independent] Let f (x) be a smooth function of x. It is often advantageous to know when a function and its derivative are independent. Show that f (x) and f 0 (x) are independent so long as f is not a solution of the differential equation f f 00 − (f 0 )2 = 0 . higher-order odes, problem set 5 1.1.2 Copyright ©Wayne Hacker 2007. All rights reserved. 4 Verifying solutions to higher-order equations Problem 13. Verify that y = c1 + c2 cos x + c3 sin x is a three-parameter family of solutions of the ODE: y 000 + y 0 = 0 on the interval (−∞, ∞). Then find a member of the family satisfying the initial conditions: y(π) = 0, y 0 (π) = 2, and y 00 (π) = −1. Problem 14. Find the largest possible interval about x = 0 for which theorem 5.1 guarantees a unique solution to the following IVP: y 00 + (tan x) y = ex , y(0) = 1, y 0 (0) = 0 . Problem 15. Determine the values of λ for which the family of curves y = c1 cos(λx) + c2 sin(λx) form a nontrivial solution to the boundary-value problem (BVP): y 00 + λy = 0 y(0) = 0, y(π) = 0 . Such a problem where only certain parameters will satisfy boundary conditions is known as an eigenvalue problem. Problem 16. For the ODE x2 y 00 + xy 0 + y = 0, show that {cos(ln x), sin(ln x)} is a fundamental solution set on the interval (−∞, ∞). Problem 17. [a nonlinear equation] a) Verify that y1 (x) = 1 and y2 (x) = ln x are solutions of the nonlinear differential equation y 00 + (y 0 )2 = 0 . (b) Do the solutions satisfy the superposition principle? In particular, is y1 + y2 also a solution to the ODE? higher-order odes, problem set 5 1.2 Copyright ©Wayne Hacker 2007. All rights reserved. 5 Differential operators dy + 8y = x + 3 in the form L(y) = f (x), dx d where L is a differential operator with constant coefficients and D = dx . Solution: All we have to do is change notation: Problem 1. Rewrite the differential equation 4 4Dy + 8y = x + 3 factor from the right −−−−−−−−−−−−−−−→ compress notation −−−−−−−−−−−−−→ (4D + 8) y = x + 3 L(y) = 1 − sin x, where L = 4D + 8 Problem 2. Rewrite the differential equation y (4) − 2y 00 + y = e−3x + e2x in the form d . L(y) = f (x), where L is a differential operator with constant coefficients and D = dx Solution: All we have to do is change notation: D4 y − 2D2 y + y = e−3x + e2x factor from the right −−−−−−−−−−−−−−−→ (D4 − 2D2 + 1)y = e−3x + e2x compress notation −−−−−−−−−−−−−→ L(y) = e−3x + e2x , where L = D4 − 2D2 + 1 Problem 3. Factor the differential operator 2D2 − 3D − 2. Solution: 2D2 − 3D − 2 = (2D + 1)(D − 2) Problem 4. Factor the differential operator D3 + 4D. Solution: D3 + 4D = D(D2 + 4). This won’t factor any farther over the reals. Problem 5. Factor the differential operator D3 + 4D2 + 3D. Solution: D3 + 4D2 + 3D = D(D2 + 4D + 3) = D(D + 1)(D + 3). Problem 6. Factor the differential operator D4 − 8D2 + 16. Solution: D4 − 8D2 + 16 = (D2 )2 − 2 · 4 · (D2 ) + 42 = (D2 − 4)2 = ((D − 2)(D + 2))2 = (D − 2)2 (D + 2)2 Problem 7. Show that D4 annihilates y = 10x3 − 2x. higher-order odes, problem set 5 Copyright ©Wayne Hacker 2007. All rights reserved. 6 Solution: D4 (10x3 − 2x) = D3 (D(10x3 ) − 2Dx) = D3 (30x2 − 2) = D2 (30Dx2 − 2D(1)) = D2 ((30 · 2x − 0) = 60D2 x = 60 · 2D(1) = 120 · 0 = 0 Problem 8. Show that (D − 2)(D + 5) annihilates y = 4e2x . Solution: (D − 2)(D + 5)(4e2x ) = 4(D + 5)(D − 2)e2x (We’ve used the fact that differential operators with constant coefficients commute) = 4(D + 5)(De2x − 2e2x ) = 4(D + 5)(2e2x − 2e2x ) = 4(D + 5)(0) = 0 Remember that since Dn+1 annihilates xn for every n ∈ {0, 1, 2, . . .}, Dn+1 must likewise annihilate any nth-degree polynomial Pn (x): D n+1 Pn (x) = D n+1 n X k=0 k ak x = n X k=0 ak D n+1 k x = n X ak · 0 = 0 k=0 Problem 9. Find a differential operator that annihilates the given function x3 (1 − 5x). Solution: Write: x3 (1 − 5x) = x3 − 5x4 . This is a 4th-degree polynomial, so D5 annihilates it. Problem 10. Find a differential operator that annihilates the given function 8x−sin x+ 10 cos(5x). Solution: Look at this term by term: D2 (8x) = 0 (D2 + 1) sin x = − sin x + sin x = 0 (D2 + 52 ) cos(5x) = −52 cos(5x) + 52 cos(5x) = 0 The operators D2 + 1 and D2 + 52 don’t factor over the real numbers. Thus: (D2 + 52 )(D2 + 1)D2 annihilates 8x − sin x + 10 cos(5x) and is the operator of lowest possible order that does so: if we removed any of the “prime” factors, then the operator would no longer annihilate the function. Copyright ©Wayne Hacker 2007. All rights reserved. 7 higher-order odes, problem set 5 Problem 11. Find linearly independent functions that are annihilated by the given operator D2 + 4D. Solution: First, we’ll factor the operator: D2 + 4D = D(D + 4) This is a second-order operator, so we’ll look for two independent solutions: one that is annihilated by D, the other by D + 4. Dy1 = 0 ⇒ (D + 4)y2 = 0 y1 = constant = 1 ⇒ y2 = e−4x 1 and e−4x are linearly independent. Problem 12. Find linearly independent functions that are annihilated by the given operator D2 (D − 5)(D − 7). Solution: This operator is already factored into prime factors. Each of these factors will annihilate one of the independent functions {y1 , y2 , y3 , y4 }. The functions y1 and y2 are independent functions that are both annihilated by D2 ; we’ll refer to them collectively as y1,2 . D2 y1,2 = 0 (D − 5)y3 = 0 (D − 7)y4 = 0 ⇒ ⇒ ⇒ y1,2 = ax + b y3 = e5x y4 = e7x ⇒ y1 = 1; y2 = x We have a 4th-order operator, so we have 4 independent solutions: {1, x, e5x , e7x } 1.3 Second-order linear homogeneous ODEs Problem 1. Find the general solution to the ODE y 00 − y = 0. Problem 2. Find the general solution to the ODE y 00 + y = 0. Problem 3. Find the general solution to the ODE y 00 − y 0 − 12y = 0. Problem 4. Find the general solution to the ODE y 00 − 2y 0 + 5y = 0. higher-order odes, problem set 5 1.4 Copyright ©Wayne Hacker 2007. All rights reserved. 8 Reduction of order for second-order ODEs Problem 1. Given that y1 (x) = xe−x is a solution to the ODE y 00 + 2y 0 + y = 0, find a second linearly independent solution to the ODE (do not use the guess method: y2 = emx ). Problem 2. Consider the ODE y 00 +2y 0 +y = 0. Since this is a constant-coefficient linear homogeneous ODE the standard approach is to seek a solution of the form: y = emx . The resulting characteristic equation is m2 + 2m + 1 = (m + 1)2 = 0 ⇒ m = 1 (a double root). However, this method only gives one solution. One short cut is to differentiate w.r.t. the parameter. However, a more straight-forward method, interns of understanding, is to use the method of reduction of order to find the second solution. Find the second independent solution using this method and verify its independence. Problem 3. Given that y1 (x) = x sin(ln x) is a solution to the ODE x2 y 00 − xy 0 + 2y = 0, find a second linearly independent solution to the ODE. Problem 4. Given that y1 (x) = 1 is a solution to the homogeneous component of the ODE y 00 +y 0 = 1, find a second linearly independent solution to the ODE and a particular solution. Write down the general solution. higher-order odes, problem set 5 1.5 Copyright ©Wayne Hacker 2007. All rights reserved. 9 Higher-order linear homogeneous ODEs Problem 1. Solve for the general solution to the ODE y 000 + 3y 00 + 2y 0 = 0 . Problem 2. Solve for the general solution to the ODE y 000 + 6y 00 + 11y 0 + 6y = 0 . Problem 3. Solve for the general solution to the ODE y 000 + 2y 00 − y 0 − 2y = 0 . Problem 4. Solve for the general solution to the ODE y (4) − y 000 + y 00 − y 0 = 0 . higher-order odes, problem set 5 1.6 Copyright ©Wayne Hacker 2007. All rights reserved. 10 Cauchy-Euler equations Problem 1. Solve the IVP x2 y 00 − 3xy 0 + 4y = 0, y(1) = 5, y 0 (1) = 3. Problem 2. Solve the IVP x2 y 00 − 4xy 0 + 6y = 0, y(−2) = 8, y 0 (−2) = 0 by transforming the equation onto the interval (0, ∞). Problem 3. Find the general solution to the inhomogeneous ODE x2 y 00 − 4xy 0 + 6y = 2 ln x by using the substitution x = et , or equivalently, t = ln x. higher-order odes, problem set 5 1.7 Copyright ©Wayne Hacker 2007. All rights reserved. 11 The method of undetermined coefficients Note: When solving inhomogeneous equations always solve for the homogeneous solution first, since the relation between the fundamental set to the homogeneous equation and the form of the inhomogeneous term determines what method we us to seek a particular solution to the inhomogeneous equation. If the inhomogeneous term is not independent from the fundamental set, then you must modify your particular solution appropriately. Problem 1. Solve for the general solution to the inhomogeneous ODE ay 00 + by = c , where a, b ∈ R+ and c ∈ R. Problem 2. Solve for the general solution to the inhomogeneous ODE y 00 + y 0 − 6y = 2x . Problem 3. Solve for the general solution to the inhomogeneous ODE y 00 + 4y = x2 − 3 sin(2x) . higher-order odes, problem set 5 1.8 Copyright ©Wayne Hacker 2007. All rights reserved. 12 Variation of parameters Problem 1. Find the general solution to the linear inhomogeneous equation y 00 + 2y 0 + y = e−x ln x on (0, ∞). Problem 2. Find the general solution to the linear inhomogeneous equation 3 1 2 00 0 2 x y + xy + x − y = x 2 on (0, ∞). 4 Problem 3. Find the general solution to the linear inhomogeneous equation x2 y 00 + xy 0 + y = sec(ln x) on (0, ∞). higher-order odes, problem set 5 1.9 1.9.1 Copyright ©Wayne Hacker 2007. All rights reserved. 13 Applications: harmonic oscillating systems Harmonic motion w/o friction and external forcing Problem 1. A 20-kg mass is attached to one end of a hanging spring. The other end of the spring is attached to the ceiling of your physics lab. The frequency of the simpleharmonic motion is observed to be 2/π, what is the spring constant? You may ignore the mass of the spring compared to the 20-kg mass. Problem 2. When a 400 newton weight is attached to a massless spring, the spring is stretched by 2 meters. A mass of 50 kg is attached to the end of the spring and released from the equilibrium position with an downward velocity 10 m/s. If we take the x-axis in the downward direction, what is the equation of motion? Problem 3. Use the alternative form of the solution x(t) = A sin (ωt + φ) to the ODE x¨ + ω 2 x = 0 to show that the time interval between two successive maxima of x(t) is 2π/ω. higher-order odes, problem set 5 1.9.2 Copyright ©Wayne Hacker 2007. All rights reserved. 14 Damped harmonic motion w/o external forcing We can add damping to a mass-spring system by immersing it in a viscous fluid. If we assume that the attached mass m is much greater than the mass of the spring mspring and the resistance is proportional to the velocity of the mass, then the governing equation with the initial conditions are of the form: x¨ + 2λx˙ + ω 2 x = 0, x(0) = x0 , x(0) ˙ = v0 , where λ is the constant of proportionality for the frictional term. If we let x(t) = emt , then the characteristic equation is m2 + 2λm + ω 2 = 0 with roots: √ m = −λ ± λ2 − ω 2 . Fact: The solution to the IVP in the case of underdamped motion (i.e., λ2 − ω 2 < 0) is of the form: √ ω 2 − λ2 t + φ . x(t) = Ae−λt sin This is known as the alternative form. The quasi-period is defined to be the roots of the equation, which can be shown to be 2π . Tperiod = √ 2 ω − λ2 higher-order odes, problem set 5 Copyright ©Wayne Hacker 2007. All rights reserved. 15 Problem 4. Solve the IVP x¨ + 2x˙ + 10x = 0, x(0) = −2, x(0) ˙ = 0. Express the solution in the alternative form. Problem 5. For the case of underdamped motion, show that the ratio between any two 2πλ consecutive maxima is exp √ω2 −λ2 . Problem 6. In the case of over-damped motion (i.e., λ2 −ω 2 > 0), show that the solution can be expressed in the form √ x(t) = Ae−λt sinh ω 2 − λ2 t + φ . higher-order odes, problem set 5 1.9.3 Copyright ©Wayne Hacker 2007. All rights reserved. 16 Forced harmonic motion w/o damping Problem 7. Solve the following initial value problem analytically x¨ + ω 2 x = AF sin(ΩF t), x(0) = 0, x(0) ˙ = v0 , where ω, ΩF , and AF are all positive constants with ω 6= ΩF (no resonance). Problem 8. Graph the solution to the initial value problem for v0 = ±1, ±2 over the time interval [0, 4π]. x¨ + x = sin(2t), x(0) = 0, x(0) ˙ = v0 . (i) How does the sign of v0 affect the solution? (ii) How does the magnitude of v0 affect the solution? (iii) Compare you numerical results here to the previous solution. higher-order odes, problem set 5 Copyright ©Wayne Hacker 2007. All rights reserved. 17 Problem 9. (near-resonance) In the notes it was shown that the solution to the IVP x¨ + ω 2 x = AF sin(ΩF t), x(0) = 0, x(0) ˙ = 0, with ω 6= ΩF was given by AF x(t) = ω AF x(t) = 2 ω " ΩF sin(ω t) − ω sin(ΩF t) Ω2F − ω 2 , which can be recast as ΩF ω sin(ω t) − sin(ΩF t) ΩF 2 −1 ω # . (1.1) We can study the case of near resonance by letting ΩF = (1 + )ω in the above solution and examine the solution for small (i.e., 0 < 1). Graph the solution for AF = 1, ω = 1, and = 1/10, 1/100, 1/1000 over the time interval [0, 10π]. What trend do you notice in the solutions x(t; ) as → 0+ ? Can you justify observation by using a Taylor expansion of the solution x(t; ) about = 0 with t fixed? Problem 10. (A variable frequency) Solve the initial value problem x¨ + x = sin(ΩF (t) · t), x(0) = 0, x(0) ˙ =0 √ over the time interval [0, 40π] for ΩF (t) = 1/ t. Since the sine function has a variable argument, you will need to solve the equation numerically using either Mathematica or Matlab. Also graph the forcing function sin(ΩF (t) · t) on the same graph. Is there any striking features that you notice between the response curve (i.e., the solution curve) and the forcing curve?
© Copyright 2024