Spectral Analysis, Seismology, Time Frequency Analysis, Signal Processing, Data Analysis.
over each taper for k in range(K): # Apply Welch Periodogram on each tapered data f, Pxx = welch(tapered_data[k, :], fs, nperseg=nperseg, noverlap=noverlap) # Accumulate the power spectrum average_power_spectrum[:len(f)] += 10 * np.log10(Pxx) # Average the power spectrum