power_spectrum binary statt rate

This commit is contained in:
Diana 2024-10-23 13:31:58 +02:00
parent e8dc6d6aa1
commit 940a461978

View File

@ -273,7 +273,7 @@ def power_spectrum(stimulus):
# computes firing rates # computes firing rates
rate = firing_rate(binary, dt = dt) rate = firing_rate(binary, dt = dt)
# creates power spectrum # creates power spectrum
freq, power = welch(rate, fs = 1/dt, nperseg = 2**16, noverlap = 2**15) freq, power = welch(binary, fs = 1/dt, nperseg = 2**16, noverlap = 2**15)
return freq, power return freq, power
def prepare_harmonic(frequency, category, num_harmonics, color): def prepare_harmonic(frequency, category, num_harmonics, color):