From 17d49af3f1fc91441455e6825f05b604d19b8e82 Mon Sep 17 00:00:00 2001 From: Nao Pross Date: Tue, 14 Dec 2021 19:44:41 +0100 Subject: Attempt CAZAC correlation codes in 16-QAM (not working) --- flowgraphs/qam_sim.grc | 231 +++++++++++++++++++++++++++++++++++++------------ flowgraphs/qam_sim.py | 199 +++++++++++++++++++++++++++++++++++++----- 2 files changed, 357 insertions(+), 73 deletions(-) diff --git a/flowgraphs/qam_sim.grc b/flowgraphs/qam_sim.grc index 02b174c..f1f9cc9 100644 --- a/flowgraphs/qam_sim.grc +++ b/flowgraphs/qam_sim.grc @@ -8,7 +8,7 @@ options: description: '' gen_cmake: 'On' gen_linking: dynamic - generate_options: no_gui + generate_options: qt_gui hier_block_src_path: '.:' id: qam_sim max_nouts: '0' @@ -32,19 +32,78 @@ options: state: enabled blocks: -- name: access_code_symbols +- name: access_code id: variable parameters: comment: '' - value: '[(-0.31622776601683794-0.31622776601683794j), (-0.31622776601683794+0.31622776601683794j), - (-0.9486832980505138-0.9486832980505138j), (-0.9486832980505138-0.31622776601683794j)]' + value: np.exp(1j * (np.pi * 13 * np.arange(20) ** 20 / 20)) states: bus_sink: false bus_source: false bus_structure: null - coordinate: [792, 1092.0] + coordinate: [32, 348.0] rotation: 0 state: true +- name: access_code_sps + id: variable + parameters: + comment: '' + value: np.asarray(np.matmul(np.matrix(access_code).T, np.matrix(np.ones(sps)))).ravel() + states: + bus_sink: false + bus_source: false + bus_structure: null + coordinate: [424, 196.0] + rotation: 0 + state: disabled +- name: access_code_sps + id: variable + parameters: + comment: '' + value: '[(1+0j), (1+0j), (1+0j), (1+0j), (-0.4539904997395467+0.8910065241883679j), + (-0.4539904997395467+0.8910065241883679j), (-0.4539904997395467+0.8910065241883679j), + (-0.4539904997395467+0.8910065241883679j), (0.30901699445433095+0.9510565162693603j), + (0.30901699445433095+0.9510565162693603j), (0.30901699445433095+0.9510565162693603j), + (0.30901699445433095+0.9510565162693603j), (-0.45399010831290065+0.8910067236300974j), + (-0.45399010831290065+0.8910067236300974j), (-0.45399010831290065+0.8910067236300974j), + (-0.45399010831290065+0.8910067236300974j), (0.30910023286329796+0.9510294664435246j), + (0.30910023286329796+0.9510294664435246j), (0.30910023286329796+0.9510294664435246j), + (0.30910023286329796+0.9510294664435246j), (0.6958765427916561+0.7181614283657487j), + (0.6958765427916561+0.7181614283657487j), (0.6958765427916561+0.7181614283657487j), + (0.6958765427916561+0.7181614283657487j), (0.6995785030650676+0.7145557487343022j), + (0.6995785030650676+0.7145557487343022j), (0.6995785030650676+0.7145557487343022j), + (0.6995785030650676+0.7145557487343022j), (-0.8755764602527171+0.4830795609952073j), + (-0.8755764602527171+0.4830795609952073j), (-0.8755764602527171+0.4830795609952073j), + (-0.8755764602527171+0.4830795609952073j), (-0.8320971948858552-0.5546298389584634j), + (-0.8320971948858552-0.5546298389584634j), (-0.8320971948858552-0.5546298389584634j), + (-0.8320971948858552-0.5546298389584634j), (0.13124680853227685+0.9913497239874998j), + (0.13124680853227685+0.9913497239874998j), (0.13124680853227685+0.9913497239874998j), + (0.13124680853227685+0.9913497239874998j), (0.9219105024602736-0.38740292390926256j), + (0.9219105024602736-0.38740292390926256j), (0.9219105024602736-0.38740292390926256j), + (0.9219105024602736-0.38740292390926256j), (0.9881969645782696-0.15318863925988194j), + (0.9881969645782696-0.15318863925988194j), (0.9881969645782696-0.15318863925988194j), + (0.9881969645782696-0.15318863925988194j), (-0.5511352729247617-0.8344159100461522j), + (-0.5511352729247617-0.8344159100461522j), (-0.5511352729247617-0.8344159100461522j), + (-0.5511352729247617-0.8344159100461522j), (-0.24665339430031927+0.9691037628036181j), + (-0.24665339430031927+0.9691037628036181j), (-0.24665339430031927+0.9691037628036181j), + (-0.24665339430031927+0.9691037628036181j), (-0.44951814008601154-0.8932712027898431j), + (-0.44951814008601154-0.8932712027898431j), (-0.44951814008601154-0.8932712027898431j), + (-0.44951814008601154-0.8932712027898431j), (-0.23300791854304934-0.9724748376673998j), + (-0.23300791854304934-0.9724748376673998j), (-0.23300791854304934-0.9724748376673998j), + (-0.23300791854304934-0.9724748376673998j), (1+0j), (1+0j), (1+0j), (1+0j), + (0.8559628469680924+0.5170373338650489j), (0.8559628469680924+0.5170373338650489j), + (0.8559628469680924+0.5170373338650489j), (0.8559628469680924+0.5170373338650489j), + (-0.8354473464520049-0.5495704971214372j), (-0.8354473464520049-0.5495704971214372j), + (-0.8354473464520049-0.5495704971214372j), (-0.8354473464520049-0.5495704971214372j), + (-0.999246638643493+0.03880921488101298j), (-0.999246638643493+0.03880921488101298j), + (-0.999246638643493+0.03880921488101298j), (-0.999246638643493+0.03880921488101298j)]' + states: + bus_sink: false + bus_source: false + bus_structure: null + coordinate: [256, 196.0] + rotation: 0 + state: disabled - name: carrier_freq id: variable parameters: @@ -54,7 +113,7 @@ blocks: bus_sink: false bus_source: false bus_structure: null - coordinate: [472, 572.0] + coordinate: [472, 756.0] rotation: 0 state: true - name: excess_bw @@ -66,7 +125,7 @@ blocks: bus_sink: false bus_source: false bus_structure: null - coordinate: [336, 1156.0] + coordinate: [336, 1332.0] rotation: 0 state: true - name: frame @@ -80,7 +139,7 @@ blocks: bus_sink: false bus_source: false bus_structure: null - coordinate: [48, 612.0] + coordinate: [48, 796.0] rotation: 0 state: disabled - name: frame_len @@ -92,7 +151,7 @@ blocks: bus_sink: false bus_source: false bus_structure: null - coordinate: [184, 252.0] + coordinate: [184, 436.0] rotation: 0 state: true - name: nfilts @@ -104,7 +163,7 @@ blocks: bus_sink: false bus_source: false bus_structure: null - coordinate: [256, 1156.0] + coordinate: [256, 1332.0] rotation: 0 state: true - name: qam_const @@ -122,7 +181,7 @@ blocks: bus_sink: false bus_source: false bus_structure: null - coordinate: [472, 356.0] + coordinate: [472, 532.0] rotation: 0 state: true - name: rrc_taps @@ -134,7 +193,7 @@ blocks: bus_sink: false bus_source: false bus_structure: null - coordinate: [432, 1156.0] + coordinate: [432, 1332.0] rotation: 0 state: true - name: samp_rate @@ -170,7 +229,7 @@ blocks: bus_sink: false bus_source: false bus_structure: null - coordinate: [48, 252.0] + coordinate: [48, 436.0] rotation: 0 state: true - name: analog_random_source_x_0 @@ -190,7 +249,7 @@ blocks: bus_sink: false bus_source: false bus_structure: null - coordinate: [48, 420.0] + coordinate: [40, 604.0] rotation: 0 state: true - name: blocks_complex_to_mag_0 @@ -206,9 +265,9 @@ blocks: bus_sink: false bus_source: false bus_structure: null - coordinate: [1088, 1080.0] + coordinate: [1088, 1272.0] rotation: 0 - state: disabled + state: enabled - name: blocks_null_source_0 id: blocks_null_source parameters: @@ -225,7 +284,7 @@ blocks: bus_sink: false bus_source: false bus_structure: null - coordinate: [48, 528.0] + coordinate: [48, 712.0] rotation: 0 state: true - name: blocks_repack_bits_bb_0 @@ -245,7 +304,7 @@ blocks: bus_sink: false bus_source: false bus_structure: null - coordinate: [768, 1300.0] + coordinate: [768, 1476.0] rotation: 0 state: true - name: blocks_stream_mux_0 @@ -264,9 +323,28 @@ blocks: bus_sink: false bus_source: false bus_structure: null - coordinate: [272, 416.0] + coordinate: [272, 600.0] rotation: 0 state: true +- name: blocks_stream_mux_1 + id: blocks_stream_mux + parameters: + affinity: '' + alias: '' + comment: '' + lengths: '[len(access_code)*4, len(testvec)//4*sps]' + maxoutbuf: '0' + minoutbuf: '0' + num_inputs: '2' + type: complex + vlen: '1' + states: + bus_sink: false + bus_source: false + bus_structure: null + coordinate: [728, 528.0] + rotation: 0 + state: disabled - name: blocks_stream_to_vector_0 id: blocks_stream_to_vector parameters: @@ -282,7 +360,7 @@ blocks: bus_sink: false bus_source: false bus_structure: null - coordinate: [992, 1312.0] + coordinate: [1000, 1488.0] rotation: 0 state: true - name: blocks_tagged_stream_align_0 @@ -300,7 +378,7 @@ blocks: bus_sink: false bus_source: false bus_structure: null - coordinate: [496, 1308.0] + coordinate: [496, 1484.0] rotation: 0 state: true - name: blocks_throttle_0 @@ -319,7 +397,7 @@ blocks: bus_sink: false bus_source: false bus_structure: null - coordinate: [744, 444.0] + coordinate: [864, 628.0] rotation: 0 state: true - name: blocks_vector_source_x_0 @@ -339,9 +417,29 @@ blocks: bus_sink: false bus_source: false bus_structure: null - coordinate: [48, 324.0] + coordinate: [48, 508.0] rotation: 0 state: true +- name: blocks_vector_source_x_1 + id: blocks_vector_source_x + parameters: + affinity: '' + alias: '' + comment: '' + maxoutbuf: '0' + minoutbuf: '0' + repeat: 'True' + tags: '[]' + type: complex + vector: access_code_sps + vlen: '1' + states: + bus_sink: false + bus_source: false + bus_structure: null + coordinate: [256, 300.0] + rotation: 0 + state: disabled - name: channels_channel_model_0 id: channels_channel_model parameters: @@ -360,7 +458,7 @@ blocks: bus_sink: false bus_source: false bus_structure: null - coordinate: [312, 716.0] + coordinate: [312, 900.0] rotation: 0 state: enabled - name: channels_selective_fading_model_0 @@ -373,7 +471,7 @@ blocks: alias: '' comment: '' delays: (0,1.8) - fDTs: ((4*carrier_freq)/(3*10e8))/samp_rate + fDTs: ((2*carrier_freq)/(3*10e8))/samp_rate mags: (1,0.12) maxoutbuf: '0' minoutbuf: '0' @@ -383,9 +481,9 @@ blocks: bus_sink: false bus_source: false bus_structure: null - coordinate: [600, 572.0] + coordinate: [600, 748.0] rotation: 0 - state: disabled + state: enabled - name: digital_constellation_decoder_cb_0 id: digital_constellation_decoder_cb parameters: @@ -399,7 +497,7 @@ blocks: bus_sink: false bus_source: false bus_structure: null - coordinate: [256, 1308.0] + coordinate: [256, 1484.0] rotation: 0 state: true - name: digital_constellation_modulator_0 @@ -420,7 +518,7 @@ blocks: bus_sink: false bus_source: false bus_structure: null - coordinate: [472, 420.0] + coordinate: [472, 604.0] rotation: 0 state: true - name: digital_corr_est_cc_0 @@ -429,18 +527,18 @@ blocks: affinity: '' alias: '' comment: '' - mark_delay: len(access_code_symbols) // 2 + mark_delay: len(access_code) // 2 maxoutbuf: '0' minoutbuf: '0' sps: '1' - symbols: access_code_symbols + symbols: access_code threshold: '0.9' threshold_method: digital.THRESHOLD_ABSOLUTE states: bus_sink: false bus_source: false bus_structure: null - coordinate: [792, 972.0] + coordinate: [792, 1148.0] rotation: 0 state: true - name: digital_lms_dd_equalizer_cc_0 @@ -453,13 +551,13 @@ blocks: maxoutbuf: '0' minoutbuf: '0' mu: 5e-3 - num_taps: '2' + num_taps: '5' sps: '1' states: bus_sink: false bus_source: false bus_structure: null - coordinate: [520, 980.0] + coordinate: [520, 1156.0] rotation: 0 state: true - name: digital_pfb_clock_sync_xxx_0 @@ -482,7 +580,7 @@ blocks: bus_sink: false bus_source: false bus_structure: null - coordinate: [248, 1004.0] + coordinate: [248, 1180.0] rotation: 0 state: true - name: fadingui_ber_0 @@ -498,7 +596,7 @@ blocks: bus_sink: false bus_source: false bus_structure: null - coordinate: [1184, 1300.0] + coordinate: [1200, 1468.0] rotation: 0 state: true - name: fadingui_datasource_0 @@ -515,7 +613,7 @@ blocks: bus_sink: false bus_source: false bus_structure: null - coordinate: [208, 612.0] + coordinate: [208, 788.0] rotation: 0 state: disabled - name: fadingui_multipath_fading_0 @@ -533,9 +631,9 @@ blocks: bus_sink: false bus_source: false bus_structure: null - coordinate: [648, 740.0] + coordinate: [648, 924.0] rotation: 0 - state: enabled + state: disabled - name: fadingui_netsink_0 id: fadingui_netsink parameters: @@ -549,7 +647,7 @@ blocks: bus_sink: false bus_source: false bus_structure: null - coordinate: [984, 620.0] + coordinate: [984, 796.0] rotation: 0 state: enabled - name: fadingui_netsink_1 @@ -565,7 +663,7 @@ blocks: bus_sink: false bus_source: false bus_structure: null - coordinate: [800, 916.0] + coordinate: [800, 1100.0] rotation: 0 state: enabled - name: fadingui_netsink_3 @@ -581,7 +679,7 @@ blocks: bus_sink: false bus_source: false bus_structure: null - coordinate: [1280, 916.0] + coordinate: [1280, 1100.0] rotation: 0 state: enabled - name: fadingui_netsink_4 @@ -597,7 +695,7 @@ blocks: bus_sink: false bus_source: false bus_structure: null - coordinate: [528, 916.0] + coordinate: [528, 1100.0] rotation: 0 state: enabled - name: fadingui_phasecorrection_0 @@ -613,7 +711,7 @@ blocks: bus_sink: false bus_source: false bus_structure: null - coordinate: [1032, 988.0] + coordinate: [1032, 1164.0] rotation: 0 state: enabled - name: import_0 @@ -718,9 +816,9 @@ blocks: bus_sink: false bus_source: false bus_structure: null - coordinate: [1304, 800.0] + coordinate: [1304, 984.0] rotation: 0 - state: disabled + state: enabled - name: qtgui_time_sink_x_0 id: qtgui_time_sink_x parameters: @@ -815,7 +913,30 @@ blocks: bus_sink: false bus_source: false bus_structure: null - coordinate: [1272, 1060.0] + coordinate: [1320, 1252.0] + rotation: 0 + state: enabled +- name: root_raised_cosine_filter_0 + id: root_raised_cosine_filter + parameters: + affinity: '' + alias: '' + alpha: excess_bw + comment: '' + decim: '1' + gain: '10' + interp: '1' + maxoutbuf: '0' + minoutbuf: '0' + ntaps: int(samp_rate) + samp_rate: samp_rate + sym_rate: 1 / sps + type: fir_filter_ccf + states: + bus_sink: false + bus_source: false + bus_structure: null + coordinate: [472, 380.0] rotation: 0 state: disabled - name: virtual_sink_0 @@ -828,7 +949,7 @@ blocks: bus_sink: false bus_source: false bus_structure: null - coordinate: [952, 444.0] + coordinate: [1048, 628.0] rotation: 0 state: true - name: virtual_sink_1 @@ -841,7 +962,7 @@ blocks: bus_sink: false bus_source: false bus_structure: null - coordinate: [1000, 756.0] + coordinate: [1000, 940.0] rotation: 0 state: true - name: virtual_sink_2 @@ -854,7 +975,7 @@ blocks: bus_sink: false bus_source: false bus_structure: null - coordinate: [1336, 988.0] + coordinate: [1336, 1164.0] rotation: 0 state: true - name: virtual_source_0 @@ -867,7 +988,7 @@ blocks: bus_sink: false bus_source: false bus_structure: null - coordinate: [48, 1308.0] + coordinate: [48, 1484.0] rotation: 0 state: true - name: virtual_source_1 @@ -880,7 +1001,7 @@ blocks: bus_sink: false bus_source: false bus_structure: null - coordinate: [48, 1052.0] + coordinate: [48, 1228.0] rotation: 0 state: true - name: virtual_source_2 @@ -893,7 +1014,7 @@ blocks: bus_sink: false bus_source: false bus_structure: null - coordinate: [56, 756.0] + coordinate: [56, 940.0] rotation: 0 state: true @@ -903,15 +1024,18 @@ connections: - [blocks_null_source_0, '0', blocks_stream_mux_0, '2'] - [blocks_repack_bits_bb_0, '0', blocks_stream_to_vector_0, '0'] - [blocks_stream_mux_0, '0', digital_constellation_modulator_0, '0'] +- [blocks_stream_mux_1, '0', blocks_throttle_0, '0'] - [blocks_stream_to_vector_0, '0', fadingui_ber_0, '0'] - [blocks_tagged_stream_align_0, '0', blocks_repack_bits_bb_0, '0'] - [blocks_throttle_0, '0', virtual_sink_0, '0'] - [blocks_vector_source_x_0, '0', blocks_stream_mux_0, '0'] +- [blocks_vector_source_x_1, '0', root_raised_cosine_filter_0, '0'] - [channels_channel_model_0, '0', channels_selective_fading_model_0, '0'] - [channels_channel_model_0, '0', fadingui_multipath_fading_0, '0'] - [channels_selective_fading_model_0, '0', fadingui_netsink_0, '0'] - [channels_selective_fading_model_0, '0', virtual_sink_1, '0'] - [digital_constellation_decoder_cb_0, '0', blocks_tagged_stream_align_0, '0'] +- [digital_constellation_modulator_0, '0', blocks_stream_mux_1, '1'] - [digital_constellation_modulator_0, '0', blocks_throttle_0, '0'] - [digital_corr_est_cc_0, '0', fadingui_phasecorrection_0, '0'] - [digital_corr_est_cc_0, '1', blocks_complex_to_mag_0, '0'] @@ -927,6 +1051,7 @@ connections: - [fadingui_phasecorrection_0, '0', fadingui_netsink_3, '0'] - [fadingui_phasecorrection_0, '0', qtgui_const_sink_x_0, '2'] - [fadingui_phasecorrection_0, '0', virtual_sink_2, '0'] +- [root_raised_cosine_filter_0, '0', blocks_stream_mux_1, '0'] - [virtual_source_0, '0', digital_constellation_decoder_cb_0, '0'] - [virtual_source_1, '0', digital_pfb_clock_sync_xxx_0, '0'] - [virtual_source_2, '0', channels_channel_model_0, '0'] diff --git a/flowgraphs/qam_sim.py b/flowgraphs/qam_sim.py index 66045c2..bd8589b 100755 --- a/flowgraphs/qam_sim.py +++ b/flowgraphs/qam_sim.py @@ -9,10 +9,25 @@ # Author: Naoki Sean Pross, Sara Cinzia 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 gnuradio import qtgui +from gnuradio.filter import firdes +import sip from gnuradio import blocks import numpy from gnuradio import channels -from gnuradio.filter import firdes from gnuradio import digital from gnuradio import gr import sys @@ -23,11 +38,40 @@ from gnuradio import eng_notation import fadingui import numpy as np +from gnuradio import qtgui -class qam_sim(gr.top_block): +class qam_sim(gr.top_block, Qt.QWidget): def __init__(self): gr.top_block.__init__(self, "QAM Sim") + Qt.QWidget.__init__(self) + self.setWindowTitle("QAM Sim") + 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", "qam_sim") + + 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 @@ -41,21 +85,107 @@ class qam_sim(gr.top_block): self.qam_const = qam_const = digital.constellation_16qam().base() self.frame_len = frame_len = len(testvec) +4 self.carrier_freq = carrier_freq = 2.4e9 - self.access_code_symbols = access_code_symbols = [(-0.31622776601683794-0.31622776601683794j), (-0.31622776601683794+0.31622776601683794j), (-0.9486832980505138-0.9486832980505138j), (-0.9486832980505138-0.31622776601683794j)] + self.access_code = access_code = np.exp(1j * (np.pi * 13 * np.arange(20) ** 20 / 20)) ################################################## # Blocks ################################################## + self.qtgui_time_sink_x_0 = qtgui.time_sink_f( + 1024, #size + samp_rate / sps, #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(0, 10) + + 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 1', 'Signal 2', '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_const_sink_x_0 = qtgui.const_sink_c( + 1024, #size + "", #name + 3 #number of inputs + ) + self.qtgui_const_sink_x_0.set_update_time(0.10) + self.qtgui_const_sink_x_0.set_y_axis(-2, 2) + self.qtgui_const_sink_x_0.set_x_axis(-2, 2) + self.qtgui_const_sink_x_0.set_trigger_mode(qtgui.TRIG_MODE_FREE, qtgui.TRIG_SLOPE_POS, 0.0, 0, "") + self.qtgui_const_sink_x_0.enable_autoscale(False) + self.qtgui_const_sink_x_0.enable_grid(False) + self.qtgui_const_sink_x_0.enable_axis_labels(True) + + + labels = ['', '', '', '', '', + '', '', '', '', ''] + widths = [1, 1, 1, 1, 1, + 1, 1, 1, 1, 1] + colors = ["blue", "red", "red", "red", "red", + "red", "red", "red", "red", "red"] + styles = [0, 0, 0, 0, 0, + 0, 0, 0, 0, 0] + markers = [0, 0, 0, 0, 0, + 0, 0, 0, 0, 0] + 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(3): + if len(labels[i]) == 0: + self.qtgui_const_sink_x_0.set_line_label(i, "Data {0}".format(i)) + else: + self.qtgui_const_sink_x_0.set_line_label(i, labels[i]) + self.qtgui_const_sink_x_0.set_line_width(i, widths[i]) + self.qtgui_const_sink_x_0.set_line_color(i, colors[i]) + self.qtgui_const_sink_x_0.set_line_style(i, styles[i]) + self.qtgui_const_sink_x_0.set_line_marker(i, markers[i]) + self.qtgui_const_sink_x_0.set_line_alpha(i, alphas[i]) + + self._qtgui_const_sink_x_0_win = sip.wrapinstance(self.qtgui_const_sink_x_0.pyqwidget(), Qt.QWidget) + self.top_grid_layout.addWidget(self._qtgui_const_sink_x_0_win) self.fadingui_phasecorrection_0 = fadingui.phasecorrection(frame_len) self.fadingui_netsink_4 = fadingui.netsink(address='udp://localhost:31417', dtype="complex", vlen=1) self.fadingui_netsink_3 = fadingui.netsink(address='udp://localhost:31419', dtype="complex", vlen=1) self.fadingui_netsink_1 = fadingui.netsink(address='udp://localhost:31418', dtype="complex", vlen=1) self.fadingui_netsink_0 = fadingui.netsink(address='udp://localhost:31416', dtype="complex", vlen=1) - self.fadingui_multipath_fading_0 = fadingui.multipath_fading(amplitudes=[0.12], delays=[1.8], los =True) self.fadingui_ber_0 = fadingui.ber(vgl=testvec + list(np.zeros(4)), vlen=frame_len,address='udp://localhost:31415') self.digital_pfb_clock_sync_xxx_0 = digital.pfb_clock_sync_ccf(sps, 2 * np.pi / 100, rrc_taps, 32, 16, 1.5, 1) - self.digital_lms_dd_equalizer_cc_0 = digital.lms_dd_equalizer_cc(2, 5e-3, 1, qam_const) - self.digital_corr_est_cc_0 = digital.corr_est_cc(access_code_symbols, 1, len(access_code_symbols) // 2, 0.9, digital.THRESHOLD_ABSOLUTE) + self.digital_lms_dd_equalizer_cc_0 = digital.lms_dd_equalizer_cc(5, 5e-3, 1, qam_const) + self.digital_corr_est_cc_0 = digital.corr_est_cc(access_code, 1, len(access_code) // 2, 0.9, digital.THRESHOLD_ABSOLUTE) self.digital_constellation_modulator_0 = digital.generic_mod( constellation=qam_const, differential=False, @@ -65,6 +195,7 @@ class qam_sim(gr.top_block): verbose=False, log=False) self.digital_constellation_decoder_cb_0 = digital.constellation_decoder_cb(qam_const) + self.channels_selective_fading_model_0 = channels.selective_fading_model( 8, ((2*carrier_freq)/(3*10e8))/samp_rate, False, 4.0, 21, (0,1.8), (1,0.12), 8 ) self.channels_channel_model_0 = channels.channel_model( noise_voltage=100e-3, frequency_offset=0, @@ -79,6 +210,7 @@ class qam_sim(gr.top_block): self.blocks_stream_mux_0 = blocks.stream_mux(gr.sizeof_char*1, [len(testvec), 4, 0]) self.blocks_repack_bits_bb_0 = blocks.repack_bits_bb(2, 8, "", False, gr.GR_LSB_FIRST) self.blocks_null_source_0 = blocks.null_source(gr.sizeof_char*1) + self.blocks_complex_to_mag_0 = blocks.complex_to_mag(1) self.analog_random_source_x_0 = blocks.vector_source_b(list(map(int, numpy.random.randint(0, 255, 400))), True) @@ -87,6 +219,7 @@ class qam_sim(gr.top_block): # Connections ################################################## self.connect((self.analog_random_source_x_0, 0), (self.blocks_stream_mux_0, 1)) + self.connect((self.blocks_complex_to_mag_0, 0), (self.qtgui_time_sink_x_0, 0)) self.connect((self.blocks_null_source_0, 0), (self.blocks_stream_mux_0, 2)) self.connect((self.blocks_repack_bits_bb_0, 0), (self.blocks_stream_to_vector_0, 0)) self.connect((self.blocks_stream_mux_0, 0), (self.digital_constellation_modulator_0, 0)) @@ -94,20 +227,29 @@ class qam_sim(gr.top_block): self.connect((self.blocks_tagged_stream_align_0, 0), (self.blocks_repack_bits_bb_0, 0)) self.connect((self.blocks_throttle_0, 0), (self.channels_channel_model_0, 0)) self.connect((self.blocks_vector_source_x_0, 0), (self.blocks_stream_mux_0, 0)) - self.connect((self.channels_channel_model_0, 0), (self.fadingui_multipath_fading_0, 0)) + self.connect((self.channels_channel_model_0, 0), (self.channels_selective_fading_model_0, 0)) + self.connect((self.channels_selective_fading_model_0, 0), (self.digital_pfb_clock_sync_xxx_0, 0)) + self.connect((self.channels_selective_fading_model_0, 0), (self.fadingui_netsink_0, 0)) self.connect((self.digital_constellation_decoder_cb_0, 0), (self.blocks_tagged_stream_align_0, 0)) self.connect((self.digital_constellation_modulator_0, 0), (self.blocks_throttle_0, 0)) + self.connect((self.digital_corr_est_cc_0, 1), (self.blocks_complex_to_mag_0, 0)) self.connect((self.digital_corr_est_cc_0, 0), (self.fadingui_phasecorrection_0, 0)) self.connect((self.digital_lms_dd_equalizer_cc_0, 0), (self.digital_corr_est_cc_0, 0)) self.connect((self.digital_lms_dd_equalizer_cc_0, 0), (self.fadingui_netsink_1, 0)) + self.connect((self.digital_lms_dd_equalizer_cc_0, 0), (self.qtgui_const_sink_x_0, 1)) self.connect((self.digital_pfb_clock_sync_xxx_0, 0), (self.digital_lms_dd_equalizer_cc_0, 0)) self.connect((self.digital_pfb_clock_sync_xxx_0, 0), (self.fadingui_netsink_4, 0)) - self.connect((self.fadingui_multipath_fading_0, 0), (self.digital_pfb_clock_sync_xxx_0, 0)) - self.connect((self.fadingui_multipath_fading_0, 0), (self.fadingui_netsink_0, 0)) + self.connect((self.digital_pfb_clock_sync_xxx_0, 0), (self.qtgui_const_sink_x_0, 0)) self.connect((self.fadingui_phasecorrection_0, 0), (self.digital_constellation_decoder_cb_0, 0)) self.connect((self.fadingui_phasecorrection_0, 0), (self.fadingui_netsink_3, 0)) + self.connect((self.fadingui_phasecorrection_0, 0), (self.qtgui_const_sink_x_0, 2)) + def closeEvent(self, event): + self.settings = Qt.QSettings("GNU Radio", "qam_sim") + self.settings.setValue("geometry", self.saveGeometry()) + event.accept() + def get_testvec(self): return self.testvec @@ -122,6 +264,7 @@ class qam_sim(gr.top_block): def set_sps(self, sps): self.sps = sps self.set_rrc_taps(firdes.root_raised_cosine(self.nfilts, self.nfilts, 1.0/float(self.sps), self.excess_bw, 45*self.nfilts)) + self.qtgui_time_sink_x_0.set_samp_rate(self.samp_rate / self.sps) def get_nfilts(self): return self.nfilts @@ -143,6 +286,8 @@ class qam_sim(gr.top_block): def set_samp_rate(self, samp_rate): self.samp_rate = samp_rate self.blocks_throttle_0.set_sample_rate(self.samp_rate) + self.channels_selective_fading_model_0.set_fDTs(((2*self.carrier_freq)/(3*10e8))/self.samp_rate) + self.qtgui_time_sink_x_0.set_samp_rate(self.samp_rate / self.sps) def get_rrc_taps(self): return self.rrc_taps @@ -168,13 +313,14 @@ class qam_sim(gr.top_block): def set_carrier_freq(self, carrier_freq): self.carrier_freq = carrier_freq + self.channels_selective_fading_model_0.set_fDTs(((2*self.carrier_freq)/(3*10e8))/self.samp_rate) - def get_access_code_symbols(self): - return self.access_code_symbols + def get_access_code(self): + return self.access_code - def set_access_code_symbols(self, access_code_symbols): - self.access_code_symbols = access_code_symbols - self.digital_corr_est_cc_0.set_mark_delay(len(self.access_code_symbols) // 2) + def set_access_code(self, access_code): + self.access_code = access_code + self.digital_corr_est_cc_0.set_mark_delay(len(self.access_code) // 2) @@ -183,21 +329,34 @@ class qam_sim(gr.top_block): def main(top_block_cls=qam_sim, 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() - def sig_handler(sig=None, frame=None): - tb.stop() - tb.wait() + tb.start() + + tb.show() - sys.exit(0) + def sig_handler(sig=None, frame=None): + Qt.QApplication.quit() signal.signal(signal.SIGINT, sig_handler) signal.signal(signal.SIGTERM, sig_handler) - tb.start() + timer = Qt.QTimer() + timer.start(500) + timer.timeout.connect(lambda: None) - tb.wait() + def quitting(): + tb.stop() + tb.wait() + qapp.aboutToQuit.connect(quitting) + qapp.exec_() if __name__ == '__main__': main() -- cgit v1.2.1