I am trying to implement adaptive hopf oscillator.
i have to implement these dynamics eqution
dx/dt = (mu - r^2)x - wy +KF(t)
dy/dt = (mu - r^2)y + wx
where , mu = 1 and initial condition of r(0)= 1 (r = root(x^2 + y^2))
F(t) is time periodic perturbation means any input signal and K is coupling strength , w is intrinsic frequency of oscillator
learning rule for this adaptation is
dw/dt = -KF(t) * [y/(root (x^2+y^2))]