显然，这是一个求导即可解决的问题。如果将时间记作t，将旋转角度记作θ，那么θ就是t的函数：

对上式求导，就得到角速度v对时间t的函数：

1. 寻找原函数

>>> import numpy as np
>>> from scipy import interpolate
>>> _x = np.linspace(192,199,8)
>>> _y = np.array([7.085,10.497,14.019,17.683,21.487,25.403,29.402,33.467])
>>> f = interpolate.interp1d(_x, _y, kind='cubic')

>>> f(193.3)
array(11.54026091)
>>> f(193.35)
array(11.71517478)
>>> f(193.351)
array(11.71867619)

2. 求解导函数

>>> def get_derivative(f, delta=1e-10):
        """导函数生成器"""
        def derivative(x):
            """导函数"""
            return (f(x+delta)-f(x))/delta
        return derivative

>>> fd = get_derivative(f)

>>> fd(193.5)
3.590830143540664
>>> fd(195.8)
3.8961776076555026