text: #!/usr/bin/env python2
  # -*- coding: utf-8 -*-
  ##################################################
  # GNU Radio Python Flow Graph
  # Title: Power spectral density and average power
  # Author: Alexandros-Apostolos A. Boulogeorgos
  # Generated: Tue Oct  8 06:33:48 2019
  ##################################################
  
  if __name__ == '__main__':
      import ctypes
      import sys
      if sys.platform.startswith('linux'):
          try:
              x11 = ctypes.cdll.LoadLibrary('libX11.so')
              x11.XInitThreads()
          except:
              print "Warning: failed to XInitThreads()"
  
  from PyQt4 import Qt
  from gnuradio import blocks
  from gnuradio import digital
  from gnuradio import eng_notation
  from gnuradio import filter
  from gnuradio import gr
  from gnuradio import qtgui
  from gnuradio.eng_option import eng_option
  from gnuradio.filter import firdes
  from optparse import OptionParser
  import sip
  import sys
  
  
  class power_spectral_density_and_average_power(gr.top_block, Qt.QWidget):
  
      def __init__(self):
          gr.top_block.__init__(self, "Power spectral density and average power")
          Qt.QWidget.__init__(self)
          self.setWindowTitle("Power spectral density and average power")
          try:
              self.setWindowIcon(Qt.QIcon.fromTheme('gnuradio-grc'))
          except:
              pass
          self.top_scroll_layout = Qt.QVBoxLayout()
          self.setLayout(self.top_scroll_layout)
          self.top_scroll = Qt.QScrollArea()
          self.top_scroll.setFrameStyle(Qt.QFrame.NoFrame)
          self.top_scroll_layout.addWidget(self.top_scroll)
          self.top_scroll.setWidgetResizable(True)
          self.top_widget = Qt.QWidget()
          self.top_scroll.setWidget(self.top_widget)
          self.top_layout = Qt.QVBoxLayout(self.top_widget)
          self.top_grid_layout = Qt.QGridLayout()
          self.top_layout.addLayout(self.top_grid_layout)
  
          self.settings = Qt.QSettings("GNU Radio", "power_spectral_density_and_average_power")
          self.restoreGeometry(self.settings.value("geometry").toByteArray())
  
          ##################################################
          # Variables
          ##################################################
          self.sr = sr = 32000
          self.samp_rate = samp_rate = 512000
          self.interpolation = interpolation = 16
          self.glfsrSourceDegree = glfsrSourceDegree = 22
  
          ##################################################
          # Blocks
          ##################################################
          self.qtgui_time_sink_x_0 = qtgui.time_sink_f(
          	1024, #size
          	samp_rate, #samp_rate
          	"", #name
          	1 #number of inputs
          )
          self.qtgui_time_sink_x_0.set_update_time(0.10)
          self.qtgui_time_sink_x_0.set_y_axis(-1, 1)
          
          self.qtgui_time_sink_x_0.set_y_label('Amplitude', "")
          
          self.qtgui_time_sink_x_0.enable_tags(-1, True)
          self.qtgui_time_sink_x_0.set_trigger_mode(qtgui.TRIG_MODE_FREE, qtgui.TRIG_SLOPE_POS, 0.0, 0, 0, "")
          self.qtgui_time_sink_x_0.enable_autoscale(False)
          self.qtgui_time_sink_x_0.enable_grid(False)
          self.qtgui_time_sink_x_0.enable_axis_labels(True)
          self.qtgui_time_sink_x_0.enable_control_panel(False)
          
          if not True:
            self.qtgui_time_sink_x_0.disable_legend()
          
          labels = ['', '', '', '', '',
                    '', '', '', '', '']
          widths = [1, 1, 1, 1, 1,
                    1, 1, 1, 1, 1]
          colors = ["blue", "red", "green", "black", "cyan",
                    "magenta", "yellow", "dark red", "dark green", "blue"]
          styles = [1, 1, 1, 1, 1,
                    1, 1, 1, 1, 1]
          markers = [-1, -1, -1, -1, -1,
                     -1, -1, -1, -1, -1]
          alphas = [1.0, 1.0, 1.0, 1.0, 1.0,
                    1.0, 1.0, 1.0, 1.0, 1.0]
          
          for i in xrange(1):
              if len(labels[i]) == 0:
                  self.qtgui_time_sink_x_0.set_line_label(i, "Data {0}".format(i))
              else:
                  self.qtgui_time_sink_x_0.set_line_label(i, labels[i])
              self.qtgui_time_sink_x_0.set_line_width(i, widths[i])
              self.qtgui_time_sink_x_0.set_line_color(i, colors[i])
              self.qtgui_time_sink_x_0.set_line_style(i, styles[i])
              self.qtgui_time_sink_x_0.set_line_marker(i, markers[i])
              self.qtgui_time_sink_x_0.set_line_alpha(i, alphas[i])
          
          self._qtgui_time_sink_x_0_win = sip.wrapinstance(self.qtgui_time_sink_x_0.pyqwidget(), Qt.QWidget)
          self.top_grid_layout.addWidget(self._qtgui_time_sink_x_0_win, 0,0,1,1)
          self.qtgui_number_sink_0 = qtgui.number_sink(
              gr.sizeof_float,
              0,
              qtgui.NUM_GRAPH_VERT,
              1
          )
          self.qtgui_number_sink_0.set_update_time(0.10)
          self.qtgui_number_sink_0.set_title("Average Power")
          
          labels = ['', '', '', '', '',
                    '', '', '', '', '']
          units = ['', '', '', '', '',
                   '', '', '', '', '']
          colors = [("black", "black"), ("black", "black"), ("black", "black"), ("black", "black"), ("black", "black"),
                    ("black", "black"), ("black", "black"), ("black", "black"), ("black", "black"), ("black", "black")]
          factor = [1, 1, 1, 1, 1,
                    1, 1, 1, 1, 1]
          for i in xrange(1):
              self.qtgui_number_sink_0.set_min(i, -1)
              self.qtgui_number_sink_0.set_max(i, 1)
              self.qtgui_number_sink_0.set_color(i, colors[i][0], colors[i][1])
              if len(labels[i]) == 0:
                  self.qtgui_number_sink_0.set_label(i, "Data {0}".format(i))
              else:
                  self.qtgui_number_sink_0.set_label(i, labels[i])
              self.qtgui_number_sink_0.set_unit(i, units[i])
              self.qtgui_number_sink_0.set_factor(i, factor[i])
          
          self.qtgui_number_sink_0.enable_autoscale(True)
          self._qtgui_number_sink_0_win = sip.wrapinstance(self.qtgui_number_sink_0.pyqwidget(), Qt.QWidget)
          self.top_grid_layout.addWidget(self._qtgui_number_sink_0_win, 0,1,2,1)
          self.qtgui_freq_sink_x_0 = qtgui.freq_sink_f(
          	1024, #size
          	firdes.WIN_BLACKMAN_hARRIS, #wintype
          	0, #fc
          	samp_rate, #bw
          	"", #name
          	1 #number of inputs
          )
          self.qtgui_freq_sink_x_0.set_update_time(0.10)
          self.qtgui_freq_sink_x_0.set_y_axis(-140, 10)
          self.qtgui_freq_sink_x_0.set_y_label('Power', 'dB')
          self.qtgui_freq_sink_x_0.set_trigger_mode(qtgui.TRIG_MODE_FREE, 0.0, 0, "")
          self.qtgui_freq_sink_x_0.enable_autoscale(False)
          self.qtgui_freq_sink_x_0.enable_grid(False)
          self.qtgui_freq_sink_x_0.set_fft_average(1.0)
          self.qtgui_freq_sink_x_0.enable_axis_labels(True)
          self.qtgui_freq_sink_x_0.enable_control_panel(False)
          
          if not True:
            self.qtgui_freq_sink_x_0.disable_legend()
          
          if "float" == "float" or "float" == "msg_float":
            self.qtgui_freq_sink_x_0.set_plot_pos_half(not True)
          
          labels = ['', '', '', '', '',
                    '', '', '', '', '']
          widths = [1, 1, 1, 1, 1,
                    1, 1, 1, 1, 1]
          colors = ["blue", "red", "green", "black", "cyan",
                    "magenta", "yellow", "dark red", "dark green", "dark blue"]
          alphas = [1.0, 1.0, 1.0, 1.0, 1.0,
                    1.0, 1.0, 1.0, 1.0, 1.0]
          for i in xrange(1):
              if len(labels[i]) == 0:
                  self.qtgui_freq_sink_x_0.set_line_label(i, "Data {0}".format(i))
              else:
                  self.qtgui_freq_sink_x_0.set_line_label(i, labels[i])
              self.qtgui_freq_sink_x_0.set_line_width(i, widths[i])
              self.qtgui_freq_sink_x_0.set_line_color(i, colors[i])
              self.qtgui_freq_sink_x_0.set_line_alpha(i, alphas[i])
          
          self._qtgui_freq_sink_x_0_win = sip.wrapinstance(self.qtgui_freq_sink_x_0.pyqwidget(), Qt.QWidget)
          self.top_grid_layout.addWidget(self._qtgui_freq_sink_x_0_win, 1,0,1,1)
          self.low_pass_filter_1 = filter.fir_filter_fff(1, firdes.low_pass(
          	1, samp_rate, 10, 10, firdes.WIN_HAMMING, 6.76))
          self.low_pass_filter_0 = filter.fir_filter_fff(1, firdes.low_pass(
          	1, samp_rate, 64000, 3200, firdes.WIN_HAMMING, 6.76))
          self.interp_fir_filter_xxx_0 = filter.interp_fir_filter_fff(interpolation, ((1, 1, 1, 1, 1, 0, 0, 0, 0, 0,1, 1, 1, 1, 1, 0, 0, 0, 0, 0)))
          self.interp_fir_filter_xxx_0.declare_sample_delay(0)
          self.digital_glfsr_source_x_0 = digital.glfsr_source_f(glfsrSourceDegree, True, 0, 1)
          self.blocks_throttle_0 = blocks.throttle(gr.sizeof_float*1, sr,True)
          self.blocks_multiply_xx_0 = blocks.multiply_vff(1)
  
          ##################################################
          # Connections
          ##################################################
          self.connect((self.blocks_multiply_xx_0, 0), (self.low_pass_filter_1, 0))    
          self.connect((self.blocks_throttle_0, 0), (self.interp_fir_filter_xxx_0, 0))    
          self.connect((self.digital_glfsr_source_x_0, 0), (self.blocks_throttle_0, 0))    
          self.connect((self.interp_fir_filter_xxx_0, 0), (self.low_pass_filter_0, 0))    
          self.connect((self.low_pass_filter_0, 0), (self.blocks_multiply_xx_0, 0))    
          self.connect((self.low_pass_filter_0, 0), (self.blocks_multiply_xx_0, 1))    
          self.connect((self.low_pass_filter_0, 0), (self.qtgui_freq_sink_x_0, 0))    
          self.connect((self.low_pass_filter_0, 0), (self.qtgui_time_sink_x_0, 0))    
          self.connect((self.low_pass_filter_1, 0), (self.qtgui_number_sink_0, 0))    
  
      def closeEvent(self, event):
          self.settings = Qt.QSettings("GNU Radio", "power_spectral_density_and_average_power")
          self.settings.setValue("geometry", self.saveGeometry())
          event.accept()
  
      def get_sr(self):
          return self.sr
  
      def set_sr(self, sr):
          self.sr = sr
          self.blocks_throttle_0.set_sample_rate(self.sr)
  
      def get_samp_rate(self):
          return self.samp_rate
  
      def set_samp_rate(self, samp_rate):
          self.samp_rate = samp_rate
          self.qtgui_time_sink_x_0.set_samp_rate(self.samp_rate)
          self.qtgui_freq_sink_x_0.set_frequency_range(0, self.samp_rate)
          self.low_pass_filter_1.set_taps(firdes.low_pass(1, self.samp_rate, 10, 10, firdes.WIN_HAMMING, 6.76))
          self.low_pass_filter_0.set_taps(firdes.low_pass(1, self.samp_rate, 64000, 3200, firdes.WIN_HAMMING, 6.76))
  
      def get_interpolation(self):
          return self.interpolation
  
      def set_interpolation(self, interpolation):
          self.interpolation = interpolation
  
      def get_glfsrSourceDegree(self):
          return self.glfsrSourceDegree
  
      def set_glfsrSourceDegree(self, glfsrSourceDegree):
          self.glfsrSourceDegree = glfsrSourceDegree
  
  
  def main(top_block_cls=power_spectral_density_and_average_power, options=None):
  
      from distutils.version import StrictVersion
      if StrictVersion(Qt.qVersion()) >= StrictVersion("4.5.0"):
          style = gr.prefs().get_string('qtgui', 'style', 'raster')
          Qt.QApplication.setGraphicsSystem(style)
      qapp = Qt.QApplication(sys.argv)
  
      tb = top_block_cls()
      tb.start()
      tb.show()
  
      def quitting():
          tb.stop()
          tb.wait()
      qapp.connect(qapp, Qt.SIGNAL("aboutToQuit()"), quitting)
      qapp.exec_()
  
  
  if __name__ == '__main__':
      main()
  