Signal Processing, Physics, Mathematics.
theta theta = signal.firwin(61, [8, 9], pass_zero=False, fs=fs) yt = signal.convolve(F, theta, mode='same') zt = np.fft.fft(yt) PSt = np.abs(zt) ** 2 # Filter alpha alpha = signal.firwin(61, [12, 13], pass_zero=False, fs=fs) ya = signal.convolve(F, alpha, mode='same') za = np.fft.fft(ya) PSa =