f2 = (arctan(-x + 2) + pi / 2) / pi
L4 = int(len(f2) / 8)
sf2 = f2[L4:-L4]
c = correlate(f1 - mean(f1), sf2 - mean(f1), 'same')
print('peak correlation occurs at:', x[argmax(c)])
plt.plot(x, c)
plt.show()
Взаимная корреляция непериодической функции с NumPy