[Python: signal cross correlation (signal/phase comparison)] #cross #correlation #signals #numpy #xcorr #np

sig_x_corr.py

import numpy as np

t = np.linspace(0, 4*np.pi, 100)
y1 = np.sin(t)
randomizer = lambda x: x + np.random.choice([-1, 1])*np.random.rand()
# y2 = np.apply_along_axis(randomizer, 0, y1) # applies once on whole array instead of per element of array
vRandomizer = np.vectorize(randomizer)
y2 = vRandomizer(y1)

corr = np.corrcoef(y1, y2)[0, 1] # choose a value from the left --> right diagonal; either [0, 1] or [1, 0]