21.04.2015 Views

Lecture Notes in Differential Equations - Bruce E. Shapiro

Lecture Notes in Differential Equations - Bruce E. Shapiro

Lecture Notes in Differential Equations - Bruce E. Shapiro

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

214 LESSON 25. HARMONIC OSCILLATIONS<br />

Figure 25.1: Harmonic oscillations with ω = 1 and <strong>in</strong>itial conditions y(0) =<br />

1, y ′ (0) = 0. (a) simple harmonic motion; (b) under-damped system with<br />

b = 0.25; (c) critically damped system with b = 2.0; (d) under-damped<br />

system with b = 0.025.<br />

Forced Oscillations<br />

F<strong>in</strong>ally, it is possible to imag<strong>in</strong>g add<strong>in</strong>g a motor to the spr<strong>in</strong>g that produces<br />

a force f(t). The result<strong>in</strong>g is system, <strong>in</strong>clud<strong>in</strong>g damp<strong>in</strong>g, called the forced<br />

harmonic oscillator:<br />

my ′′ = −Cy ′′ − ky + f(t) (25.17)<br />

If we def<strong>in</strong>e a force per unit mass F (t) = f(t)/m and ω and the parameter<br />

b as before, this becomes, <strong>in</strong> standard form,<br />

y ′′ + by ′ + ω 2 y = F (t) (25.18)<br />

This is the general second-order l<strong>in</strong>ear equation with constant coefficients<br />

that are positive. Thus any second order l<strong>in</strong>ear differential equation with<br />

positive constant coefficients describes a harmonic oscillator of some sort.<br />

The particular solution is<br />

∫ ∫<br />

y P = e r2t e (r1−r2)u e −r1s F (s)ds (25.19)<br />

t<br />

where r 1 and r 2 are the roots of the characteristic equation. For example,<br />

suppose the system driven by a force function<br />

Then<br />

u<br />

F (t) = F 0 s<strong>in</strong> αt (25.20)<br />

∫ ∫<br />

y P =F 0 e r2t e (r1−r2)u e −r1s s<strong>in</strong> αs dsdu (25.21)<br />

t<br />

u<br />

= F ∫<br />

0<br />

α 2 + r1<br />

2 e r2t e −r2u (α cos αu + r 1 s<strong>in</strong> αu)du (25.22)<br />

t

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!