Signal Processing, Physics, Mathematics.
np.abs(za) ** 2 # Filter beta beta = signal.firwin(61, [19, 20], pass_zero=False, fs=fs) yb = signal.convolve(F, beta, mode='same') zb = np.fft.fft(yb) PSb = np.abs(zb) ** 2 # RPR calculation numd = np.sum(PSd, axis=0) numt = np.sum(PSt, axis=0) numa = np.sum(PSa, axis=0) numb = np.sum(PSb, axis=0)