Report for course EGME 511 (Advanced Mechanical Vibration). California state university, Fullerton, CA. Van der Pol diﬀerential equation is given by
In this analysis, we will consider the case only for positive . The above equation will be solved numerically using Matlab's ODE45 for diﬀerent initial conditions, and the phase portrait (velocity vs. displacement) is plotted to show the limit cycle for diﬀerent initial conditions.
To use ODE45, one must ﬁrst convert the above second order ODE to two ODE's, each of which is ﬁrst order. Letting results in
The system of equations to be solved by ODE45 is the following
Subject to initial conditions and . In the Matlab implementation below, the values of and the initial conditions are deﬁned at the top of the code. This needs to be modiﬁed to change the initial conditions before running the program again.
The program was run for a number of diﬀerent initial conditions (diﬀerent and diﬀerent ). In all cases, was used. It is noticed that one limit cycle is reached in all cases. Both the and the phase plot where shown. The following are output from 3 diﬀerent runs made. The title on the plots show the initial conditions used.
A non-linear second order ODE was solved numerically using Matlab's ode45. The solution to the Van Der Pol was found to contain a limit cycle in the phase portrait when starting from any initial conditions.