Section 3.4. Second Order Nonhomogeneous Equations y + p(x)y + q(x)y = f(x) (N) The corresponding homogeneous equation y + p(x)y + q(x)y = 0 (H) is called the reduced equation of (N). 1
General Results THEOREM 1: If z = z 1 (x) and z = z 2 (x) are solutions of equation (N), then y(x) = z 1 (x) z 2 (x) is a solution of equation (H). 2
THEOREM 2: Let y = y 1 (x) and y = y 2 (x) be linearly independent solutions of the reduced equation (H) and let z = z(x) be a particular solution of (N). Then y(x) = C 1 y 1 (x) + C 2 y 2 (x) + z(x) is the general solution of (N). 3
The general solution of (N) consists of the general solution of the reduced equation (H) plus a particular solution of (N): y = C 1 y 1 (x) + C 2 y 2 (x) }{{} + z(x). }{{} gen soln (H) + part soln (N) 4
To find the general solution of (N) you need to find: (i) a fundamental set of solutions y 1, y 2 of the reduced equation (H), and (ii) a particular solution z of (N). 5
THEOREM: (Superposition Principle) If z = z f (x) and z = z g (x) are particular solutions of y + p(x)y + q(x)y = f(x) and y + p(x)y + q(x)y = g(x) respectively, then z(x) = z f (x) + z g (x) is a particular solution of y + p(x)y + q(x)y = f(x) + g(x). 6
Variation of Parameters Let y = y 1 (x) and y = y 2 (x) be independent solutions of the reduced equation (H) and let W(x) = W(x) = y 1 y 2 y 2y 1 be their Wronskian. Set z(x) = y 1 (x)u(x) + y 2 (x)v(x) where u and v are to be determined so that z is a solution of (N). 7
u (x) = y 2(x)f(x), W(x) u(x) = y 2 (x)f(x) W(x) dx; v (x) = y 1(x)f(x) W(x) v(x) = y 1 (x)f(x) W(x) dx and z(x) = y 1 (x) y 2 (x)f(x) W(x) dx+y 2 (x) y 1 (x)f(x) W(x) dx 8
Examples: 1. {y 1 (x) = x 2, y 2 (x) = x 4 } is a fundamental set of solutions of the reduced equation of y 5 x y + 8 x2y = 4x3 Find a particular solution z of the equation. 9
y 1 = x 2, y 2 = x 4 W[y 1, y 2 ] = 10
2. Find the general solution of y 4y + 4y = e2x x 11
3. Find the general solution of y + y 6y = 3e 2x 12
Section 3.5. Undetermined Coefficients aka Guessing Find a particular solution of y + ay + by = f(x) NOTE: This method can be used only when: 1. constant coefficients 2. f is an exponential function 13
If f(x) = ce rx set z(x) = Ae rx Example: Find a particular solution of y 5y + 6y = 7e 4x. Set z = Ae 4x where A is to be determined. 14
Answer: z = 1 6 e 4x. The general solution of the differential equation is: y = C 1 e 2x + C 2 e 3x + 1 6 e 4x. 15
If f(x) = c cos βx, d sin βx, or c cos βx + d sin βx, set z(x) = A cos βx + B sin βx where A, B are to be determined 16
Example: Find a particular solution of y 2y + y = 3cos 2x. Set z = Acos 2x + B sin 2x 17
Answer: z = 9 25 cos 2x 12 25 sin 2x. The general solution of the differential equation is: y = C 1 e x +C 2 xe x 9 25 cos 2x 12 25 sin 2x. 18
Example: Find a particular solution of y 2y +5y = 2cos 3x 4sin 3x+e 2x Set z = Acos 3x + B sin 3x + Ce 2x where A, B, C are to be determined. 19
y 2y +5y = 2cos 3x 4sin 3x+e 2x 20
Answer: z = 8 13 cos 3x + 1 13 sin 3x + 1 5 e2x. 21
If f(x) = ce αx cos βx, de αx sin βx or ce αx cos βx + de αx sin βx set z(x) = Ae αx cos βx + Be αx sin βx where A, B are to be determined. 22
Example: Find a particular solution of y + 9y = 4e x sin 2x. Set z = Ae x cos 2x + Be x sin 2x Answer: z = 4 13 ex cos 2x+ 6 13 ex sin 2x. 23
Example: Find a particular solution z of y + y 6y = 3e 2x. 24
A BIG Difficulty: The trial solution z is a solution of the reduced equation. Example 3, Section 3.4: z = 3 5 x e2x 25
Example: Find a particular solution z of y 6y + 9y = 5e 3x. 26
Example: Find a particular solution z of y + 4y = 2cos 2x. 27
Example: Find a particular solution z of y 5y + 6y = 4e 2x + 3 28
Example: Find a particular solution z of y 3y = 4e 3x 2 29
Answers: z 1 = 3 5 x e2x z 2 = 5 2 x2 e 3x z 3 = 1 2 xsin 2x z 4 = 4xe 2x + 1 2 z 5 = 4 3 x e3x 2 3 x 30
The Method of Undetermined Coefficients A. Applies only to equations of the form y + ay + by = f(x) where a, b are constants and f is an exponential function. c.f Variation of Parameters 31
B. Basic Case: If: f(x) = ae rx set z = Ae rx. f(x) = ccos βx, dsin βx, or ccos βx + dsin βx, set z = Acos βx + B sin βx. f(x) = ce αx cos βx, de αx sin βx or ce αx cos βx + de αx sin βx, set z = Ae αx cos βx + Be αx sin βx. 32
BUT: If z satisfies the reduced equation, try xz; if xz also satisfies the reduced equation, then x 2 z will give a particular solution. 33
C. General Case: If f(x) = p(x)e rx where p is a polynomial of degree n, then set z = P(x)e rx where P is a polynomial of degree n with undetermined coefficients. 34
Example: Find a particular solution of y y 6y = (2x 2 1)e 2x. Set z = (Ax 2 + Bx + C)e 2x. Answer: z = ( 1 2 x2 4 3 x 16 9 ) e 2x. 35
If f(x) = p(x)cos βx + q(x)sin βx where p, q are polynomials of degree n, then set z = P(x)cos βx + Q(x)sin βx where P, Q are polynomials of degree n with undetermined coefficients. 36
Example: y 2y 3y = 3 cos x + (x 2)sin x. Set z = (Ax+B)cos x+(cx+d)sin x. Answer: z = ( 1 10 x 47 ) ( 1 cos x 50 5 x 2 ) sin x. 25 37
If f(x) = p(x)e αx cos βx + q(x)e αx sin βx where p, q are polynomials of degree n, then set z = P(x)e αx cos βx + Q(x)e αx sin βx where P, Q are polynomials of degree n with undetermined coefficients. 38
Example: y + 4y = 2x e x cos x. Set z = (Ax+B)e x cos x+(cx+d)e x sin x Answer: z = 1 25 (10x 7)ex cos x+ 1 25 (5x 1)ex sin x. 39
BUT: Warning!!! If any part of z satisfies the reduced equation, try xz; if any part of xz also satisfies the reduced equation, then x 2 z will give a particular solution. 40
Examples: 1. Give the form of a particular solution of y 4y 5y = 2cos 3x 5e 5x + 4. Answer: z = Acos 3x + B sin 3x + Cxe 5x + D 41
Answer: z = Ax + B + Cx 2 e 4x 42 2. Give the form of a particular solution of y + 8y + 16y = 2x 1 + 7e 4x.
C cos 2x + Dsin 2x + Ee 2x 43 3. Give the form of a particular solution of y + y = 4sin x cos 2x + 2e 2x. Answer: z = Axcos x + Bxsin x +
4. Give the form of the general solution of y + 9y = 4cos 2x + 3sin 2x Answer: y = C 1 cos 3x + C 2 sin 3x + Acos 2x + B sin 2x 44
5. Give the form of a particular solution of y + 9y = 4cos 3x + 3sin 2x Answer: z = Axcos 3x + Bxsin 3x + C cos 2x + Dsin 2x 45
6. Give the form of the general solution of y + 4y + 4y = 4xe 2x + 3 Answer: y = C 1 e 2x + C 2 xe 2x + (Ax 3 + Bx 2 )e 2x + C 46
7. Give the form of the general solution of y + 4y + 4y = 4e 2x sin 2x + 3x Answer: y = C 1 e 2x + C 2 xe 2x + Ae 2x cos 2x + Be 2x sin 2x + Cx + D 47
8. Give the form of the general solution of y + 4y = 4sin 2x + 3 Answer: y = C 1 e 4x +C 2 +Acos 2x+B sin 2x+Cx 48
9. Give the form of a particular solution of y + 2y + 10y = 2e 3x sin x + 4e 3x Answer: z = Ae 3x cos x + Be 3x sin x + Ce 3x 49
10. Give the form of the general solution of y + 2y + 10y = 2e x sin 3x + 2e x Answer: y = C 1 e x cos 3x+C 2 e x sin 3x+ Axe x cos 3x + Bxe x sin 3x + Ce x 50
11. Give the form of a particular solution of y 2y 8y = 2cos 3x (3x+1)e 2x 4 Answer: z = Acos 3x+B sin 3x+(Cx 2 +Dx)e 2x +E 51
12. Give the form of a particular solution of y 2y 8y = 2cos 3x 3xe 2x 3x Answer: z = Acos 3x+B sin 3x+(Cx 2 +Dx)e 2x +Ex + F 52
13. Find the general solution of y 4y + 4y = 4e 2x + e2x x Answer: y = C 1 e 2x + C 2 xe 2x 2x 2 e 2x + xe 2x ln x 53
Summary: 1. Variation of parameters: Can be applied to any linear nonhomogeneous equations, but requires a fundamental set of solutions of the reduced equation. 54
2. Undetermined coefficients: Is limited to linear nonhomogeneous equations with constant coefficients, and f must be an exponential function, f(x) = ae rx, f(x) = ccos βx + dsin βx, f(x) = ce αx cos βx + de αx sin βx, or p(x)f(x) p a polynomial. 55
In cases where both methods are applicable, the method of undetermined coefficients is usually more efficient and, hence, the preferable method. 56
Section 3.6. Vibrating Mechanical Systems 57
I. Undamped, Free Vibrations Hooke s Law: The restoring force of a spring is proportional to the displacement: F = ky, k > 0. Newton s Second Law: Force equals mass times acceleration: F = ma = m d2 y dt 2. Mathematical model: m d2 y dt 2 = ky 58
which can be written where ω = k/m. d 2 y dt 2 + ω2 y = 0 ω/2π is called the natural frequency of the system.
The general solution of this equation is: y = C 1 sin ωt + C 2 cos ωt which can be written y = Asin(ωt + φ). A is called the amplitude, φ is called the phase shift. 59
y A t A
II. Damped, Free Vibrations: A resistance force R (e.g., friction) proportional to the velocity v = y and acting in a direction opposite to the motion: R = cy with c > 0. Force Equation: F = ky(t) cy (t). Newton s Second Law: F = ma = my 60
Mathematical Model: my (t) = ky(t) cy (t) or y + c m y + k m y = 0 (c, k, m constant)
Characteristic equation: r 2 + c m r + k m = 0. Roots r = c ± c 2 4km. 2m There are three cases to consider: c 2 4km < 0, c 2 4km > 0, c 2 4km = 0. 61
Case 1: c 2 4km < 0. Complex roots: (Underdamped) r 1 = c 2m + iω, r 2 = c 2m iω where ω = 4km c 2 2m. General solution: y = e ( c/2m)t (C 1 cos ωt + C 2 sin ωt) 62
or y(t) = A e ( c/2m)t sin(ωt + φ 0 ) where A and φ 0 are constants, NOTE: The motion is oscillatory AND y(t) 0 as t. 63
Underdamped Case: y t 64
Case 2: c 2 4km > 0. Two distinct real roots: (Overdamped) r 1 = c + c 2 4km, r 2 = c c 2 4km. 2m 2m General solution: y(t) = y = C 1 e r 1t + C 2 e r 2t. The motion is nonoscillatory. 65
NOTE: Since c 2 4km < c 2 = c, r 1 and r 2 are both negative and y(t) 0 as t. 66
Case 3: c 2 4km = 0. One real root: (Critically Damped) r 1 = c 2m, General solution: y(t) = y = C 1 e (c/2m) t +C 2 t e (c/2m) t. The motion is nonoscillatory and y(t) 0 as t. 67
Overdamped and Critically Damped Cases: y y t t y t 68
III. Forced Free Vibrations Apply an external force G to an undamped, freely vibrating system Force Equation: F = ky + G. Mathematical Model: my = ky+g or y + k m y = G m, a nonhomogeneous equation. 69
A periodic external force: G = α cos γt, α, γ > 0 const. Force Equation: F = ky + α cos γt Mathematical Model: y + k m y = α m cos γt y + ω 2 y = α m cos γt where ω = k/m. 70
ω/2π is called the natural frequency of the system, γ/2π is called the applied frequency. 71
Case 1: γ ω. y + ω 2 y = α m cos γt General solution, reduced equation: y = Asin(ωt + φ 0 ). Form of particular solution (undetermined coefficients): z = Acos γt + B sin γt. A particular solution: z = α/m ω 2 γ 2 cos γt. 72
General solution: y = Asin(ωt + φ 0 ) + α/m ω 2 cos γt. γ2
ω/γ rational: periodic motion 1 ω/γ irrational: not periodic y 1 t 73
Case 2: γ = ω. y + ω 2 y = α m cos ωt General solution, reduced equation: y = Asin(ωt + φ 0 ). Form of particular solution (undetermined coefficients): z = A tcos γt + B tsin γt. A particular solution: α 2ωm t sin ωt. 74
General solution: y = Asin(ωt + φ 0 ) + α 2ωm t sin ωt. Resonance: Unbounded oscillation!
Resonance: y = Asin(ωt + φ 0 ) + α 2ωm t sin ωt. y t 75
IV. Forced Damped Vibrations Apply an external force G to a damped, freely vibrating system Force Equation: F = ky cy + G. Mathematical Model: my = ky cy + G or y + c m y + k m y = G m, a nonhomogeneous equation. 76
A periodic external force: G = α cos γt, α, γ > 0 const. Force Equation: F = ky cy + αcos γt Mathematical Model: y + c m y + k m y = α m cos γt or y + α y + β y = α m cos γt where α = c/m, β = k/m. 77
General solution, reduced equation: Case I: distinct real roots: Y c (t) = C 1 e r 1t + C 2 e r 2t, r 1 < 0, r 2 < 0 Note: lim Y c (t) = 0. t 78
Case II: real, equal roots: Y c (t) = C 1 e rt + C 2 t e rt, r 1 = r 2 = r < 0 Note: lim Y c (t) = 0. t Case III: complex roots: Y c (t) = C 1 e rt cos ωt + C 2 e rt sin ωt, r 1 = r + ωi, r 2 = r ωi < 0, r < 0. Note: lim Y c (t) = 0. t 79
Particular solution of (N) y + α y + β y = α m cos γt will have the form: Z(t) = Acos γt + B sin γt. General solution of (N): y(t) = Y c (t) + Z(t) Note: lim y(t) = Z(t) t 80
Y c (t), the general solution of the reduced equation, is called the transient solution. Z(t) the particular solution of (N), is called the steady state solution. Example: y + 2y + 5y = cos t General solution y = C 1 e t cos 2t + C 2 e t sin 2t+ 1 5 cos t + 1 10 sin t 81
Steady state and transient solutions. y 3 2 1 1 5 x 1 y 3 2 1 1 5 x 1 82