0
我试图通过RTL-SDR获得特定频率的功率。我正在修改我在网上找到的FFT示例。这里缩写代码(删除多余的东西):在特定频率下获得功率,SDR和FFT
import matplotlib.mlab as mlab
import rtlsdr
NFFT=1024
dwell = 0.016
sample_rate = 2.4e6
offset = 200e3
freq = 100e6
sdr = rtlsdr.RtlSdr(0)
sdr.set_sample_rate(sample_rate)
sdr.set_manual_gain_enabled(1)
sdr.set_gain(22.9)
sdr.freq_correction = 0
numsamples = next_2_to_pow(int(dwell * sample_rate))
freq = freq - offset # avoid dc spike
sdr.set_center_freq(freq)
samples = sdr.read_samples(numsamples)
powers, freqs = mlab.psd(samples, NFFT=NFFT, Fs=sample_rate/1e6, window=hamming(NFFT))
bin_offset = int(offset/(sample_rate/NFFT))
freq = int(freq + offset)
pwr = float("{:.2f}".format(10 * math.log10(powers[ int(len(powers)/2) + bin_offset ])))
这根本行不通。即使在我注入信号时,频率上的功率也保持不变(本底噪声电平)。
我有两个理论为什么这不起作用。 1)RTL-SDR返回I/Q数据,这种方法没有考虑到(?)2)我对FFT的理解听起来不够完善。
这是什么?我该如何解决它?