Implement CrossSpectrum for unevenly sampled data?

[dhuppenkothen]

Maybe this would be a good idea. :)

[SamiraBhattacharya]

Hi, after some research, I have come up with this code. Do give me feedback on this.

from stingray import Lightcurve, Crossspectrum

# Generate some fake data
t1 = np.sort(np.random.uniform(0, 10, 100))
y1 = np.sin(2*np.pi*t1) + np.random.normal(0, 0.1, len(t1))
lc1 = Lightcurve(t1, y1)

t2 = np.sort(np.random.uniform(0, 10, 100))
y2 = np.sin(4*np.pi*t2) + np.random.normal(0, 0.1, len(t2))
lc2 = Lightcurve(t2, y2)

# Compute the cross spectrum
cs = Crossspectrum(lc1, lc2)
freq, power = cs.compute_crossspectrum(normalization='leahy')

# Plot the cross spectrum
import matplotlib.pyplot as plt
plt.plot(freq, power)
plt.xlabel('Frequency')
plt.ylabel('Power')
plt.show()

[matteobachetti]

Addressed by @pupperemeritus' #737 🥳