Spectral Analysis, Seismology, Time Frequency Analysis, Signal Processing, Data Analysis.
average_power_spectrum[:len(f)] /= K # Visualize fig, ax = plt.subplots(1, 1, figsize=(10, 6), sharex=True) ax.plot(f, average_power_spectrum[:len(f)]) ax.set_xlabel('Frequency (Hz)') ax.set_ylabel('Power/Frequency (dB/Hz)') ax.set_title('Multitaper Power Spectrum using Welch\'s Method')