fit_harmonic_oscillator

gala.dynamics.fit_harmonic_oscillator(orbit, omega0=[1.0, 1.0, 1.0])

Fit the toy harmonic oscillator potential to the sum of the energy residuals relative to the mean energy by minimizing the function

\[\begin{split}f(\boldsymbol{\omega}) = \sum_i (\frac{1}{2}v_i^2 + \Phi_{\rm sho}(x_i\,|\,\boldsymbol{\omega}) - <E>)^2\end{split}\]

Parameters:

orbit : CartesianOrbit

omega0 : array_like (optional)

Initial frequency guess.

Returns:

omegas : float

Best-fit harmonic oscillator frequencies.