NR IQ 比较工具

已于 2022-04-03 11:16:35 修改
阅读量4.9k 收藏
点赞数
分类专栏： NR/LTE 文章标签： NR IQ 比较工具
于 2022-04-02 23:21:08 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/xnxygxp/article/details/123931158
版权

文章目录

需求说明

将生成的iq bin文件跟ref bin文件一一对比，允许精度误差为3bit。
需要满足多cc的比较，并且将输出的IQ 错误点，按slot，symb，ant，I_value，Q_value输出到文件中。

C语言实现示例代码

#include <stdio.h>
#include <stdint.h>
#include <string.h>
#include <unistd.h>

//macro
//C:/Users/xueping.guo/Desktop/testVec_tdd_case1case94_passLog/dspVec/case94
//C:/Users/xueping.guo/Desktop/testVec/caseD20
#define TEST_CASE_BASE_PATH            "C:/Users/xueping.guo/Desktop/testVec/case159"
#define DBG_INFO                       printf
#define NR_DL_MAX_RB_NUM              (273)
#define NR_SC_NUM_PER_RB              (12)
#define NR_MAX_SYMB_NUM_PER_SUBFRAM   (14)
#define NR_MAX_ANT_NUM                (4)  //maybe larger
#define NR_MAX_SLOT_PER_RADIO         (20)
#define IQ_DATA_SIZE_PER_ANT           (NR_DL_MAX_RB_NUM * NR_SC_NUM_PER_RB * sizeof(uint32_t))


//declaration function
int32_t loadFileData(uint8_t *dstBufP, char* srcFileName);

//global variable
uint8_t phyIqDataArry[IQ_DATA_SIZE_PER_ANT * NR_MAX_ANT_NUM * NR_MAX_SYMB_NUM_PER_SUBFRAM * NR_MAX_SLOT_PER_RADIO];
uint8_t refIqDataArry[IQ_DATA_SIZE_PER_ANT * NR_MAX_ANT_NUM * NR_MAX_SYMB_NUM_PER_SUBFRAM * NR_MAX_SLOT_PER_RADIO];


int main(int argc, char* argv[])
{
    uint8_t* phyIqDataBufP = phyIqDataArry;
    uint8_t* refIqDataBufP = refIqDataArry;
    uint32_t phyIqDataSize;
    uint32_t refIqDataSize;
    char srcFileNameArry[100];
    char logFileNameArry[100];
    uint16_t realAntNum = 0;
    uint8_t totalCcNum = 0;

    //test case base address
    char* testCaseBaseAddrP = argv[1];
    
    if (NULL == testCaseBaseAddrP)
    {
        DBG_INFO("[ERROR:] Test case base addr is null, please input correct path!!!\n");
        return 0;
    }
    else
    {
        DBG_INFO("[INFO:] TestCase base address is: %s\n", testCaseBaseAddrP);
    }

    //calc total cc num
    memset(srcFileNameArry, 0, sizeof(srcFileNameArry));
    sprintf(srcFileNameArry, "%s%s", testCaseBaseAddrP, "/out");
    //cc0
    if(access(srcFileNameArry, F_OK) == 0)
    {
        totalCcNum++;
    }

    memset(srcFileNameArry, 0, sizeof(srcFileNameArry));
    sprintf(srcFileNameArry, "%s%s", testCaseBaseAddrP, "/out1");
    //cc1
    if(access(srcFileNameArry, F_OK) == 0)
    {
        totalCcNum++;
    }

    memset(srcFileNameArry, 0, sizeof(srcFileNameArry));
    sprintf(srcFileNameArry, "%s%s", testCaseBaseAddrP, "/out2");
    //cc2
    if(access(srcFileNameArry, F_OK) == 0)
    {
        totalCcNum++;
    }

    memset(srcFileNameArry, 0, sizeof(srcFileNameArry));
    sprintf(srcFileNameArry, "%s%s", testCaseBaseAddrP, "/out3");
    //cc3
    if(access(srcFileNameArry, F_OK) == 0)
    {
        totalCcNum++;
    }
    
    for (uint8_t ccIdx = 0; ccIdx < totalCcNum; ++ccIdx)
    {
        //initialize flag per ccIdx
        uint8_t iqPassFlag = 1;

        if (0 == ccIdx)
        {
            //read phy iq data from file dlSymbPortData4g.bin
            memset(srcFileNameArry, 0, sizeof(srcFileNameArry));
            sprintf(srcFileNameArry, "%s%s", testCaseBaseAddrP, "/out/dlSymbPortData.bin");
            phyIqDataSize = (uint32_t)loadFileData(phyIqDataBufP, srcFileNameArry);

            //read ref iq data
            memset(srcFileNameArry, 0, sizeof(srcFileNameArry));
            sprintf(srcFileNameArry, "%s%s", testCaseBaseAddrP, "/ref/dlSymbPortData.bin");
            refIqDataSize = (uint32_t)loadFileData(refIqDataBufP, srcFileNameArry);
        }
        else if (1 == ccIdx)
        {
            memset(srcFileNameArry, 0, sizeof(srcFileNameArry));
            sprintf(srcFileNameArry, "%s%s", testCaseBaseAddrP, "/out1/dlSymbPortData.bin");
            phyIqDataSize = (uint32_t)loadFileData(phyIqDataBufP, srcFileNameArry);

            //read ref iq data
            memset(srcFileNameArry, 0, sizeof(srcFileNameArry));
            sprintf(srcFileNameArry, "%s%s", testCaseBaseAddrP, "/ref1/dlSymbPortData.bin");
            refIqDataSize = (uint32_t)loadFileData(refIqDataBufP, srcFileNameArry);
        }
        else if (2 == ccIdx)
        {
            memset(srcFileNameArry, 0, sizeof(srcFileNameArry));
            sprintf(srcFileNameArry, "%s%s", testCaseBaseAddrP, "/out2/dlSymbPortData.bin");
            phyIqDataSize = (uint32_t)loadFileData(phyIqDataBufP, srcFileNameArry);

            //read ref iq data
            memset(srcFileNameArry, 0, sizeof(srcFileNameArry));
            sprintf(srcFileNameArry, "%s%s", testCaseBaseAddrP, "/ref2/dlSymbPortData.bin");
            refIqDataSize = (uint32_t)loadFileData(refIqDataBufP, srcFileNameArry);
        }
        else if (3 == ccIdx)
        {
            memset(srcFileNameArry, 0, sizeof(srcFileNameArry));
            sprintf(srcFileNameArry, "%s%s", testCaseBaseAddrP, "/out3/dlSymbPortData.bin");
            phyIqDataSize = (uint32_t)loadFileData(phyIqDataBufP, srcFileNameArry);

            //read ref iq data
            memset(srcFileNameArry, 0, sizeof(srcFileNameArry));
            sprintf(srcFileNameArry, "%s%s", testCaseBaseAddrP, "/ref3/dlSymbPortData.bin");
            refIqDataSize = (uint32_t)loadFileData(refIqDataBufP, srcFileNameArry);
        }
        else
        {
            DBG_INFO("[ERROR:] ccIdx is invalid, ccIdx: %u\n", ccIdx);
        }

        //open compare log file
        memset(logFileNameArry, 0, sizeof(logFileNameArry));
        sprintf(logFileNameArry, "%s%s%u%s", testCaseBaseAddrP, "/iqCmpLogResult_cc", ccIdx, ".log");
        
        //open log file and clear dirty data
        FILE* logFileOpenP = fopen(logFileNameArry, "w+");

        if (phyIqDataSize != refIqDataSize)
        {
            fprintf(logFileOpenP, "[ERROR:] compare failed because of the different length of data file!!!\n");
            fprintf(logFileOpenP, "phyIqDataSize = %u, refIqDataSize = %u\n", phyIqDataSize, refIqDataSize);
            fclose(logFileOpenP);
            printf("[ERROR:] finished the compare because of the different length!!!\n");
            return 0;
        }

        //get reality ant num
        //273RB*12RE*4BYTE*1ANT*14SYM*20TTI = 3669120
        //273RB*12RE*4BYTE*2ANT*14SYM*20TTI = 7338240
        if (phyIqDataSize == 3669120)
        {
            realAntNum = 1;
            DBG_INFO("numTxPhyAnt = %u\n", realAntNum);
        }
        else if (7338240 == phyIqDataSize)
        {
            realAntNum = 2;
            DBG_INFO("numTxPhyAnt = %u\n", realAntNum);
        }
        else if (14676480 == phyIqDataSize)
        {
            realAntNum = 4;
            DBG_INFO("numTxPhyAnt = %u\n", realAntNum);
        }
        else
        {
            iqPassFlag = 0;
            DBG_INFO("cc%u IQ data file size is wrong, not support now!!!\n", ccIdx);
            return 0;
        }

        //subframe loop
        uint32_t iqDataSizeBytePerAnt = NR_DL_MAX_RB_NUM * NR_SC_NUM_PER_RB * sizeof(uint32_t);
        uint32_t totalNumSc           = NR_DL_MAX_RB_NUM * NR_SC_NUM_PER_RB;
        for (uint32_t slotIdx = 0; slotIdx < NR_MAX_SLOT_PER_RADIO; ++slotIdx)
        {
            for (uint32_t symIdx = 0; symIdx < NR_MAX_SYMB_NUM_PER_SUBFRAM; ++symIdx)
            {
                for (uint32_t antIdx = 0; antIdx < realAntNum; ++antIdx)
                {

                    uint16_t* phyIq_16bit_p = (uint16_t*)(phyIqDataBufP + slotIdx * iqDataSizeBytePerAnt * realAntNum * 14 + symIdx * iqDataSizeBytePerAnt * realAntNum + iqDataSizeBytePerAnt * antIdx);
                    uint16_t* refIq_16bit_p = (uint16_t*)(refIqDataBufP + slotIdx * iqDataSizeBytePerAnt * realAntNum * 14 + symIdx * iqDataSizeBytePerAnt * realAntNum + iqDataSizeBytePerAnt * antIdx);

                    for (uint32_t scIdx = 0; scIdx < totalNumSc; ++scIdx)
                    {
                        uint32_t  startLocation = slotIdx * iqDataSizeBytePerAnt * realAntNum * 14 + symIdx * iqDataSizeBytePerAnt * realAntNum + antIdx * iqDataSizeBytePerAnt + scIdx * 4;
                        uint16_t absVal_1 = ((phyIq_16bit_p[2 * scIdx] - refIq_16bit_p[2 * scIdx]) > 0 )? (phyIq_16bit_p[2 * scIdx] - refIq_16bit_p[2 * scIdx]) : (refIq_16bit_p[2 * scIdx] - phyIq_16bit_p[2 * scIdx]);
                        uint16_t absVal_2 = ((phyIq_16bit_p[2 * scIdx + 1] - refIq_16bit_p[2 * scIdx + 1]) > 0 )? (phyIq_16bit_p[2 * scIdx + 1] - refIq_16bit_p[2 * scIdx + 1]) : (refIq_16bit_p[2 * scIdx + 1] - phyIq_16bit_p[2 * scIdx + 1]);

                        if ((absVal_1 <= 3) && (absVal_2 <= 3))
                        {
                            continue;
                        }
                        else
                        {
                            iqPassFlag = 0;
                            //printf("[ERROR:] compare fail, slotIdx = %u, symIdx = %u, antIdx = %u, scIdx = %u\n", slotIdx, symIdx, antIdx, scIdx);
                            fprintf(logFileOpenP, "[ERROR:] cc%u compare fail, slotIdx = %u, symIdx = %u, antIdx = %u, scIdx = %u,  phy_I = %04X,  ref_I = %04X,  phy_Q = %04X,  ref_Q = %04X, location = %08X\n", ccIdx, slotIdx, symIdx, antIdx, scIdx,
                                                                                                                            phyIq_16bit_p[2 * scIdx], refIq_16bit_p[2 * scIdx],
                                                                                                                            phyIq_16bit_p[2 * scIdx + 1], refIq_16bit_p[2 * scIdx + 1],
                                                                                                                            startLocation);
                        }
                    }
                }
            }
        }

        if (iqPassFlag)
        {
            DBG_INFO("cc%u IQ compare pass!!!\n", ccIdx);
            fprintf(logFileOpenP, "cc%u IQ compare pass!!!", ccIdx);
        }
        else
        {
            DBG_INFO("cc%u IQ compare fail!!!\n", ccIdx);
            fprintf(logFileOpenP, "cc%u IQ compare fail!!!", ccIdx);
        }
    }

	//close log file
	fclose(logFileOpenP);
    return 0;
}

int32_t loadFileData(uint8_t *dstBufP, char* srcFileName)
{
	int32_t   size = -1;
	FILE*     srcFileP;
	size_t    readSize = 0;
	
	srcFileP = fopen(srcFileName, "rb");
	if (NULL == srcFileP)
	{
		DBG_INFO("[ERROR]: The file: %s can not be opened.\n", srcFileName);
		fclose(srcFileP);
		return 0;
	}
	
	//get file size and then read data
	fseek(srcFileP, 0, SEEK_END);
	size = ftell(srcFileP);
	fseek(srcFileP, 0, SEEK_SET);
	
	readSize = fread((void*)dstBufP, (size_t)1, (size_t)size, srcFileP);
	
	if (0 == readSize)
	{
		DBG_INFO("[ERROR]: Can not read the file %s.\n", srcFileName);
	}
	fclose(srcFileP);
	
	return size;
}
SmileToLifeForever
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
基于I/Q数据的5G控制信道盲检
weixin_70923796的博客
09-17 346
CRC校验不通过，则不存在DCI，继续盲检，直至遍历所有期待的PDCCH位置。为了研究真实环境5G下行I/Q解码流程，构建“5G核心网+5G基站+测试终端+I/Q数字采集系统”的测试环境，采集基站下行空口信号I/Q数据，并基于3GPP第15版本标准，开展5G下行控制信道PDCCH盲检流程研究，包括候选PDCCH时频域位置确定、信道估计与相关检测，PDCCH比特序列解码流程（含解扰，解速率匹配，极化码译码，循环冗余校验，信息解析等环节），展示了关键层三参数、协议流程、算法原理以及PDCCH解码星座图。
“IQ波形文件生成器”C++实现和Python实现的对比
Sagittarius_Warrior的博客
07-05 2253
本文主要借助我在项目中编写的一个小软件，来对一窥C++实现和Python实现在“编码效率”、“运行效率”、“内存占用”等方面的差异。当然，这个对比只是针对这类项目，不具普遍性。 一、问题     将两个几M大小的文本格式txt文件（一个存储的是float型的I路数据，另一个存储的是float型的Q路数据）合成一个二进制格式的bin文件（波形文件）。以本次实验的“bt_wave_i.tx
NR及LTE中的IQ数据与信息、比特率、码元、波特率之间的关系
weixin_42764060的博客
06-14 1461
信息：对信源进行数字编码后的数据，基本单位是bit。比特率：信息的速率称为比特率(bit/s、bps)，通常用Rb表示。
5G NR物理层各个步骤各项参数梳理记录
ddatalent的博客
02-22 6893
本文主要是对的一个翻译和梳理，同时引用自己项目的各项参数作为一个例子，以供日后复习，也便于大家一起学习。
IQ基带信号
Go Forward Forever!
12-31 3863
射频信号可以下变频得到较窄带宽的基带IQ信号，反之，较窄带的基带IQ信号可以上变频成射频信号发送。在IQ信号层面可以进行多样的调制、解调处理。
[4G&5G专题-9]：前传接口 CPRI与OBSAI的那些事
热门推荐
文火冰糖（王文兵）的博客
01-07 1万+
第1章. 什么是CPRI协议概述; 第2章. CPRI协议的分层架构; 第3章. CPRI协议速率与SFP的关系; 第4章. CPRI协议的管理面以及与以太网的关系 第5章. CPRI协议的用户面以及与IQ的关系 第6章. CPRI协议的同步面与时钟同步 第7章. CPRI协议的功能实现 第8章CPRI协议的不足 第9章. OBSAI...
5G Toolbox? Getting Started Guide_5Gtoolbox_5gtoolbox使用_5G_5Gma
09-11
这个工具箱提供了丰富的函数和模型，涵盖了5G新无线电（NR）的标准规范，包括物理层（PHY）和高层协议栈。下面我们将深入探讨5G Toolbox的核心功能、使用方法以及在5G通信中的应用。 1. **5G NR物理层模型**：5G ...
罗德与施瓦茨的R&STD-SCDMA测量方案全攻略
03-04
这些软件工具专门用于评估和优化基站性能，提供详细的基站发射机和接收机测试功能，确保网络质量。 总之，罗德与施瓦茨的R&STD-SCDMA测量方案涵盖了从信号生成到分析的全过程，为TD-SCDMA网络的开发、测试和优化...
AgilentN5182A.rar_agilent_matlab agilent_ofdm_ofdm agilent_psk i
07-14
- **脚本文件**：这些可能是.m文件，包含MATLAB代码，用于生成PSK、FSK和OFDM信号的参数设置、符号生成、IQ调制等过程。 - **配置文件**：可能包含设备配置参数，如频率、功率、调制类型等，以确保信号与Agilent N...
电信设备-多协议接口数字基带信道模拟器.zip
09-18
它是一种能够模拟真实通信环境中各种信道特性的工具，如衰落、干扰、噪声以及多径传播等。通过这种模拟，工程师可以在实验室环境下测试和优化通信系统的性能，而无需实际部署到现场。本资料重点讨论了多协议接口的...
基于FPGA的OFDM调制解调系统的仿真-源码
09-30
7. **软件/插件开发**：在FPGA实现OFDM系统时，可能需要使用专门的开发工具，如Xilinx的Vivado或Intel的Quartus II。此外，MATLAB或Simulink可用于设计验证和仿真，而ModelSim等仿真软件可以帮助调试代码。 8. **...
公开的基带IQ信号文件格式
JXISH001的博客
03-30 2574
交换基带IQ信号文件格式 目的 通信信号处理时常要记录IQ信号，需要一种开放、方便、高效的交流方法。 方法 用有意义的文件名，相同文件名保存2个文件，json保存采集配置，bin保存I16二进制交织排列的IQ数据。 例如： filename.json filename.bin 代码 json读取 参考本人另一文章《》https://editor.csdn.net/md/?articleId=115328488 ...
第一次通过iq 数据 txt 描述文本 生成 bin 文件
HongDaYu的搬砖之路
01-06 923
#include <iostream> #include <boost/filesystem/operations.hpp> #include <boost/filesystem/path.hpp> #include <boost/filesystem/convenience.hpp> #include <boost/program_options.hpp> #include <vector> #include <string&g
NR PHY -- 学习随笔
weixin_44519974的博客
08-22 465
NR 物理层协议学习随笔记
倍加福ENA58IL-R-ProfiNET手册(译文)
07-11
本手册描述了配备PROFINET接口的Pepperl+Fuchs绝对旋转编码器如何集成到PROFINET网络中。 本手册适用于以下绝对旋转编码器类型： ■■Exx58N-...PN... ■Exx58N-。。。请注意。。。 ■■ENA58IL-...B17 ■ENA58IL-。。。B17
框架完整论文答辩PPT模板 (2)
07-11
【作品名称】：框架完整论文答辩PPT模板 (2) 【适用人群】：适用于希望学习不同技术领域的小白或进阶学习者。可作为毕设项目、课程设计、大作业、工程实训或初期项目立项。
手持式DR设备，全球前8强生产商排名及市场份额调研数据.pdf
07-11
QYResearch是全球知名的大型咨询公司，行业涵盖各高科技行业产业链细分市场，横跨如半导体产业链（半导体设备及零部件、半导体材料、集成电路、制造、封测、分立器件、传感器、光电器件）、光伏产业链（设备、硅料/硅片、电池片、组件、辅料支架、逆变器、电站终端）、新能源汽车产业链（动力电池及材料、电驱电控、汽车半导体/电子、整车、充电桩）、通信产业链（通信系统设备、终端设备、电子元器件、射频前端、光模块、4G/5G/6G、宽带、IoT、数字经济、AI）、先进材料产业链（金属材料、高分子材料、陶瓷材料、纳米材料等）、机械制造产业链（数控机床、工程机械、电气机械、3C自动化、工业机器人、激光、工控、无人机）、食品药品、医疗器械、农业等。 邮箱：market@qyresearch.com
大学生创业计划书(16)三套材料.doc
最新发布
07-11
大学生创业计划书(16)三套材料.doc
大学生创业计划书(4)三份资料.docx
07-11
大学生创业计划书(4)三份资料.docx
tkinter fram
08-12
您好！关于Tkinter框架，我可以为您提供一些基本信息和帮助。Tkinter是Python的标准GUI工具包，用于创建图形用户界面（Graphical User Interface，简称GUI）。它提供了一组用于创建窗口、按钮、标签、文本框等GUI元素的类和方法。 您可以使用Tkinter创建各种交互式应用程序，从简单的工具到复杂的桌面应用程序都可以实现。它具有跨平台的特性，可在Windows、Mac和Linux等操作系统上运行。 要开始使用Tkinter，您需要导入Tkinter模块，并创建一个Tk对象作为主窗口。然后，您可以使用各种Tkinter类创建和配置GUI元素，并使用布局管理器来安排它们在窗口中的位置。 以下是一个简单的Tkinter示例代码，演示了如何创建一个简单的窗口： ```python import tkinter as tk # 创建主窗口 root = tk.Tk() # 创建标签 label = tk.Label(root, text="Hello, Tkinter!") label.pack() # 运行主循环 root.mainloop() ``` 这段代码创建了一个带有“Hello, Tkinter!”文本的标签，并将其放置在主窗口中。然后，通过调用`mainloop()`方法运行程序的主循环，使窗口保持打开状态。 这只是Tkinter的基础知识，您可以进一步探索Tkinter的各种功能和特性，如按钮、文本框、事件处理等。希望这可以帮助到您！如果您有任何进一步的问题，请随时提问。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值