摘要
使用gnuradio观测理想信道和随机信道的眼图;
平台信息
- GNU Radio Companion:3.10.8.0 (Python 3.10.13)
原理简介
数字基带传输系统和数字带通传输系统
- 数字基带传输系统
数字基带传输系统是指将未经调制的基带信号直接送上信道的系统。基带信号是指信源（信息源，也称发终端）发出的没有经过调制（进行频谱搬移和变换）的原始电信号，其特点是频率较低，信号频谱从零频附近开始，具有低通形式。
在某些具有低通特性的有线信道中，特别是在传输距离不太远的情况下，基带信号可以不经过载波调制而直接进行传输，这样的传输系统，称为数字基带传输系统。 - 数字带通传输系统
与之相对，数字带通传输系统是指包括调制和解调过程的数字传输系统。在这个系统中，基带信号被调制到指定的频带上再送上信道，在接收端则采用相应的解调方式原恢复出基带信号。
数字带通传输系统的主要特点是信号频谱通常是带通型的，适合于在带通型信道中传输。
总的来说，数字基带传输系统和数字带通传输系统的主要区别在于信号是否经过调制。数字基带传输系统传输的是未经调制的基带信号，而数字带通传输系统传输的是经过调制的信号。
| 数字通信系统及其信道的结构 |
|---|
 |
理想信道和随机信道
[https://zhuanlan.zhihu.com/p/664291541]
[https://zhuanlan.zhihu.com/p/626116095]
[https://www.zhihu.com/question/49575216]
[https://www.gaussianwaves.com/2020/08/rician-flat-fading-channel-simulation/]
信道是任何一种通信系统中必不可少的组成部分。任何一个通信系统都可以视为由发送，信道与接收三部分组成。信道通常指以传输媒介为基础的信号通道。
信号在信道中传输，可能遇到的影响主要有信道加性噪声 、 信号幅度衰减和相位失真 、 信道特性的非线性 、带宽限制和多径失真等。
- 理想信道是指信号传输无失真的信道.
- 随机信道:
- 恒参信道:传输媒质性质稳定从而使得信道特性稳定,如高斯信道(AWGN).
- 随参信道:传输媒质性质的随机变化从而使得信道特性也随机变化,如Nakagami-m信道等.
高斯信道通常用于描述无障碍环境中的传统无线通信系统，所以不考虑多径效应的影响，相当于恒参信道。
信号经过多条路径到达接收机；到达接收机的多径信号：如果同相（相位差是2π的整数倍时），叠加后信号会增强；如果反相（相位差是π的奇数倍时），叠加后信号会减弱。由此造成接收信号的幅度变化，这就是多径衰落。
眼图
[https://www.zhihu.com/question/50907793]
眼图是一种用于观察数字信号在传输过程中随时间变化的图形。通过眼图，本研究可以直观地了解到数字信号的许多重要指标，如幅度、时间、噪声等。以下是一些主要的观察点：
- 眼高：表示信号的平均幅度。
- 眼交叉比：眼图中“眼睛”张开的大小程度。
- 眼幅度：信号的最大和最小幅度。
- 眼电平：逻辑1和逻辑0的电平。
- 眼信噪比：表示信号的噪声水平。
- 眼垂直张开度和眼水平张开度：分别表示信号在垂直和水平方向的张开程度。
- 质量因素：衡量信号质量的一个综合指标。
除此之外，还可以观察以下参数： - 眼宽：表示数字波形的平均时间宽度。
- 眼上升时间和眼下降时间：分别表示数字波形的上升和下降速度。
- 上冲和下冲：表示数字波形在达到峰值或谷值时的超过设定电压的程度。
通过观察这些参数，本研究可以对数字信号的质量进行定性和定量的分析，进而对信号传输系统进行调整和优化。
眼图（Eye Diagram）是用余辉方式累积叠加显示采集到的串行信号的比特位的结果，叠加后的图形形状看起来和眼睛很像，故名眼图。眼图是一系列数字信号在示波器上累积而显示的图形，它包含了丰富的信息，从眼图上可以观察出码间串扰和噪声的影响，体现了数字信号整体的特征，从而估计系统优劣程度。
眼图实际上就是数字信号的一系列不同二进制码按一定的规律在示波器屏幕上累积后的显示，简单地说，由于示波器具有余辉功能，只要将捕获的所有波形按每三个比特分别地叠加累积，从而就形成了眼图。
8 种比特跳变组合及其层叠形成的眼图：
| 8种比特跳变组合及其层叠形成的眼图 |
|---|
 |
由于理想情况下数字信号是矩形，没有上升时间和下降时间，因而如上的眼图是不太像眼睛的。眼图是一种信号显示模式。这种模式不需要重复的波形，有助于确定噪声、抖动、 失真、信号强度以及许多其他测量结果。它能够给出系统性能的总体统计视图， 因为它查看的是比特流中每个比特组合的叠加结果。这种模式需要的是同步时钟信号触发。每当发生触发事件时（允许示波器重新准备触发），示波器会对数据进行采样，并且在屏幕上构建所有可能的 1 和 0 组合。全频时钟和分频时钟均可用于触发，但是，如果码型的长度是时钟分频比的偶数倍，则眼图会缺少组合，因而不完整。此外，如果将数据用作自身的触发信号，则眼图看起来可能完整，但示波器只会在数据码型的上升沿上触发。想要得到精确的眼图测试结果，本研究应当避免采用这种方式。眼图反映的是链路是传输的所有数字信号的整体特性。
- 实时波形能够反映波形的细节信息，如观察上升/下降边沿、过冲、单调性等。
- 眼图能够体现信号的整体特征。
- 实时波形很好，可以说明信号品质没有问题吗？不一定，只能代表某些比特。
- 眼图很好，可以说明信号品质没有问题吗？当然可以。代表整体。
| 眼图波形的采样过程 | 眼图的指标 | 有串扰和没有串扰的受扰信号眼图 | 现实中的眼图 |
|---|---|---|---|
 |
 |
 |
 |
实现
仿真图
| 理想信道 | 随机信道 |
|---|---|
 |
 |
仿真图生成的代码
- 理想信道
#!/usr/bin/env python3
# -*- coding: utf-8 -*-
#
# SPDX-License-Identifier: GPL-3.0
#
# GNU Radio Python Flow Graph
# Title: final-base-yummy
# Author: qsbye
# GNU Radio version: 3.10.8.0
from PyQt5 import Qt
from gnuradio import qtgui
from gnuradio import blocks
import numpy
from gnuradio import channels
from gnuradio.filter import firdes
from gnuradio import digital
from gnuradio import gr
from gnuradio.fft import window
import sys
import signal
from PyQt5 import Qt
from argparse import ArgumentParser
from gnuradio.eng_arg import eng_float, intx
from gnuradio import eng_notation
from gnuradio.qtgui import Range, RangeWidget
from PyQt5 import QtCore
import sip
class final_base_yummy(gr.top_block, Qt.QWidget):
def __init__(self):
gr.top_block.__init__(self, "final-base-yummy", catch_exceptions=True)
Qt.QWidget.__init__(self)
self.setWindowTitle("final-base-yummy")
qtgui.util.check_set_qss()
try:
self.setWindowIcon(Qt.QIcon.fromTheme('gnuradio-grc'))
except BaseException as exc:
print(f"Qt GUI: Could not set Icon: {str(exc)}", file=sys.stderr)
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", "final_base_yummy")
try:
geometry = self.settings.value("geometry")
if geometry:
self.restoreGeometry(geometry)
except BaseException as exc:
print(f"Qt GUI: Could not restore geometry: {str(exc)}", file=sys.stderr)
##################################################
# Variables
##################################################
self.sps = sps = 4
self.qpsk = qpsk = digital.constellation_rect([0.707+0.707j, -0.707+0.707j, -0.707-0.707j, 0.707-0.707j], [0, 1, 2, 3],
4, 2, 2, 1, 1).base()
self.nfilts = nfilts = 32
self.variable_adaptive_algorithm_0 = variable_adaptive_algorithm_0 = digital.adaptive_algorithm_cma( qpsk, .0001, 4).base()
self.samp_rate = samp_rate = 32000
self.rrc_taps = rrc_taps = firdes.root_raised_cosine(nfilts, nfilts, 1.0/float(sps), 0.35, 11*sps*nfilts)
self.phase_bw = phase_bw = 0.0628
self.excess_bw = excess_bw = 0.35
self.epsilon = epsilon = 1.0
self.delay = delay = 0
##################################################
# Blocks
##################################################
self._epsilon_range = Range(0.5, 3.0, 0.1, 1.0, 200)
self._epsilon_win = RangeWidget(self._epsilon_range, self.set_epsilon, "信道时延", "counter_slider", float, QtCore.Qt.Horizontal)
self.top_grid_layout.addWidget(self._epsilon_win, 0, 0, 1, 1)
for r in range(0, 1):
self.top_grid_layout.setRowStretch(r, 1)
for c in range(0, 1):
self.top_grid_layout.setColumnStretch(c, 1)
self.qtgui_eye_sink_x_0 = qtgui.eye_sink_c(
2048, #size
samp_rate, #samp_rate
2, #number of inputs
None
)
self.qtgui_eye_sink_x_0.set_update_time(0.10)
self.qtgui_eye_sink_x_0.set_samp_per_symbol(sps)
self.qtgui_eye_sink_x_0.set_y_axis(-1.5, 1.5)
self.qtgui_eye_sink_x_0.set_y_label('Amplitude', "")
self.qtgui_eye_sink_x_0.enable_tags(True)
self.qtgui_eye_sink_x_0.set_trigger_mode(qtgui.TRIG_MODE_FREE, qtgui.TRIG_SLOPE_POS, 0.0, 0, 0, "")
self.qtgui_eye_sink_x_0.enable_autoscale(False)
self.qtgui_eye_sink_x_0.enable_grid(False)
self.qtgui_eye_sink_x_0.enable_axis_labels(True)
self.qtgui_eye_sink_x_0.enable_control_panel(False)
labels = ['Re{receiver input}', 'Im{receiver input}', 'Re{after sync.}', 'Im{after sync.}', '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', 'blue', 'blue', 'blue', '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, 0, 1, 1, 1,
1, 1, 1, 1, 1]
markers = [-1, -1, -1, -1, -1,
-1, -1, -1, -1, -1]
for i in range(4):
if len(labels[i]) == 0:
if (i % 2 == 0):
self.qtgui_eye_sink_x_0.set_line_label(i, "Eye [Re{{Data {0}}}]".format(round(i/2)))
else:
self.qtgui_eye_sink_x_0.set_line_label(i, "Eye [Im{{Data {0}}}]".format(round((i-1)/2)))
else:
self.qtgui_eye_sink_x_0.set_line_label(i, labels[i])
self.qtgui_eye_sink_x_0.set_line_width(i, widths[i])
self.qtgui_eye_sink_x_0.set_line_color(i, colors[i])
self.qtgui_eye_sink_x_0.set_line_style(i, styles[i])
self.qtgui_eye_sink_x_0.set_line_marker(i, markers[i])
self.qtgui_eye_sink_x_0.set_line_alpha(i, alphas[i])
self._qtgui_eye_sink_x_0_win = sip.wrapinstance(self.qtgui_eye_sink_x_0.qwidget(), Qt.QWidget)
self.top_grid_layout.addWidget(self._qtgui_eye_sink_x_0_win, 1, 0, 2, 2)
for r in range(1, 3):
self.top_grid_layout.setRowStretch(r, 1)
for c in range(0, 2):
self.top_grid_layout.setColumnStretch(c, 1)
self.digital_pfb_clock_sync_xxx_0 = digital.pfb_clock_sync_ccf(sps, phase_bw, rrc_taps, nfilts, (nfilts/2), 1.5, 2)
self._delay_range = Range(0, 200, 1, 0, 200)
self._delay_win = RangeWidget(self._delay_range, self.set_delay, "Delay", "counter_slider", float, QtCore.Qt.Horizontal)
self.top_grid_layout.addWidget(self._delay_win, 0, 3, 1, 1)
for r in range(0, 1):
self.top_grid_layout.setRowStretch(r, 1)
for c in range(3, 4):
self.top_grid_layout.setColumnStretch(c, 1)
self.channels_channel_model_0 = channels.channel_model(
noise_voltage=0,
frequency_offset=0,
epsilon=1.0,
taps=[1],
noise_seed=0,
block_tags=False)
self.blocks_throttle2_0 = blocks.throttle( gr.sizeof_gr_complex*1, samp_rate, True, 0 if "auto" == "auto" else max( int(float(0.1) * samp_rate) if "auto" == "time" else int(0.1), 1) )
self.blocks_interleaved_char_to_complex_0 = blocks.interleaved_char_to_complex(False,1.0)
self.blocks_file_sink_0_0_0 = blocks.file_sink(gr.sizeof_char*1, 'final-yummy-random', False)
self.blocks_file_sink_0_0_0.set_unbuffered(True)
self.blocks_file_sink_0_0 = blocks.file_sink(gr.sizeof_char*1, 'final-yummy', False)
self.blocks_file_sink_0_0.set_unbuffered(True)
self.blocks_complex_to_interleaved_char_0 = blocks.complex_to_interleaved_char(False, 1.0)
self.analog_random_source_x_0 = blocks.vector_source_b(list(map(int, numpy.random.randint(0, 2, 1000))), True)
##################################################
# Connections
##################################################
self.connect((self.analog_random_source_x_0, 0), (self.blocks_file_sink_0_0_0, 0))
self.connect((self.analog_random_source_x_0, 0), (self.blocks_interleaved_char_to_complex_0, 0))
self.connect((self.blocks_complex_to_interleaved_char_0, 0), (self.blocks_file_sink_0_0, 0))
self.connect((self.blocks_interleaved_char_to_complex_0, 0), (self.blocks_throttle2_0, 0))
self.connect((self.blocks_throttle2_0, 0), (self.channels_channel_model_0, 0))
self.connect((self.channels_channel_model_0, 0), (self.blocks_complex_to_interleaved_char_0, 0))
self.connect((self.channels_channel_model_0, 0), (self.digital_pfb_clock_sync_xxx_0, 0))
self.connect((self.channels_channel_model_0, 0), (self.qtgui_eye_sink_x_0, 0))
self.connect((self.digital_pfb_clock_sync_xxx_0, 0), (self.qtgui_eye_sink_x_0, 1))
def closeEvent(self, event):
self.settings = Qt.QSettings("GNU Radio", "final_base_yummy")
self.settings.setValue("geometry", self.saveGeometry())
self.stop()
self.wait()
event.accept()
def get_sps(self):
return self.sps
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), 0.35, 11*self.sps*self.nfilts))
self.qtgui_eye_sink_x_0.set_samp_per_symbol(self.sps)
def get_qpsk(self):
return self.qpsk
def set_qpsk(self, qpsk):
self.qpsk = qpsk
def get_nfilts(self):
return self.nfilts
def set_nfilts(self, nfilts):
self.nfilts = nfilts
self.set_rrc_taps(firdes.root_raised_cosine(self.nfilts, self.nfilts, 1.0/float(self.sps), 0.35, 11*self.sps*self.nfilts))
def get_variable_adaptive_algorithm_0(self):
return self.variable_adaptive_algorithm_0
def set_variable_adaptive_algorithm_0(self, variable_adaptive_algorithm_0):
self.variable_adaptive_algorithm_0 = variable_adaptive_algorithm_0
def get_samp_rate(self):
return self.samp_rate
def set_samp_rate(self, samp_rate):
self.samp_rate = samp_rate
self.blocks_throttle2_0.set_sample_rate(self.samp_rate)
self.qtgui_eye_sink_x_0.set_samp_rate(self.samp_rate)
def get_rrc_taps(self):
return self.rrc_taps
def set_rrc_taps(self, rrc_taps):
self.rrc_taps = rrc_taps
self.digital_pfb_clock_sync_xxx_0.update_taps(self.rrc_taps)
def get_phase_bw(self):
return self.phase_bw
def set_phase_bw(self, phase_bw):
self.phase_bw = phase_bw
self.digital_pfb_clock_sync_xxx_0.set_loop_bandwidth(self.phase_bw)
def get_excess_bw(self):
return self.excess_bw
def set_excess_bw(self, excess_bw):
self.excess_bw = excess_bw
def get_epsilon(self):
return self.epsilon
def set_epsilon(self, epsilon):
self.epsilon = epsilon
def get_delay(self):
return self.delay
def set_delay(self, delay):
self.delay = delay
def main(top_block_cls=final_base_yummy, options=None):
if gr.enable_realtime_scheduling() != gr.RT_OK:
gr.logger("realtime").warning("Error: failed to enable real-time scheduling.")
qapp = Qt.QApplication(sys.argv)
tb = top_block_cls()
tb.start()
tb.show()
def sig_handler(sig=None, frame=None):
tb.stop()
tb.wait()
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)
qapp.exec_()
if __name__ == '__main__':
main()
- 随机信道
#!/usr/bin/env python3
# -*- coding: utf-8 -*-
#
# SPDX-License-Identifier: GPL-3.0
#
# GNU Radio Python Flow Graph
# Title: Eye Sink Example
# GNU Radio version: 3.10.8.0
from PyQt5 import Qt
from gnuradio import qtgui
from gnuradio import blocks
import numpy
from gnuradio import channels
from gnuradio.filter import firdes
from gnuradio import digital
from gnuradio import gr
from gnuradio.fft import window
import sys
import signal
from PyQt5 import Qt
from argparse import ArgumentParser
from gnuradio.eng_arg import eng_float, intx
from gnuradio import eng_notation
from gnuradio.qtgui import Range, RangeWidget
from PyQt5 import QtCore
import sip
class pyqt_eye_sink(gr.top_block, Qt.QWidget):
def __init__(self):
gr.top_block.__init__(self, "Eye Sink Example", catch_exceptions=True)
Qt.QWidget.__init__(self)
self.setWindowTitle("Eye Sink Example")
qtgui.util.check_set_qss()
try:
self.setWindowIcon(Qt.QIcon.fromTheme('gnuradio-grc'))
except BaseException as exc:
print(f"Qt GUI: Could not set Icon: {str(exc)}", file=sys.stderr)
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", "pyqt_eye_sink")
try:
geometry = self.settings.value("geometry")
if geometry:
self.restoreGeometry(geometry)
except BaseException as exc:
print(f"Qt GUI: Could not restore geometry: {str(exc)}", file=sys.stderr)
##################################################
# Variables
##################################################
self.symbol_rate = symbol_rate = 100000
self.sps = sps = 8
self.nfilts = nfilts = 32
self.alpha = alpha = 0.35
self.time_offset = time_offset = 1.00
self.taps = taps = [1.0 + 0.0j, ]
self.samp_rate = samp_rate = symbol_rate*sps
self.rrc_taps = rrc_taps = firdes.root_raised_cosine(nfilts, nfilts, 1.0/float(sps), alpha, 45*nfilts)
self.qpsk = qpsk = digital.constellation_rect([0.707+0.707j, -0.707+0.707j, -0.707-0.707j, 0.707-0.707j], [0, 1, 2, 3],
4, 2, 2, 1, 1).base()
self.qam_64 = qam_64 = digital.constellation_rect([(-7+7j),(-5+7j),(-3+7j),(-1+7j),(+1+7j),(+3+7j),(+5+7j),(+7+7j),(-7+5j),(-5+5j),(-3+5j),(-1+5j),(+1+5j),(+3+5j),(+5+5j),(+7+5j),(-7+3j),(-5+3j),(-3+3j),(-1+3j),(+1+3j),(+3+3j),(+5+3j),(+7+3j),(-7+1j),(-5+1j),(-3+1j),(-1+1j),(+1+1j),(+3+1j),(+5+1j),(+7+1j),(-7-1j),(-5-1j),(-3-1j),(-1-1j),(+1-1j),(+3-1j),(+5-1j),(+7-1j),(-7-3j),(-5-3j),(-3-3j),(-1-3j),(+1-3j),(+3-3j),(+5-3j),(+7-3j),(-7-5j),(-5-5j),(-3-5j),(-1-5j),(+1-5j),(+3-5j),(+5-5j),(+7-5j),(-7-7j),(-5-7j),(-3-7j),(-1-7j),(+1-7j),(+3-7j),(+5-7j),(+7-7j)]
, [26, 27,25, 24, 8, 9, 11, 10,30,31,29,28,12, 13,15,14,22,23,21,20,4,5,7,6,18,19,17,16,0,1,3,2,50,51,49,48,32,33,35,34,54,55,53,52,36,37,39,38,62,63,61,60,44,45,47,46,58,59,57,56,40,41,43,42],
4, 2, 2, 1, 1).base()
self.qam_16 = qam_16 = digital.constellation_rect(digital.qam_16()[0], digital.qam_16()[1],
4, 2, 2, 1, 1).base()
self.noise_volt = noise_volt = 0.0001
self.bpsk = bpsk = digital.constellation_rect([1, -1], [0, 1],
4, 2, 0, 1, 1).base()
##################################################
# Blocks
##################################################
self._time_offset_range = Range(0.999, 1.001, 0.0001, 1.00, 200)
self._time_offset_win = RangeWidget(self._time_offset_range, self.set_time_offset, "Timing Offset", "counter_slider", float, QtCore.Qt.Horizontal)
self.top_grid_layout.addWidget(self._time_offset_win, 0, 2, 1, 1)
for r in range(0, 1):
self.top_grid_layout.setRowStretch(r, 1)
for c in range(2, 3):
self.top_grid_layout.setColumnStretch(c, 1)
self._noise_volt_range = Range(0, 1, 0.01, 0.0001, 200)
self._noise_volt_win = RangeWidget(self._noise_volt_range, self.set_noise_volt, "Noise Voltage", "counter_slider", float, QtCore.Qt.Horizontal)
self.top_grid_layout.addWidget(self._noise_volt_win, 0, 0, 1, 1)
for r in range(0, 1):
self.top_grid_layout.setRowStretch(r, 1)
for c in range(0, 1):
self.top_grid_layout.setColumnStretch(c, 1)
self._alpha_range = Range(0.0001, 1, 0.01, 0.35, 200)
self._alpha_win = RangeWidget(self._alpha_range, self.set_alpha, "Roll off", "counter_slider", float, QtCore.Qt.Horizontal)
self.top_grid_layout.addWidget(self._alpha_win, 0, 1, 1, 1)
for r in range(0, 1):
self.top_grid_layout.setRowStretch(r, 1)
for c in range(1, 2):
self.top_grid_layout.setColumnStretch(c, 1)
self.qtgui_eye_sink_x_0 = qtgui.eye_sink_c(
2048, #size
samp_rate, #samp_rate
2, #number of inputs
None
)
self.qtgui_eye_sink_x_0.set_update_time(0.10)
self.qtgui_eye_sink_x_0.set_samp_per_symbol(sps)
self.qtgui_eye_sink_x_0.set_y_axis(-1.5, 1.5)
self.qtgui_eye_sink_x_0.set_y_label('Amplitude', "")
self.qtgui_eye_sink_x_0.enable_tags(True)
self.qtgui_eye_sink_x_0.set_trigger_mode(qtgui.TRIG_MODE_FREE, qtgui.TRIG_SLOPE_POS, 0.0, 0, 0, "")
self.qtgui_eye_sink_x_0.enable_autoscale(False)
self.qtgui_eye_sink_x_0.enable_grid(False)
self.qtgui_eye_sink_x_0.enable_axis_labels(True)
self.qtgui_eye_sink_x_0.enable_control_panel(False)
labels = ['Re{receiver input}', 'Im{receiver input}', 'Re{after sync.}', 'Im{after sync.}', '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', 'blue', 'blue', 'blue', '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, 0, 1, 1, 1,
1, 1, 1, 1, 1]
markers = [-1, -1, -1, -1, -1,
-1, -1, -1, -1, -1]
for i in range(4):
if len(labels[i]) == 0:
if (i % 2 == 0):
self.qtgui_eye_sink_x_0.set_line_label(i, "Eye [Re{{Data {0}}}]".format(round(i/2)))
else:
self.qtgui_eye_sink_x_0.set_line_label(i, "Eye [Im{{Data {0}}}]".format(round((i-1)/2)))
else:
self.qtgui_eye_sink_x_0.set_line_label(i, labels[i])
self.qtgui_eye_sink_x_0.set_line_width(i, widths[i])
self.qtgui_eye_sink_x_0.set_line_color(i, colors[i])
self.qtgui_eye_sink_x_0.set_line_style(i, styles[i])
self.qtgui_eye_sink_x_0.set_line_marker(i, markers[i])
self.qtgui_eye_sink_x_0.set_line_alpha(i, alphas[i])
self._qtgui_eye_sink_x_0_win = sip.wrapinstance(self.qtgui_eye_sink_x_0.qwidget(), Qt.QWidget)
self.top_grid_layout.addWidget(self._qtgui_eye_sink_x_0_win, 1, 0, 2, 2)
for r in range(1, 3):
self.top_grid_layout.setRowStretch(r, 1)
for c in range(0, 2):
self.top_grid_layout.setColumnStretch(c, 1)
self.qtgui_const_sink_x_0_0 = qtgui.const_sink_c(
2048, #size
'RC, synced', #name
1, #number of inputs
None # parent
)
self.qtgui_const_sink_x_0_0.set_update_time(0.10)
self.qtgui_const_sink_x_0_0.set_y_axis((-2), 2)
self.qtgui_const_sink_x_0_0.set_x_axis((-2), 2)
self.qtgui_const_sink_x_0_0.set_trigger_mode(qtgui.TRIG_MODE_FREE, qtgui.TRIG_SLOPE_POS, 0.0, 0, "")
self.qtgui_const_sink_x_0_0.enable_autoscale(False)
self.qtgui_const_sink_x_0_0.enable_grid(False)
self.qtgui_const_sink_x_0_0.enable_axis_labels(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", "red"]
styles = [0, 0, 0, 0, 0,
0, 0, 0, 0, 0]
markers = [0, 0, 0, 0, 0,
0, 0, 0, 0, 0]
alphas = [0.5, 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_const_sink_x_0_0.set_line_label(i, "Data {0}".format(i))
else:
self.qtgui_const_sink_x_0_0.set_line_label(i, labels[i])
self.qtgui_const_sink_x_0_0.set_line_width(i, widths[i])
self.qtgui_const_sink_x_0_0.set_line_color(i, colors[i])
self.qtgui_const_sink_x_0_0.set_line_style(i, styles[i])
self.qtgui_const_sink_x_0_0.set_line_marker(i, markers[i])
self.qtgui_const_sink_x_0_0.set_line_alpha(i, alphas[i])
self._qtgui_const_sink_x_0_0_win = sip.wrapinstance(self.qtgui_const_sink_x_0_0.qwidget(), Qt.QWidget)
self.top_grid_layout.addWidget(self._qtgui_const_sink_x_0_0_win, 2, 2, 1, 1)
for r in range(2, 3):
self.top_grid_layout.setRowStretch(r, 1)
for c in range(2, 3):
self.top_grid_layout.setColumnStretch(c, 1)
self.qtgui_const_sink_x_0 = qtgui.const_sink_c(
2048, #size
'RRC @ receiver input', #name
1, #number of inputs
None # parent
)
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", "green", "black", "cyan",
"magenta", "yellow", "dark red", "dark green", "red"]
styles = [0, 0, 0, 0, 0,
0, 0, 0, 0, 0]
markers = [0, 0, 0, 0, 0,
0, 0, 0, 0, 0]
alphas = [0.5, 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_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.qwidget(), Qt.QWidget)
self.top_grid_layout.addWidget(self._qtgui_const_sink_x_0_win, 1, 2, 1, 1)
for r in range(1, 2):
self.top_grid_layout.setRowStretch(r, 1)
for c in range(2, 3):
self.top_grid_layout.setColumnStretch(c, 1)
self.digital_pfb_clock_sync_xxx_0_0 = digital.pfb_clock_sync_ccf((sps*1.001), (6.28/100.0), rrc_taps, nfilts, (nfilts/2), 1.5, 1)
self.digital_pfb_clock_sync_xxx_0 = digital.pfb_clock_sync_ccf((sps*1.001), (6.28/100.0), rrc_taps, nfilts, (nfilts/2), 1.5, sps)
self.digital_constellation_modulator_0 = digital.generic_mod(
constellation=qpsk,
differential=False,
samples_per_symbol=sps,
pre_diff_code=True,
excess_bw=alpha,
verbose=False,
log=False,
truncate=False)
self.channels_channel_model_0 = channels.channel_model(
noise_voltage=noise_volt,
frequency_offset=0,
epsilon=time_offset,
taps=taps,
noise_seed=0,
block_tags=False)
self.blocks_throttle_0 = blocks.throttle(gr.sizeof_gr_complex*1, samp_rate,True)
self.analog_random_source_x_0 = blocks.vector_source_b(list(map(int, numpy.random.randint(0, 256, 10000))), True)
##################################################
# Connections
##################################################
self.connect((self.analog_random_source_x_0, 0), (self.digital_constellation_modulator_0, 0))
self.connect((self.blocks_throttle_0, 0), (self.channels_channel_model_0, 0))
self.connect((self.channels_channel_model_0, 0), (self.digital_pfb_clock_sync_xxx_0, 0))
self.connect((self.channels_channel_model_0, 0), (self.digital_pfb_clock_sync_xxx_0_0, 0))
self.connect((self.channels_channel_model_0, 0), (self.qtgui_const_sink_x_0, 0))
self.connect((self.channels_channel_model_0, 0), (self.qtgui_eye_sink_x_0, 0))
self.connect((self.digital_constellation_modulator_0, 0), (self.blocks_throttle_0, 0))
self.connect((self.digital_pfb_clock_sync_xxx_0, 0), (self.qtgui_eye_sink_x_0, 1))
self.connect((self.digital_pfb_clock_sync_xxx_0_0, 0), (self.qtgui_const_sink_x_0_0, 0))
def closeEvent(self, event):
self.settings = Qt.QSettings("GNU Radio", "pyqt_eye_sink")
self.settings.setValue("geometry", self.saveGeometry())
self.stop()
self.wait()
event.accept()
def get_symbol_rate(self):
return self.symbol_rate
def set_symbol_rate(self, symbol_rate):
self.symbol_rate = symbol_rate
self.set_samp_rate(self.symbol_rate*self.sps)
def get_sps(self):
return self.sps
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.alpha, 45*self.nfilts))
self.set_samp_rate(self.symbol_rate*self.sps)
self.qtgui_eye_sink_x_0.set_samp_per_symbol(self.sps)
def get_nfilts(self):
return self.nfilts
def set_nfilts(self, nfilts):
self.nfilts = nfilts
self.set_rrc_taps(firdes.root_raised_cosine(self.nfilts, self.nfilts, 1.0/float(self.sps), self.alpha, 45*self.nfilts))
def get_alpha(self):
return self.alpha
def set_alpha(self, alpha):
self.alpha = alpha
self.set_rrc_taps(firdes.root_raised_cosine(self.nfilts, self.nfilts, 1.0/float(self.sps), self.alpha, 45*self.nfilts))
def get_time_offset(self):
return self.time_offset
def set_time_offset(self, time_offset):
self.time_offset = time_offset
self.channels_channel_model_0.set_timing_offset(self.time_offset)
def get_taps(self):
return self.taps
def set_taps(self, taps):
self.taps = taps
self.channels_channel_model_0.set_taps(self.taps)
def get_samp_rate(self):
return self.samp_rate
def set_samp_rate(self, samp_rate):
self.samp_rate = samp_rate
self.blocks_throttle_0.set_sample_rate(self.samp_rate)
self.qtgui_eye_sink_x_0.set_samp_rate(self.samp_rate)
def get_rrc_taps(self):
return self.rrc_taps
def set_rrc_taps(self, rrc_taps):
self.rrc_taps = rrc_taps
self.digital_pfb_clock_sync_xxx_0.update_taps(self.rrc_taps)
self.digital_pfb_clock_sync_xxx_0_0.update_taps(self.rrc_taps)
def get_qpsk(self):
return self.qpsk
def set_qpsk(self, qpsk):
self.qpsk = qpsk
def get_qam_64(self):
return self.qam_64
def set_qam_64(self, qam_64):
self.qam_64 = qam_64
def get_qam_16(self):
return self.qam_16
def set_qam_16(self, qam_16):
self.qam_16 = qam_16
def get_noise_volt(self):
return self.noise_volt
def set_noise_volt(self, noise_volt):
self.noise_volt = noise_volt
self.channels_channel_model_0.set_noise_voltage(self.noise_volt)
def get_bpsk(self):
return self.bpsk
def set_bpsk(self, bpsk):
self.bpsk = bpsk
def main(top_block_cls=pyqt_eye_sink, options=None):
qapp = Qt.QApplication(sys.argv)
tb = top_block_cls()
tb.start()
tb.show()
def sig_handler(sig=None, frame=None):
tb.stop()
tb.wait()
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)
qapp.exec_()
if __name__ == '__main__':
main()
效果
| 理想信道眼图 | 随机信道眼图 |
|---|---|
 |
 |