Page 151 - 6734

import matplotlib.pyplot as plt

def deriv(y,t):
m=1.0; j=1.0; c=0.1
return [y[1], (-j*y[0]-c*y[1])/m] # повертає
значення функцій dy/dt та dy'/dt
t = np.linspace(0.0, 10.0, 100) # час
yinit = np.array([1.0, 0.0]) # початкові умови для
`y` та `y'`
y = odeint(deriv, yinit, t) # інтегруємо систему диф.
рівнянь, повертає двовимірний масив зі значеннями `y`
та `y'`
plt.plot(t, y[:,0], 'k-', t, y[:,1], 'k--') #
переміщення і швидкість
plt.xlabel("t");plt.ylabel("y,
y'");plt.grid();plt.show()

Рисунок 25 – Функції переміщення y (-) і швидкості y’ (–)

scipy.interpolate – інтерполяція
Інтерполяція – це спосіб знаходження проміжних значень
величини за її відомим дискретним набором значень [9, 45, 53]. Для
інтерполяції і апроксимації сплайнами застосовують функції з
150