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

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

前言

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 word
    for (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 items
    int 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 数据格式大致如下面形式:

目录
相关文章
|
7月前
|
数据可视化 物联网 Python
GNU Radio简介及流程图搭建
GNU Radio简介及流程图搭建
202 0
|
7月前
|
数据可视化 关系型数据库 编译器
【C/C++ 单线程性能分析工具 Gprof】 GNU的C/C++ 性能分析工具 Gprof 使用全面指南
【C/C++ 单线程性能分析工具 Gprof】 GNU的C/C++ 性能分析工具 Gprof 使用全面指南
1099 2
|
7月前
解决GNU Radio+USRP实现OFDM收发在接收端QPSK星座图映射无“抖动”问题
解决GNU Radio+USRP实现OFDM收发在接收端QPSK星座图映射无“抖动”问题
191 0
|
7月前
|
Python
GNU Radio创建Zadoff-Chu序列python OOT块
GNU Radio创建Zadoff-Chu序列python OOT块
92 0
|
C++
warning: non-static data member initializers only available with -std=c++11 or -std=gnu++11
warning: non-static data member initializers only available with -std=c++11 or -std=gnu++11
201 0
|
IDE Linux 开发工具
Linux入门学习教程:GNU C及将Vim打造成C/C++的半自动化IDE
C语言在Linux系统中的重要性自然是无与伦比、不可替代,所以我写Linux江湖系列不可能不提C语言。C语言是我的启蒙语言,感谢C语言带领我进入了程序世界。虽然现在不靠它吃饭,但是仍免不了经常和它打交道,特别是在Linux系统下。
1301 0