The differential equation of the undamped oscillator is
with
.
Here is one such:
![[Graphics:Images/index_gr_3.gif]](Images/index_gr_3.gif)
To make this second order diffeq into a system of two first order diffeq's, you put
and then replace
with
and replace
with
:
![[Graphics:Images/index_gr_10.gif]](Images/index_gr_10.gif)
|
|
|
|
Clean this up:
![[Graphics:Images/index_gr_13.gif]](Images/index_gr_13.gif)
|
|
|
|
Now look at:
![[Graphics:Images/index_gr_16.gif]](Images/index_gr_16.gif)
![[Graphics:Images/index_gr_17.gif]](Images/index_gr_17.gif)
Throw in a trajectory:
![[Graphics:Images/index_gr_18.gif]](Images/index_gr_18.gif)
![[Graphics:Images/index_gr_19.gif]](Images/index_gr_19.gif)
Look at the corresponding solution
:
![[Graphics:Images/index_gr_21.gif]](Images/index_gr_21.gif)
![[Graphics:Images/index_gr_22.gif]](Images/index_gr_22.gif)
Another:
![[Graphics:Images/index_gr_23.gif]](Images/index_gr_23.gif)
![[Graphics:Images/index_gr_24.gif]](Images/index_gr_24.gif)
And the corresponding solution
:
![[Graphics:Images/index_gr_26.gif]](Images/index_gr_26.gif)
![[Graphics:Images/index_gr_27.gif]](Images/index_gr_27.gif)
Compare the trajectories on the same axes:
![[Graphics:Images/index_gr_28.gif]](Images/index_gr_28.gif)
![[Graphics:Images/index_gr_29.gif]](Images/index_gr_29.gif)
Compare the solution curves
:
![[Graphics:Images/index_gr_31.gif]](Images/index_gr_31.gif)
![[Graphics:Images/index_gr_32.gif]](Images/index_gr_32.gif)
Play with other starting points and then describe how solutions of the undamped oscillator
look.
Why do you think folks call these things undamped oscillators?
Here is the diffeq of an undamped oscillator as studied above:
![[Graphics:Images/index_gr_34.gif]](Images/index_gr_34.gif)
To get a damped version of this oscillator, you throw in a shock absorber term
with
positive but not too large:
![[Graphics:Images/index_gr_38.gif]](Images/index_gr_38.gif)
To make this second order differential equation into a system of two first order differential equations, you put
and then replace
with
and replace
with
:
![[Graphics:Images/index_gr_45.gif]](Images/index_gr_45.gif)
|
|
|
|
Clean this up:
![[Graphics:Images/index_gr_48.gif]](Images/index_gr_48.gif)
|
|
|
|
Now look at:
![[Graphics:Images/index_gr_51.gif]](Images/index_gr_51.gif)
![[Graphics:Images/index_gr_52.gif]](Images/index_gr_52.gif)
Throw in a trajectory:
![[Graphics:Images/index_gr_53.gif]](Images/index_gr_53.gif)
![[Graphics:Images/index_gr_54.gif]](Images/index_gr_54.gif)
Look at the corresponding solution
:
![[Graphics:Images/index_gr_56.gif]](Images/index_gr_56.gif)
![[Graphics:Images/index_gr_57.gif]](Images/index_gr_57.gif)
Another:
![[Graphics:Images/index_gr_58.gif]](Images/index_gr_58.gif)
![[Graphics:Images/index_gr_59.gif]](Images/index_gr_59.gif)
Look at the corresponding solution
:
![[Graphics:Images/index_gr_61.gif]](Images/index_gr_61.gif)
![[Graphics:Images/index_gr_62.gif]](Images/index_gr_62.gif)
Compare the trajectories on the same axes:
![[Graphics:Images/index_gr_63.gif]](Images/index_gr_63.gif)
![[Graphics:Images/index_gr_64.gif]](Images/index_gr_64.gif)
Look at the two solution curves:
![[Graphics:Images/index_gr_65.gif]](Images/index_gr_65.gif)
![[Graphics:Images/index_gr_66.gif]](Images/index_gr_66.gif)
Play with other starting points and then describe how solutions of the damped oscillator
with
small and positive, and
positive look.
Why do you think folks call these things damped oscillators?
You can get an undamped oscillator this way:
![[Graphics:Images/index_gr_70.gif]](Images/index_gr_70.gif)
You can get a damped oscillator by making
positive and small:
![[Graphics:Images/index_gr_73.gif]](Images/index_gr_73.gif)
You can get a negatively damped (self-excited) oscillator by making
negative and small:
![[Graphics:Images/index_gr_76.gif]](Images/index_gr_76.gif)
Here's how you can get a van der Pol oscillator:
![[Graphics:Images/index_gr_78.gif]](Images/index_gr_78.gif)
When
, do you expect the van der Pol oscillator to behave like an undamped, a damped or a negatively damped oscillator?
When
or
, do you expect the van der Pol oscillator to behave like an undamped, a damped or a negatively damped oscillator?
Another look at a van der Pol oscillator:
![[Graphics:Images/index_gr_83.gif]](Images/index_gr_83.gif)
To make this second order differential equation into a system of two first order differential equations, you put
and then replace
by
and replace
by
:
![[Graphics:Images/index_gr_90.gif]](Images/index_gr_90.gif)
|
|
|
|
Clean this up:
![[Graphics:Images/index_gr_93.gif]](Images/index_gr_93.gif)
|
|
|
|
Now look at:
![[Graphics:Images/index_gr_96.gif]](Images/index_gr_96.gif)
![[Graphics:Images/index_gr_98.gif]](Images/index_gr_98.gif)
Cowabunga.
Throw in a trajectory:
![[Graphics:Images/index_gr_99.gif]](Images/index_gr_99.gif)
![[Graphics:Images/index_gr_100.gif]](Images/index_gr_100.gif)
Look at the corresponding solution
:
![[Graphics:Images/index_gr_102.gif]](Images/index_gr_102.gif)
![[Graphics:Images/index_gr_103.gif]](Images/index_gr_103.gif)
Just like heartbeats in E.R.
Another:
![[Graphics:Images/index_gr_104.gif]](Images/index_gr_104.gif)
![[Graphics:Images/index_gr_105.gif]](Images/index_gr_105.gif)
And the corresponding solution
:
![[Graphics:Images/index_gr_107.gif]](Images/index_gr_107.gif)
![[Graphics:Images/index_gr_108.gif]](Images/index_gr_108.gif)
Compare the two trajectories on the same axes:
![[Graphics:Images/index_gr_109.gif]](Images/index_gr_109.gif)
![[Graphics:Images/index_gr_110.gif]](Images/index_gr_110.gif)
Folks like to call this a "limit cycle."
See the two solution curves:
![[Graphics:Images/index_gr_111.gif]](Images/index_gr_111.gif)
![[Graphics:Images/index_gr_112.gif]](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.