CN113672541A - 一种基于gpu的pcm/fm遥测信号非相干解调实现方法 - Google Patents

一种基于gpu的pcm/fm遥测信号非相干解调实现方法 Download PDF

Info

Publication number
CN113672541A
CN113672541A CN202110790954.7A CN202110790954A CN113672541A CN 113672541 A CN113672541 A CN 113672541A CN 202110790954 A CN202110790954 A CN 202110790954A CN 113672541 A CN113672541 A CN 113672541A
Authority
CN
China
Prior art keywords
data
gpu
signal
module
pcm
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202110790954.7A
Other languages
English (en)
Other versions
CN113672541B (zh
Inventor
焦义文
孙宽飞
马宏
杨文革
吴涛
李冬
李超
高泽夫
柴恬怡
李远鹏
滕飞
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Peoples Liberation Army Strategic Support Force Aerospace Engineering University
Original Assignee
Peoples Liberation Army Strategic Support Force Aerospace Engineering University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Peoples Liberation Army Strategic Support Force Aerospace Engineering University filed Critical Peoples Liberation Army Strategic Support Force Aerospace Engineering University
Priority to CN202110790954.7A priority Critical patent/CN113672541B/zh
Publication of CN113672541A publication Critical patent/CN113672541A/zh
Application granted granted Critical
Publication of CN113672541B publication Critical patent/CN113672541B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F13/00—Interconnection of, or transfer of information or other signals between, memories, input/output devices or central processing units
    • G06F13/38—Information transfer, e.g. on bus
    • G06F13/42—Bus transfer protocol, e.g. handshake; Synchronisation
    • G06F13/4204—Bus transfer protocol, e.g. handshake; Synchronisation on a parallel bus
    • G06F13/4221—Bus transfer protocol, e.g. handshake; Synchronisation on a parallel bus being an input/output bus, e.g. ISA bus, EISA bus, PCI bus, SCSI bus
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F7/00—Methods or arrangements for processing data by operating upon the order or content of the data handled
    • G06F7/38—Methods or arrangements for performing computations using exclusively denominational number representation, e.g. using binary, ternary, decimal representation
    • G06F7/48—Methods or arrangements for performing computations using exclusively denominational number representation, e.g. using binary, ternary, decimal representation using non-contact-making devices, e.g. tube, solid state device; using unspecified devices
    • G06F7/57—Arithmetic logic units [ALU], i.e. arrangements or devices for performing two or more of the operations covered by groups G06F7/483 – G06F7/556 or for performing logical operations
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D10/00—Energy efficient computing, e.g. low power processors, power management or thermal management

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Computing Systems (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Digital Transmission Methods That Use Modulated Carrier Waves (AREA)

Abstract

本发明提供基于GPU的PCM/FM遥测信号并行非相干解调实现方法，所述方法基于具有GPU的遥测解调系统实现，所述具有GPU的遥测解调系统包括信号采集模块、高速接口模块、核心运算模块和系统控制模块；所述信号采集模块对信号进行数字采样及调节信号幅度；所述高速接口模块接收所述信号采集模块输出的数据信号，并将所述数据信号发送到系统控制模块；所述系统控制模块将分段后的数据分配给所述核心运算模块；所述核心运算模块包括多个GPU；所述方法包括：并行数字下变频、并行FIR滤波、并行鉴频解调和多普勒频偏估计。本发明的方案，可靠性更高，扩展性更好；与传统基于GPU的PCM/FM遥测信号解调器相比，本发明并行性更好，可实现信号的实时解调。

Description

一种基于GPU的PCM/FM遥测信号非相干解调实现方法
技术领域
本发明涉及通信技术领域，尤其涉及一种基于GPU的PCM/FM遥测信号非相干解调实现方法。
背景技术
飞行器遥测是导弹、火箭、卫星等航天器试验和运行过程中必不可少的重要支持系统，能够实时监测航天器内部工作状态、电气性能、环境参数等重要信息，为航天器性能检测、效能评估及故障分析提供依据。现如今随着航天事业的深入发展，全球各个航天大国对宇宙的探索和对宇宙资源的争夺愈加强烈，这就使得遥测技术突显出更为重要的作用。
脉冲编码调制/调频(PCM/FM)技术具有较强的抗尾焰效应能力、抗噪声性能强、发射机功率高等特点，成为国内外航空航天遥测领域长期采用的一种主流体制。PCM/FM体制在遥测系统中运用多年，调频信号本身具有较好的抗干扰能力，传统解调方式为非相干解调，非相干解调具有结构简单的解调设计电路。
经过几十年的发展，遥测系统已经从早期的功能分立的专用设备逐步发展成为多功能数字遥测系统，具有数字化程度高、多功能一体化等优点，较好地满足了现有任务的需要。面向未来任务日益多样化，试验环境日趋复杂化，对遥测系统的灵活适应性提出了更高要求。传统基于软件无线电设计理念的遥测系统侧重于功能的“软件实现”，尽管近年来软件化程度越来越高，但其灵活适应性仍然受到限制，主要体现在:①系统控制和执行设备仍然相互依赖，硬件资源不可切割和按需分配，性能升级与功能扩展需要对硬件环境进行重新部署，这种单纯通过扩展设备量和增加系统复杂度的方式，不适应未来日益增多的试验任务需求；②面向既定任务和功能需求的定制化研发模式，难以快速响应用户的即时需求，不适应未来复杂多变的试验任务需求。因此，迫切需要寻求对传统软件无线电遥测系统的设计理念和研发模式的转型，进一步提高遥测系统的灵活适应性。
传统的PCM/FM遥测解调多采用基于FPGA的硬件板卡实现，核心解调算法由FPGA完成。常用的基于FPGA的PCM/FM遥测信号解调器由信号处理模块、接口模块和数据采集模块组成，系统的结构框图如图1所示。其中数据采集模块将模拟信号进行数字采样，信号处理模块完成信号的下变频、滤波、鉴频解调，接口模块负责将解调数据输出。
图1所示的结构以基于FPGA的信号处理板卡为核心，外加一些时钟管理和电源管理等模块，该板卡主要完成遥测信号的解调以及数据格式转换，但存在以下不足：
(1)硬件系统结构复杂，嵌入式系统开发周期长、成本高，升级维护难度大；(2)软硬件之间高度的耦合性使得系统专用性强，一旦设计制造完成，功能扩展和修改困难大；(3)系统开发门槛高，完成系统的整体开发和调试需要开发人员克服硬件系统、软件语言、处理算法等多方面困难。
发明内容
为解决上述技术问题，本发明提出了一种基于GPU的PCM/FM遥测信号非相干解调实现方法，所述方法，用于解决现有技术中系统控制和执行设备仍然相互依赖，单纯通过扩展设备量和增加系统复杂度的方式，不适应未来日益增多的试验任务需求；面向既定任务和功能需求的定制化研发模式，难以快速响应用户的即时需求，不适应未来复杂多变的试验任务需求；硬件系统结构复杂，功能难于扩展，软硬件之间耦合度高的技术问题。
根据本发明的第一方面，提供一种基于GPU的PCM/FM遥测信号非相干解调实现方法，
所述具有GPU的遥测解调系统包括信号采集模块、高速接口模块、核心运算模块和系统控制模块；所述信号采集模块包括ADC模块、10MHz频标，用于对信号进行数字采样及调节信号幅度；所述高速接口模块包括高速以太网卡、存储服务器及处理模块，所述高速接口模块用于接收所述信号采集模块输出的数据信号，并将所述数据信号发送到系统控制模块；所述系统控制模块包括CPU及缓存单元，所述缓存单元用于缓存所述数据信号，CPU用于调度任务，对所述缓存单元中的数据分段，并将分段后的数据分配给所述核心运算模块；所述核心运算模块包括多个GPU，所述核心运算模块以并行方式处理所述系统控制模块发送的分段后的数据，GPU之间通过Nvlink总线互联；CPU和GPU之间的数据交互以及CPU对GPU的系统调度通过PCIE总线完成；
所述方法包括以下步骤：
步骤S601：所述信号采集模块包括ADC模块和数据分发服务(DDS，DataDistribution Servic)中间件，在所述数据采集模块中，模拟数据经过幅度调节后由模数转换模块进行采样；采样后的数字信号利用DDS中间件进行发布，然后通过万兆网络发送到所述高速接口模块的万兆网卡；
步骤S602：所述高速接口模块通过16x PCIE总线将采集到的数据输入到DDS中间件的缓存区，并在系统控制模块的调度下开始数据的循环缓存；
步骤S603：在所述系统控制模块的调度下，对DDS中间件缓冲区内的数据分段，GPU按段订阅数据，获取PCM/FM信号；
步骤S604：各GPU收到调度指令并获得数据后，将数据转化为浮点数；
步骤S605：各GPU利用数字控制振荡器(NCO)产生的信号进行并行下变频以降低频率，并根据估计的多普勒频偏误差对该数字控制振荡器的信号作实时消除误差处理；
步骤S606：对并行下变频后的数据进行并行时域滤波运算，以消除高频信号的干扰；
步骤S607：将滤波后的数据进行并行叉积鉴频运算，鉴频结果即为所需的调制信号f(t)；将得到的鉴频结果分为两路，一路用于解调数据的输出，即对滤波后的数据进行累加操作并执行帧同步后作为输出数据；另一路用于估计多普勒频偏误差；
步骤S608：根据鉴频运算得到的鉴频结果估计多普勒频偏误差，并将此误差作为所述估计的多普勒频偏误差，反馈给NCO，以供下次下变频运算时使用；
步骤S609：基于解调数据接口要求对同步后的数据进行格式化输出。
根据本发明的第二方面，提供一种计算机可读存储介质，所述存储介质中存储有多条指令；所述多条指令，用于由处理器加载并执行如前所述的基于GPU的PCM-FM遥测信号非相干解调实现方法。
根据本发明第三方面，提供一种计算机可读存储介质，所述存储介质中存储有多条指令；所述多条指令，用于由处理器加载并执行如前所述的基于GPU的PCM/FM遥测信号并行非相干解调实现方法。
根据本发明的上述方案，本发明对传统软件无线电遥测系统的设计思路进行改进，采用CPU+GPU的异构并行计算平台，以GPU作为核心的数据处理器件，用CPU完成任务的分配和调度，设计了基于GPU的PCM/FM遥测信号解调方法及装置，利用图形处理单元(GPU)的高灵活性和高效并行数据处理能，进一步提高遥测系统的灵活适应性，本发明可以有效降低遥测解调器的开发难度；商用器件和软件运算可有效提高系统对各种不同接口遥测解调系统的适应能力，可根据性能需求灵活配置GPU，克服现有技术不足。具有重构灵活、扩展性好等优势。利用GPU并行运算提高运算效率。可靠性更高，扩展性更好，设计实现更简单，可满足测控系统中PCM/FM遥测信号的实时解调需求。与传统基于FPGA的PCM/FM遥测信号解调器相比，本发明可靠性更高，扩展性更好；与传统基于GPU的PCM/FM遥测信号解调器相比，本发明并行性更好，可实现信号的实时解调。
基于GPU的PCM/FM遥测信号非相干解调方法充分利用GPU强大的并行浮点运算能力和灵活的可重构特性，重新设计系统架构并映射和优化解调算法。由于适应了软件无线电发展的趋势，具有更高的可重构特性，该系统将在我国未来的遥测系统中具有广泛的应用前景。本发明提出和设计的PCM/FM遥测信号解调方法具有如下明显优势：
(1)系统开发周期短。OpenCL和CUDA等基于GPU的开发框架使得GPU上的编程较FPGA开发门槛更低，这些开发框架里包含大量的接口函数和库文件，可以像开发普通的C程序一样灵活调用这些库函数，这一特性大大节省了GPU程序的开发周期。
(2)系统重构简单方便。GPU核函数的运行是在CPU的调度下顺序执行的，而数字信号处理的流程可以被分成多个部分通过不同的核函数执行完成。数字信号处理流程的修改可通过仅修改这些核函数完成，而不必更改底层硬件。这就使得GPU数字信号处理流程的修改和重置与全硬件实现相比更加快速简便。
(3)浮点运算能力强，运算精度高。GPU本身具备超强的双精度和单精度并行浮点运算能力，同时GPU具数十GB的片上存储空间，这些优异的性能为本地信号的高精度实现生成提供了更多可供选择的形式，为高精度数字信号处理提供了可能。
(4)扩展性好，升级维护方便，更加易于增加新的功能。基于CUDA开发的数字信号处理系统基于同一架构的硬件平台而实现了所有软件无线电的算法，数据重用性更好。该系统经过一次ADC便可将采样数据反复用于信道综合、信号同步、频率相位提取、数据解调等软件无线电功能的实现，终端功能的扩展仅需要增加软件模块，而不需要重新设计硬件平台。
上述说明仅是本发明技术方案的概述，为了能够更清楚了解本发明的技术手段，并可依照说明书的内容予以实施，以下以本发明的较佳实施例并配合附图详细说明如后。
附图说明
构成本发明的一部分的附图用来提供对本发明的进一步理解，本发明提供如下附图进行说明。在附图中：
图1为现有技术中基于FPGA的PCM/FM遥测解调系统结构示意图；
图2为本发明一个实施方式的基于GPU的PCM/FM遥测解调系统结构示意图；
图3为本发明一个实施方式的基于GPU的PCM/FM遥测信号非相干解调实现方法流程示意图；
图4为本发明PCM/FM遥测信号解调实现原理示意图；
图5为本发明一个实施方式的基于GPU的PCM/FM遥测信号非相干鉴频解调示意图；
图6为本发明一个实施方式的并行非相干鉴频方法流程示意图；
图7为本发明一个实施方式的帧头EB90的调制信号序列；
图8为本发明一个实施方式的PCM/FM调制信号时域图；
图9为本发明一个实施方式的PCM/FM调制信号频谱图；
图10为本发明一个实施方式的基于GPU的正交下变频后频谱图；
图11为本发明一个实施方式的FIR低通滤波后频谱图；
图12为本发明一个实施方式的基于GPU的PCM/FM遥测信号解调出的信号幅度；
图13为本发明一个实施方式的基于GPU的PCM/FM遥测信号非相干解调输出信号。
具体实施方式
首先结合图2-3说明本发明一个实施方式的基于GPU的PCM/FM遥测信号非相干解调方法流程。所述基于GPU的PCM/FM遥测信号非相干解调方法，基于具有GPU的遥测解调系统实现，如图2所示，所述具有GPU的遥测解调系统包括信号采集模块、高速接口模块、核心运算模块和系统控制模块；所述信号采集模块包括ADC模块、10MHz频标，用于对信号进行数字采样及调节信号幅度；所述高速接口模块包括高速以太网卡、存储服务器及处理模块，所述高速接口模块用于接收所述信号采集模块输出的数据信号，并将所述数据信号发送到系统控制模块；所述系统控制模块包括CPU及缓存单元，所述缓存单元用于缓存所述数据信号，CPU用于调度任务，对所述缓存单元中的数据分段，并将分段后的数据分配给所述核心运算模块；所述核心运算模块包括多个GPU，所述核心运算模块以并行方式处理所述系统控制模块发送的分段后的数据，进行并行数字下变频、并行FIR滤波、并行鉴频、多普勒频偏估计以及位同步，GPU之间通过Nvlink总线互联，CPU和GPU之间的数据交互以及CPU对GPU的系统调度通过PCIE总线完成。
本实施例中，所述具有GPU的遥测解调系统采用CPU+GPU的异构并行计算平台，所述具有GPU的遥测解调系统包括信号采集模块、高速接口模块、核心运算模块和系统控制模块。所述数据采集模块由ADC模块、10MHz频标等部分构成，用于完成模拟中频信号的数字采样，同时辅以信号的幅度调节功能。所述系统控制模块是系统管控中心，由PC主机和缓存组成，用于分发数据、各模块的管理调度以及逻辑控制计算。所述核心运算模块由多块GPU计算显卡组成，是具有GPU的遥测解调系统的核心模块，此模块上进行并行运算，进行并行数字下变频、并行FIR滤波、并行鉴频、多普勒频偏估计以及位同步。
本实施例中，多个GPU之间通过高速Nvlink总线互联，完成中间数据的交互，CPU和GPU之间的数据交互以及CPU对GPU的系统调度通过高速PCIE总线完成。各个模块之间通过高速以太网连接，上述实现方式提高了灵活性和可扩展性。
从调频特征来看待PCM/FM信号，其解调可以采用鉴频方式。但是遥测传输信道复杂，遥测信号在传输过程中受多径衰落的影响，难以提取相干载波，所以相干解调鉴频难于实施。在工程上常采用非相干鉴频解调方法，这种方法结构简单，容易实现，且对因多径衰落造成的幅度和相位方面的失真有较强的抗干扰能力，PCM/FM遥测信号解调原理如图4所示，首先使用AD转换模块将模拟信号转换成数字信号。然后数字下变频模块(Digital DownConverter，简称DDC)，将中频输入信号下变频到零中频信号，低通滤波器滤除高频干扰信号。之后进行内插、重采样以减小解调计算量，加快解调速度。利用I路、Q路信号进行非相干叉积鉴频，将得到的鉴频结果分为两路，一路用于估计多普勒频偏，然后将频率偏差反馈到NCO进行频率修正；另一路用于解调数据的输出。
本实施例中，将数据并行与任务并行相结合，由系统控制模块中的数据预处理模块将缓存单元中的数据分解为多个片段，每个片段的端部均与其相邻的片段端部具有部分数据重叠，用以保持数据的连续性和依赖性。例如，将每一片段的头部与其相邻的片段的尾部作若干字节重叠，将每一片段的尾部与其相邻的片段的头部作若干字节重叠。再对每个片段按流水方式进行任务分解，所述核心运算模块包括并行下变频模块、并行FIR低通滤波模块和并行鉴频模块，由一个核函数各自完成一个任务，多个核函数实现任务的并行处理。本实施例的方式，在遥测解调中，既解决了数据连续性问题，又解决了并行处理问题。
如图3、图5-图6所示，图6为基于CUDA(Compute Unified Device Architecture,统一计算设备架构)架构的并行非相干鉴频流程，数据需要经过下变频、FIR低通滤波和鉴频操作。数据的处理包括：(1)在数据预处理阶段，根据数据特性和计算资源数量将数据分解为M段；(2)对每一个数据段，在数据段内，根据频率生成NCO，执行并行下变频，将中心频率搬移到基带；(3)对下变频后的数据执行并行滤波；(4)对滤波后的数据进行差分鉴频；并利用鉴频结果计算多普勒频偏，并根据多普勒频偏结果对NCO进行多普勒误差补偿；(5)对M个数据段产生的鉴频结果进行同步处理，输出处理结果。
本实施例中，利用数据分段将数据和流程最大限度解耦合，提高系统的并行性，用空间复杂度换取数据处理时间复杂度，实现处理数据的并行化。基于GPU的PCM/FM遥测信号并行非相干解调方法包括：
步骤S601：所述信号采集模块包括ADC模块和数据分发服务(DDS，DataDistribution Servic)中间件，在所述数据采集模块中，模拟数据经过幅度调节后由模数转换模块进行采样；采样后的数字信号利用DDS中间件进行发布，然后通过万兆网络发送到所述高速接口模块的万兆网卡；
步骤S602：所述高速接口模块通过16x PCIE总线将采集到的数据输入到DDS中间件的缓存区，并在系统控制模块的调度下开始数据的循环缓存；
步骤S603：在所述系统控制模块的调度下，对DDS中间件缓冲区内的数据分段，GPU按段订阅数据，获取PCM/FM信号；
本实施例中，获取的PCM/FM信号为：
Figure BDA0003161063250000091
其中，t为时间，SPCM/FM(t)为接收到的PCM/FM信号，Ac为幅度，ωc为载波频率，kf为调制指数，f(τ)为调制信号，θ为初始相位；
步骤S604：各GPU收到调度指令并获得数据后，将数据转化为浮点数；
步骤S605：各GPU利用数字控制振荡器(NCO)产生的信号进行并行下变频以降低频率，并根据估计的多普勒频偏误差对该数字控制振荡器的信号作实时消除误差处理；
本实施例中，进行数字下变频以降低频率，即将接收到的PCM/FM信号分别与NCO产生的正交本振的信号cos(ωct)和sin(ωct)相乘，得到：
Figure BDA0003161063250000101
Figure BDA0003161063250000102
步骤S606：对并行下变频后的数据进行并行时域滤波运算，以消除高频信号的干扰；
本实施例中，滤波后获得同相分量与正交分量，
Figure BDA0003161063250000103
Figure BDA0003161063250000104
其中，I(t)为同相分量，Q(t)为正交分量。
步骤S607：将滤波后的数据进行并行叉积鉴频运算，鉴频结果即为所需的调制信号f(t)；将得到的鉴频结果分为两路，一路用于解调数据的输出，即对滤波后的数据进行累加操作并执行帧同步后作为输出数据；另一路用于估计多普勒频偏误差；
本实施例中，所述鉴频运算是通过I路、Q路信号进行非相干叉积鉴频，
鉴频的目的是计算所需的调制信号f(t)，
Figure BDA0003161063250000105
Figure BDA0003161063250000111
本实施例中，频率鉴别的目的即为找到所需的调制信号f(t)，由于调频信号的振幅是确定的，所以：
I2(t)+Q2(t)＝1 (公式4)
因此，公式4转化为：
f(t)＝I'(t)Q(t)-Q'(t)I(t) (公式5)
采用数字化鉴频，由数学运算可得：
Figure BDA0003161063250000112
其中，I'(n)为同相分量I(n)的导数，I(n)为数字化形式的同相分量，I(n-1)为延迟一位的同相分量，Q'(n)为数字化形式的正交分量Q(n)的导数，Q(n-1)为延迟一位的正交分量。
数字化形式的同相分量I(n)、数字化形式的正相分量Q(n)，两路基带信号的表达式为：
Figure BDA0003161063250000113
其中，T为采样周期；
然后对相位进行差分，则公式5化简为：
Figure BDA0003161063250000121
由公式8可知，当2πkfTf(n)小于预设量级时，“≈”关系成立。2πkf和f(n)的值是固定不变的，因此，在公式8中能够调整的值只有采样周期T。所以要采样频率fs足够大，本实施例中，fs＝56M，能保证非相干鉴频解调的准确性，就能得到调制信号f(t)。
步骤S608：根据鉴频运算得到的鉴频结果估计多普勒频偏误差，并将此误差作为所述估计的多普勒频偏误差，反馈给NCO，以供下次下变频运算时使用；
步骤S609：基于解调数据接口要求对同步后的数据进行格式化输出。
进一步地，又一实施例用以说明本申请的有效性。
本实施例的实验环境为：1)通用计算机：型号为W580-G20；2)CPU：2颗英特尔Xeon(至强)E5-2680处理器，主频2.4GHz；3)操作系统：CentOS Linux release 7.9.2009；4)GPU：NVIDA Tesla K20c；5)软件环境：NVIDIA显卡驱动455.23，CUDA11.1。
首先以随机方式仿真产生调制信号，帧头采用EB90、码型为NRZ-L，生成的数据序列如图7所示；把信号加到高斯白噪声信道中，以这些随机数为种子进行扩展。通过对信号调频，实现数字正交调制的中频信号。生成的PCM/FM调制信号如图8所示，载波频率14MHz、码速率2Mbit/s，采样率56MHz，频谱图如图9所示。
正交下变频成基带信号，下变频后的频谱图以及FIR滤波后的频谱图分别如图10和图11所示，从图中可以看出下变频后信号变成零中频、FIR低通滤波很好的滤除干扰。图12展示的是解调后的信号幅度，进行判决后的解调信号如图13所示，可以看出解调信号的波形和原输入信号的波形一致，因此可知本文设计可以正确的解调出原信号。
对1s调制信号进行解调，利用Matlab软件进行完全串行的解调用时295s；而基于GPU的实验平台进行解调用时367ms，加速比804。从实验结果看，基于GPU可以实现PCM/FM信号的实时解调。
本发明实施例进一步给出一种计算机可读存储介质，所述存储介质中存储有多条指令；所述多条指令，用于由处理器加载并执行如前所述的基于GPU的PCM/FM遥测信号并行非相干解调方法。
本发明实施例进一步给出一种计算机可读存储介质，所述存储介质中存储有多条指令；所述多条指令，用于由处理器加载并执行如前所述的基于GPU的PCM/FM遥测信号非相干解调方法。
需要说明的是，在不冲突的情况下，本发明中的实施例及实施例中的特征可以相互组合。
在本发明所提供的几个实施例中，应该理解到，所揭露的系统，装置和方法，可以通过其它的方式实现。例如，以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如，多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外，在本发明各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用硬件加软件功能单元的形式实现。
上述以软件功能单元的形式实现的集成的单元，可以存储在一个计算机可读取存储介质中。上述软件功能单元存储在一个存储介质中，包括若干指令用以使得一台计算机装置(可以是个人计算机，实体机服务器，或者网络云服务器等，需安装Windows、Linux或者Windows Server操作系统)执行本发明各个实施例所述方法的部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(Read-Only Memory，ROM)、随机存取存储器(RandomAccess Memory，RAM)、磁碟或者光盘等各种可以存储程序代码的介质。
以上所述，仅是本发明的较佳实施例而已，并非对本发明作任何形式上的限制，依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰，均仍属于本发明技术方案的范围内。

Claims (8)

1.一种基于GPU的PCM-FM遥测信号非相干解调实现方法，其特征在于，所述具有GPU的遥测解调系统包括信号采集模块、高速接口模块、核心运算模块和系统控制模块；所述信号采集模块包括ADC模块、10MHz频标，用于对信号进行数字采样及调节信号幅度；所述高速接口模块包括高速以太网卡、存储服务器及处理模块，所述高速接口模块用于接收所述信号采集模块输出的数据信号，并将所述数据信号发送到系统控制模块；所述系统控制模块包括CPU及缓存单元，所述缓存单元用于缓存所述数据信号，CPU用于调度任务，对所述缓存单元中的数据分段，并将分段后的数据分配给所述核心运算模块；所述核心运算模块包括多个GPU，所述核心运算模块以并行方式处理所述系统控制模块发送的分段后的数据，GPU之间通过Nvlink总线互联；CPU和GPU之间的数据交互以及CPU对GPU的系统调度通过PCIE总线完成；
所述方法包括以下步骤：
步骤S601：所述信号采集模块包括ADC模块和数据分发服务中间件，在所述数据采集模块中，模拟数据经过幅度调节后由模数转换模块进行采样；采样后的数字信号利用DDS中间件进行发布，然后通过万兆网络发送到所述高速接口模块的万兆网卡；
步骤S602：所述高速接口模块通过16x PCIE总线将采集到的数据输入到DDS中间件的缓存区，并在系统控制模块的调度下开始数据的循环缓存；
步骤S603：在所述系统控制模块的调度下，对DDS中间件缓冲区内的数据分段，GPU按段订阅数据，获取PCM/FM信号；
步骤S604：各GPU收到调度指令并获得数据后，将数据转化为浮点数；
步骤S605：各GPU利用数字控制振荡器(NCO)产生的信号进行并行下变频以降低频率，并根据估计的多普勒频偏误差对该数字控制振荡器的信号作实时消除误差处理；
步骤S606：对并行下变频后的数据进行并行时域滤波运算，以消除高频信号的干扰；
步骤S607：将滤波后的数据进行并行叉积鉴频运算，鉴频结果即为所需的调制信号f(t)；将得到的鉴频结果分为两路，一路用于解调数据的输出，即对滤波后的数据进行累加操作并执行帧同步后作为输出数据；另一路用于估计多普勒频偏误差；
步骤S608：根据鉴频运算得到的鉴频结果估计多普勒频偏误差，并将此误差作为所述估计的多普勒频偏误差，反馈给NCO，以供下次下变频运算时使用；
步骤S609：基于解调数据接口要求对同步后的数据进行格式化输出。
2.如权利要求1所述的基于GPU的PCM/FM遥测信号非相干解调实现方法，其特征在于，将数据并行与任务并行相结合，由数据预处理模块将缓存单元中的数据分解为多个片段，每个片段的端部均与其相邻的片段端部具有部分数据重叠。
3.如权利要求1所述的基于GPU的PCM/FM遥测信号非相干解调实现方法，其特征在于，步骤S603中，获取的PCM/FM信号为：
Figure FDA0003161063240000021
其中，t为时间，SPCM/FM(t)为接收到的PCM/FM信号，Ac为幅度，ωc为载波频率，kf为调制指数，f(τ)为调制信号，θ为初始相位。
4.如权利要求3所述的基于GPU的PCM/FM遥测信号非相干解调实现方法，其特征在于，步骤S605中，下变频，即将接收到的PCM/FM信号分别与NCO产生的正交本振的信号cos(ωct)和sin(ωct)相乘，得到：
Figure FDA0003161063240000022
Figure FDA0003161063240000023
5.如权利要求4所述的基于GPU的PCM/FM遥测信号非相干解调实现方法，其特征在于，步骤S606中，滤波后获得同相分量与正交分量，
Figure FDA0003161063240000031
Figure FDA0003161063240000032
其中，I(t)为同相分量，Q(t)为正交分量。
6.如权利要求5所述的基于GPU的PCM/FM遥测信号非相干解调实现方法，其特征在于，步骤S607中，
所述鉴频运算是通过I路、Q路信号进行非相干叉积鉴频，
鉴频的目的是计算所需的调制信号f(t)，
Figure FDA0003161063240000033
当2πkfTf(n)小于预设量级时，“≈”关系成立，T为采样周期，f(n)为调制信号。
7.一种计算机可读存储介质，所述存储介质中存储有多条指令；所述多条指令，用于由处理器加载并执行如权利要求1-6中任一项所述的基于GPU的PCM/FM遥测信号非相干解调实现方法。
8.一种计算机可读存储介质，所述存储介质中存储有多条指令；所述多条指令，用于由处理器加载并执行如权利要求1-6中任一项所述的基于GPU的PCM/FM遥测信号并行非相干解调实现方法。
CN202110790954.7A 2021-07-13 2021-07-13 一种基于gpu的pcm/fm遥测信号非相干解调实现方法 Active CN113672541B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110790954.7A CN113672541B (zh) 2021-07-13 2021-07-13 一种基于gpu的pcm/fm遥测信号非相干解调实现方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110790954.7A CN113672541B (zh) 2021-07-13 2021-07-13 一种基于gpu的pcm/fm遥测信号非相干解调实现方法

Publications (2)

Publication Number Publication Date
CN113672541A true CN113672541A (zh) 2021-11-19
CN113672541B CN113672541B (zh) 2022-08-12

Family

ID=78539056

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110790954.7A Active CN113672541B (zh) 2021-07-13 2021-07-13 一种基于gpu的pcm/fm遥测信号非相干解调实现方法

Country Status (1)

Country Link
CN (1) CN113672541B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116015322A (zh) * 2023-03-28 2023-04-25 北京中天星控科技开发有限公司 一种多体制机载遥测发射系统及方法

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6570441B1 (en) * 1999-11-05 2003-05-27 Nec Corporation Incoherent demodulator and method of incoherently demodulating an IF signal
US20100194635A1 (en) * 2009-01-30 2010-08-05 Analog Devices, Inc Method and apparatus for software gps receiver
CN105117170A (zh) * 2015-08-24 2015-12-02 浪潮(北京)电子信息产业有限公司 一种计算机系统架构
CN108226921A (zh) * 2017-12-28 2018-06-29 四川九洲空管科技有限责任公司 基于cpu和gpu架构的二次雷达信息处理系统
US20180364361A1 (en) * 2017-06-15 2018-12-20 Microsoft Technology Licensing, Llc Cloud-offloaded position calculation with on-device acquisition
CN109495410A (zh) * 2018-09-28 2019-03-19 西南电子技术研究所（中国电子科技集团公司第十研究所） 高动态pcm/fm信号载波频率精确估计方法
CN111683111A (zh) * 2020-04-17 2020-09-18 中国人民解放军战略支援部队航天工程大学 一种基于gpu的干涉测量多相信道化基带转换系统
CN111682880A (zh) * 2020-04-17 2020-09-18 中国人民解放军战略支援部队航天工程大学 一种基于gpu的流式架构宽带信号数字下变频系统
CN111935047A (zh) * 2020-06-30 2020-11-13 西安空间无线电技术研究所 一种低信噪比损失的fm信号解调方法及系统

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6570441B1 (en) * 1999-11-05 2003-05-27 Nec Corporation Incoherent demodulator and method of incoherently demodulating an IF signal
US20100194635A1 (en) * 2009-01-30 2010-08-05 Analog Devices, Inc Method and apparatus for software gps receiver
CN105117170A (zh) * 2015-08-24 2015-12-02 浪潮(北京)电子信息产业有限公司 一种计算机系统架构
US20180364361A1 (en) * 2017-06-15 2018-12-20 Microsoft Technology Licensing, Llc Cloud-offloaded position calculation with on-device acquisition
CN108226921A (zh) * 2017-12-28 2018-06-29 四川九洲空管科技有限责任公司 基于cpu和gpu架构的二次雷达信息处理系统
CN109495410A (zh) * 2018-09-28 2019-03-19 西南电子技术研究所（中国电子科技集团公司第十研究所） 高动态pcm/fm信号载波频率精确估计方法
CN111683111A (zh) * 2020-04-17 2020-09-18 中国人民解放军战略支援部队航天工程大学 一种基于gpu的干涉测量多相信道化基带转换系统
CN111682880A (zh) * 2020-04-17 2020-09-18 中国人民解放军战略支援部队航天工程大学 一种基于gpu的流式架构宽带信号数字下变频系统
CN111935047A (zh) * 2020-06-30 2020-11-13 西安空间无线电技术研究所 一种低信噪比损失的fm信号解调方法及系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
刘燕都等: "基于CUDA的数字调相信号并行解调程序设计", 《通信技术》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116015322A (zh) * 2023-03-28 2023-04-25 北京中天星控科技开发有限公司 一种多体制机载遥测发射系统及方法

Also Published As

Publication number Publication date
CN113672541B (zh) 2022-08-12

Similar Documents

Publication Publication Date Title
CN107315168B (zh) 一种软件化雷达信号数据处理系统及方法
CN108226921B (zh) 基于cpu和gpu架构的二次雷达信息处理系统
CN102325058B (zh) 一种变频系统群时延测试方法
CN111239774B (zh) 一种卫星导航数字信号生成方法及服务器
CN113672541B (zh) 一种基于gpu的pcm/fm遥测信号非相干解调实现方法
CN102147473A (zh) 一种同频多系统卫星导航信号生成系统
CN103716084A (zh) 遥感卫星数据接收系统重构方法
CN103389500A (zh) 一种卫星导航系统数字中频信号生成装置
CN212694060U (zh) 一种卫星导航接收板卡、装置和系统
CN108008362A (zh) 一种中频雷达目标回波及干扰产生模块
CN104111465A (zh) 连续波干扰信号的估计器、估计方法、消除器及消除方法
CN114050950B (zh) 一种基于gpu的pcm/fm多符号检测位同步方法
CN101762818B (zh) 基于码存储的gps软件接收机基带信号实时跟踪方法
CN113708916B (zh) 一种基于gpu的pcm/fm遥测信号并行多符号检测方法
CN106526513A (zh) 基于异构双核的磁共振接收机
CN109218238B (zh) 一种基于希尔伯特变换的实信号多普勒频移方法
CN101848038A (zh) 基于软件无线电滤波的无线传输装置及其方法
CN104090993B (zh) 一种甚长基线干涉测量相关处理实现方法
CN115905813B (zh) 一种基于vpx架构的可重构信号处理板的重构方法
CN114755630B (zh) 一种基于soc的调频连续波雷达
CN111130545B (zh) 一种数模混合微系统dac/adc单元回环测试系统
Yin et al. Design and Implementation of Multi-rate Digital Signal Processing Software Based on FPGA
CN102045078B (zh) 一种基于fpga的软件接收机系统及实现方法
CN220289718U (zh) 一种用于检测相位的装置及系统
CN113553291A (zh) 数据采样的系统、方法、装置、设备及存储介质

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant