GNU Radio之OFDM Carrier Allocator底层C++实现

文章目录

  • 前言
  • 一、OFDM Carrier Allocator 简介
  • 二、底层 C++ 实现
    • 1、make 函数
    • 2、ofdm_carrier_allocator_cvc_impl 函数
    • 3、calculate_output_stream_length 函数
    • 4、work 函数
    • 5、~ofdm_carrier_allocator_cvc_impl 函数
  • 三、OFDM 数据格式


前言

OFDM Carrier Allocator 是 OFDM 子载波分配模块,也即串并转换模块。该模块的作用是给每个子载波分配相应的值,数据相应地实现串并转换。本文记录其底层 C++ 代码实现。


一、OFDM Carrier Allocator 简介

OFDM Carrier Allocator:OFDM 子载波分配模块,也即串并转换模块。该模块的作用是给每个子载波分配相应的值,数据相应地实现串并转换。

双击 OFDM Carrier Allocator 这个模块,会出现下面这样的属性图:
在这里插入图片描述
官方示例中这里设总载波个数为 64,载波类型包括以下几种:数据子载波(Occupied carriers)、导频子载波(Pilot carriers)、虚拟子载波(Virtual carriers, VC)和直流子载波(DC carrier)。 其中数据子载波可以用来携带数据和同步字,导频载波用于发射特定符号数据,空载波的值默认为 0。在试验中发送端 OFDM 子载波的编号规则是:为了保持左右对称,将 64 个子载波编号为 [-32,31]。其中,最左边 6 个子载波 [-32,-27] 和最右边的 5 个子载波 [26,31] 称为虚拟载波,即不使用的载波它的值默认为 0,其作用是为了减少带外功率和抑制 ACI;编号为[-26,-22]、[-20,-8]、[-6,-1]、[1,6]、[8,20]、[22,26] 都是数据子载波;编号为 -21、-7、7、21 为导频子载波;编号为 0 为直流子载波,作用是降低 OFDM 信号的峰值信噪比以及硬件电路的复杂度;

由此可知,64 个子载波中有效子载波只有 48 个这样我们便可以计算出承载一帧有效数据(包括帧头和负载共 448 个符号)所需要的 OFDM 符号数(将 64 个子载波组成的并行数据称为一个 OFDM 符号),即 448/48=9.33 个 OFDM 符号,也就是说传输一帧有效数据需要 10 个 OFDM 符号。另外,为了接收端解调方便,实验中的 OFDM 信号需要加入 2 个符号的同步字置于一帧有效数据的开头。所以,最后传输的一帧数据中包含了 12 个 OFDM 字符, 2 个同步字和 10 个有效数据

二、底层 C++ 实现

根据维基百科官方资料可以看到其实现的 C++ 代码所在的源文件:OFDM_Carrier_Allocator
在这里插入图片描述
这里我们重点查看 ofdm_carrier_allocator_cvc_impl.cc 文件的内容,使用 source insight 我们可以看到这 5 个实现函数:
在这里插入图片描述
当你在 GNU Radio Companion (GRC) 图中使用 OFDM Carrier Allocator 模块时,模块的生命周期会遵循特定的调用顺序,涉及你提到的函数。这些函数的调用顺序反映了模块从创建到执行再到销毁的整个过程。以下是一个高层次的概览:

    1. make()
    • 作用:这是一个静态工厂方法,用于创建 OFDM Carrier Allocator 模块的实例。在 GRC 生成的 Python 代码中,当流图初始化并创建模块实例时,会首先调用此方法
    • 调用时机:流图初始化时。
    1. ofdm_carrier_allocator_cvc_impl()
    • 作用:这是构造函数,由 make() 方法内部调用,用于初始化模块实例。它设置模块的初始状态,包括配置参数和内部变量的初始化。
    • 调用时机:紧随 make() 方法之后。
    1. calculate_output_stream_length(const gr_vector_int& ninput_items)
    • 作用:这个方法计算输出流的长度,基于输入项的数量。它用于确保流控制逻辑正确处理输入和输出之间的大小关系,尽管它可能不在所有模块中都显式实现。
    • 调用时机:这个方法在处理数据(即调用 work() 方法)之前被调用,但具体是否被调用以及调用时机取决于 GNU Radio 框架的版本和具体实现。
    1. work(int noutput_items, gr_vector_int& ninput_items, gr_vector_const_void_star& input_items, gr_vector_void_star& output_items)
    • 作用:这是核心处理函数,负责实际的数据处理逻辑,包括读取输入数据,执行 OFDM 载波分配,以及输出处理后的数据
    • 调用时机:流图运行时,反复调用。每次处理一批数据时都会调用。
    1. ~ofdm_carrier_allocator_cvc_impl()
    • 作用:析构函数,用于清理资源,如内部分配的内存等。
    • 调用时机:模块销毁时调用,比如流图停止运行并进行清理时。

make() 开始,通过构造函数初始化,到 work() 方法的重复调用处理数据,再到析构函数的清理工作,这一系列函数定义了 OFDM Carrier Allocator 模块在 GNU Radio 流图中的生命周期和行为。

1、make 函数

ofdm_carrier_allocator_cvc::sptr ofdm_carrier_allocator_cvc::make(int fft_len,const std::vector<std::vector<int>>& occupied_carriers,const std::vector<std::vector<int>>& pilot_carriers,const std::vector<std::vector<gr_complex>>& pilot_symbols,const std::vector<std::vector<gr_complex>>& sync_words,const std::string& len_tag_key,const bool output_is_shifted)
{return gnuradio::get_initial_sptr(new ofdm_carrier_allocator_cvc_impl(fft_len,occupied_carriers,pilot_carriers,pilot_symbols,sync_words,len_tag_key,output_is_shifted));
}

这段代码展示了 GNU Radio 中一个典型的工厂方法模式实现,用于创建和初始化 ofdm_carrier_allocator_cvc 模块的一个实例。这个方法属于 ofdm_carrier_allocator_cvc 类,是一个静态成员函数,意味着它可以不通过对象实例而直接调用。这个方法的目的是简化对象的创建和管理,同时隐藏具体实现类 ofdm_carrier_allocator_cvc_impl 的细节。

  • 函数参数
    • fft_len:整数,指定FFT(快速傅立叶变换)的长度。
    • occupied_carriers:一个二维整数向量,指定了哪些载波被用于传输数据。
    • pilot_carriers:一个二维整数向量,指定了哪些载波被用作导频。
    • pilot_symbols:一个二维 gr_complex 向量,包含导频载波上的符号。
    • sync_words:一个二维 gr_complex 向量,包含用于帧同步的同步字。
    • len_tag_key:字符串,指定了长度标签的键,用于流标签系统中。
    • output_is_shifted:布尔值,指示输出是否应该被频移。
  • 函数实现
    函数体内部,make 方法通过调用 new 来创建 ofdm_carrier_allocator_cvc_impl 的一个新实例,将所有接收到的参数传递给其构造函数。然后,它使用 gnuradio::get_initial_sptr 函数将裸指针包装为智能指针(gnuradio::get_initial_sptr 是 GNU Radio中用于创建智能指针的辅助函数,保证了对象的安全管理和自动内存管理)。
  • 智能指针
    返回值是 ofdm_carrier_allocator_cvc::sptr,即指向 ofdm_carrier_allocator_cvc 的智能指针。在 GNU Radio 中,智能指针主要用于管理模块实例的生命周期,避免内存泄露。智能指针在对象不再被使用时会自动释放关联的资源。

2、ofdm_carrier_allocator_cvc_impl 函数

这段代码负责在 OFDM 系统中配置和准备用于传输的帧,包括数据载波、导频载波及其符号,以及同步字的分配。通过精确控制这些参数,可以优化信号的传输效率和质量,同时简化接收端的信号处理流程。

ofdm_carrier_allocator_cvc_impl::ofdm_carrier_allocator_cvc_impl(// 参数说明:int fft_len,												// FFT(快速傅里叶变换)的长度。const std::vector<std::vector<int>>& occupied_carriers,		// 被占用的载波,即用于传输数据的频率位置。const std::vector<std::vector<int>>& pilot_carriers,		// 导频载波,用于辅助接收端的信道估计。const std::vector<std::vector<gr_complex>>& pilot_symbols,	// 导频符号,已知的信号,用于帮助接收端进行信道估计。const std::vector<std::vector<gr_complex>>& sync_words,		// 同步字,用于帧同步。const std::string& len_tag_key,								// 长度标签键,用于标记流的长度。const bool output_is_shifted)								// 指示输出是否被移位的布尔值。// 构造函数初始化列表 tagged_stream_block 表明该组件是一个标记流块,用于处理带有特定标记的数据流。: tagged_stream_block("ofdm_carrier_allocator_cvc",io_signature::make(1, 1, sizeof(gr_complex)),io_signature::make(1, 1, sizeof(gr_complex) * fft_len),len_tag_key),// 成员变量初始化:d_fft_len(fft_len),						// FFT长度d_occupied_carriers(occupied_carriers),	// 被占用的载波d_pilot_carriers(pilot_carriers),			// 导频载波d_pilot_symbols(pilot_symbols),			// 导频符号d_sync_words(sync_words),					// 同步字d_symbols_per_set(0),						// 记录在一个OFDM符号集中,所有被占用的载波上分配的数据和导频符号的总数;"集"指的是一个OFDM帧d_output_is_shifted(output_is_shifted)	// 指示输出是否被移位的布尔值。
{// Sanity checks// If that is is null, the input is wrong -> force user to use ((),) in python// 检查被占用的载波是否为空,如果是,则抛出异常。if (d_occupied_carriers.empty()) {throw std::invalid_argument("Occupied carriers must be of type vector of vector i.e. ((),).");}// 检查每个数据载波索引是否在有效范围内(即非负且小于FFT长度),并对于output_is_shifted为真的情况,调整载波索引for (unsigned i = 0; i < d_occupied_carriers.size(); i++) {for (unsigned j = 0; j < d_occupied_carriers[i].size(); j++) {if (occupied_carriers[i][j] < 0) {d_occupied_carriers[i][j] += d_fft_len;}if (d_occupied_carriers[i][j] > d_fft_len || d_occupied_carriers[i][j] < 0) {throw std::invalid_argument("data carrier index out of bounds");}if (d_output_is_shifted) {d_occupied_carriers[i][j] =(d_occupied_carriers[i][j] + fft_len / 2) % fft_len;}}}// 检查  导频载波是否为空,如果是,则抛出异常。if (d_pilot_carriers.empty()) {throw std::invalid_argument("Pilot carriers must be of type vector of vector i.e. ((),).");}// 检查每个导频载波索引是否在有效范围内(即非负且小于FFT长度),并对于output_is_shifted为真的情况,调整载波索引for (unsigned i = 0; i < d_pilot_carriers.size(); i++) {for (unsigned j = 0; j < d_pilot_carriers[i].size(); j++) {if (d_pilot_carriers[i][j] < 0) {d_pilot_carriers[i][j] += d_fft_len;}if (d_pilot_carriers[i][j] > d_fft_len || d_pilot_carriers[i][j] < 0) {throw std::invalid_argument("pilot carrier index out of bounds");}if (d_output_is_shifted) {d_pilot_carriers[i][j] = (d_pilot_carriers[i][j] + fft_len / 2) % fft_len;}}}// 检查导频符号是否为空,如果是,则抛出异常。if (d_pilot_symbols.empty()) {throw std::invalid_argument("Pilot symbols must be of type vector of vector i.e. ((),).");}// 确保每个导频载波的大小与相应的导频符号向量的大小匹配。for (unsigned i = 0; i < std::max(d_pilot_carriers.size(), d_pilot_symbols.size());i++) {if (d_pilot_carriers[i % d_pilot_carriers.size()].size() !=d_pilot_symbols[i % d_pilot_symbols.size()].size()) {throw std::invalid_argument("pilot_carriers do not match pilot_symbols");}}// 验证同步字的长度与 FFT 长度相匹配。for (unsigned i = 0; i < d_sync_words.size(); i++) {if (d_sync_words[i].size() != (unsigned)d_fft_len) {throw std::invalid_argument("sync words must be fft length");}}// 计算每组被占用的符号总数。for (unsigned i = 0; i < d_occupied_carriers.size(); i++) {d_symbols_per_set += d_occupied_carriers[i].size();}set_tag_propagation_policy(TPP_DONT);	// 置标记传播策略为 TPP_DONT,不传播任何标记。set_relative_rate((uint64_t)d_symbols_per_set, (uint64_t)d_occupied_carriers.size());	// 设置相对速率,影响处理数据的速率。
}

3、calculate_output_stream_length 函数

这段代码是 ofdm_carrier_allocator_cvc_impl 类的一个成员函数,用于计算输出流的长度,基于输入项的数量。这个函数对于确保 OFDM 调制过程中输入和输出流的长度匹配至关重要。

int ofdm_carrier_allocator_cvc_impl::calculate_output_stream_length(const gr_vector_int& ninput_items) 	// 一个整数向量,其中每个元素代表对应输入流中的项目数目; 对于大多数基于流的处理模块,这个向量可能只有一个元素,因为它们通常只处理一个输入流。
{int nin = ninput_items[0];			// 表示输入流中的项目数int nout = (nin / d_symbols_per_set) * d_occupied_carriers.size();	// 计算的输出项的基本数目int k = 0;// 处理余下的输入项; 使用一个循环来处理不足以形成一个完整 OFDM 符号集的余下输入项。for (int i = 0; i < nin % d_symbols_per_set; k++) {	// i 变量用于累加当前处理的输入项数目,而 k 变量用于迭代被占用的载波集合。nout++;i += d_occupied_carriers[k % d_occupied_carriers.size()].size();	// 确保了即使在处理最后一批输入项时,也能正确考虑每个被占用载波集合的大小}return nout + d_sync_words.size();	// 这个长度包括了基于输入项的计算结果,以及为同步字预留的额外输出项
}

这个函数通过精确计算在给定输入项数目的情况下,需要生成多少个输出项来确保 OFDM 调制的正确性。它不仅考虑了每个 OFDM 符号集内的符号总数和被占用载波的集合大小,还考虑了同步字的需求。通过这种方式,ofdm_carrier_allocator_cvc_impl 能够确保输出流长度恰当地匹配输入流长度和 OFDM 调制过程的需求。

示例中 64 个子载波编号为 [-32,31],将其变为 [0, 63]
在这里插入图片描述

4、work 函数

work 函数负责在 OFDM 调制过程中的载波分配工作,包括同步字、数据符号和导频符号的映射。通过这种方式,确保了 OFDM 信号能够被正确地组装并准备好发送。这个过程不仅涉及到信号处理的实际操作,还包括流标签的管理,以确保流中的元数据被正确处理和传递。

int ofdm_carrier_allocator_cvc_impl::work(int noutput_items,gr_vector_int& ninput_items,gr_vector_const_void_star& input_items,gr_vector_void_star& output_items)
{// 输入输出数据类型转换,以便于处理复数数据const gr_complex* in = (const gr_complex*)input_items[0];gr_complex* out = (gr_complex*)output_items[0];std::vector<tag_t> tags;memset((void*)out, 0x00, sizeof(gr_complex) * d_fft_len * noutput_items);	// 初始化输出缓冲区// Copy Sync wordfor (unsigned i = 0; i < d_sync_words.size(); i++) {memcpy((void*)out, (void*)&d_sync_words[i][0], sizeof(gr_complex) * d_fft_len);	// 将预定义的同步字复制到输出流的开始部分,这些同步字用于接收端的帧同步。out += d_fft_len;}// Copy data symbols// 将输入数据符号映射到被占用的 OFDM 载波上。这一过程考虑了载波的集合配置和当前处理的数据符号集合。long n_ofdm_symbols = 0; // Number of output itemsint curr_set = 0;int symbols_to_allocate = d_occupied_carriers[0].size();int symbols_allocated = 0;for (int i = 0; i < ninput_items[0]; i++) {if (symbols_allocated == 0) {// Copy all tags associated with these input symbols onto this OFDM symbol// 复制与输入符号相关联的所有标签到输出 OFDM 符号上,保持流标签系统中的信息。get_tags_in_range(tags,0,nitems_read(0) + i,nitems_read(0) + std::min(i + symbols_to_allocate, (int)ninput_items[0]));for (unsigned t = 0; t < tags.size(); t++) {add_item_tag(0,nitems_written(0) + n_ofdm_symbols +(n_ofdm_symbols == 0 ? 0 : d_sync_words.size()),tags[t].key,tags[t].value);}n_ofdm_symbols++;}out[(n_ofdm_symbols - 1) * d_fft_len +d_occupied_carriers[curr_set][symbols_allocated]] = in[i];symbols_allocated++;if (symbols_allocated == symbols_to_allocate) {curr_set = (curr_set + 1) % d_occupied_carriers.size();symbols_to_allocate = d_occupied_carriers[curr_set].size();symbols_allocated = 0;}}// Copy pilot symbols// 在数据符号映射完成后,将导频符号映射到指定的导频载波上。这些导频符号对于接收端的信道估计至关重要。for (int i = 0; i < n_ofdm_symbols; i++) {for (unsigned k = 0; k < d_pilot_carriers[i % d_pilot_carriers.size()].size();k++) {out[i * d_fft_len + d_pilot_carriers[i % d_pilot_carriers.size()][k]] =d_pilot_symbols[i % d_pilot_symbols.size()][k];}}return n_ofdm_symbols + d_sync_words.size();	// 返回处理完成的OFDM符号数量,包括数据符号和同步字。
}

5、~ofdm_carrier_allocator_cvc_impl 函数

析构函数,用于清理资源,如内部分配的内存等。

ofdm_carrier_allocator_cvc_impl::~ofdm_carrier_allocator_cvc_impl() {}

三、OFDM 数据格式

经过子载波分配模块和 IFFT 后的 OFDM 数据格式大致如下面形式:
在这里插入图片描述


我的qq:2442391036,欢迎交流!


本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.hqwc.cn/news/571886.html

如若内容造成侵权/违法违规/事实不符,请联系编程知识网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

Docker进阶:Docker Swarm —弹性伸缩调整服务的副本数量

Docker进阶&#xff1a;Docker Swarm —弹性伸缩调整服务的副本数量 1、 创建一个Nginx服务&#xff08;Manager节点&#xff09;2、查看服务状态&#xff08;Manager节点&#xff09;3、测试访问&#xff08;Worker节点&#xff09;4、查看服务日志&#xff08;Manager节点&am…

mysql 高阶语句 与视图

目录 一 前言 二 msql 高阶语句使用方法 &#xff08;一&#xff09; 查询并排序&#xff08;order by&#xff09; 1&#xff0c;排序方式 2&#xff0c;查询指定列 并排序 3, 条件判断&#xff08;过滤指定行&#xff09; 再查询指定列 并排序 4&#xff0c;查…

IDEA的Scala环境搭建

目录 前言 Scala的概述 Scala环境的搭建 一、配置Windows的JAVA环境 二、配置Windows的Scala环境 编写一个Scala程序 前言 学习Scala最好先掌握Java基础及高级部分知识&#xff0c;文章正文中会提到Scala与Java的联系&#xff0c;简单来讲Scala好比是Java的加强版&#x…

IDEA2023版本创建spring boot项目时,Java版本无法选择Java8问题解决

先简单说下出现本问题的原因&#xff1a; spring boot3.0发布时提到未来Java17将会成为主流版本&#xff0c;所有的Java EE Api都需要迁移到Jakarta EE上来。而spring boot3.0及以上版本已经不支持Java8了&#xff0c;支持Java17及以上版本。同时官方支持项目初始化的 Spring B…

智慧公厕的全域感知、全网协同、全业务融合和全场景智慧赋能

公共厕所是城市的重要组成部分&#xff0c;为市民提供基本的生活服务。然而&#xff0c;传统的公厕管理模式存在诸多问题&#xff0c;如排队等候时间长、卫生状况差、空气质量差等&#xff0c;严重影响市民的出行和生活质量。为了解决这些问题&#xff0c;智慧公厕应运而生&…

以智慧公厕建设助推城市的高质量发展

近年来&#xff0c;随着城市化进程的加快&#xff0c;城市基础设施建设日益完善&#xff0c;其中智慧公厕的建设成为了城市高质量发展的重要组成部分。智慧公厕以其智能化管理、数字化使用和信息化运行的特点&#xff0c;将公共厕所的管理水平提升到了一个全新的高度&#xff0…

YOLOv9改进策略:IoU优化 | Wasserstein Distance Loss,助力小目标涨点

&#x1f4a1;&#x1f4a1;&#x1f4a1;本文独家改进&#xff1a;基于Wasserstein距离的小目标检测评估方法 Wasserstein Distance Loss | 亲测在多个数据集能够实现涨点&#xff0c;对小目标、遮挡物性能提升明显 &#x1f4a1;&#x1f4a1;&#x1f4a1;MS COCO和PASC…

element-ui实现各种证件照上传预览下载组件封装,图片上传回显及长宽自定义功能单个图片上传功能附带源码

element-ui实现证件照上传预览下载组件封装 效果&#xff1a; 参数说明 我只写了两个参数&#xff0c;后续有需求再对其组件进行丰富~ 参数说明fileListProp用来存储上传后后端返回的图片UR了uploadUrl图片上传返回的URL后端接口地址widthProp图片上传框的宽度heightProp图片…

Bash and a Tough Math Puzzle 线段树维护区间gcd

还是一道很不错的题目&#xff0c;很容易想到用一棵线段树来维护区间gcd 注意用倍数来剪枝就好了&#xff0c;很是一到很好的题目的 #include<iostream> #include<vector> using namespace std; const int N 5e510; int n,q; struct Segment{int l,r;int d; }tr[…

智慧体育场馆的优势都有哪些?

体育场馆作为体育产业和事业发展的重要载体&#xff0c;全民对健康和运动的需求越来越大&#xff0c;体育馆的需求也更大。而以前的体育场馆管理不仅人工成本高&#xff0c;人民的使用和消费也不方便。因此智慧体育馆的出现大大降低了运营人力成本及现金管理风险&#xff0c;大…

20240320-2-线性回归+逻辑回归

线性回归于逻辑回归面试题 1. 简单介绍一下线性回归。 **线性回归(Linear Regression)是利用称为线性回归方程的最小平方函数对一个或多个自变量和因变量之间关系进行建模的一种回归分析。**这种函数是一个或多个称为回归系数的模型参数的线性组合。只有一个自变量的情况称为简…

day35 贪心算法 part04● 860.柠檬水找零 ● 406.根据身高重建队列 ● 452. 用最少数量的箭引爆气球

一遍过&#xff0c;就是有10的时候尽量把10用掉&#xff0c;20是用不掉的&#xff0c;如果发现不够找回了&#xff0c;就return false。 class Solution { public:bool lemonadeChange(vector<int>& bills) {unordered_map<int,int> mp;for(int i0;i<bill…