Except for the external driving source, this circuit is similar to the circuits considered in a previous post. To write state equations, we model the neon-bulb as a nonlinear resistor in parallel with an small linear inductor and a small linear capacitor.

When the sinusoidal source is turned off, the oscillator oscillates freely with a frequency $\omega_1$. The addition of the sinusoidal voltage source add another frequency $\omega_2$ to the system. One can think of the Poincaré map (generated by sampling the flow at frequency $\omega_2$) as a map on the circle, the circle being the limit cycle. In fact, for certain parameters, we see behavior similar to those obtained from the circle map.

The complete state equations of the circuit is given by:

\[ \begin{array}{lcl} \frac{dv}{dt} & = &

-\frac{1}{RC}v - \frac{1}{C}i + \frac{E}{RC} \\

\frac{di}{dt} & = & \frac{1}{L}v - \frac{f(i)}{L}

- \frac{E_0 \sin(\omega t)}{L}

\end{array}

\]

By adjusting the capacitor $C$, we change the free running frequency of the oscillator, i.e. changing the ratio between the two frequency. We observe frequency locking where the oscillator oscillates at a submultiple of the driving frequency for a range of capacitor values. If we vary the driving frequency, we observe the Devil's staircase, similar to the one observed in the circle map.

**References:**

Van der Pol, B. and Van der Mark, J., “Frequency demultiplication”, Nature, 120, 363-364, 1927.

Kennedy, M. P. & Chua, L. O., “Van der Pol and chaos,” IEEE Trans. Circuits Syst. CAS-33, 974–980, 1986.

Kennedy, M. P., Krieg, K. R. & Chua, L. O., “The Devil’s staircase: The electrical engineer’s fractal,” IEEE Trans. Circuits Syst. CAS-36, 1133–1139, 1989.