pyqtgraph_plot_sweep.py

#!/usr/bin/env python

#####################################################################
## This example makes use of sweep_device.py to perform
## a single 100 MHz sweep monitoring and plot
##
## To perform a more sophisticated sweep setup with many entries
## for different frequency location, look into SweepEntry() of util.py
#####################################################################

# import required libraries
from pyqtgraph.Qt import QtGui, QtCore
import pyqtgraph as pg
import sys
import numpy as np
from pyrf.devices.thinkrf import WSA
from pyrf.sweep_device import SweepDevice

# Constants for configuration
RFE_MODE = 'SH'
START_FREQ = 2350e6
STOP_FREQ = 3050e6
SPP = 16384
PPB = 1
RBW = 125e6 / (SPP * PPB)  # 125 MHz is the sampling rate
ATTENUATION = 0
GAIN = 'HIGH'

# Setup the plotting window
class MainApplication(pg.GraphicsWindow):

    def __init__(self, dut):
        super(MainApplication, self).__init__()
        self.dut = dut

    # use ';' to do SCPI command on the run
    def keyPressEvent(self, event):
        if event.text() == ';':
            cmd, ok = QtGui.QInputDialog.getText(win, 'Enter SCPI Command',
                        'Enter SCPI Command:')
            if ok:
                if '?' not in cmd:
                    dut.scpiset(cmd)

# initialize an RTSA (aka WSA) device handle
dut = WSA()
win = MainApplication(dut)

# get device's IP and connect
if len(sys.argv) > 1:
    ip = sys.argv[1]
else:
    ip, ok = QtGui.QInputDialog.getText(win, 'Open Device',
                'Enter a hostname or IP address:')
dut.connect(ip)

# initialize RTSA configurations
dut.reset()
dut.request_read_perm()
dut.psfm_gain(GAIN)

# setup a sweep device
sweepdev = SweepDevice(dut)

# initialize the plot
win.resize(1000, 600)
win.setWindowTitle("PYRF FFT Plot Example - " + ip)
fft_plot = win.addPlot(title="Power Vs. Frequency")
fft_plot.enableAutoRange('xy', True)
curve = fft_plot.plot(pen='g')
fft_plot.setLabel('bottom', text='Frequency', units='Hz', unitPrefix=None)

# initialize the y-axis of the plot
plot_ymin = -130
plot_ymax = 10
fft_plot.setYRange(plot_ymin, plot_ymax)
fft_plot.setLabel('left', text='Power', units='dBm', unitPrefix=None)

# Get data and plot
def update():
    global dut, curve, fft_plot, fstart, fstop, print_once

    # get sweep spectral data of a given RBW and plot
    fstart, fstop, spectra_data = sweepdev.capture_power_spectrum(START_FREQ,
                                STOP_FREQ,
                                RBW,
                                {'attenuator':ATTENUATION},
                                mode = RFE_MODE)
    freq_range = np.linspace(fstart , fstop, len(spectra_data))
    curve.setData(freq_range, spectra_data, pen='g')

timer = QtCore.QTimer()
timer.timeout.connect(update)
timer.start(0)


## Start Qt event loop unless running in interactive mode or using pyside.
if __name__ == '__main__':
    import sys
    if (sys.flags.interactive != 1) or not hasattr(QtCore, 'PYQT_VERSION'):
        QtGui.QApplication.instance().exec_()