def simulate_data(N):

    t = np.arange(N)
    a = 0.1 * 2*np.pi  # base frequency
    b = 0.01 * 2*np.pi # frequency of the frequency change
    c = 8              # magnitude of the frequency change
    f = np.sin( a * (t + c*np.sin(b*t)) )
    y = f + 0.1*np.random.randn(N)

    return (y, f)
Returns the angular units.
dataset/ETHPy150Open bayespy/bayespy/bayespy/demos/lssm_tvd.py/simulate_data--crossover with--dataset/ETHPy150Open django/django/django/contrib/gis/db/backends/base/models.py/SpatialRefSysMixin.angular_units