Previous Next

Van der Pol Oscillator Simulation (No Noise)

[x_out] = vanderpol_sim( muxin_outstep)

Differential Equation
We use  x_1 (t) and  x_2 (t) to denote the oscillator position and velocity as a function of time. The ordinary differential equation for the Van der Pol oscillator with no noise satisfies the differential equation  \[
      x_1 '(t) & = & x_2 (t)
      x_2 '(t) & = & \mu [ 1 - x_1(t)^2 ] x_2 (t) - x_1(t)

Is a scalar specifying the value of  \mu is the differential equation above.

is a column vector with two elements specifying the initial value for  x(t) \in \B{R}^2  . To be specific,  x_1 (0) is equal to xi(1) and  x_2(0) is equal to xi(2) .

is an integer scalar that specifies the number of time points at which the approximate solution to the ODE is returned.

is a scalar that specifies the difference in time between points at which the solution of the ODE is approximated.

is a matrix with row size two and column size n_out that contains the approximation solution to the ODE. To be specific, for k = 1 , ... , n_out , x_out(i,k) is an approximation for  x_i [ (k-1) \Delta t ] .

A fourth-order Runge method with step size step is used to approximate the solution of the ODE.

The file vanderpol_sim_ok.m is an example and test of vanderpol_sim. It returns true, if the test passes, and false otherwise.
Input File: example/nonlinear/vanderpol_sim.m