I am currently applying a short-term Fourier transform to some electrophysiological signal (variable 'my_signal', in mV/m^2). I would like to get the spectral amplitude in the same unit (mV/m^2/Hz). I did:
f, t, Zxx = sg.stft(my_signal, fs = sf, window = 'boxcar', nperseg = window_size_step, noverlap = (window_size_step-hop_size),return_onesided = True)
amp_spectrum = np.abs(Zxx)
I have a doubt on the unit of variable 'amp_spectrum'. Do I need to apply some additional scaling?