CN114859296B - 具有空间扫描特性的高频地波雷达信号处理方法 - Google Patents

具有空间扫描特性的高频地波雷达信号处理方法 Download PDF

Info

Publication number
CN114859296B
CN114859296B CN202210444733.9A CN202210444733A CN114859296B CN 114859296 B CN114859296 B CN 114859296B CN 202210444733 A CN202210444733 A CN 202210444733A CN 114859296 B CN114859296 B CN 114859296B
Authority
CN
China
Prior art keywords
transmitting
signal
azimuth
distance
signals
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
Application number
CN202210444733.9A
Other languages
English (en)
Other versions
CN114859296A (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.)
Harbin Institute of Technology
Original Assignee
Harbin Institute of Technology
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 Harbin Institute of Technology filed Critical Harbin Institute of Technology
Priority to CN202210444733.9A priority Critical patent/CN114859296B/zh
Publication of CN114859296A publication Critical patent/CN114859296A/zh
Application granted granted Critical
Publication of CN114859296B publication Critical patent/CN114859296B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • 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
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Abstract

具有空间扫描特性的高频地波雷达信号处理方法，属于高频地波雷达信号技术领域。本发明针对现有高频地波雷达信号不具备发射波束的空间扫描功能的问题。主要包括对雷达信号的设计方法：设定包含Q个方位个数的发射波束方位集合，生成Q路发射方位信号，每路发射方位信号中发射方位信号的个数与发射天线阵列中发射单元个数相同；设定每个信号调制周期包括P个脉冲的待发射时间编码信号；使P＝Q，按序号将每路发射方位信号中的每个发射方位信号与对应脉冲序号的待发射时间编码信号相乘获得Q路具有方位信息的空间‑时间‑编码信号，每路空间‑时间‑编码信号中包括M个空间‑时间‑编码信号。本发明设计的雷达信号具备波束空间扫描能力。

Description

具有空间扫描特性的高频地波雷达信号处理方法
技术领域
本发明涉及具有空间扫描特性的高频地波雷达信号处理方法，属于高频地波雷达信号技术领域。
背景技术
高频地波雷达工作于短波频段(3MHz-30MHz)的频率低端，是一种可海洋遥感和观测海洋观测对象的雷达系统。由于高频信号波长较长，经典的高频地波雷达在发射电磁波时，一般不进行空间扫描，而是采用以对数周期天线为主的天线，直接向大范围空间辐射，此时发射形成的观测波束宽度较宽，不具备发射波束的空间扫描功能。
由于发射波束无空间扫描功能，经典高频地波雷达的雷达信号设计及其处理自然而然地与发射空间方位无关。在设计时，仅考虑信号在时间-频率域完成即可。进入二十一世纪后，单站高频地波雷达面临着强电离层杂波问题的困扰。为了抑制不同方向电离层杂波影响，高频地波雷达的信号设计方法从时间上的信号正交特性入手，发展了正交编码信号体制的高频地波雷达信号，从而使雷达信号设计方法及其处理发展为时间-频率-相位联合域的同时多维信号优化问题。
随着对短波频段电离层杂波、海杂波空时特性研究的不断深入，挖掘各类信号的空时特性，并基于此发展新的信号设计技术，逐渐成为新型高频地波雷达的重要发展方向。相应地，空时杂波抑制方法、空时联合信号处理方法层出不穷，也对高频地波雷达提出了空间信息最优化利用的需求。相比大范围空间辐射模式，在实际的大多数应用中，越来越需要对特定的海域，或观测对象进行持续关注，使观测模式向具有方向特性转变。这些转变和需求都迫切使得高频地波雷达具有发射波束的空间扫描功能，扩展现有信号设计的维度，从时间-频率-相位联合域发展成为空间-时间-频率-相位联合域。因此，空间信息与雷达时间信号结合的有效性问题就成为了新型高频地波雷达的核心问题，急需发展新的空间-时间雷达信号形式，以在实现空间扫描的同时，有效提取观测对象信息。
发明内容
针对现有方法设计的高频地波雷达信号不具备发射波束的空间扫描功能的问题，本发明提供一种具有空间扫描特性的高频地波雷达信号处理方法。
本发明的一种具有空间扫描特性的高频地波雷达信号处理方法，包括具有空间扫描特性的高频地波雷达信号的设计方法，具体包括：
步骤一：设定包含Q个方位个数的发射波束方位集合，生成Q路发射方位信号，每路发射方位信号中发射方位信号的个数与发射天线阵列中发射单元个数相同；
步骤二：设定每个信号调制周期包括P个脉冲的待发射时间编码信号；其中Q、M和P均为正整数；
步骤三：使P＝Q，按序号将每路发射方位信号中的每个发射方位信号与对应脉冲序号的待发射时间编码信号相乘，获得Q路具有方位信息的空间-时间-编码信号，Q路具有方位信息的空间-时间-编码信号发射至观测空间进行观测；每路空间-时间-编码信号中包括M个空间-时间-编码信号。
根据本发明的具有空间扫描特性的高频地波雷达信号处理方法，步骤一中，生成的Q路发射方位信号中第q路发射方位信号表示为：
Figure BDA0003616188010000021
式中βm为发射方位信号，θq为发射波束方位，m为发射天线阵列中发射单元的序号，d为发射单元间距，c为光速，f0是发射电磁波频率；q为发射方位序号，0≤q≤Q-1；
待发射时间编码信号表示为：
Figure BDA0003616188010000022
式中s(t)为待发射时间编码信号，t为时间，p为信号调制周期中脉冲序号，I为每个脉冲宽度内的编码个数，i为脉冲宽度内编码序号，且0≤i≤I；
Figure BDA0003616188010000023
为编码，rect(·)为门函数：
Figure BDA0003616188010000024
tb为每个编码的时间长度，T为脉冲重复周期。
根据本发明的具有空间扫描特性的高频地波雷达信号处理方法，步骤二中，每路空间-时间-编码信号中，对应第m个发射单元的空间-时间-编码信号为：
Figure BDA0003616188010000025
式中sm(t)为第m个发射单元的空间-时间-编码信号，sp(t)为待发射时间编码信号s(t)的第p个脉冲，p＝q。
根据本发明的具有空间扫描特性的高频地波雷达信号处理方法，还包括对具有方位信息的空间-时间-编码信号的回波信号进行处理的方法，具体包括：
步骤四：Q路具有方位信息的空间-时间-编码信号被观测对象反射后，获得Q个方位的观测对象回波信号；
步骤五：对Q个方位的观测对象回波信号进行脉冲压缩处理，得到发射方位-接收通道-距离频谱矩阵；所述接收通道的个数N为接收天线阵列中接收单元的个数；N为正整数；
步骤六：将设定个数信号调制周期的发射方位-接收通道-距离频谱矩阵按照信号调制周期序号排列，获得排序后发射方位-接收通道-周期-距离频谱矩阵；
步骤七：对所有发射方位中每个发射方位的接收通道-周期-距离频谱使用与发射方位相同指向的波束加权系数进行数字波束形成处理，得到波束-周期-距离频谱；
步骤八：对波束-周期-距离频谱中的每一列离散化距离单元信号，进行多普勒分析处理，得到具有空间扫描特性的波束-距离-多普勒频谱矩阵；
步骤九：对发射波束方位集合中的每一个方位，重复步骤二至步骤八，完成具有凝视效果的观测空间扫描。
根据本发明的具有空间扫描特性的高频地波雷达信号处理方法，步骤四中，观测对象回波信号表示为：
Figure BDA0003616188010000031
式中
Figure BDA0003616188010000032
为第q个发射波束方位下第n个接收通道的观测对象回波信号，n＝0，1，2，……，N-1；τ为观测对象时间延迟，l为信号调制周期计数；θp为观测对象回波信号方位；
Figure BDA0003616188010000033
R为观测对象距离发射天线阵列中心的距离，R＞＞N·d；
根据θp＝θq，观测对象回波信号的表达式中，
Figure BDA0003616188010000034
的计算结果为常数Aq。
根据本发明的具有空间扫描特性的高频地波雷达信号处理方法，步骤五中获得发射方位-接收通道-距离频谱的方法包括：
第q个发射波束方位下的观测对象回波信号
Figure BDA0003616188010000035
进行脉冲压缩处理，得到单个信号调制周期当前方位的接收通道-距离频谱
Figure BDA0003616188010000036
Figure BDA0003616188010000041
式中r为观测对象相对高频地波雷达的离散化距离序号；
由Q个发射方位的观测对象回波信号获得单个信号调制周期的Q个接收通道-距离频谱，形成发射方位-接收通道-距离频谱矩阵，所述矩阵维度为Q×N×rmax，其中rmax为观测对象相对高频地波雷达的最大离散化距离序号。
根据本发明的具有空间扫描特性的高频地波雷达信号处理方法，步骤六中，排序后发射方位-接收通道-周期-距离频谱矩阵W为：
Figure BDA0003616188010000042
式中
Figure BDA0003616188010000043
为第q个发射方位第n个接收通道的设定个数信号调制周期的距离频谱：
Figure BDA0003616188010000044
式中x为离散化距离单元信号，lmax为最大信号调制周期序号，
Figure BDA0003616188010000045
中x的行号l为信号调制周期序号，列号为离散化距离r的序号。
根据本发明的具有空间扫描特性的高频地波雷达信号处理方法，步骤七中，得到波束-周期-距离频谱的方法包括：
Figure BDA0003616188010000046
式中Bq为波束-周期-距离频谱，
Figure BDA0003616188010000047
为与θq相关的第n个接收通道波束加权系数，所有
Figure BDA0003616188010000048
组成波束加权系数集合w(θq)：
Figure BDA0003616188010000049
式中an来自长度为N的窗函数。
根据本发明的具有空间扫描特性的高频地波雷达信号处理方法，所述窗函数包括汉宁窗函数和切比雪夫窗函数；
当窗函数为汉宁窗函数，则：
Figure BDA0003616188010000051
根据本发明的具有空间扫描特性的高频地波雷达信号处理方法，步骤八中，对波束-周期-距离频谱Bq中的每一列离散化距离单元信号进行多普勒分析处理，得到第q个发射波束方位的距离-多普勒频谱矩阵Dq：
Figure BDA0003616188010000052
式中y为多普勒信号，y的行号dmax为多普勒单元序号，y的列号为离散化距离r的序号；取dmax＝lmax，则：
第r列的多普勒信号yd(r)为：
Figure BDA0003616188010000053
式中wl来自给定lmax的窗函数。
本发明的有益效果：为实现高频地波雷达空间波束扫描，本发明方法在传统时间-编码信号设计方法基础上，基于高频地波雷达的空间波束扫描需求，设计空间-时间域复合的编码信号，并进一步实现了信号处理的全过程，可输出具有空间凝视特性的编码信号的波束-距离-多普勒频谱，完成了扫描方位内回波信息的无失真提取。
本发明能够使得编码形式的高频地波雷达观测信号具备空间扫描特性，可使得高频地波雷达观测波束在空间中自由扫描，既可以实施有目的地长时关注特殊区域或观测对象，还可以跟随观测对象运动持续改变。同时由于信号具备空间特性，因此可以增强雷达系统的抗空间干扰能力，有利于观测对象检测和航迹跟踪，具有实施简单方便、可自适应的改变辐射信号的空间特性等特点。
本发明可用于具有空间扫描功能或空间凝视特性的高频地波雷达或高频雷达，提升其波束空间扫描能力。
附图说明
图1是本发明所述具有空间扫描特性的高频地波雷达信号的设计方法的流程图；
图2是本发明所述具有空间扫描特性的高频地波雷达信号处理方法的整体流程图；
图3是本发明具体实施例步骤五中，距离频谱示意图；
图4是本发明具体实施例步骤九中，距离-多普勒频谱示意图；
图5是本发明具体实施例步骤九中，波束方位为-30°时的波束-距离-多普勒频谱结果示意图；
图6是本发明具体实施例步骤九中，波束方位为-10°时的波束-距离-多普勒频谱结果示意图；
图7是本发明具体实施例步骤九中，波束方位为0°时的波束-距离-多普勒频谱结果示意图；
图8是本发明具体实施例步骤九中，波束方位为10°时的波束-距离-多普勒频谱结果示意图；
图9是本发明具体实施例步骤九中，波束方位为30°时的波束-距离-多普勒频谱结果示意图；
图10是本发明具体实施例步骤九中，多普勒分析处理采用的汉宁窗函数时间域波形图；
图11是本发明具体实施例步骤九中，观测对象能量随波束方向变化曲线示意图；
图12是本发明具体实施例步骤十中，当发射方位变化集合中均为固定方向时，观测对象信号能量在波束-距离-多普勒矩阵集合中的响应示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动的前提下所获得的所有其它实施例，都属于本发明保护的范围。
需要说明的是，在不冲突的情况下，本发明中的实施例及实施例中的特征可以相互组合。
下面结合附图和具体实施例对本发明作进一步说明，但不作为本发明的限定。
具体实施方式一、结合图1所示，本发明提供了一种具有空间扫描特性的高频地波雷达信号处理方法，包括具有空间扫描特性的高频地波雷达信号的设计方法，具体包括：
步骤一：设定包含Q个方位个数的发射波束方位集合，生成Q路发射方位信号，每路发射方位信号中发射方位信号的个数与发射天线阵列中发射单元个数相同；
步骤二：设定每个信号调制周期包括P个脉冲的待发射时间编码信号；其中Q、M和P均为正整数；
步骤三：使P＝Q，按序号将每路发射方位信号中的每个发射方位信号与对应脉冲序号的待发射时间编码信号相乘，获得Q路具有方位信息的空间-时间-编码信号，Q路具有方位信息的空间-时间-编码信号发射至观测空间进行观测；每路空间-时间-编码信号中包括M个空间-时间-编码信号。
进一步，步骤一中，生成的Q路发射方位信号中第q路发射方位信号表示为：
Figure BDA0003616188010000071
式中βm为发射方位信号，θq为发射波束方位，{θq}为事先指定的发射波束方位集合，m为发射天线阵列中发射单元的序号，d为发射单元间距，c为光速，为3×108；f0是发射电磁波频率；q为发射方位序号，0≤q≤Q-1，Q可以为1或任意整数个，Q路发射方位信号与发射天线阵列中发射单元一一对应；
其中发射天线阵列采用均匀直线形式。
待发射时间编码信号表示为：
Figure BDA0003616188010000072
式中s(t)为待发射时间编码信号，t为时间，p为信号调制周期中脉冲序号，I为每个脉冲宽度内的编码个数，i为脉冲宽度内编码序号，且0≤i≤I；
Figure BDA0003616188010000073
为编码，rect(·)为门函数：
Figure BDA0003616188010000074
tb为每个编码的时间长度，T为脉冲重复周期。
根据
Figure BDA0003616188010000075
编码集合
Figure BDA0003616188010000076
中编码采用具备正交特性的编码，总个数为P×I个，以Frank编码为例，其编码形式为：
Figure BDA0003616188010000077
其中φi+pP＝2π·i·p/(P×I),i＝0…I-1,p＝0…P-1。
再进一步，步骤二中，每路空间-时间-编码信号中，对应第m个发射单元的空间-时间-编码信号为：
Figure BDA0003616188010000078
式中sm(t)为第m个发射单元上的空间-时间-编码信号，sp(t)为待发射时间编码信号s(t)的第p个脉冲，此时，p＝q。
显然，发射脉冲上集成复合了事先设计的发射方位信息，使得空间信息与时间-编码信号有效结合，也使得辐射出的多路空间-时间-编码信号在脉冲序列上实现了波束指向。这些波束指向既可以随时间固定，产生持续关注效果，也可以随脉冲间不断变化，从而产生发射波束扫描效果。
再进一步，结合图2所示，本实施方式还包括对具有方位信息的空间-时间-编码信号的回波信号进行处理的方法，具体包括：
步骤四：Q路具有方位信息的空间-时间-编码信号被观测对象反射后，获得Q个方位的观测对象回波信号；所述观测对象回波信号中包含了观测对象方位、距离时延、径向速度信息；
步骤五：对Q个方位的观测对象回波信号进行脉冲压缩处理，得到发射方位-接收通道-距离频谱矩阵；所述接收通道的个数N为接收天线阵列中接收单元的个数；N为正整数；
步骤六：将设定个数信号调制周期的发射方位-接收通道-距离频谱矩阵按照信号调制周期序号排列，获得排序后发射方位-接收通道-周期-距离频谱矩阵；
步骤七：对所有发射方位中每个发射方位的接收通道-周期-距离频谱使用与发射方位相同指向的波束加权系数进行数字波束形成处理，汇总所有发射方位的结果，得到波束-周期-距离频谱；
步骤八：对波束-周期-距离频谱中的每一列离散化距离单元信号，进行多普勒分析处理，得到具有空间扫描特性的波束-距离-多普勒频谱矩阵；处理完成所有接收波束后，处理结果可以组成波束-距离-多普勒三维矩阵，观测对象信息位于该矩阵中，呈现为多维尖峰特性；
步骤九：对发射波束方位集合中的每一个方位，重复步骤二至步骤八，完成具有凝视效果的观测空间扫描。至此，实现了空间扫描特性的编码信号设计及其处理过程，提取了所有发射方位方向的观测对象信息。
再进一步，步骤四中，观测对象回波信号表示为：
Figure BDA0003616188010000081
式中
Figure BDA0003616188010000082
为第q个发射波束方位下第n个接收通道的观测对象回波信号，n＝0，1，2，……，N-1；τ为观测对象时间延迟，l为信号调制周期计数，一个信号调制周期的时间长度Tm为脉冲个数与脉冲重复周期的乘积，即Tm＝PT；θp为观测对象回波信号方位；0≤t；
Figure BDA0003616188010000091
R为观测对象距离发射天线阵列中心的距离，R＞＞N·d；显然，若观测对象运动，则观测对象距离R随时间变化；接收天线阵列中接收单元与发射单元间距相同。
显然，当回波信号方位与发射方向相同时，根据θp＝θq，观测对象回波信号的表达式中，与发射方位θq相关的求和项
Figure BDA0003616188010000092
的计算结果为常数Aq。
再进一步，步骤五中获得发射方位-接收通道-距离频谱的方法包括：
第q个发射波束方位下的观测对象回波信号
Figure BDA0003616188010000093
进行脉冲压缩处理，得到单个信号调制周期当前方位的接收通道-距离频谱
Figure BDA0003616188010000094
Figure BDA0003616188010000095
式中r为观测对象相对高频地波雷达的离散化距离序号，由多通道回波信号
Figure BDA0003616188010000096
中的时延τ导出；
由Q个发射方位的观测对象回波信号获得单个信号调制周期的Q个接收通道-距离频谱，形成发射方位-接收通道-距离频谱矩阵，所述矩阵维度为Q×N×rmax，其中rmax为观测对象相对高频地波雷达的最大离散化距离序号。
再进一步，步骤六中，排序后发射方位-接收通道-周期-距离频谱矩阵W为：
Figure BDA0003616188010000097
式中
Figure BDA0003616188010000098
为第q个发射方位第n个接收通道的设定个数信号调制周期的距离频谱：
Figure BDA0003616188010000099
式中x为离散化距离单元信号，lmax为最大信号调制周期序号，
Figure BDA00036161880100000910
中x的行号l为信号调制周期序号，列号为离散化距离r的序号。
再进一步，步骤七中，得到波束-周期-距离频谱的方法包括：
Figure BDA0003616188010000101
式中Bq为波束-周期-距离频谱，
Figure BDA0003616188010000102
为与θq相关的第n个接收通道波束加权系数，所有
Figure BDA0003616188010000103
组成波束加权系数集合w(θq)：
Figure BDA0003616188010000104
式中an来自长度为N的窗函数。
再进一步，所述窗函数包括汉宁窗函数和切比雪夫窗函数；
当窗函数为汉宁窗函数，则：
Figure BDA0003616188010000105
显然，该接收波束的周期-距离信号Bq指向了发射方向θq，包含了该方向的所有回波。
依次取出其他发射方位的多通道-周期-距离频谱矩阵，并按照上述方法处理。处理后，接收波束指向与事先设计的发射方位一一对应，涵盖了所有方向的观测对象回波，成为了波束-周期-距离频谱，其维度为Q×lmax×rmax。
再进一步，步骤八中，对波束-周期-距离频谱Bq中的每一列离散化距离单元信号进行多普勒分析处理，得到第q个发射波束方位的距离-多普勒频谱矩阵Dq：
Figure BDA0003616188010000106
式中y为多普勒信号，y的行号dmax为多普勒单元序号，y的列号为离散化距离r的序号；大多数情况下，取dmax＝lmax，dmax＝pmax则：
第r列的多普勒信号yd(r)为：
Figure BDA0003616188010000107
式中wl来自给定lmax的窗函数。所述窗函数包括汉宁窗函数、海明窗函数等窗函数，可以依据多普勒分析需求自行选择。
如此，所有接收波束的距离-多普勒频谱Dq组成了波束-距离-多普勒三维矩阵ΓARD。在矩阵ΓARD中，观测对象可在其对应方位的波束、距离、径向速度的波束q、距离单元r、多普勒单元d处呈现为多维尖峰。取观测对象所在方位为0度波束的距离-多普勒频谱，结果如说明书附图12所示，其中颜色越深表明其幅度越高。同时，噪声和杂波在该矩阵中表现为面状等其他非尖峰特性，显然观测对象尖峰所在的距离、速度就是观测对象本身具有的信息，因此可以使用恒虚警检测发现该观测对象与噪声杂波的区别，并估计其真实方位、距离、多普勒速度等信息。这样就得到了空间扫描特性的编码信号及其处理过程，得到的结果体现了随着时间变化其方位的波束-距离-多普勒三维矩阵。
本实施方式中，若所需观测方向为固定方向，则更改发射方位集合中元素为固定值，重复步骤二至步骤八，就得到了空间凝视特性的编码信号及其处理过程，得到的结果体现了随着时间变化，始终观察固定方位的波束-距离-多普勒三维矩阵结果。
显然，本发明提及的编码信号及其处理方法的空间扫描特性表现在发射波束方位θq的变化上。当θq始终指向某一方位时，高频地波雷达观测产生长时关注某方向的效果；当θq随着观测对象运动引起的角度变化而相应变化时，高频地波雷达观测则产生波束跟随效果。
具体实施例：
具有空间扫描特性的高频地波雷达信号处理方法，包括以下步骤：
步骤一：设发射电磁波频率为10MHz，发射天线阵列单元间距为电磁波波长的二分之一，即间距为15米，且天线阵列中发射单元个数为8，发射波束方位集合{θp}中方位个数为5，设定为-30度、-10度、0度、+10度、+20度，得到的发射方位信号共有5组，每组中包含信号8个。
步骤二：设定待发射的时间编码信号中脉冲重复周期T为2ms，脉冲个数为5，编码时间长度为50us，单组编码集合长度L为5，此时编码集合总个数为25个。
步骤三：将多路发射方位信号与时间编码信号相乘，获得具有方位信息的多路空间-时间-编码信号。
步骤四：设单站高频地波雷达的接收天线阵列天线单元个数为8个，且天线单元间距也为15米，在接收天线阵列的每个单元上，相比发射信号，得到多通道观测对象回波信号。
作为示例，设定观测对象距离为150km、径向速度为15m/s、方位为0度，则相对发射信号，单站高频地波雷达的观测对象时间延迟为1ms，一个调制周期时长为10ms。
步骤五：按照脉冲上的发射方位属性，选取波束方位相同的脉冲，此时5个波束方位，构成5组沿时间推进的脉冲序列，然后进行脉冲压缩处理，获得该发射方位的多通道-距离频谱。此时，在时间采样间隔为50us条件下，距离频谱结果如图3所示，此时距离单元大小为7.5km，显然，观测对象尖峰位于第20个距离单元，即150km处。
步骤六：计算5个发射方位的多通道观测对象回波信号，获得单个调制周期的5个多通道-距离频谱，形成发射方位-多通道-距离矩阵，其维度为5×8×40。
步骤七：将多个调制周期的发射方位-多通道-距离频谱按照调制周期序号排列，调制周期为256个，重排为发射方位-多通道-周期-距离频谱矩阵，其维度5×8×40×256。
步骤八：某个发射方位的多通道-周期-距离频谱矩阵，使用与发射方位相同指向的波束加权系数进行数字波束形成处理，得到接收波束的周期-距离信号。此时，接收波束指向与发射方位相同，即每个发射方位产生一个接收波束。汇总所有发射方位的结果，得到波束-周期-距离频谱，其维度为5×40×256。
步骤九：针对每一个波束-周期-距离频谱中的每一列离散化距离单元的信号，进行多普勒分析处理，可以得波束-距离-多普勒频谱矩阵Dq，其维度为5×40×256。结果如图4至图11所示。其中多普勒分析处理采用的窗函数为256点的汉宁窗，其时间域波形如图10所示。
可以看到，当发射波束与观测对象方位对准时，观测对象反射能量最大，表明具有空间扫描特性的编码信号设计合理，并正确实现了扫描处理。
步骤十，当发射波束方位集合中元素有且仅有一个方向0度时，重复步骤二至步骤九，就得到了始终指向0度的空间扫描结果，此时波束-距离-多普勒频谱仅有指向0度的距离-多普勒频谱，如图12所示。若发射波束方位不变化，则波束-距离-多普勒结果也仅包含指向0度的距离-多普勒频谱，当观测对象进入该方向时，可以获得最大反射能量。表明了空间扫描时的凝视观测效果。
综上所述，本发明提出的具有发射波束空间扫描功能的编码信号及其处理方法，使得高频地波雷达能够在实现发射波束空间扫描的同时，实现正确地提取观测对象回波信息。具有实施简单方便、可自适应的改变权值等特点。
虽然在本文中参照了特定的实施方式来描述本发明，但是应该理解的是，这些实施例仅仅是本发明的原理和应用的示例。因此应该理解的是，可以对示例性的实施例进行许多修改，并且可以设计出其他的布置，只要不偏离所附权利要求所限定的本发明的精神和范围。应该理解的是，可以通过不同于原始权利要求所描述的方式来结合不同的从属权利要求和本文中所述的特征。还可以理解的是，结合单独实施例所描述的特征可以使用在其它所述实施例中。

Claims (9)

1.一种具有空间扫描特性的高频地波雷达信号处理方法，其特征在于包括具有空间扫描特性的高频地波雷达信号的设计方法，具体包括：
步骤一：设定包含Q个方位个数的发射波束方位集合，生成Q路发射方位信号，每路发射方位信号中发射方位信号的个数与发射天线阵列中发射单元个数相同；
步骤二：设定每个信号调制周期包括P个脉冲的待发射时间编码信号；其中Q、M和P均为正整数；
步骤三：使P＝Q，按序号将每路发射方位信号中的每个发射方位信号与对应脉冲序号的待发射时间编码信号相乘，获得Q路具有方位信息的空间-时间-编码信号，Q路具有方位信息的空间-时间-编码信号发射至观测空间进行观测；每路空间-时间-编码信号中包括M个空间-时间-编码信号；
步骤一中，生成的Q路发射方位信号中第q路发射方位信号表示为：
Figure FDA0004112608570000011
式中βm为发射方位信号，θq为发射波束方位，m为发射天线阵列中发射单元的序号，d为发射单元间距，c为光速，f0是发射电磁波频率；q为发射方位序号，0≤q≤Q-1；
步骤二中，待发射时间编码信号表示为：
Figure FDA0004112608570000012
式中s(t)为待发射时间编码信号，t为时间，p为信号调制周期中脉冲序号，I为每个脉冲宽度内的编码个数，i为脉冲宽度内编码序号，且0≤i≤l；
Figure FDA0004112608570000013
为编码，rect(·)为门函数：
Figure FDA0004112608570000014
tb为每个编码的时间长度，T为脉冲重复周期。
2.根据权利要求1所述的具有空间扫描特性的高频地波雷达信号处理方法，其特征在于，步骤三中，每路空间-时间-编码信号中，对应第m个发射单元的空间-时间-编码信号为：
Figure FDA0004112608570000015
式中sm(t)为第m个发射单元的空间-时间-编码信号，sp(t)为待发射时间编码信号s(t)的第p个脉冲，p＝q。
3.根据权利要求2所述的具有空间扫描特性的高频地波雷达信号处理方法，其特征在于，还包括对具有方位信息的空间-时间-编码信号的回波信号进行处理的方法，具体包括：
步骤四：Q路具有方位信息的空间-时间-编码信号被观测对象反射后，获得Q个方位的观测对象回波信号；
步骤五：对Q个方位的观测对象回波信号进行脉冲压缩处理，得到发射方位-接收通道-距离频谱矩阵；所述接收通道的个数N为接收天线阵列中接收单元的个数；N为正整数；
步骤六：将设定个数信号调制周期的发射方位-接收通道-距离频谱矩阵按照信号调制周期序号排列，获得排序后发射方位-接收通道-周期-距离频谱矩阵；
步骤七：对所有发射方位中每个发射方位的接收通道-周期-距离频谱使用与发射方位相同指向的波束加权系数进行数字波束形成处理，得到波束-周期-距离频谱；
步骤八：对波束-周期-距离频谱中的每一列离散化距离单元信号，进行多普勒分析处理，得到具有空间扫描特性的波束-距离-多普勒频谱矩阵；
步骤九：对发射波束方位集合中的每一个方位，重复步骤二至步骤八，完成具有凝视效果的观测空间扫描。
4.根据权利要求3所述的具有空间扫描特性的高频地波雷达信号处理方法，其特征在于，步骤四中，观测对象回波信号表示为：
Figure FDA0004112608570000021
式中
Figure FDA0004112608570000022
为第q个发射波束方位下第n个接收通道的观测对象回波信号，n＝0，1，2，……，N-1；τ为观测对象时间延迟，l为信号调制周期计数；θp为观测对象回波信号方位；
Figure FDA0004112608570000023
R为观测观测对象距离发射天线阵列中心的距离，R＞＞N·d；
根据θp＝θq，观测对象回波信号的表达式中，
Figure FDA0004112608570000024
的计算结果为常数Aq。
5.根据权利要求4所述的具有空间扫描特性的高频地波雷达信号处理方法，其特征在于，步骤五中获得发射方位-接收通道-距离频谱的方法包括：
第q个发射波束方位下的观测对象回波信号
Figure FDA0004112608570000025
进行脉冲压缩处理，得到单个信号调制周期当前方位的接收通道-距离频谱
Figure FDA0004112608570000031
Figure FDA0004112608570000032
式中r为观测对象相对高频地波雷达的离散化距离序号；
由Q个发射方位的观测对象回波信号获得单个信号调制周期的Q个接收通道-距离频谱，形成发射方位-接收通道-距离频谱矩阵，所述矩阵维度为Q×N×rmax，其中rmax为观测对象相对高频地波雷达的最大离散化距离序号。
6.根据权利要求5所述的具有空间扫描特性的高频地波雷达信号处理方法，其特征在于，步骤六中，排序后发射方位-接收通道-周期-距离频谱矩阵W为：
Figure FDA0004112608570000033
式中
Figure FDA0004112608570000034
为第q个发射方位第n个接收通道的设定个数信号调制周期的距离频谱：
Figure FDA0004112608570000035
式中x为离散化距离单元信号，lmax为最大信号调制周期序号，
Figure FDA0004112608570000036
中x的行号l为信号调制周期序号，列号为离散化距离r的序号。
7.根据权利要求6所述的具有空间扫描特性的高频地波雷达信号处理方法，其特征在于，步骤七中，得到波束-周期-距离频谱的方法包括：
Figure FDA0004112608570000037
式中Bq为波束-周期-距离频谱，
Figure FDA0004112608570000038
为与θq相关的第n个接收通道波束加权系数，所有
Figure FDA0004112608570000039
组成波束加权系数集合w(θq)：
Figure FDA00041126085700000310
式中an来自长度为N的窗函数。
8.根据权利要求7所述的具有空间扫描特性的高频地波雷达信号处理方法，其特征在于，所述窗函数包括汉宁窗函数和切比雪夫窗函数；
当窗函数为汉宁窗函数，则：
Figure FDA0004112608570000041
9.根据权利要求8所述的具有空间扫描特性的高频地波雷达信号处理方法，其特征在于，步骤八中，对波束-周期-距离频谱Bq中的每一列离散化距离单元信号进行多普勒分析处理，得到第q个发射波束方位的距离-多普勒频谱矩阵Dq：
Figure FDA0004112608570000042
式中y为多普勒信号，y的行号dmax为多普勒单元序号，y的列号为离散化距离r的序号；取dmax＝lmax，则：
第r列的多普勒信号yd(r)为：
Figure FDA0004112608570000043
式中wl来自给定lmax的窗函数。
CN202210444733.9A 2022-04-26 2022-04-26 具有空间扫描特性的高频地波雷达信号处理方法 Active CN114859296B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210444733.9A CN114859296B (zh) 2022-04-26 2022-04-26 具有空间扫描特性的高频地波雷达信号处理方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210444733.9A CN114859296B (zh) 2022-04-26 2022-04-26 具有空间扫描特性的高频地波雷达信号处理方法

Publications (2)

Publication Number Publication Date
CN114859296A CN114859296A (zh) 2022-08-05
CN114859296B true CN114859296B (zh) 2023-04-14

Family

ID=82633954

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210444733.9A Active CN114859296B (zh) 2022-04-26 2022-04-26 具有空间扫描特性的高频地波雷达信号处理方法

Country Status (1)

Country Link
CN (1) CN114859296B (zh)

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB0327729D0 (en) * 2001-07-26 2003-12-31 Codar Ocean Sensors Ltd Multi-station HF FMCW radar frequency sharing with GPS time modulation multiplexing
CN1893415A (zh) * 2005-06-29 2007-01-10 美国博通公司 使用复数信号交换的块编码模块和块编码方法
CN101295019A (zh) * 2008-05-23 2008-10-29 北京航空航天大学 一种线性调频连续波合成孔径雷达视频信号的仿真方法
CN101507135A (zh) * 2006-08-21 2009-08-12 皇家飞利浦电子股份有限公司 用于多站点和多波束传输的空间-时间/空间-频率编码
US9154172B1 (en) * 2013-12-31 2015-10-06 Hrl Laboratories, Llc Time encoded circuits and methods and a time encoder based beamformer for use in receiving and transmitting applications
CN105445701A (zh) * 2015-11-11 2016-03-30 西安电子科技大学 Ddma-mimo雷达目标的单脉冲角度估计方法
CN107942327A (zh) * 2017-10-20 2018-04-20 西安电子科技大学 基于脉冲相位编码的单通道hrws‑sar成像方法
CN110412568A (zh) * 2019-07-05 2019-11-05 西安电子科技大学 基于扩展方位角相位编码的距离模糊杂波抑制方法
CN112180367A (zh) * 2020-09-09 2021-01-05 中国科学院空天信息创新研究院 多输入多输出合成孔径雷达的回波分离方法、装置及终端
RU2746799C1 (ru) * 2020-07-14 2021-04-21 Федеральное Государственное Казенное Военное Образовательное Учреждение Высшего Образования "Тихоокеанское Высшее Военно-Морское Училище Имени С.О. Макарова" Министерства Обороны Российской Федерации (Г. Владивосток) Статистический способ имитационной защиты загоризонтной радиолокационной станции поверхностной волны
CN113033754A (zh) * 2020-12-08 2021-06-25 中国海洋大学 基于协同情景进化的高频地波雷达目标跟踪算法测评方法
CN113572509A (zh) * 2021-06-23 2021-10-29 南京理工大学 基于时间调制阵列的mimo雷达系统及波束形成方法

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB0327729D0 (en) * 2001-07-26 2003-12-31 Codar Ocean Sensors Ltd Multi-station HF FMCW radar frequency sharing with GPS time modulation multiplexing
CN1893415A (zh) * 2005-06-29 2007-01-10 美国博通公司 使用复数信号交换的块编码模块和块编码方法
CN101507135A (zh) * 2006-08-21 2009-08-12 皇家飞利浦电子股份有限公司 用于多站点和多波束传输的空间-时间/空间-频率编码
CN101295019A (zh) * 2008-05-23 2008-10-29 北京航空航天大学 一种线性调频连续波合成孔径雷达视频信号的仿真方法
US9154172B1 (en) * 2013-12-31 2015-10-06 Hrl Laboratories, Llc Time encoded circuits and methods and a time encoder based beamformer for use in receiving and transmitting applications
CN105445701A (zh) * 2015-11-11 2016-03-30 西安电子科技大学 Ddma-mimo雷达目标的单脉冲角度估计方法
CN107942327A (zh) * 2017-10-20 2018-04-20 西安电子科技大学 基于脉冲相位编码的单通道hrws‑sar成像方法
CN110412568A (zh) * 2019-07-05 2019-11-05 西安电子科技大学 基于扩展方位角相位编码的距离模糊杂波抑制方法
RU2746799C1 (ru) * 2020-07-14 2021-04-21 Федеральное Государственное Казенное Военное Образовательное Учреждение Высшего Образования "Тихоокеанское Высшее Военно-Морское Училище Имени С.О. Макарова" Министерства Обороны Российской Федерации (Г. Владивосток) Статистический способ имитационной защиты загоризонтной радиолокационной станции поверхностной волны
CN112180367A (zh) * 2020-09-09 2021-01-05 中国科学院空天信息创新研究院 多输入多输出合成孔径雷达的回波分离方法、装置及终端
CN113033754A (zh) * 2020-12-08 2021-06-25 中国海洋大学 基于协同情景进化的高频地波雷达目标跟踪算法测评方法
CN113572509A (zh) * 2021-06-23 2021-10-29 南京理工大学 基于时间调制阵列的mimo雷达系统及波束形成方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
张忠传 ; 牟善祥 ; 李成国 ; 赵红梅 ; .宽带相控阵雷达信号空时自适应仿真分析.系统仿真学报.2007,(第17期),4068-4078. *

Also Published As

Publication number Publication date
CN114859296A (zh) 2022-08-05

Similar Documents

Publication Publication Date Title
CN111693983B (zh) 互补波形构建方法及模块、认知雷达系统及波形发射方法
CN114152916A (zh) 一种基于先导信息的mimo雷达波形正交方法
CN111830482B (zh) 基于捷变ofdm的fda雷达目标定位方法
CN106338723A (zh) 一种基于互质脉冲重复间隔的空时自适应处理方法及装置
Rabaste et al. Signal waveforms and range/angle coupling in coherent colocated MIMO radar
CN103308892B (zh) 一种多扫频雷达信号的产生及压缩方法
CN111736126A (zh) 一种基于低成本安防雷达的信号处理方法
CN111736141A (zh) 基于循环编码阵列的雷达通信一体化设计方法
Yang et al. DDMA MIMO radar system for low, slow, and small target detection
Lu et al. Co-located MIMO radar target detection in cluttered and noisy environment based on 2D block sparse recovery
CN112965034B (zh) 一种改善天波雷达慢时间相位编码信号多普勒容限的方法
CN108490425B (zh) 一种双基地mimo雷达的测角方法
CN114859296B (zh) 具有空间扫描特性的高频地波雷达信号处理方法
CN111580040A (zh) 双基地展开互质阵列mimo雷达dod和doa降维估计方法
CN114764136A (zh) 基于多时间尺度耦合网络的雷达抗干扰波形产生方法
CN114325595A (zh) 联合发射阵列与脉冲编码的机载mimo-sar距离模糊分离方法
CN115407335A (zh) 基于OFDM-chirp信号和DBF处理的MIMO-SAR解模糊方法及装置
CN115685080A (zh) 一种雷达时间-多普勒双域联合编码波形设计与处理方法
CN114624694A (zh) 一种基于非等间隔调制ddma-mimo的雷达信号处理方法
CN113204018A (zh) 一种基于lfmcw的mimo雷达快速信号处理方法
CN113093146A (zh) 基于量子跳跃逃逸机制的mimo雷达正交波形设计方法
CN111025292A (zh) 一种四发mimo雷达发射波形设计方法
Bezoušek et al. MIMO radar signals with better correlation characteristics
CN113204017B (zh) 一种基于lfmcw的mimo雷达信号处理方法
CN113376590B (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