From 62e4fe1cac74a1055c20212f246d3547e976bd4f Mon Sep 17 00:00:00 2001 From: sara Date: Tue, 12 Oct 2021 16:56:07 +0200 Subject: Simulations Versuch von AM und SSB --- tests/Simulation/AM&DSB_SC/Simulation.py | 439 +++++++++++++++++++++++++++++++ 1 file changed, 439 insertions(+) create mode 100755 tests/Simulation/AM&DSB_SC/Simulation.py (limited to 'tests/Simulation/AM&DSB_SC/Simulation.py') diff --git a/tests/Simulation/AM&DSB_SC/Simulation.py b/tests/Simulation/AM&DSB_SC/Simulation.py new file mode 100755 index 0000000..86fa17d --- /dev/null +++ b/tests/Simulation/AM&DSB_SC/Simulation.py @@ -0,0 +1,439 @@ +#!/usr/bin/env python3 +# -*- coding: utf-8 -*- + +# +# SPDX-License-Identifier: GPL-3.0 +# +# GNU Radio Python Flow Graph +# Title: SSB Simulation +# Author: Sara Halter +# GNU Radio version: 3.8.2.0 + +from distutils.version import StrictVersion + +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 PyQt5 import Qt +from PyQt5.QtCore import QObject, pyqtSlot +from gnuradio import qtgui +from gnuradio.filter import firdes +import sip +from gnuradio import analog +from gnuradio import blocks +from gnuradio import filter +from gnuradio import fosphor +from gnuradio.fft import window +from gnuradio import gr +import sys +import signal +from argparse import ArgumentParser +from gnuradio.eng_arg import eng_float, intx +from gnuradio import eng_notation +from gnuradio.qtgui import Range, RangeWidget + +from gnuradio import qtgui + +class Simulation(gr.top_block, Qt.QWidget): + + def __init__(self): + gr.top_block.__init__(self, "SSB Simulation ") + Qt.QWidget.__init__(self) + self.setWindowTitle("SSB Simulation ") + qtgui.util.check_set_qss() + 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", "Simulation") + + try: + if StrictVersion(Qt.qVersion()) < StrictVersion("5.0.0"): + self.restoreGeometry(self.settings.value("geometry").toByteArray()) + else: + self.restoreGeometry(self.settings.value("geometry")) + except: + pass + + ################################################## + # Variables + ################################################## + self.samp_rate = samp_rate = 256000 + self.audio_rate = audio_rate = 32000 + self.tuning = tuning = 51500 + self.samp_rate_1 = samp_rate_1 = 192000 + self.reverse = reverse = 1 + self.decim = decim = samp_rate/audio_rate + self.carrier_freq = carrier_freq = 16000 + + ################################################## + # Blocks + ################################################## + self._tuning_range = Range(48000, 58000, 100, 51500, 200) + self._tuning_win = RangeWidget(self._tuning_range, self.set_tuning, 'tuning', "counter_slider", float) + self.top_grid_layout.addWidget(self._tuning_win) + # Create the options list + self._reverse_options = (-1, 1, ) + # Create the labels list + self._reverse_labels = ('Upper', 'Lower', ) + # Create the combo box + self._reverse_tool_bar = Qt.QToolBar(self) + self._reverse_tool_bar.addWidget(Qt.QLabel('Sideband' + ": ")) + self._reverse_combo_box = Qt.QComboBox() + self._reverse_tool_bar.addWidget(self._reverse_combo_box) + for _label in self._reverse_labels: self._reverse_combo_box.addItem(_label) + self._reverse_callback = lambda i: Qt.QMetaObject.invokeMethod(self._reverse_combo_box, "setCurrentIndex", Qt.Q_ARG("int", self._reverse_options.index(i))) + self._reverse_callback(self.reverse) + self._reverse_combo_box.currentIndexChanged.connect( + lambda i: self.set_reverse(self._reverse_options[i])) + # Create the radio buttons + self.top_grid_layout.addWidget(self._reverse_tool_bar) + self.qtgui_time_sink_x_0 = qtgui.time_sink_f( + 1024, #size + audio_rate, #samp_rate + "Modulated", #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(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) + self.qtgui_time_sink_x_0.enable_stem_plot(False) + + + labels = ['Signal ', 'Träger', 'Signal 3', 'Signal 4', 'Signal 5', + 'Signal 6', 'Signal 7', 'Signal 8', 'Signal 9', 'Signal 10'] + 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] + styles = [1, 1, 1, 1, 1, + 1, 1, 1, 1, 1] + markers = [-1, -1, -1, -1, -1, + -1, -1, -1, -1, -1] + + + for i in range(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) + self.qtgui_freq_sink_x_0_0_0 = qtgui.freq_sink_c( + 2048, #size + firdes.WIN_BLACKMAN_hARRIS, #wintype + 0, #fc + samp_rate_1, #bw + "", #name + 1 + ) + self.qtgui_freq_sink_x_0_0_0.set_update_time(0.10) + self.qtgui_freq_sink_x_0_0_0.set_y_axis(-140, 10) + self.qtgui_freq_sink_x_0_0_0.set_y_label('Relative Gain', 'dB') + self.qtgui_freq_sink_x_0_0_0.set_trigger_mode(qtgui.TRIG_MODE_FREE, 0.0, 0, "") + self.qtgui_freq_sink_x_0_0_0.enable_autoscale(False) + self.qtgui_freq_sink_x_0_0_0.enable_grid(False) + self.qtgui_freq_sink_x_0_0_0.set_fft_average(1.0) + self.qtgui_freq_sink_x_0_0_0.enable_axis_labels(True) + self.qtgui_freq_sink_x_0_0_0.enable_control_panel(False) + + + + 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 range(1): + if len(labels[i]) == 0: + self.qtgui_freq_sink_x_0_0_0.set_line_label(i, "Data {0}".format(i)) + else: + self.qtgui_freq_sink_x_0_0_0.set_line_label(i, labels[i]) + self.qtgui_freq_sink_x_0_0_0.set_line_width(i, widths[i]) + self.qtgui_freq_sink_x_0_0_0.set_line_color(i, colors[i]) + self.qtgui_freq_sink_x_0_0_0.set_line_alpha(i, alphas[i]) + + self._qtgui_freq_sink_x_0_0_0_win = sip.wrapinstance(self.qtgui_freq_sink_x_0_0_0.pyqwidget(), Qt.QWidget) + self.top_grid_layout.addWidget(self._qtgui_freq_sink_x_0_0_0_win) + self.qtgui_freq_sink_x_0_0 = qtgui.freq_sink_c( + 2048, #size + firdes.WIN_BLACKMAN_hARRIS, #wintype + 0, #fc + samp_rate_1, #bw + "", #name + 1 + ) + self.qtgui_freq_sink_x_0_0.set_update_time(0.10) + self.qtgui_freq_sink_x_0_0.set_y_axis(-140, 10) + self.qtgui_freq_sink_x_0_0.set_y_label('Relative Gain', 'dB') + self.qtgui_freq_sink_x_0_0.set_trigger_mode(qtgui.TRIG_MODE_FREE, 0.0, 0, "") + self.qtgui_freq_sink_x_0_0.enable_autoscale(False) + self.qtgui_freq_sink_x_0_0.enable_grid(False) + self.qtgui_freq_sink_x_0_0.set_fft_average(1.0) + self.qtgui_freq_sink_x_0_0.enable_axis_labels(True) + self.qtgui_freq_sink_x_0_0.enable_control_panel(False) + + + + 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 range(1): + if len(labels[i]) == 0: + self.qtgui_freq_sink_x_0_0.set_line_label(i, "Data {0}".format(i)) + else: + self.qtgui_freq_sink_x_0_0.set_line_label(i, labels[i]) + self.qtgui_freq_sink_x_0_0.set_line_width(i, widths[i]) + self.qtgui_freq_sink_x_0_0.set_line_color(i, colors[i]) + self.qtgui_freq_sink_x_0_0.set_line_alpha(i, alphas[i]) + + self._qtgui_freq_sink_x_0_0_win = sip.wrapinstance(self.qtgui_freq_sink_x_0_0.pyqwidget(), Qt.QWidget) + self.top_grid_layout.addWidget(self._qtgui_freq_sink_x_0_0_win) + self.qtgui_freq_sink_x_0 = qtgui.freq_sink_f( + 1024, #size + firdes.WIN_BLACKMAN_hARRIS, #wintype + 0, #fc + samp_rate, #bw + "", #name + 1 + ) + 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('Relative Gain', '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) + + + 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 range(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) + self.freq_xlating_fir_filter_xxx_0 = filter.freq_xlating_fir_filter_ccc(8, firdes.low_pass(1.0,samp_rate,3000,100), 0, samp_rate) + self.fosphor_glfw_sink_c_0 = fosphor.glfw_sink_c() + self.fosphor_glfw_sink_c_0.set_fft_window(firdes.WIN_BLACKMAN_hARRIS) + self.fosphor_glfw_sink_c_0.set_frequency_range(0, samp_rate) + self.blocks_throttle_0 = blocks.throttle(gr.sizeof_float*1, 48000,True) + self.blocks_repeat_0 = blocks.repeat(gr.sizeof_float*1, 4) + self.blocks_multiply_xx_0_0_0 = blocks.multiply_vff(1) + self.blocks_multiply_xx_0_0 = blocks.multiply_vff(1) + self.blocks_multiply_xx_0 = blocks.multiply_vcc(1) + self.blocks_multiply_const_vxx_0_0_1 = blocks.multiply_const_ff(10) + self.blocks_multiply_const_vxx_0_0_0 = blocks.multiply_const_ff(1) + self.blocks_multiply_const_vxx_0_0 = blocks.multiply_const_ff(100e-3) + self.blocks_multiply_const_vxx_0 = blocks.multiply_const_cc(100e-6) + self.blocks_float_to_complex_0 = blocks.float_to_complex(1) + self.blocks_complex_to_float_0 = blocks.complex_to_float(1) + self.blocks_add_xx_0 = blocks.add_vff(1) + self.band_pass_filter_0 = filter.fir_filter_ccf( + 1, + firdes.band_pass( + 1, + samp_rate_1, + 16300, + 19000, + 200, + firdes.WIN_HAMMING, + 6.76)) + self.analog_sig_source_x_1 = analog.sig_source_c(samp_rate_1, analog.GR_SIN_WAVE, 16000, 1, 0, 0) + self.analog_sig_source_x_0_0_0 = analog.sig_source_f(audio_rate, analog.GR_SIN_WAVE, 1500, 1, 0, 0) + self.analog_sig_source_x_0_0 = analog.sig_source_f(audio_rate, analog.GR_COS_WAVE, 1500, 1, 0, 0) + self.analog_sig_source_x_0 = analog.sig_source_c(samp_rate, analog.GR_TRI_WAVE, 5000, 1, 0, 0) + self.analog_const_source_x_0 = analog.sig_source_f(0, analog.GR_CONST_WAVE, 0, 0, 0) + + + + ################################################## + # Connections + ################################################## + self.connect((self.analog_const_source_x_0, 0), (self.blocks_float_to_complex_0, 0)) + self.connect((self.analog_sig_source_x_0, 0), (self.blocks_multiply_const_vxx_0, 0)) + self.connect((self.analog_sig_source_x_0_0, 0), (self.blocks_multiply_xx_0_0, 0)) + self.connect((self.analog_sig_source_x_0_0_0, 0), (self.blocks_multiply_xx_0_0_0, 1)) + self.connect((self.analog_sig_source_x_1, 0), (self.blocks_multiply_xx_0, 0)) + self.connect((self.band_pass_filter_0, 0), (self.fosphor_glfw_sink_c_0, 0)) + self.connect((self.band_pass_filter_0, 0), (self.qtgui_freq_sink_x_0_0_0, 0)) + self.connect((self.blocks_add_xx_0, 0), (self.blocks_multiply_const_vxx_0_0, 0)) + self.connect((self.blocks_complex_to_float_0, 0), (self.blocks_multiply_xx_0_0, 1)) + self.connect((self.blocks_complex_to_float_0, 1), (self.blocks_multiply_xx_0_0_0, 0)) + self.connect((self.blocks_float_to_complex_0, 0), (self.blocks_multiply_xx_0, 1)) + self.connect((self.blocks_multiply_const_vxx_0, 0), (self.freq_xlating_fir_filter_xxx_0, 0)) + self.connect((self.blocks_multiply_const_vxx_0_0, 0), (self.blocks_throttle_0, 0)) + self.connect((self.blocks_multiply_const_vxx_0_0, 0), (self.qtgui_freq_sink_x_0, 0)) + self.connect((self.blocks_multiply_const_vxx_0_0, 0), (self.qtgui_time_sink_x_0, 0)) + self.connect((self.blocks_multiply_const_vxx_0_0_0, 0), (self.blocks_add_xx_0, 1)) + self.connect((self.blocks_multiply_const_vxx_0_0_1, 0), (self.blocks_repeat_0, 0)) + self.connect((self.blocks_multiply_xx_0, 0), (self.band_pass_filter_0, 0)) + self.connect((self.blocks_multiply_xx_0, 0), (self.qtgui_freq_sink_x_0_0, 0)) + self.connect((self.blocks_multiply_xx_0_0, 0), (self.blocks_add_xx_0, 0)) + self.connect((self.blocks_multiply_xx_0_0_0, 0), (self.blocks_multiply_const_vxx_0_0_0, 0)) + self.connect((self.blocks_repeat_0, 0), (self.blocks_float_to_complex_0, 1)) + self.connect((self.blocks_throttle_0, 0), (self.blocks_multiply_const_vxx_0_0_1, 0)) + self.connect((self.freq_xlating_fir_filter_xxx_0, 0), (self.blocks_complex_to_float_0, 0)) + + + def closeEvent(self, event): + self.settings = Qt.QSettings("GNU Radio", "Simulation") + self.settings.setValue("geometry", self.saveGeometry()) + event.accept() + + def get_samp_rate(self): + return self.samp_rate + + def set_samp_rate(self, samp_rate): + self.samp_rate = samp_rate + self.set_decim(self.samp_rate/self.audio_rate) + self.analog_sig_source_x_0.set_sampling_freq(self.samp_rate) + self.freq_xlating_fir_filter_xxx_0.set_taps(firdes.low_pass(1.0,self.samp_rate,3000,100)) + self.qtgui_freq_sink_x_0.set_frequency_range(0, self.samp_rate) + self.fosphor_glfw_sink_c_0.set_frequency_range(0, self.samp_rate) + + def get_audio_rate(self): + return self.audio_rate + + def set_audio_rate(self, audio_rate): + self.audio_rate = audio_rate + self.set_decim(self.samp_rate/self.audio_rate) + self.analog_sig_source_x_0_0.set_sampling_freq(self.audio_rate) + self.analog_sig_source_x_0_0_0.set_sampling_freq(self.audio_rate) + self.qtgui_time_sink_x_0.set_samp_rate(self.audio_rate) + + def get_tuning(self): + return self.tuning + + def set_tuning(self, tuning): + self.tuning = tuning + + def get_samp_rate_1(self): + return self.samp_rate_1 + + def set_samp_rate_1(self, samp_rate_1): + self.samp_rate_1 = samp_rate_1 + self.analog_sig_source_x_1.set_sampling_freq(self.samp_rate_1) + self.band_pass_filter_0.set_taps(firdes.band_pass(1, self.samp_rate_1, 16300, 19000, 200, firdes.WIN_HAMMING, 6.76)) + self.qtgui_freq_sink_x_0_0.set_frequency_range(0, self.samp_rate_1) + self.qtgui_freq_sink_x_0_0_0.set_frequency_range(0, self.samp_rate_1) + + def get_reverse(self): + return self.reverse + + def set_reverse(self, reverse): + self.reverse = reverse + self._reverse_callback(self.reverse) + + def get_decim(self): + return self.decim + + def set_decim(self, decim): + self.decim = decim + + def get_carrier_freq(self): + return self.carrier_freq + + def set_carrier_freq(self, carrier_freq): + self.carrier_freq = carrier_freq + + + + + +def main(top_block_cls=Simulation, options=None): + if gr.enable_realtime_scheduling() != gr.RT_OK: + print("Error: failed to enable real-time scheduling.") + + if StrictVersion("4.5.0") <= StrictVersion(Qt.qVersion()) < StrictVersion("5.0.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 sig_handler(sig=None, frame=None): + Qt.QApplication.quit() + + signal.signal(signal.SIGINT, sig_handler) + signal.signal(signal.SIGTERM, sig_handler) + + timer = Qt.QTimer() + timer.start(500) + timer.timeout.connect(lambda: None) + + def quitting(): + tb.stop() + tb.wait() + + qapp.aboutToQuit.connect(quitting) + qapp.exec_() + +if __name__ == '__main__': + main() -- cgit v1.2.1