CN116736237B - 一种非连续谱波形设计方法、装置、设备和介质 - Google Patents
一种非连续谱波形设计方法、装置、设备和介质 Download PDFInfo
- Publication number
- CN116736237B CN116736237B CN202311023172.6A CN202311023172A CN116736237B CN 116736237 B CN116736237 B CN 116736237B CN 202311023172 A CN202311023172 A CN 202311023172A CN 116736237 B CN116736237 B CN 116736237B
- Authority
- CN
- China
- Prior art keywords
- time domain
- optimization problem
- carrier signal
- modulation symbol
- waveform
- 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.)
- Active
Links
- 238000001228 spectrum Methods 0.000 title claims abstract description 84
- 238000000034 method Methods 0.000 title claims abstract description 82
- 238000013461 design Methods 0.000 title claims abstract description 38
- 238000005457 optimization Methods 0.000 claims abstract description 63
- 108010076504 Protein Sorting Signals Proteins 0.000 claims abstract description 53
- 230000008569 process Effects 0.000 claims description 19
- 238000004590 computer program Methods 0.000 claims description 14
- 230000006870 function Effects 0.000 claims description 13
- 239000011159 matrix material Substances 0.000 claims description 12
- 238000005070 sampling Methods 0.000 claims description 9
- 239000000969 carrier Substances 0.000 claims description 8
- 238000012546 transfer Methods 0.000 claims description 7
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000001514 detection method Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000008094 contradictory effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000001556 precipitation Methods 0.000 description 1
- 238000007493 shaping process Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000001308 synthesis method Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/28—Details of pulse systems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical 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
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Radar, Positioning & Navigation (AREA)
- Computational Mathematics (AREA)
- Algebra (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Remote Sensing (AREA)
- Pure & Applied Mathematics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本申请属于雷达波形设计技术领域，涉及一种非连续谱波形设计方法、装置、设备和介质。方法包括：获取非连续谱信号的基本参数，并生成初始多载波信号；根据初始多载波信号的连续表示与离散表示，得到初始多载波信号的时域波形，并转换到频域，得到阻带能量；以最小化阻带能量为目标建立优化问题，并转化为实值优化问题；定义拉格朗日函数，初始化拉格朗日乘子；根据拉格朗日乘子，计算调制符号序列；根据调制符号序列，计算时域信号序列；根据调制符号序列和时域信号序列，判断满足输出条件时，输出时域信号序列对应的时域波形为实值优化问题的解。采用本方法能够在干扰密集条件下实现非连续谱波形设计。
Description
技术领域
本申请涉及雷达波形设计技术领域，特别是涉及一种非连续谱波形设计方法、装置、设备和介质。
背景技术
雷达作为一种非接触式传感器，已经在许多领域中得到广泛应用。在自动驾驶领域，许多关键技术均与雷达息息相关，如盲区检测、变道辅助和防碰撞预警等功能。在空中交通管制中，机场场面监视雷达用来监视跑道以及停机坪上的飞机和车辆，使管制员全面了解和掌握机场上各类目标的分布和活动情况。在气象观测中，雷达可以用来探测大气现象，分辨降水类型、分布、移动或演变。在物联网场景下，雷达可以安装于智能家居中，用于监测物体、人体状态。
随着各种各样电子设备的增多，雷达工作面对的电磁环境愈加复杂且恶劣，不同电子设备占据着不同的频段。当干扰在整个工作频段内密集分布时，频谱资源紧张，雷达难以找到足够大的工作带宽。
非连续谱波形应运而生，通过对发射信号的频谱进行设计，可以避开工作频带内的干扰，利用多个离散、非连续的频带保证工作带宽。
现有的技术中，非连续谱波形设计方法有：基于频谱塑形的波形合成方法以获得给定频谱模板的恒模相位编码信号，基于MM算法解决恒模约束下的频谱模板匹配问题以获得非连续谱波形。
然而，当干扰密集时，频谱资源少，要寻找较大的可用带宽较为困难，现有方法难以进行非连续谱波形设计。
发明内容
基于此，有必要针对上述技术问题，提供一种非连续谱波形设计方法、装置、设备和介质，能够在干扰密集条件下实现非连续谱波形设计。
一种非连续谱波形设计方法，包括：
获取非连续谱信号的基本参数，并根据所述基本参数，生成初始多载波信号；
根据所述初始多载波信号的连续表示与离散表示，得到初始多载波信号的时域波形；将所述初始多载波信号的时域波形转换到频域，并得到阻带能量；
以最小化阻带能量为目标建立优化问题，并将所述优化问题转化为实值优化问题；
定义拉格朗日函数，初始化拉格朗日乘子；根据所述拉格朗日乘子，计算调制符号序列；根据所述调制符号序列，计算时域信号序列；
根据所述调制符号序列和所述时域信号序列，判断满足输出条件时，输出所述时域信号序列对应的时域波形为所述实值优化问题的解。
在一个实施例中，判断不满足输出条件时：
根据所述时域信号序列，更新拉格朗日乘子；根据更新后的拉格朗日乘子，更新所述调制符号序列；根据更新后的调制符号序列，更新所述时域信号序列；
根据更新后的调制符号序列和更新后的时域信号序列进行判断，直至满足输出条件。
在一个实施例中，根据所述初始多载波信号的连续表示与离散表示，得到初始多载波信号的时域波形，包括：
所述初始多载波信号的连续表示为：
；
；
式中，为初始多载波信号的波形的连续表示，/>为第/>个子载波，/>为第/>个子载波上的加权调制波形，/>为虚数符号，/>为初始多载波信号的起始频率，/>为子载波之间的间隔，/>为初始多载波信号的脉冲宽度；/>为子载波上的调制信号，共/>个，/>为宽度为/>幅度为1的矩形脉冲；
所述初始多载波信号的离散表示，即多个子载波信号的离散表示为：
；
；
式中，为初始多载波信号的波形的离散表示，/>为转置，/>为采样间隔，/>为每个脉冲的采样点，/>；
所述初始多载波信号的时域波形为：
；
定义子载波上调制信号所组成的序列以及由多个子载波信号的离散表示/>所构成的多载波信号矩阵分别为：
；
；
则：
；
式中，分别为/>在第/>个采样点的值，/>为调制符号序列到时域信号序列的转移矩阵，/>。
在一个实施例中，以最小化阻带能量为目标建立优化问题，并将所述优化问题转化为实值优化问题，包括：
考虑最小化阻带能量，在恒模约束的情况下，建立优化问题：
；
；
；
式中，为阻带能量；
将所述优化问题转化为实值优化问题：
；
；
；
其中：
；
；
；
；
；
式中，调制符号序列，/>为时域信号序列，/>为转移矩阵/>的实值矩阵，/>为取矩阵/>的实部，/>为取矩阵/>的虚部。
在一个实施例中，根据所述拉格朗日乘子，计算调制符号序列，包括：
；
简化为：
；
其中：
；
；
则解为：
；
式中，为惩罚因子，/>为拉格朗日乘子。
在一个实施例中，根据所述调制符号序列，计算时域信号序列，包括：
；
简化为：
；
；
其中：
；
则解为：
。
在一个实施例中，根据所述时域信号序列，更新拉格朗日乘子，包括：
。
一种非连续谱波形设计装置，包括：
获取模块，用于获取非连续谱信号的基本参数，并根据所述基本参数，生成初始多载波信号；
转换模块，用于根据所述初始多载波信号的连续表示与离散表示，得到初始多载波信号的时域波形；将所述初始多载波信号的时域波形转换到频域，并得到阻带能量；
建模模块，用于以最小化阻带能量为目标建立优化问题，并将所述优化问题转化为实值优化问题；
计算模块，用于定义拉格朗日函数，初始化拉格朗日乘子；根据所述拉格朗日乘子，计算调制符号序列；根据所述调制符号序列，计算时域信号序列；
输出模块，用于根据所述调制符号序列和所述时域信号序列，判断满足输出条件时，输出所述时域信号序列对应的时域波形为所述实值优化问题的解。
一种计算机设备，包括存储器和处理器，所述存储器存储有计算机程序，所述处理器执行所述计算机程序时实现以下步骤：
获取非连续谱信号的基本参数，并根据所述基本参数，生成初始多载波信号；
根据所述初始多载波信号的连续表示与离散表示，得到初始多载波信号的时域波形；将所述初始多载波信号的时域波形转换到频域，并得到阻带能量；
以最小化阻带能量为目标建立优化问题，并将所述优化问题转化为实值优化问题；
定义拉格朗日函数，初始化拉格朗日乘子；根据所述拉格朗日乘子，计算调制符号序列；根据所述调制符号序列，计算时域信号序列；
根据所述调制符号序列和所述时域信号序列，判断满足输出条件时，输出所述时域信号序列对应的时域波形为所述实值优化问题的解。
一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现以下步骤：
获取非连续谱信号的基本参数，并根据所述基本参数，生成初始多载波信号；
根据所述初始多载波信号的连续表示与离散表示，得到初始多载波信号的时域波形；将所述初始多载波信号的时域波形转换到频域，并得到阻带能量；
以最小化阻带能量为目标建立优化问题，并将所述优化问题转化为实值优化问题；
定义拉格朗日函数，初始化拉格朗日乘子；根据所述拉格朗日乘子，计算调制符号序列；根据所述调制符号序列，计算时域信号序列；
根据所述调制符号序列和所述时域信号序列，判断满足输出条件时，输出所述时域信号序列对应的时域波形为所述实值优化问题的解。
上述非连续谱波形设计方法、装置、设备和介质，生成初始多载波信号，定义拉格朗日函数，建立离散谱波形设计的优化模型，并转化为实值优化问题，采用交替向量乘子法（Alternating Direction Method of Multipliers，ADMM），将实值优化问题分为多个子问题，并进行迭代求解。本申请提出了一种频域离散的非连续谱波形也就是离散谱波形（在频谱上可以分为通带和阻带两个部分，在阻带内其频谱幅度较低，而在通带内其频谱幅度较高且接***坦）的设计方法，解决了频域干扰密集场景下雷达探测的带宽需求与紧张的频谱资源之间的矛盾，通过将信号频谱离散化，保证雷达的工作带宽，进一步增强了雷达在频谱密集场景下的工作能力，实现电磁频谱拥堵环境下频谱的高效利用，可以用于雷达平台有多部电子设备协同工作的场景。
附图说明
图1为一个实施例中一种非连续谱波形设计方法的应用场景图；
图2为一个实施例中一种非连续谱波形设计方法的流程示意图；
图3为另一个实施例中一种非连续谱波形设计方法的流程示意图；
图4为一个实施例中干扰密集场景下的雷达信号频谱；
图5为一个实施例中一种非连续谱波形设计装置的结构框图；
图6为一个实施例中计算机设备的内部结构图。
具体实施方式
为了使本申请的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本申请进行进一步详细说明。应当理解，此处描述的具体实施例仅仅用以解释本申请，并不用于限定本申请。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。
需要说明，本申请实施例中所有方向性指示（诸如上、下、左、右、前、后……）仅用于解释在某一特定姿态（如附图所示）下各部件之间的相对位置关系、运动情况等，如果该特定姿态发生改变时，则该方向性指示也相应地随之改变。
另外，在本申请中如涉及“第一”、“第二”等的描述仅用于描述目的，而不能理解为指示或暗示其相对重要性或者隐含指明所指示的技术特征的数量。由此，限定有“第一”“第二”的特征可以明示或者隐含地包括至少一个该特征。在本申请的描述中，“多组”的含义是至少两组，例如两组，三组等，除非另有明确具体的限定。
在本申请中，除非另有明确的规定和限定，术语“连接”、“固定”等应做广义理解，例如，“固定”可以是固定连接，也可以是可拆卸连接，或成一体；可以是机械连接，也可以是电连接，还可以是物理连接或无线通信连接；可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部的连通或两个元件的相互作用关系，除非另有明确的限定。对于本领域的普通技术人员而言，可以根据具体情况理解上述术语在本申请中的具体含义。
另外，本申请各个实施例之间的技术方案可以相互结合，但是必须是以本领域普通技术人员能够实现为基础，当技术方案的结合出现相互矛盾或无法实现时应当认为这种技术方案的结合不存在，也不在本申请要求的保护范围之内。
本申请提供的方法，可以应用于如图1所示的应用环境中。其中，终端102通过网络与服务器104进行通信，终端102可以包括但不限于是各种个人计算机、笔记本电脑、智能手机、平板电脑和便携式可穿戴设备，服务器104可以是各类门户网站、工作***后台对应的服务器等。
本申请提供了一种非连续谱波形设计方法，如图2所示，在一个实施例中，以该方法应用于图1中的终端为例进行说明，包括：
步骤202，获取非连续谱信号的基本参数，并根据基本参数，生成初始多载波信号。
在本步骤中，基本参数包括：信号起止频率、阻带起止频率以及信号脉宽；其中，信号起止频率以及阻带起止频率需要根据实际情况设定。至于如何根据基本参数生成初始多载波信号属于现有技术，在此不再赘述。
步骤204，根据初始多载波信号的连续表示与离散表示，得到初始多载波信号的时域波形；将初始多载波信号的时域波形转换到频域，并得到阻带能量。
具体地：
初始多载波信号的连续表示（也就是秒的发射信号波形）为：
；
；
式中，为初始多载波信号的波形的连续表示，/>为第/>个子载波，，也就是说，共有N个子载波，/>为第/>个子载波上的加权调制波形，/>为虚数符号，/>为初始多载波信号的起始频率，/>为子载波之间的间隔，/>为初始多载波信号的脉冲宽度；/>为子载波上的调制信号，共/>个，/>，C为复数集，/>为宽度为/>幅度为1的矩形脉冲；
初始多载波信号的离散表示（也就是通过采样将模拟信号变为数字信号，对子载波进行离散化表示，即多个子载波信号的离散表示）为：
；
；
式中，为初始多载波信号的波形的离散表示，/>为转置，/>为采样间隔，/>为每个脉冲的采样点，/>；
初始多载波信号的时域波形为：
；
定义子载波上调制信号所组成的序列以及由多个子载波信号的离散表示/>所构成的多载波信号矩阵分别为：
；
；
则：
；
式中，分别为/>在第/>个采样点的值，/>为调制符号序列到时域信号序列的转移矩阵，/>，/>。
采用傅里叶变换，将初始多载波信号的时域波形转换到频域，并得到阻带能量：
设的傅里叶变换为/>，则基带信号（雷达发射的信号）频率中，通带频率集合/>和阻带频率集合/>内的发射能量分别为：/>和/>。
的傅里叶变换为/>，/>的傅里叶变换为/>，则：
；
；
定义为对应于/>的离散频率集，/>为对应于/>的离散频率集。/>的每行对应一个离散频率，定义/>为一个从/>修改的矩阵，将/>中所有行复制到/>，并将对应于/>的行置零。类似地，通过将/>中对应于/>的行置零来创建矩阵/>。则阻带和通带上的能量可以分别表示为/>和/>，其中，H为共轭转置。
步骤206，以最小化阻带能量为目标建立优化问题，并将优化问题转化为实值优化问题。
具体地：
考虑最小化阻带能量，在恒模约束的情况下，建立优化问题：
；
；
；
令：
；
则优化问题可以表示为：
；
；
；
式中，为阻带能量，/>；
将优化问题转化为实值优化问题：
；
；
；
其中：
；
；
；
；
；
式中，为调制符号序列，/>为时域信号序列，/>为转移矩阵/>的实值矩阵，即调制符号与时域信号的转移矩阵，/>为取矩阵/>的实部，/>为取矩阵/>的虚部。
步骤208，定义拉格朗日函数，初始化拉格朗日乘子；根据拉格朗日乘子，计算调制符号序列；根据调制符号序列，计算时域信号序列。
具体地：
定义拉格朗日函数：
；
式中，为拉格朗日乘子，/>为惩罚因子。
根据拉格朗日乘子，计算调制符号序列，即给定，更新/>：
；
忽略与无关的项后，简化为：
；
其中：
；
；
则实值优化问题的解为：
；
式中，为惩罚因子，/>为拉格朗日乘子。
根据调制符号序列，计算时域信号序列，即给定，更新/>：
；
忽略与无关的项后，简化为：
；
；
其中：
；
则实值优化问题的解为：
。
步骤210，根据调制符号序列和时域信号序列，判断满足输出条件时，输出时域信号序列对应的时域波形为所述实值优化问题的解。
在本步骤中，满足输出条件具体为：满足收敛条件或达到预设的最大迭代次数，其中，收敛条件为：，/>为预先设定的公差参数。
上述非连续谱波形设计方法，生成初始多载波信号，定义拉格朗日函数，建立离散谱波形设计的优化模型，并转化为实值优化问题，采用交替向量乘子法（AlternatingDirection Method of Multipliers，ADMM），将实值优化问题分为多个子问题，并进行迭代求解。本申请提出了一种频域离散的非连续谱波形也就是离散谱波形（在频谱上可以分为通带和阻带两个部分，在阻带内其频谱幅度较低，而在通带内其频谱幅度较高且接***坦）的设计方法，解决了频域干扰密集场景下雷达探测的带宽需求与紧张的频谱资源之间的矛盾，通过将信号频谱离散化，保证雷达的工作带宽，进一步增强了雷达在频谱密集场景下的工作能力，实现电磁频谱拥堵环境下频谱的高效利用，可以用于雷达平台有多部电子设备协同工作的场景。
应该理解的是，虽然图2的流程图中的各个步骤按照箭头的指示依次显示，但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明，这些步骤的执行并没有严格的顺序限制，这些步骤可以以其它的顺序执行。而且，图2中的至少一部分步骤可以包括多个子步骤或者多个阶段，这些子步骤或者阶段并不必然是在同一时刻执行完成，而是可以在不同的时刻执行，这些子步骤或者阶段的执行顺序也不必然是依次进行，而是可以与其它步骤或者其它步骤的子步骤或者阶段的至少一部分轮流或者交替地执行。
如图3所示，在一个实施例中，判断不满足输出条件时：根据时域信号序列，更新拉格朗日乘子；根据更新后的拉格朗日乘子，更新调制符号序列；根据更新后的调制符号序列，更新时域信号序列；根据更新后的调制符号序列和更新后的时域信号序列进行判断，直至满足输出条件。
具体地：
根据时域信号序列，更新拉格朗日乘子，即给定，更新：
。
也就是说，其具体流程为：
（1）输入：；
（2）初始化：，/>，/>；
（3）for,其中/>为预设的最大迭代次数；
（4）更新：/>；
（5）更新：/>；
（6）更新：/>；
（7）end当或/>时；
（8）从中获得波形序列/>。
在本实施例中，首先输入非连续谱信号的基本参数：信号起止频率、阻带起止频率以及信号脉宽；接着进行初始化，根据基本参数，生成一个初始的多载波信号；然后建立优化问题，并定义拉格朗日函数；判定是否满足收敛条件，若满足收敛条件，则输出时域信号；若不满足收敛条件，则依次更新子载波上的调制信号序列、时域信号序列、拉格朗日乘子，通过迭代的方法确定，直至收敛。
如图4所示，在干扰密集场景下，采用本方法设计波形，得到雷达信号频谱，其中，实线为雷达信号频谱，虚线为干扰频谱。由图4中可以看出，在干扰较为密集的场景下，波形可以在干扰所占用频谱位置处形成缺口，以避开同平台其他电子设备所占用的频谱。
本申请还提供了一种非连续谱波形设计装置，如图5所示，在一个实施例中，包括：获取模块502、转换模块504、建模模块506、计算模块508和输出模块510，其中：
获取模块502，用于获取非连续谱信号的基本参数，并根据基本参数，生成初始多载波信号；
转换模块504，用于根据初始多载波信号的连续表示与离散表示，得到初始多载波信号的时域波形；将初始多载波信号的时域波形转换到频域，并得到阻带能量；
建模模块506，用于以最小化阻带能量为目标建立优化问题，并将优化问题转化为实值优化问题；
计算模块508，用于定义拉格朗日函数，初始化拉格朗日乘子；根据拉格朗日乘子，计算调制符号序列；根据调制符号序列，计算时域信号序列；
输出模块510，用于根据调制符号序列和时域信号序列，判断满足输出条件时，输出时域信号序列对应的时域波形为所述实值优化问题的解。
关于一种非连续谱波形设计装置的具体限定可以参见上文中对于一种非连续谱波形设计方法的限定，在此不再赘述。上述装置中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算机设备中的处理器中，也可以以软件形式存储于计算机设备中的存储器中，以便于处理器调用执行以上各个模块对应的操作。
在一个实施例中，提供了一种计算机设备，该计算机设备可以是终端，其内部结构图可以如图6所示。该计算机设备包括通过***总线连接的处理器、存储器、网络接口、显示屏和输入装置。其中，该计算机设备的处理器用于提供计算和控制能力。该计算机设备的存储器包括非易失性存储介质、内存储器。该非易失性存储介质存储有操作***和计算机程序。该内存储器为非易失性存储介质中的操作***和计算机程序的运行提供环境。该计算机设备的网络接口用于与外部的终端通过网络连接通信。该计算机程序被处理器执行时以实现一种非连续谱波形设计方法。该计算机设备的显示屏可以是液晶显示屏或者电子墨水显示屏，该计算机设备的输入装置可以是显示屏上覆盖的触摸层，也可以是计算机设备外壳上设置的按键、轨迹球或触控板，还可以是外接的键盘、触控板或鼠标等。
本领域技术人员可以理解，图6中示出的结构，仅仅是与本申请方案相关的部分结构的框图，并不构成对本申请方案所应用于其上的计算机设备的限定，具体的计算机设备可以包括比图中所示更多或更少的部件，或者组合某些部件，或者具有不同的部件布置。
在一个实施例中，提供了一种计算机设备，包括存储器和处理器，该存储器存储有计算机程序，该处理器执行计算机程序时实现上述实施例中方法的步骤。
在一个实施例中，提供了一种计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现上述实施例中方法的步骤。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一非易失性计算机可读取存储介质中，该计算机程序在执行时，可包括如上述各方法的实施例的流程。其中，本申请所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用，均可包括非易失性和/或易失性存储器。非易失性存储器可包括只读存储器（ROM）、可编程ROM（PROM）、电可编程ROM（EPROM）、电可擦除可编程ROM（EEPROM）或闪存。易失性存储器可包括随机存取存储器（RAM）或者外部高速缓冲存储器。作为说明而非局限，RAM以多种形式可得，诸如静态RAM（SRAM）、动态RAM（DRAM）、同步DRAM（SDRAM）、双数据率SDRAM（DDRSDRAM）、增强型SDRAM（ESDRAM）、同步链路（Synchlink） DRAM（SLDRAM）、存储器总线（Rambus）直接RAM（RDRAM）、直接存储器总线动态RAM（DRDRAM）、以及存储器总线动态RAM（RDRAM）等。
以上实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本申请的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对本申请专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本申请构思的前提下，还可以做出若干变形和改进，这些都属于本申请的保护范围。因此，本申请专利的保护范围应以所附权利要求为准。
Claims (9)
1.一种非连续谱波形设计方法，其特征在于，包括：
获取非连续谱信号的基本参数，并根据所述基本参数，生成初始多载波信号；
根据所述初始多载波信号的连续表示与离散表示，得到初始多载波信号的时域波形；将所述初始多载波信号的时域波形转换到频域，并得到阻带能量；
以最小化阻带能量为目标建立优化问题，并将所述优化问题转化为实值优化问题；
定义拉格朗日函数，初始化拉格朗日乘子；根据所述拉格朗日乘子，计算调制符号序列；根据所述调制符号序列，计算时域信号序列；
根据所述调制符号序列和所述时域信号序列，判断满足输出条件时，输出所述时域信号序列对应的时域波形为所述实值优化问题的解；
判断不满足输出条件时：
根据所述时域信号序列，更新拉格朗日乘子；根据更新后的拉格朗日乘子，更新所述调制符号序列；根据更新后的调制符号序列，更新所述时域信号序列；
根据更新后的调制符号序列和更新后的时域信号序列进行判断，直至满足输出条件。
2.根据权利要求1所述的非连续谱波形设计方法，其特征在于，根据所述初始多载波信号的连续表示与离散表示，得到初始多载波信号的时域波形，包括：
所述初始多载波信号的连续表示为：
，
，
式中，为初始多载波信号的波形的连续表示，/>为第/>个子载波，/>为第/>个子载波上的加权调制波形，/>为虚数符号，/>为初始多载波信号的起始频率，/>为子载波之间的间隔，/>为初始多载波信号的脉冲宽度；/>为子载波上的调制信号，共/>个，/>为宽度为/>幅度为1的矩形脉冲；
所述初始多载波信号的离散表示，即多个子载波信号的离散表示为：
，
，
式中，为初始多载波信号的波形的离散表示，/>为转置，/>为采样间隔，/>为每个脉冲的采样点，/>=0,1,…,/>-1；
所述初始多载波信号的时域波形为：
，
定义子载波上的调制信号所组成的序列以及多个子载波信号的离散表示/>所构成的多载波信号矩阵分别为：
，
，
则：
，
式中，分别为/>在第0,1,…,/>-1个采样点的值，/>为调制符号序列到时域信号序列的转移矩阵，/>。
3.根据权利要求2所述的非连续谱波形设计方法，其特征在于，以最小化阻带能量为目标建立优化问题，并将所述优化问题转化为实值优化问题，包括：
考虑最小化阻带能量，在恒模约束的情况下，建立优化问题：
，
，
，
式中，为阻带能量；
将所述优化问题转化为实值优化问题：
，
，
，
其中：
，
，
，
，
，
式中，为调制符号序列，/>为时域信号序列，/>为转移矩阵/>的实值矩阵，/>为取矩阵/>的实部，/>为取矩阵/>的虚部。
4.根据权利要求3所述的非连续谱波形设计方法，其特征在于，根据所述拉格朗日乘子，计算调制符号序列，包括：
，
简化为：
，
其中：
，
，
则解为：
，
式中，为惩罚因子，/>为拉格朗日乘子。
5.根据权利要求4所述的非连续谱波形设计方法，其特征在于，根据所述调制符号序列，计算时域信号序列，包括：
，
简化为：
，
，
其中：
，
则解为：
。
6.根据权利要求5所述的非连续谱波形设计方法，其特征在于，根据所述时域信号序列，更新拉格朗日乘子，包括：
。
7.一种非连续谱波形设计装置，其特征在于，包括：
获取模块，用于获取非连续谱信号的基本参数，并根据所述基本参数，生成初始多载波信号；
转换模块，用于根据所述初始多载波信号的连续表示与离散表示，得到初始多载波信号的时域波形；将所述初始多载波信号的时域波形转换到频域，并得到阻带能量；
建模模块，用于以最小化阻带能量为目标建立优化问题，并将所述优化问题转化为实值优化问题；
计算模块，用于定义拉格朗日函数，初始化拉格朗日乘子；根据所述拉格朗日乘子，计算调制符号序列；根据所述调制符号序列，计算时域信号序列；
输出模块，用于根据所述调制符号序列和所述时域信号序列，判断满足输出条件时，输出所述时域信号序列对应的时域波形为所述实值优化问题的解；判断不满足输出条件时：根据所述时域信号序列，更新拉格朗日乘子；根据更新后的拉格朗日乘子，更新所述调制符号序列；根据更新后的调制符号序列，更新所述时域信号序列；根据更新后的调制符号序列和更新后的时域信号序列进行判断，直至满足输出条件。
8.一种计算机设备，包括存储器和处理器，所述存储器存储有计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至6中任一项所述方法的步骤。
9.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现权利要求1至6中任一项所述的方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311023172.6A CN116736237B (zh) | 2023-08-15 | 2023-08-15 | 一种非连续谱波形设计方法、装置、设备和介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311023172.6A CN116736237B (zh) | 2023-08-15 | 2023-08-15 | 一种非连续谱波形设计方法、装置、设备和介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116736237A CN116736237A (zh) | 2023-09-12 |
| CN116736237B true CN116736237B (zh) | 2023-10-31 |
Family
ID=87911905
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311023172.6A Active CN116736237B (zh) | 2023-08-15 | 2023-08-15 | 一种非连续谱波形设计方法、装置、设备和介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116736237B (zh) |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7376688B1 (en) * | 2001-01-09 | 2008-05-20 | Urbain A. von der Embse | Wavelet multi-resolution waveforms |
| US8908992B1 (en) * | 2011-04-12 | 2014-12-09 | Stc.Unm | System and methods of regularized optimization for matrix factorization and image and video reconstruction |
| CN104898113A (zh) * | 2015-06-19 | 2015-09-09 | 哈尔滨工业大学 | 一种多输入多输出雷达波形设计方法 |
| CN105137422A (zh) * | 2015-09-09 | 2015-12-09 | 哈尔滨工业大学 | 一种非连续谱的连续调相信号设计方法 |
| CN106961406A (zh) * | 2016-01-11 | 2017-07-18 | 中兴通讯股份有限公司 | 多载波***的数据调制、解调方法、帧生成方法及节点 |
| CN113489573A (zh) * | 2021-09-07 | 2021-10-08 | 中国人民解放军国防科技大学 | 双基地雷达通信一体化波形设计方法、计算机设备及介质 |
| CN114884787A (zh) * | 2022-04-22 | 2022-08-09 | 华中科技大学 | 适用快时变信道的波形可控多载波通信方法、装置及*** |
| CN116389212A (zh) * | 2023-04-07 | 2023-07-04 | 东南大学 | 基于admm的低复杂度多载波调制信号峰均比抑制方法 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2013101342A1 (en) * | 2011-12-29 | 2013-07-04 | Raytheon Bbn Technologies Corp. | Non-contiguous spectral-band modulator and method for non-contiguous spectral-band modulation |
| EP4026286A4 (en) * | 2019-09-05 | 2022-11-09 | ZTE Corporation | LOW PEAK POWER RATIO (PAPR) MODULATION SCHEME |
| US20210132211A1 (en) * | 2019-10-30 | 2021-05-06 | University Of Kansas | Physical waveform optimization for multiple-beam multifunction digital arrays |
-
2023
- 2023-08-15 CN CN202311023172.6A patent/CN116736237B/zh active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7376688B1 (en) * | 2001-01-09 | 2008-05-20 | Urbain A. von der Embse | Wavelet multi-resolution waveforms |
| US8908992B1 (en) * | 2011-04-12 | 2014-12-09 | Stc.Unm | System and methods of regularized optimization for matrix factorization and image and video reconstruction |
| CN104898113A (zh) * | 2015-06-19 | 2015-09-09 | 哈尔滨工业大学 | 一种多输入多输出雷达波形设计方法 |
| CN105137422A (zh) * | 2015-09-09 | 2015-12-09 | 哈尔滨工业大学 | 一种非连续谱的连续调相信号设计方法 |
| CN106961406A (zh) * | 2016-01-11 | 2017-07-18 | 中兴通讯股份有限公司 | 多载波***的数据调制、解调方法、帧生成方法及节点 |
| CN113489573A (zh) * | 2021-09-07 | 2021-10-08 | 中国人民解放军国防科技大学 | 双基地雷达通信一体化波形设计方法、计算机设备及介质 |
| CN114884787A (zh) * | 2022-04-22 | 2022-08-09 | 华中科技大学 | 适用快时变信道的波形可控多载波通信方法、装置及*** |
| CN116389212A (zh) * | 2023-04-07 | 2023-07-04 | 东南大学 | 基于admm的低复杂度多载波调制信号峰均比抑制方法 |
Non-Patent Citations (3)
| Title |
|---|
| MIMO天波雷达非连续谱正交相位编码波形设计;薛永华 等;信号处理;第31卷(第10期);第1255-1262页 * |
| Spectrally Compatible MIMO Radar Beampattern Design Under Constant Modulus Constraints;Khaled Alhujaili 等;IEEE Transactions on Aerospace and Electronic Systems;第56卷(第6期);第4749-4766页 * |
| 非连续谱雷达信号设计综述;位寅生 等;雷达学报;第11卷(第2期);第183-197页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116736237A (zh) | 2023-09-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104076351B (zh) | 一种用于高速高机动目标的相参积累检测方法 | |
| CN106990399B (zh) | 用于目标跟踪的组网雷达***功率和带宽联合分配方法 | |
| Morzfeld et al. | What the collapse of the ensemble Kalman filter tells us about particle filters | |
| Qiao et al. | Adaptive lasso echo state network based on modified Bayesian information criterion for nonlinear system modeling | |
| CN116047427B (zh) | 一种小样本雷达有源干扰识别方法 | |
| Yang et al. | Kalman filters for linear continuous‐time fractional‐order systems involving coloured noises using fractional‐order average derivative | |
| CN116736237B (zh) | 一种非连续谱波形设计方法、装置、设备和介质 | |
| Ray et al. | Scaling and nonlinear behaviour of daily mean temperature time series across India | |
| CN115830709A (zh) | 基于联邦边缘学习的动作识别方法、服务器及电子设备 | |
| CN103679715A (zh) | 一种基于非负矩阵分解的手机图像特征提取方法 | |
| Wu | Communication modulation recognition algorithm based on STFT mechanism in combination with unsupervised feature-learning network | |
| CN106842204A (zh) | 一种高分辨率聚束sar二维缺失数据成像处理方法 | |
| CN112034447B (zh) | 一种提高脉冲多普勒雷达检测性能的二相调制方法及装置 | |
| CN105425234A (zh) | 基于多任务贝叶斯压缩感知的距离-多普勒成像方法 | |
| Jain et al. | A Cloud‐Based Machine Learning Approach to Reduce Noise in ECG Arrhythmias for Smart Healthcare Services | |
| Xueyun et al. | Application of kernel PCA for foetal ECG estimation | |
| Yang et al. | Federated learning over wireless backhaul for distributed micro‐Doppler radars: Deep learning aided gradient estimation | |
| CN116381629A (zh) | 基于实时动态卷积的雷达大型目标模拟***以及方法 | |
| CN109886126A (zh) | 一种基于动态采样机制与rbf神经网络的区域车辆密度估计方法 | |
| CN114418068A (zh) | 基于多步神经网络结构的毫米波雷达稀疏点云重构方法 | |
| Li et al. | A chaotification model based on sine and cosecant functions for enhancing chaos | |
| Liu et al. | Enhanced algorithms for deploying the minimum sensors to construct a wireless sensor network having full coverage of critical square grids | |
| CN109031289B (zh) | 认知捷变频雷达波形设计方法和装置 | |
| CN112583748B (zh) | 一种信道估计方法、装置及电子设备 | |
| CN110133647A (zh) | 一种分块稀疏贝叶斯学习isar成像散射系数估计方法 |
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 |