G.6)  Damped oscillators, undamped oscillators and van der Pol's nonlinear oscillator*

G.6.a) Undamped oscillators

The differential equation of the undamped oscillator is
         [Graphics:Images/index_gr_1.gif]  with  [Graphics:Images/index_gr_2.gif].
Here is one such:

[Graphics:Images/index_gr_3.gif]
[Graphics:Images/index_gr_4.gif]

To make this second order diffeq into a system of two first order diffeq's, you put
             [Graphics:Images/index_gr_5.gif]
and then replace [Graphics:Images/index_gr_6.gif] with [Graphics:Images/index_gr_7.gif] and replace [Graphics:Images/index_gr_8.gif] with [Graphics:Images/index_gr_9.gif]:

[Graphics:Images/index_gr_10.gif]
[Graphics:Images/index_gr_11.gif]
[Graphics:Images/index_gr_12.gif]

Clean this up:

[Graphics:Images/index_gr_13.gif]
[Graphics:Images/index_gr_14.gif]
[Graphics:Images/index_gr_15.gif]

Now look at:

[Graphics:Images/index_gr_16.gif]

[Graphics:Images/index_gr_17.gif]

Throw in a trajectory:

[Graphics:Images/index_gr_18.gif]

[Graphics:Images/index_gr_19.gif]

Look at the corresponding solution [Graphics:Images/index_gr_20.gif]:

[Graphics:Images/index_gr_21.gif]

[Graphics:Images/index_gr_22.gif]

Another:

[Graphics:Images/index_gr_23.gif]

[Graphics:Images/index_gr_24.gif]

And the corresponding solution [Graphics:Images/index_gr_25.gif]:

[Graphics:Images/index_gr_26.gif]

[Graphics:Images/index_gr_27.gif]

Compare the trajectories on the same axes:

[Graphics:Images/index_gr_28.gif]

[Graphics:Images/index_gr_29.gif]

Compare the solution curves [Graphics:Images/index_gr_30.gif]:

[Graphics:Images/index_gr_31.gif]

[Graphics:Images/index_gr_32.gif]

Play with other starting points and then describe how solutions of the undamped oscillator
          [Graphics:Images/index_gr_33.gif]
look.
Why do you think folks call these things undamped oscillators?

G.6.b) Damped oscillators

Here is the diffeq of an undamped oscillator as studied above:

[Graphics:Images/index_gr_34.gif]
[Graphics:Images/index_gr_35.gif]

To get a damped version of this oscillator, you throw in a shock absorber term
            [Graphics:Images/index_gr_36.gif]
with [Graphics:Images/index_gr_37.gif] positive but not too large:

[Graphics:Images/index_gr_38.gif]
[Graphics:Images/index_gr_39.gif]

To make this second order differential equation into a system of two first order differential equations, you put
             [Graphics:Images/index_gr_40.gif]
and then replace [Graphics:Images/index_gr_41.gif] with [Graphics:Images/index_gr_42.gif] and replace [Graphics:Images/index_gr_43.gif] with [Graphics:Images/index_gr_44.gif]:

[Graphics:Images/index_gr_45.gif]
[Graphics:Images/index_gr_46.gif]
[Graphics:Images/index_gr_47.gif]

Clean this up:

[Graphics:Images/index_gr_48.gif]
[Graphics:Images/index_gr_49.gif]
[Graphics:Images/index_gr_50.gif]

Now look at:

[Graphics:Images/index_gr_51.gif]

[Graphics:Images/index_gr_52.gif]

Throw in a trajectory:

[Graphics:Images/index_gr_53.gif]

[Graphics:Images/index_gr_54.gif]

Look at the corresponding solution [Graphics:Images/index_gr_55.gif]:

[Graphics:Images/index_gr_56.gif]

[Graphics:Images/index_gr_57.gif]

Another:

[Graphics:Images/index_gr_58.gif]

[Graphics:Images/index_gr_59.gif]

Look at the corresponding solution [Graphics:Images/index_gr_60.gif]:

[Graphics:Images/index_gr_61.gif]

[Graphics:Images/index_gr_62.gif]

Compare the trajectories on the same axes:

[Graphics:Images/index_gr_63.gif]

[Graphics:Images/index_gr_64.gif]

Look at the two solution curves:

[Graphics:Images/index_gr_65.gif]

[Graphics:Images/index_gr_66.gif]

Play with other starting points and then describe how solutions of the damped oscillator
          [Graphics:Images/index_gr_67.gif]
with [Graphics:Images/index_gr_68.gif] small and positive, and [Graphics:Images/index_gr_69.gif] positive look.
Why do you think folks call these things damped oscillators?

G.6.d.i) The nonlinear van der Pol heartbeat oscillator

You can get an undamped oscillator this way:

[Graphics:Images/index_gr_70.gif]
[Graphics:Images/index_gr_71.gif]

You can get a damped oscillator by making [Graphics:Images/index_gr_72.gif] positive and small:

[Graphics:Images/index_gr_73.gif]
[Graphics:Images/index_gr_74.gif]

You can get a negatively damped (self-excited) oscillator by making [Graphics:Images/index_gr_75.gif] negative and small:

[Graphics:Images/index_gr_76.gif]
[Graphics:Images/index_gr_77.gif]

Here's how you can get a van der Pol oscillator:

[Graphics:Images/index_gr_78.gif]
[Graphics:Images/index_gr_79.gif]

When  [Graphics:Images/index_gr_80.gif], do you expect the van der Pol oscillator to behave like an undamped, a damped or a negatively damped oscillator?

When [Graphics:Images/index_gr_81.gif] or [Graphics:Images/index_gr_82.gif], do you expect the van der Pol oscillator to behave like an undamped, a damped or a negatively damped oscillator?

G.6.d.ii)

Another look at a van der Pol oscillator:

[Graphics:Images/index_gr_83.gif]
[Graphics:Images/index_gr_84.gif]

To make this second order differential equation into a system of two first order differential equations, you put
             [Graphics:Images/index_gr_85.gif]
and then replace [Graphics:Images/index_gr_86.gif] by [Graphics:Images/index_gr_87.gif] and replace [Graphics:Images/index_gr_88.gif] by [Graphics:Images/index_gr_89.gif]:

[Graphics:Images/index_gr_90.gif]
[Graphics:Images/index_gr_91.gif]
[Graphics:Images/index_gr_92.gif]

Clean this up:

[Graphics:Images/index_gr_93.gif]
[Graphics:Images/index_gr_94.gif]
[Graphics:Images/index_gr_95.gif]

Now look at:

[Graphics:Images/index_gr_96.gif]
[Graphics:Images/index_gr_97.gif]

[Graphics:Images/index_gr_98.gif]

Cowabunga.
Throw in a trajectory:

[Graphics:Images/index_gr_99.gif]

[Graphics:Images/index_gr_100.gif]

Look at the corresponding solution [Graphics:Images/index_gr_101.gif]:

[Graphics:Images/index_gr_102.gif]

[Graphics:Images/index_gr_103.gif]

Just like heartbeats in E.R.
Another:

[Graphics:Images/index_gr_104.gif]

[Graphics:Images/index_gr_105.gif]

And the corresponding solution [Graphics:Images/index_gr_106.gif]:

[Graphics:Images/index_gr_107.gif]

[Graphics:Images/index_gr_108.gif]

Compare the two trajectories on the same axes:

[Graphics:Images/index_gr_109.gif]

[Graphics:Images/index_gr_110.gif]

Folks like to call this a "limit cycle."

See the two solution curves:

[Graphics:Images/index_gr_111.gif]

[Graphics:Images/index_gr_112.gif]

Heartbeats.
Look carefully at the flow plot and try to see how the flow forces the limit cycle to happen.


Converted by Mathematica      November 4, 1999