CN112213693B - 相位编码雷达信号码元宽度的解算方法及系统 - Google Patents
相位编码雷达信号码元宽度的解算方法及系统 Download PDFInfo
- Publication number
- CN112213693B CN112213693B CN202010995393.XA CN202010995393A CN112213693B CN 112213693 B CN112213693 B CN 112213693B CN 202010995393 A CN202010995393 A CN 202010995393A CN 112213693 B CN112213693 B CN 112213693B
- Authority
- CN
- China
- Prior art keywords
- code element
- width information
- element width
- pulse
- phase
- 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
- 238000000034 method Methods 0.000 title claims abstract description 40
- 238000005070 sampling Methods 0.000 claims abstract description 27
- 230000035772 mutation Effects 0.000 claims abstract description 16
- 238000012545 processing Methods 0.000 claims description 56
- 238000006243 chemical reaction Methods 0.000 claims description 10
- 230000005855 radiation Effects 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 7
- 238000010606 normalization Methods 0.000 claims description 6
- 238000004422 calculation algorithm Methods 0.000 claims description 5
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 claims description 4
- 101100517651 Caenorhabditis elegans num-1 gene Proteins 0.000 claims description 3
- 238000004590 computer program Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims description 3
- 238000013139 quantization Methods 0.000 claims description 3
- 108091026890 Coding region Proteins 0.000 abstract description 7
- 238000012544 monitoring process Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 4
- 238000012795 verification Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000035485 pulse pressure Effects 0.000 description 1
- 238000012163 sequencing technique Methods 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
- 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/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
- G01S7/418—Theoretical aspects
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明提供一种相位编码雷达信号码元宽度的解算方法及系统，涉及电磁环境信号监测技术领域，所述系统采用该解算方法，包括步骤：对输入的采样样本，提取相位差数据，统计突变点信息，并统计各个突变点之间的距离的最小值，以此作为相位编码信号的码元宽度信息；根据码元特征，验证该始码元宽度信息是否真实有效，如果有效，则作为真实的码元宽度；如果无效则依次以距离最小值的1/2、1/3以及最小值和次小值的最大公约数作为码元宽度，利用以上方法再次验证，直至检测出有效码元宽度为止。本发明能适应相位编码雷达信号的任意随机编码序列，特别是当所有子脉冲相位均与相邻子脉冲相位有相同时，也能正确有效地估算出相位编码信号的码元宽度。
Description
技术领域
本发明涉及电磁环境信号监测技术领域，特别是涉及一种相位编码雷达信号码元宽度的解算方法及系统。
背景技术
相位编码信号是雷达脉冲脉内调制信号的一种，广泛运用于脉压雷达系统，相应地，对相位编码信号的识别与特征参数的提取，就成为电子侦察系统非常重要的一环。雷达信号中的相位编码信号主要包括二相编码信号、四相编码信号，常用码型有m序列码、巴克码、随机序列码、泰勒码等。快速而又准确地获得相位编码信号的类型、码元宽度、码元个数、码元规律等特征参数，从而提供有效的作战支援，是电子侦察的重要内容。
目前对相位编码雷达信号的码元宽度的估算，一般是提取采样样本的突变点信息，并统计各个突变点之间的距离的最小值，以此作为相位编码信号的码元宽度信息进行解算。这就要求在相位编码信号的编码序列中必须至少有一个子脉冲与其前后相邻的两个子脉冲相位都不同，对二相编码信号，表现为与该子脉冲前后相邻的两个子脉冲的相位相同。例如附图1的虚线框内所示，子脉冲10相位为“-”，子脉冲9和11相位均为“+”，才能保证得到准确的估算。实际上，相位编码雷达信号脉冲的脉间编码序列是有随机的，当编码序列不满足以上条件，例如附图2、附图3和附图4中，都没有单独的一个子脉冲与其前后相邻子脉冲都不同相的，运用这种估算方法就不准确了。
发明内容
本发明所要解决的技术问题是：为了克服现有技术中的不足，本发明提供一种相位编码雷达信号码元宽度的解算方法及系统，特别是在相位编码信号的编码序列中不存在子脉冲与前后相邻子脉冲相位均不同的情况时，也能完成对相位编码雷达信号的码元宽度的正确估算。
本发明解决其技术问题所要采用的技术方案是：一种相位编码雷达信号码元宽度的解算方法，先用传统方法估算码元宽度，记为CodeWid0，根据信号原理验证该CodeWid0是否为有效的码宽，如果有效，则作为真实的码元宽度；如果无效则依次以CodeWid0/2、CodeWid0/3以及最大公约数作为码元宽度，利用以上方法进行再次验证，直至检测出有效的码元宽度为止。本发明的算法只解算原码宽度，并不提取相位编码信息的码元规律，因此，对两相、四相或其他编码类型均适用。其步骤如下：
步骤1：对输入的采样样本，提取相位差数据，统计突变点信息，并统计各个突变点之间的距离，以突变点之间距离的最小值作为相位编码信号的初始码元宽度信息，记为IniCodeWid。
步骤2：根据码元特征，验证初始码元宽度信息IniCodeWid是否真实有效，如果是真实有效的，则以初始码元宽度信息IniCodeWid为相位编码信号的真实码元宽度信息，解算真值，否则转入步骤3。
针对传统解算方法出现的问题，用求取最大公约数的方法可以解决，但是该算法的实现相对比较耗时，而一般1/2、1/3出现的概率比较高，所以先用这两个值试探。
步骤3：以突变点之间距离的最小值的1/2作为新的初始码元宽度信息，记为CodeWid0，即CodeWid0=IniCodeWid/2，根据码元特征，验证新的初始码元宽度信息CodeWid0是否真实有效，如果是真实有效的，则以新的初始码元宽度信息CodeWid0为相位编码信号的真实码元宽度信息，解算真值，否则转入步骤4。
步骤4：以突变点之间距离的最小值的1/3作为新的初始码元宽度信息，记为CodeWid0，即CodeWid0=IniCodeWid/3，根据码元特征，验证新的初始码元宽度信息CodeWid0是否真实有效，如果是真实有效，则以此CodeWid0为相位编码信号的真实码元宽度信息，解算真值，否则转入步骤5。
步骤5：计算突变点之间的距离的最小值和次小值的最大公约数，以最大公约数作为相位编码信号的码元宽度信息，记为CodeWid1，并以此CodeWid1作为相位编码信号的真实码元宽度信息，解算真值。
进一步，步骤1具体包括以下步骤：
步骤1.1：对输入采样样本提取相位差数据，并对相位差数据做归一化处理；
步骤1.2：对步骤1.1中进行归一化处理后的相位差数据进行统计，找出采样样本的相位差基值basephasediff；
步骤1.3：以相位差基值basephasediff为基础，找出采样样本的所有突变点的峰值的信息，峰值数量记为peaknum，峰值记为Peak[i].value，峰值点序号记为Peak[i].point，其中，i取值为0~peaknum-1；
步骤1.4：分析步骤1.3获得的所有峰值的信息，剔除伪峰，记录所有有效峰值的信息，峰值个数记为realpeaknum，峰值记为RealPeak[i].value，峰值点序号记为RealPeak[i]. point，其中，i取值为0~realpeaknum-1；
步骤1.5：统计相邻有效峰值之间的距离，并将距离从小到大进行排序，排序后的有效峰值距离记为RealPeak[i].dist，其中，i取值为0~realpeaknum-2；以有效峰值最小距离RealPeak[0].dist作为相位编码信号的初始码元宽度信息，记为IniCodeWid。
具体的，所述码元特征指相邻突变点之间的距离是码元宽度的整数倍。
进一步，根据码元特征，步骤2中验证初始码元宽度信息IniCodeWid是否真实有效的具体方法为：判断各突变点之间的距离是否都是初始码元宽度信息IniCodeWid的整数倍，如果都是初始码元宽度信息IniCodeWid的整数倍，则判断所述初始码元宽度信息IniCodeWid为真实有效，则以初始码元宽度信息IniCodeWid为相位编码信号的真实码元宽度信息，解算真值。
进一步，根据码元特征，步骤3和步骤4中验证新的初始码元宽度信息CodeWid0是否真实有效的具体方法为：判断各突变点之间的距离是否都是新的初始码元宽度信息CodeWid0的整数倍，如果都是新的初始码元宽度信息CodeWid0的整数倍，则判断所述新的初始码元宽度信息CodeWid0为真实有效，则以此CodeWid0为相位编码信号的真实码元宽度信息，解算真值。
至此，可以估算出相位编码雷达信号各种编码序列的码元宽度。
一种相位编码雷达信号码元宽度的解算系统，包括主控制台、DSP处理单元、FPGA处理单元和AD转换单元，其中，
主控制台，用于下发工作控制指令到DSP处理单元和FPGA处理单元，控制DSP处理单元和FPGA处理单元进入相应的工作模式；
AD转换单元，用于实时采集来自前端天线及微波模块的模拟信号，经量化处理转换成数字脉冲信号；
FPGA处理单元，用于采集并处理来自AD转换单元的数字脉冲信号，并进行参数解算，完成对数字脉冲信号的载频、脉冲宽度、幅度、到达时间以及方位的测量和处理，根据主控制台的工作控制指令形成脉冲描述字或脉内采样样本数据，并将脉冲描述字或脉内采样样本数据发送至DSP处理单元；
DSP处理单元，内部存储有实现相位编码雷达信号码元宽度解算的计算机程序，用于根据主控制台下发的工作控制指令对FPGA处理单元送来的脉冲描述字进行分选，提取出各辐射源信号，并解算各辐射源信号的参数；或者根据主控制台下发的工作控制指令接受FPGA处理单元送来的脉内采样样本数据，进行脉内调制特征分析，并估算脉内特征参数。
进一步，所述工作模式包括分选模式和脉内分析模式，其中，分选模式指针对输入的在时域上混叠在一起的多部雷达的脉冲序列，利用各种分选算法，分离出每一部雷达各自对应的脉冲序列并提取出每一部雷达的参数的工作模式；脉内分析模式指接受FPGA处理单元送来的脉内采样样本数据，进行脉内调制特征分析，并估算脉内特征参数。工作流程是先通过分选提取出雷达参数，再对感兴趣的雷达的脉内特征进行分析。
本发明的有益效果是：
（1）能适应相位编码雷达信号的任意编码序列，特别是在编码序列中不存在子脉冲与其前后相邻子脉冲相位均不同的情况时，也能完成对相位编码雷达信号的码元宽度的正确估算；
（2）能防止因不能正确提取相位编码雷达信号脉内的码元宽度参数，而造成后端对目标属性和威胁等级的识别失败及最终导致的指控系统决策错误。
附图说明
下面结合附图和实施例对本发明作进一步说明。
图1是典型的13位巴克码编码序列示意图。
图2是随机编码序列示意图。
图3是随机编码序列示意图。
图4是随机编码序列示意图。
图5是本发明相位编码雷达信号码元宽度的解算方法的流程图。
具体实施方式
本发明充分考虑相位编码雷达信号脉内编码序列产生的随机性和电磁环境信号监测系统的实时性，基于传统的对相位编码雷达信号码元宽度的估算方法进行进一步优化，以适应对各种随机码序列码元宽度的估算。
本发明对相位编码雷达信号码元宽度的估算方法为：先用传统方法估算码元宽度，记为CodeWid0，根据码元特征验证该CodeWid0是否为有效的码宽，如果有效，则作为真实的码元宽度；如果无效则依次以CodeWid0/2、CodeWid0/3以及最大公约数作为码元宽度，利用以上方法进行再次验证，直至检测出有效的码元宽度为止。
本发明中码元特征指相邻突变点之间的距离，即相邻峰值之间的距离，是码元宽度的整数倍。根据这一码元特征，对本发明的相位编码雷达信号的码元宽度的估算方法进行详细说明。
如图5所示，本发明的相位编码雷达信号的码元宽度的估算方法，包括以下步骤：
步骤1：对输入的采样样本，提取相位差数据，统计突变点信息，并统计各个突变点之间的距离，以突变点之间距离的最小值作为相位编码信号的初始码元宽度信息，记为IniCodeWid，具体包括以下步骤：
步骤1.1：对输入采样样本提取相位差数据，并对相位差数据做归一化处理；
步骤1.2：对步骤1.1中进行归一化处理后的相位差数据进行统计，找出采样样本的相位差基值basephasediff。
步骤1.3：以相位差基值basephasediff为基础，找出采样样本的所有突变点的峰值的信息，峰值数量记为peaknum，峰值记为Peak[i].value，峰值点序号记为Peak[i].point，其中，i取值为0~peaknum-1。
步骤1.4：分析步骤1.3获得的所有峰值的信息，剔除伪峰，记录所有有效峰值的信息，峰值个数记为realpeaknum，峰值记为RealPeak[i].value，峰值点序号记为RealPeak[i]. point，其中，i取值为0~realpeaknum-1。
步骤1.5：统计相邻有效峰值之间的距离，并将距离从小到大进行排序，排序后的有效峰值距离记为RealPeak[i].dist，其中，i取值为0~realpeaknum-2；以有效峰值最小距离RealPeak[0].dist作为相位编码信号的初始码元宽度信息，记为IniCodeWid。
步骤2：根据码元特征，各峰值之间的距离应该都是码元宽度信息的的整数倍，以此作为依据，验证初始码元宽度信息IniCodeWid是否真实有效，本实施例中，步骤1中已经剔除伪峰，因此，只需对有效峰值之间的距离进行验证即可，也就是，如果每个有效峰值距离RealPeak[i].dist均是初始码元宽度信息IniCodeWid的整数倍，则判定该初始码元宽度信息IniCodeWid为真实有效，因此，以该初始码元宽度信息IniCodeWid为相位编码信号的真实码元宽度信息，解算真值，否则转入步骤3。
当有效峰值最小距离RealPeak[0].dist不满足码元宽度的码元特征时，则需选择有效峰值最小距离RealPeak[0].dist的1/2或1/3进行判断，进入后序步骤。
步骤3：以有效峰值最小距离RealPeak[0].dist的1/2，也就是初始码元宽度信息IniCodeWid的1/2，作为新的初始码元宽度信息，记为CodeWid0，即CodeWid0=IniCodeWid/2，以IniCodeWid/2作为新的初始码元宽度信息CodeWid0，根据码元特征，各峰值之间的距离应该都是码元宽度信息的的整数倍，验证该始码元宽度信息CodeWid0是否真实有效，如果每个有效峰值距离RealPeak[i].dist均是新的初始码元宽度信息CodeWid0的整数倍，则判定该新的初始码元宽度信息CodeWid0为真实有效，则以此CodeWid0为相位编码信号信号的真实码元宽度信息，解算真值，否则转入步骤4。到该步骤即可解算出图2所示情况的码元宽度真值。
步骤4，以有效峰值最小距离RealPeak[0].dist的1/3，也就是初始码元宽度信息IniCodeWid的1/3，作为新的初始码元宽度信息，记为CodeWid0，即CodeWid0=IniCodeWid/3，以IniCodeWid/3作为新的初始码元宽度信息CodeWid0，根据码元特征，各峰值之间的距离应该都是码元宽度信息的的整数倍，验证该始码元宽度信息CodeWid0是否真实有效，如果每个有效峰值距离RealPeak[i].dist均是新的初始码元宽度信息CodeWid0的整数倍，则判定该新的初始码元宽度信息CodeWid0为真实有效，则以此CodeWid0为相位编码信号信号的真实码元宽度信息，解算真值，否则转入步骤5；到该步骤即可解算出图3所示情况的码元宽度真值。
步骤5，计算突变点峰值之间的距离最小值和次小值的最大公约数，以此作为相位编码信号的码元宽度信息，记为CodeWid1，并以此CodeWid1作为信号的真实码元宽度信息，解算真值。到该步骤即可解算出图4所示情况的码元宽度真值。
至此，可以估算出相位编码雷达信号各种编码序列的码元宽度。
为了实现上述方法在实际问题中的应用，结合上述方法给出了一种相位编码雷达信号码元宽度的解算系统，为了达到对输入脉冲流的实时、全脉冲采集处理，本实施例基于AD+FPGA+DSP的硬件架构，包括主控制台、DSP处理单元、FPGA处理单元和AD转换单元，其中，主控制台，用于下发工作控制指令到DSP处理单元和FPGA处理单元，控制DSP处理单元和FPGA处理单元进入相应的工作模式；AD转换单元实时采集来自前端天线及微波模块的模拟信号，经量化处理转换成数字信号；FPGA处理单元采集并处理来自AD转换单元的数字脉冲信号，进行参数解算，完成对数字脉冲信号的载频、脉冲宽度、幅度、到达时间、方位等各参数的测量、处理，根据主控制台下发的工作控制指令可形成脉冲描述字（PDW）或脉内采样样本数据，并将脉冲描述字或脉内采样样本数据发送至DSP处理单元；DSP处理单元，内部存储有实现相位编码雷达信号码元宽度解算的计算机程序，根据主控制台下发的工作控制指令可对FPGA处理单元送来的PDW进行分选，提取出各辐射源信号，并解算各辐射源信号的详细参数；或者根据主控制台下发的工作控制指令接受FPGA处理单元送来的脉内采样样本数据，进行脉内调制特征分析，并估算脉内特征参数。本发明的估算方法在DSP的脉内分析模块中实施，用于对相位编码雷达信号的脉内特征参数之一——码元宽度进行估算。
作为优选本实施例中主控制台采用普通计算机即可；DSP处理单元采用型号为TMS320C6657的DSP；FPGA处理单元采用Xilinx公司的型号为XCKU060的FPGA；AD转换单元采用ADI公司的AD9689。
本实施例有两种工作模式：分选模式和脉内分析模式，由主控制台根据实时侦收到的雷达信号环境进行决策，并下发工作控制指令引导本模块的实时工作模式。其中，分选模式指对FPGA处理单元送来的PDW进行分选，提取出各辐射源信号，并解算各辐射源信号的详细参数；脉内分析模式指接受FPGA处理单元送来的脉内采样样本数据，进行脉内调制特征分析，并估算脉内特征参数。先通过分选提取出雷达信号，再通过脉内分析对感兴趣的雷达进行进一步的脉内特征分析。
正常工作流程为：主控制台先下发分选命令，进入分选模式，主控制台根据分选结果选择需要进行脉内分析的信号，并下发工作控制指令从分选模式切换至脉内分析模块的脉内分析模式。
具体实现步骤如下：
1、接收并解析工作控制指令。
接收来自主控制台的工作控制指令，控制FPGA处理单元和DSP处理单元进入相应的工作模式。本发明专注于脉内分析模式。如果是脉内分析命令，则转入下一步。
2、用传统方法估算码元宽度。
对输入的采样样本，提取相位差数据，统计相位差峰值信息及其峰值之间的距离信息并排序，有效峰值个数记为距离记为realpeaknum，有效峰值距离记为RealPeak[i].dist（i取值为0~realpeaknum-2），以最小距离RealPeak[0].dist作为相位编码信号的码元宽度信息，记为IniCodeWid。
3、验证IniCodeWid是否为有效的码元宽度信息。
根据码元特征，各峰值之间的距离应该都是码元宽度信息的的整数倍，以此作为依据，验证该始码元宽度信息IniCodeWid是否真实有效，如果是，则以此IniCodeWid为信号的真实码元宽度信息，解算真值，否则转入下一步。
4、以IniCodeWid/2为码元宽度信息，并验证其是否有效。
以IniCodeWid/2作为新的初始码元宽度信息CodeWid0，根据码元特征，各峰值之间的距离应该都是码元宽度信息的的整数倍，验证该始码元宽度信息CodeWid0是否真实有效，如果是，则以此CodeWid0为信号的真实码元宽度信息，解算真值，否则转入下一步。
5、以IniCodeWid/3为码元宽度信息，并验证其是否有效。
以IniCodeWid/3作为新的初始码元宽度信息CodeWid0，根据码元特征，各峰值之间的距离应该都是码元宽度信息的的整数倍，验证该始码元宽度信息CodeWid0是否真实有效，如果是，则以此CodeWid0为信号的真实码元宽度信息，解算真值，否则转入下一步。
6、以突变点之间距离的最大公约数作为码元宽度信息。
计算突变点之间的距离的最大公约数，以此作为相位编码信号的码元宽度信息，记为CodeWid1，并以此CodeWid1作为信号的真实码元宽度信息，解算真值。
经过以上步骤，可以估算出相位编码雷达信号各种编码序列的码元宽度，特别是在随机编码序列中不存在子脉冲与其前后相邻子脉冲相位均不同的情况时，也能完成对相位编码雷达信号的码元宽度的正确估算，进而防止因不能正确提取相位编码雷达信号脉内的码元宽度参数，而造成后端对目标属性和威胁等级的识别失败及最终导致的指控系统决策失误。
以上述依据本发明的理想实施例为启示，通过上述的说明内容，相关的工作人员完全可以在不偏离本发明的范围内，进行多样的变更以及修改。本项发明的技术范围并不局限于说明书上的内容，必须要根据权利要求范围来确定其技术性范围。
Claims (7)
1.一种相位编码雷达信号码元宽度的解算方法，其特征在于：包括以下步骤：
步骤1：对输入的采样样本，提取相位差数据，统计突变点信息，并统计各个突变点之间的距离，以突变点之间距离的最小值作为相位编码信号的初始码元宽度信息，记为IniCodeWid；
步骤2：根据码元特征，验证初始码元宽度信息IniCodeWid是否真实有效，如果是真实有效的，则以初始码元宽度信息IniCodeWid为相位编码信号的真实码元宽度信息，解算真值，否则转入步骤3；
步骤3：以突变点之间距离的最小值的1/2作为新的初始码元宽度信息，记为CodeWid0，根据码元特征，验证新的初始码元宽度信息CodeWid0是否真实有效，如果是真实有效的，则以新的初始码元宽度信息CodeWid0为相位编码信号的真实码元宽度信息，解算真值，否则转入步骤4；
步骤4：以突变点之间距离的最小值的1/3作为新的初始码元宽度信息，记为CodeWid0，根据码元特征，验证新的初始码元宽度信息CodeWid0是否真实有效，如果是真实有效，则以此CodeWid0为相位编码信号的真实码元宽度信息，解算真值，否则转入步骤5；
步骤5：计算突变点之间的距离的最小值和次小值的最大公约数，以最大公约数作为相位编码信号的码元宽度信息，记为CodeWid1，并以此CodeWid1作为相位编码信号的真实码元宽度信息，解算真值。
2.如权利要求1所述的相位编码雷达信号码元宽度的解算方法，其特征在于：所述步骤1具体包括以下步骤：
步骤1.1：对输入采样样本提取相位差数据，并对相位差数据做归一化处理；
步骤1.2：对步骤1.1中进行归一化处理后的相位差数据进行统计，找出采样样本的相位差基值basephasediff；
步骤1.3：以相位差基值basephasediff为基础，找出采样样本的所有突变点的峰值的信息，峰值数量记为peaknum，峰值记为Peak[i].value，峰值点序号记为Peak[i].point，其中，i取值为0~peaknum-1；
步骤1.4：分析步骤1.3获得的所有峰值的信息，剔除伪峰，记录所有有效峰值的信息，峰值个数记为realpeaknum，峰值记为RealPeak[i].value，峰值点序号记为RealPeak[i].point，其中，i取值为0~realpeaknum-1；
步骤1.5：统计相邻有效峰值之间的距离，并将距离从小到大进行排序，排序后的有效峰值距离记为RealPeak[i].dist，其中，i取值为0~realpeaknum-2；以有效峰值最小距离RealPeak[0].dist作为相位编码信号的初始码元宽度信息，记为IniCodeWid。
3.如权利要求1或2所述的相位编码雷达信号码元宽度的解算方法，其特征在于：所述码元特征指相邻突变点之间的距离是码元宽度的整数倍。
4.如权利要求3所述的相位编码雷达信号码元宽度的解算方法，其特征在于：步骤2中验证初始码元宽度信息IniCodeWid是否真实有效具体为：判断各突变点之间的距离是否都是初始码元宽度信息IniCodeWid的整数倍，如果都是初始码元宽度信息IniCodeWid的整数倍，则判断所述初始码元宽度信息IniCodeWid为真实有效，则以初始码元宽度信息IniCodeWid为相位编码信号的真实码元宽度信息，解算真值。
5.如权利要求3所述的相位编码雷达信号码元宽度的解算方法，其特征在于：步骤3和步骤4中验证新的初始码元宽度信息CodeWid0是否真实有效具体为：判断各突变点之间的距离是否都是新的初始码元宽度信息CodeWid0的整数倍，如果都是新的初始码元宽度信息CodeWid0的整数倍，则判断所述新的初始码元宽度信息CodeWid0为真实有效，则以此CodeWid0为相位编码信号的真实码元宽度信息，解算真值。
6.一种相位编码雷达信号码元宽度的解算系统，其特征在于：包括主控制台、DSP处理单元、FPGA处理单元和AD转换单元，其中，
主控制台，用于下发工作控制指令到DSP处理单元和FPGA处理单元，控制DSP处理单元和FPGA处理单元进入相应的工作模式；
AD转换单元，用于实时采集来自前端天线及微波模块的模拟信号，经量化处理转换成数字脉冲信号；
FPGA处理单元，用于采集并处理来自AD转换单元的数字脉冲信号，并进行参数解算，完成对数字脉冲信号的载频、脉冲宽度、幅度、到达时间以及方位的测量和处理，根据主控制台的工作控制指令形成脉冲描述字或脉内采样样本数据，并将脉冲描述字或脉内采样样本数据发送至DSP处理单元；
DSP处理单元，内部存储有实现相位编码雷达信号码元宽度解算的计算机程序，用于根据主控制台下发的工作控制指令对FPGA处理单元送来的脉冲描述字进行分选，提取出各辐射源信号，并解算各辐射源信号的参数；或者根据主控制台下发的工作控制指令接受FPGA处理单元送来的脉内采样样本数据，进行脉内调制特征分析，并估算脉内特征参数。
7.如权利要求6所述的相位编码雷达信号码元宽度的解算系统，其特征在于：所述工作模式包括分选模式和脉内分析模式，其中，分选模式指针对输入的在时域上混叠在一起的多部雷达的脉冲序列，利用各种分选算法，分离出每一部雷达各自对应的脉冲序列并提取出每一部雷达的参数的工作模式；脉内分析模式指接受FPGA处理单元送来的脉内采样样本数据，进行脉内调制特征分析，并估算脉内特征参数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010995393.XA CN112213693B (zh) | 2020-09-21 | 2020-09-21 | 相位编码雷达信号码元宽度的解算方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010995393.XA CN112213693B (zh) | 2020-09-21 | 2020-09-21 | 相位编码雷达信号码元宽度的解算方法及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112213693A CN112213693A (zh) | 2021-01-12 |
| CN112213693B true CN112213693B (zh) | 2024-02-06 |
Family
ID=74049807
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010995393.XA Active CN112213693B (zh) | 2020-09-21 | 2020-09-21 | 相位编码雷达信号码元宽度的解算方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112213693B (zh) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1838650A (zh) * | 2006-04-19 | 2006-09-27 | 武汉虹信通信技术有限责任公司 | 一种实现实时检测实际通信波特率的装置和方法 |
| CN104702470A (zh) * | 2015-03-31 | 2015-06-10 | 南京埃斯顿自动控制技术有限公司 | 基于fpga的波特率在线检测方法 |
| CN108761441A (zh) * | 2018-06-14 | 2018-11-06 | 中国人民解放军战略支援部队信息工程大学 | 一种雷达信号相位编码规律识别方法及装置 |
| CN208207200U (zh) * | 2018-06-04 | 2018-12-07 | 西北工业大学 | 一种基于多通道穿墙雷达数据采集的高分辨目标成像系统 |
| CN109164447A (zh) * | 2018-09-26 | 2019-01-08 | 安徽博微长安电子有限公司 | L波段低空监视雷达的电磁频谱侦察装置 |
| CN109672515A (zh) * | 2018-12-25 | 2019-04-23 | 西安思丹德信息技术有限公司 | 一种低信噪比信号的位同步方法 |
| CN110749878A (zh) * | 2019-09-24 | 2020-02-04 | 北京理工大学 | 一种基于数字信道化的脉内频率编码信号参数测量方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7414943B2 (en) * | 2005-01-06 | 2008-08-19 | Lsi Corporation | Method and apparatus for extracting information from burst cutting area of recording medium |
-
2020
- 2020-09-21 CN CN202010995393.XA patent/CN112213693B/zh active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1838650A (zh) * | 2006-04-19 | 2006-09-27 | 武汉虹信通信技术有限责任公司 | 一种实现实时检测实际通信波特率的装置和方法 |
| CN104702470A (zh) * | 2015-03-31 | 2015-06-10 | 南京埃斯顿自动控制技术有限公司 | 基于fpga的波特率在线检测方法 |
| CN208207200U (zh) * | 2018-06-04 | 2018-12-07 | 西北工业大学 | 一种基于多通道穿墙雷达数据采集的高分辨目标成像系统 |
| CN108761441A (zh) * | 2018-06-14 | 2018-11-06 | 中国人民解放军战略支援部队信息工程大学 | 一种雷达信号相位编码规律识别方法及装置 |
| CN109164447A (zh) * | 2018-09-26 | 2019-01-08 | 安徽博微长安电子有限公司 | L波段低空监视雷达的电磁频谱侦察装置 |
| CN109672515A (zh) * | 2018-12-25 | 2019-04-23 | 西安思丹德信息技术有限公司 | 一种低信噪比信号的位同步方法 |
| CN110749878A (zh) * | 2019-09-24 | 2020-02-04 | 北京理工大学 | 一种基于数字信道化的脉内频率编码信号参数测量方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112213693A (zh) | 2021-01-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108226869B (zh) | 一种重频参差的雷达信号的检测方法 | |
| CN110673109B (zh) | 一种星载大光斑激光雷达全波形数据分解方法 | |
| CN111190146B (zh) | 一种基于视觉图形特征的复杂信号分选方法 | |
| CN108549078B (zh) | 一种雷达脉冲信号跨信道合并及检测方法 | |
| CN110988834B (zh) | 一种基于自适应阈值双门限的脉冲到达时间测量方法 | |
| CN115499279B (zh) | 一种基于相位差分的脉内调制类型识别方法 | |
| CN112130142A (zh) | 一种复杂运动目标微多普勒特征提取方法及系统 | |
| CN108226890B (zh) | 基于时间向统计的机场异物雷达检测方法 | |
| CN112380992A (zh) | 一种加工过程监控数据准确性评估与优化方法及装置 | |
| CN112213693B (zh) | 相位编码雷达信号码元宽度的解算方法及系统 | |
| CN114296032A (zh) | 一种基于fpga的实时复杂信号识别与参数估计方法 | |
| KR101651541B1 (ko) | 슬라이딩 윈도우기반 ｉｆｆ 신호탐지 방법 | |
| CN113986711A (zh) | 一种时间序列数据的峰值检测方法、装置及设备 | |
| CN110632563B (zh) | 一种基于短时傅里叶变换的脉内频率编码信号参数测量方法 | |
| CN115980689A (zh) | 基于点云检测的辐射源信号分选方法、装置、设备及介质 | |
| CN113848391B (zh) | 一种脉冲信号检测和提取方法 | |
| CN113052053B (zh) | 一种坐标时间序列周期项的提取方法 | |
| CN115219991A (zh) | 一种基于希尔伯特变换的二相编码调制信号识别方法 | |
| Hodzic et al. | Spatial analysis of target signatures | |
| CN112763989B (zh) | 一种基于cdif的抖动信号分选方法 | |
| CN117310636B (zh) | 一种固定脉冲重复间隔测量方法、设备和介质 | |
| CN113687325B (zh) | 基于lp和hrrp模型的遮蔽小目标检测方法 | |
| CN118604746B (zh) | 基于分段处理和小波变换的雷达调制识别与参数估计方法 | |
| CN118311534B (zh) | 基于机器视觉检测技术的参差雷达信号分选方法与系统 | |
| Xie et al. | Deinterleaving Method for Radar Signals with Dwell and Switch Pulse Repetition Intervals in a Prior Unknown Scenario |
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 |