CN109061576B - 一种射频掩护信号下的阵列雷达信号发射接收方法 - Google Patents
一种射频掩护信号下的阵列雷达信号发射接收方法 Download PDFInfo
- Publication number
- CN109061576B CN109061576B CN201810650437.8A CN201810650437A CN109061576B CN 109061576 B CN109061576 B CN 109061576B CN 201810650437 A CN201810650437 A CN 201810650437A CN 109061576 B CN109061576 B CN 109061576B
- Authority
- CN
- China
- Prior art keywords
- signal
- radar
- party
- signals
- transmitting
- 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
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
- G01S7/285—Receivers
- G01S7/292—Extracting wanted echo-signals
-
- 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
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
本发明公开了一种射频掩护信号下的阵列雷达信号发射接收方法,其思路为:我方雷达发射信号,确定我方雷达发射信号中的干扰信号,然后构建目标函数并进行求解,得到我方雷达发射信号中的最优探测信号矩阵后与我方雷达发射信号中的干扰信号在雷达信号发射天线网络同时发射,经对方目标反射后得到干扰低的雷达探测回波信号,雷达信号接收天线网络对所述干扰低的雷达探测回波信号进行接收,并对接收到的干扰低的雷达探测回波信号进行融合积累后,得到融合积累后的雷达探测回波信号,对融合积累后的雷达探测回波信号进行脉冲压缩处理后,从而得到对方目标位置,所述对方目标位置为一种射频掩护信号下的阵列雷达信号发射接收结果。
Description
技术领域
本发明属于通信技术领域,特别涉及一种射频掩护信号下的阵列雷达信号发射接收方法,适用于对对方进行有效干扰的同时使我方获得探测能力。
背景技术
在现代的雷达对抗系统中,根据干扰能量来源的不同,可分为有源干扰和无源干扰两种;在这两种干扰方式中,有源干扰的配置更加灵活、破坏性更大,使其成为当前研究的最广泛的一种干扰形式。
雷达的有源干扰,指的是雷达通过发射特定的信号干扰对方的雷达,根据发射信号的类型分为压制性干扰和欺骗性干扰两种;在压制性干扰中,雷达发射大能量的压制性干扰信号用于干扰对方的雷达,使其失去工作能力,但此时我方需要调用其它雷达发射探测信号用于探测对方的位置,这样会造成我方雷达在时间、空间及资源上极大的浪费。
发明内容
针对上述现有技术存在的问题,本发明的目的在于提出一种射频掩护信号下的阵列雷达信号发射接收方法,该种射频掩护信号下的阵列雷达信号发射接收方法是在利用己方电磁干扰机对敌方进行干扰的同时发射相对于干扰信号正交的雷达探测信号,在对敌方进行有效干扰的同时使己方获得探测能力,同时保证己方探测行为的隐蔽性。
为达到上述技术目的,本发明采用如下技术方案予以实现。
一种射频掩护信号下的阵列雷达信号发射接收方法,包括以下步骤:
步骤1,确定我方和对方,我方包括我方雷达,对方包括对方雷达和对方目标,对方雷达置于对方目标中;我方雷达发射信号,我方雷达发射的信号包括探测信号和干扰信号;
步骤2,确定我方雷达发射信号中的干扰信号;
步骤3,根据我方雷达发射信号中的探测信号和干扰信号,构建目标函数;
步骤4,对所述目标函数进行求解,得到我方雷达发射信号中的最优探测信号矩阵;
步骤5,设定雷达信号发射天线网络和雷达信号接收天线网络,然后将我方雷达发射信号中的最优探测信号矩阵和我方雷达发射信号中的干扰信号在所述雷达信号发射天线网络同时发射,经对方目标反射后得到干扰低的雷达探测回波信号;
步骤6,雷达信号接收天线网络对所述干扰低的雷达探测回波信号进行接收,并对接收到的干扰低的雷达探测回波信号进行融合积累后,得到融合积累后的雷达探测回波信号,对融合积累后的雷达探测回波信号进行脉冲压缩处理后,从而得到对方目标位置,所述对方目标位置为一种射频掩护信号下的阵列雷达信号发射接收结果。
本发明的有益效果:
本发明方法采用射频掩护的思想,同时设计探测信号与干扰信号,探测信号与干扰信号同时工作,互不影响,能达到干扰对方的同时探测信号的目的;探测信号具有低距离旁瓣,可防止大功率目标回波淹没附近距离单元中的小功率目标回波;另外,本发明方法利用单部雷达即可达到同时干扰和探测的效果,与现有的雷达探测与干扰形式相比,节约资源和成本。
附图说明
下面结合附图和具体实施方式对本发明作进一步详细说明。
图1为本发明的一种射频掩护信号下的阵列雷达信号发射接收方法流程图;
图2a为本发明发射端为第一发射天线阵列时的处理过程流程示意图;
图2b为本发明发射端为第二发射天线阵列时的处理过程流程示意图;
图3是干扰信号与探测信号的互相关仿真示意图;
图4是探测信号的自相关仿真示意图;
图5是干扰信号与探测信号的互相关仿真示意图;
图6是探测信号的自相关仿真示意图;
图7是干扰信号的采样率为40MHz的条件下干扰信号与探测信号的互相关图;
图8是干扰信号的采样率为80MHz的条件下干扰信号与探测信号的互相关图;
图9是干扰信号长度为探测信号长度的5倍条件下干扰信号与探测信号的互相关图;
图10是干扰信号长度为探测信号长度的5倍条件下探测信号自相关图;
图11是干扰信号加窗的情况下干扰信号与探测信号的互相关图;
图12是干扰信号加窗的情况下探测信号自相关图。
具体实施方式
参照图1,为本发明的一种射频掩护信号下的阵列雷达信号发射接收方法流程图;其中所述射频掩护信号下的阵列雷达信号发射接收方法,包括以下步骤:
步骤1,确定我方和对方,我方包括我方雷达,对方包括对方雷达和对方目标,对方雷达置于对方目标中;我方雷达发射信号,且我方雷达发射的信号包括探测信号和干扰信号,其中探测信号用于探测对方目标,干扰信号用于干扰对方目标中的对方雷达。
步骤2,确定我方雷达发射信号中干扰信号的参数:根据我方雷达发射信号的天线所能承受的最大发射功率设定我方雷达发射信号中干扰信号的时宽及带宽;通过所述我方雷达发射信号中干扰信号的时宽及带宽得到的干扰信号功率应小于我方雷达发射信号的天线所能承受的最大发射功率。
具体地,本发明方法采用射频掩护的思想确定我方雷达发射信号中的干扰信号和我方雷达发射信号中的探测信号,我方雷达发射信号中的干扰信号类型确定为线性调频信号,根据线性调频信号的表达式及干扰信号的时宽及带宽,得到我方雷达发射信号中的连续干扰信号sjam(t),其表达式为:
其中,rect(·)表示矩形窗函数,T表示我方雷达发射信号中的连续干扰信号时宽,exp(·)表示指数函数,j表示虚数单位,π表示圆周率,μ表示我方雷达发射信号中的连续干扰信号调频率,μ=B/T,B表示我方雷达发射信号中的连续干扰信号带宽,t表示时间变量。
对我方雷达发射信号中的连续干扰信号sjam(t)进行采样点数为Nt的采样后得到的离散信号,记为我方雷达发射信号中的干扰信号Sjam;Nt≥1。
根据实际需要确定我方雷达发射信号中的探测信号总个数Nt以及我方雷达发射信号中每个探测信号的码元长度Ns。
步骤3,根据我方雷达发射信号中的探测信号和干扰信号,构建目标函数。
3a)设定我方雷达发射信号中的探测信号包括Na个角频率,将所述Na个角频率分别进行归一化后得到Na个归一化角频率,记为我方雷达发射探测信号时的归一化角频率方向集合fa,fam表示我方雷达发射探测信号时的第m个归一化角频率,m=1,2,…,Na,Na表示我方雷达发射探测信号时的归一化角频率方向集合fa中包括的归一化角频率总个数。
根据我方雷达发射探测信号时的归一化角频率方向集合fa,得到我方雷达发射信号中探测信号的自相关峰值旁瓣电平PSLra的计算公式如下:
其中,PSLra表示我方雷达发射信号中探测信号的自相关峰值旁瓣电平,max表示取最大值操作,k=1,2,…,Ns-1,Ns表示我方雷达发射信号中每个探测信号的码元长度,∈表示属于,|·|表示取模值操作,at(fam)表示我方雷达发射信号中探测信号在fam方向形成的导向矢量,
at(fam)=[1,exp(j2πfam),…,exp(j2π(Nt-1)fam)]T,Nt表示我方雷达发射信号中的探测信号总个数,(·)H表示取共轭转置操作,(·)T表示取转置操作,S表示我方雷达发射信号中的探测信号,si表示我方雷达发射信号中的第i个探测信号,表示(Ns-k)×k维全零矩阵,0k×k表示k×k维全零矩阵,表示(Ns-k)×(Ns-k)维全1矩阵,表示k×(Ns-k)维全零矩阵。
3b)将我方雷达发射信号中的探测信号发射后经由对方目标反射回来、并由我方雷达接收到的信号,记为我方雷达接收的探测信号,所述我方雷达接收的探测信号包括Na个角频率,将所述Na个角频率分别归一化后得到我方雷达接收探测信号时的归一化角频率集合f'a,
根据我方雷达接收探测信号时的归一化角频率集合f'a和我方雷达发射探测信号时的归一化角频率方向集合fa,得到我方雷达发射信号中的探测信号与我方雷达接收的探测信号之间的互相关距离旁瓣电平PSLrc计算公式如下:
其中,PSLrc表示我方雷达发射信号中的探测信号与我方雷达接收的探测信号之间的互相关距离旁瓣电平,p=-(Ns-1),-(Ns-2),…,0,1,2,…,Ns-1,Ns表示我方雷达发射信号中每个探测信号的码元长度,n=1,2,…,Nb,Nb表示我方雷达发射信号中的探测信号的归一化角频率组合集合中包括的归一化角频率组合总个数,at(fam)表示我方雷达发射信号中的探测信号在fam方向形成的导向矢量,at(f'am)表示我方雷达发射信号中的探测信号在f'am方向形成的导向矢量,ar(fam)表示我方雷达接收探测信号时在fam方向形成的导向矢量,ar(f'am)表示我方雷达接收探测信号时在f'am方向形成的导向矢量,
at(f'am)=[1,exp(j2πf'am),…,exp(j2π(Nt-1)f'am)]T,
ar(f'am)=[1,exp(j2πf'am),…,exp(j2π(Nr-1)f'am)]T,
ar(fam)=[1,exp(j2πfam),…,exp(j2π(Nr-1)fam)]T,Nr表示我方雷达包括的接收天线个数,Nt表示我方雷达发射信号中的探测信号总个数,表示(Ns-p)×p维全零矩阵,
3c)设定我方雷达发射信号中的干扰信号包括Na个角频率,将所述Na个角频率分别归一化后得到我方雷达发射干扰信号时的Na个归一化角频率fb,
然后根据我方雷达发射探测信号时的归一化角频率方向集合fa和我方雷达发射干扰信号时的Na个归一化角频率fb,得到我方雷达发射信号中的探测信号与干扰信号在相同方向的互相关旁瓣电平PSLjam1的计算公式如下:
其中,PSLjam1表示我方雷达发射信号中的探测信号与干扰信号在相同方向的互相关旁瓣电平,p=-(Ns-1),-(Ns-2),…,0,1,2,…,Ns-1,Ns表示我方雷达发射信号中每个探测信号的码元长度,Sjam表示我方雷达发射信号中的干扰信号,S表示我方雷达发射信号中的探测信号,at(fbm)表示我方雷达发射信号中的探测信号在fbm方向形成的导向矢量,
at(fbm)=[1,exp(j2πfbm),…,exp(j2π(Nt-1)fbm)]T,m=1,2,…,Na,Nt表示我方雷达发射诱导信号中的探测信号总个数,Nr表示我方雷达包括的接收天线总个数,(·)T表示取转置操作,表示(Ns-p)×p维全零矩阵,0p×p表示p×p维全零矩阵,表示(Ns-p)×(Ns-p)维全1矩阵,表示p×(Ns-p)维全零矩阵,(·)H表示取共轭转置操作。
3d)将我方雷达发射信号中的干扰信号发射后经由对方目标反射回来、并由我方雷达接收到的信号,记为我方雷达接收的干扰信号,我方雷达接收的干扰信号包括Na个角频率,将所述Na个角频率分别归一化后得到我方雷达接收干扰信号时的归一化角频率集合f'b,
然后根据我方雷达发射探测信号时的归一化角频率方向集合fa和我方雷达接收干扰信号时的归一化角频率集合f'b,得到我方雷达发射信号中的探测信号与接收的干扰信号在不同方向的互相关旁瓣电平PSLjam2的计算公式如下:
其中,PSLjam2表示我方雷达发射信号中的探测信号与接收的干扰信号在不同方向的互相关旁瓣电平,p=-(Ns-1),-(Ns-2),…,0,1,2,…,Ns-1,Ns表示我方雷达发射信号中每个探测信号的码元长度,at(fam)表示我方雷达发射信号中的探测信号在fam方向形成的导向矢量,at(f'bm)表示我方雷达发射信号中的探测信号在f'bm方向形成的导向矢量,ar(fam)表示我方雷达接收干扰信号时在fam方向形成的导向矢量,ar(f'bm)表示我方雷达接收干扰信号时在f'bm方向形成的导向矢量,at(f'bm)=[1,exp(j2πf'bm),…,exp(j2π(Nt-1)f'bm)]T,
ar(f'bm)=[1,exp(j2πf'bm),…,exp(j2π(Nr-1)f'bm)]T,Nt表示我方雷达发射信号中的探测信号总个数,Nr表示我方雷达包括的接收天线个数,表示(Ns-p)×p维全零矩阵,0p×p表示p×p维全零矩阵,表示(Ns-p)×(Ns-p)维全1矩阵,表示p×(Ns-p)维全零矩阵,(·)H表示取共轭转置操作,(·)T表示取转置操作,Sjam表示我方雷达发射信号中的干扰信号,S表示我方雷达发射信号中的探测信号。
3e)以同时最小化我方雷达发射信号中探测信号的自相关峰值旁瓣电平PSLra、我方雷达发射信号中探测信号与我方雷达接收的探测信号之间的互相关距离旁瓣电平PSLrc、我方雷达发射信号中探测信号与干扰信号在相同方向的互相关旁瓣电平PSLjam1和我方雷达发射信号中探测信号与接收的干扰信号在不同方向的互相关旁瓣电平PSLjam2为准则,构建目标函数:
其中,PSL表示我方雷达发射信号的峰值旁瓣电平值,min表示取最小值操作,P表示我方雷达发射信号中的探测信号相位矩阵,||·||p表示取p范数操作,α、β、η分别表示正的加权系数。
具体地,运用SQP对目标函数进行迭代求解,得到我方雷达发射信号中的探测信号最优相位矩阵SQP算法可以在短时间内找到一个局部最小解,但其结果的产生依赖于初值,因此在求解我方雷达发射信号中的探测信号最优相位矩阵时需要设置初值然后再进行优化。
4a)初始化:设置我方雷达发射信号中的探测信号总个数Nt以及我方雷达发射信号中每个探测信号的码元长度Ns;设置在求解时最为匹配的循环次数N;设置最高旁瓣标志为无穷大;令l表示第l次循环,l的初始值为1,l=1,2,…,N,N为大于1的正整数。
4b)确定第l次循环后我方雷达发射信号中的探测信号初始化相位矩阵Pl,即第l次循环后我方雷达发射信号中的探测信号初始化相位矩阵Pl中的每个元素均设置成一个随机相位值,每个相位值的取值范围均为0~2π。
4c)调用SQP算法优化第l次循环后我方雷达发射信号中的探测信号初始化相位矩阵Pl后,得到第k次优化后我方雷达发射信号中的探测信号相位矩阵进而计算得到第l次优化后我方雷达发射信号中的探测信号矩阵循环次数与优化次数相等且一一对应。
4d)根据第l次优化后我方雷达发射信号中的探测信号矩阵利用目标函数计算得到第l次循环后我方雷达发射信号中的探测信号峰值旁瓣电平值PSLl,比较第l次循环后我方雷达发射信号中的探测信号峰值旁瓣电平值PSLl与最高旁瓣标志的大小,如果第l次循环后我方雷达发射信号中的探测信号峰值旁瓣电平值PSLl小于最高旁瓣标志,则保存此次优化的结果,即将第k次优化后我方雷达发射信号中的探测信号相位矩阵作为第g个我方雷达发射信号中的探测信号优化相位矩阵Pg,将第l次优化后我方雷达发射信号中的探测信号矩阵作为第g个我方雷达发射信号中的优化探测信号矩阵g的初始值为1,并令g的值加1;然后令l的值加1,返回4b);如果第l次循环后我方雷达发射信号中的探测信号峰值旁瓣电平值PSLl大于或等于最高旁瓣标志,则忽略此次优化的结果,然后令l的值加1,并返回4b)。
直到N次循环结束后,此时得到g个我方雷达发射信号中的探测信号优化相位矩阵,[P1,P2,…,Pg],以及g个我方雷达发射信号中的优化探测信号矩阵然后根据目标函数分别计算第1个我方雷达发射信号中的探测信号峰值旁瓣电平值PSL1至第g个我方雷达发射信号中的探测信号峰值旁瓣电平值PSLg,记为我方雷达发射信号中探测信号的g个峰值旁瓣电平值,比较我方雷达发射信号中探测信号的g个峰值旁瓣电平值大小,将其中最小的峰值旁瓣电平值所对应的优化相位矩阵,记为我方雷达发射信号中的探测信号最优相位矩阵然后根据我方雷达发射信号中的探测信号最优相位矩阵计算得到我方雷达发射信号中的最优探测信号矩阵所述我方雷达发射信号中的最优探测信号矩阵包括Nt个最优探测信号。
5.1确定我方雷达包括雷达信号发射天线网络和雷达信号接收天线网络,雷达信号发射天线网络包括A个发射天线,当雷达信号发射天线网络中的A个发射天线个数与我方雷达发射信号中探测信号的波形个数取值相等时,雷达信号发射天线网络包括Nt个发射天线,雷达信号发射天线网络中的Nt个发射天线规格完全相同,进而组成第一发射天线阵列,此时A=Nt。
当雷达信号发射天线网络中的A个发射天线个数是我方雷达发射信号中的探测信号总个数2倍时,即雷达信号发射天线网络包括2Nt个发射天线,雷达信号发射天线网络中的2Nt个发射天线规格完全相同,进而组成第二发射天线阵列,此时A=2Nt;所述2Nt个发射天线分为Nt个探测信号天线和Nt个干扰信号天线,进而得到Nt个发射天线组,每个发射天线组包括1个探测信号天线和1个干扰信号天线,且每个发射天线组分别对应我方雷达发射信号中的最优探测信号矩阵中1个最优探测信号。
雷达信号接收天线网络是由Nr个规格完全相同的发射天线组成的天线阵列,由于干扰信号波束并没有携带目标信息,需要将干扰信号波束的能量降到最低,以防止干扰信号波束对探测信号造成干扰。
5.2如图2a所示,当A=Nt时:在第i个发射天线中,我方雷达发射信号中的最优探测信号矩阵中第i个最优探测信号通过调相器,得到第i个最优探测信号波束,第i个最优探测信号波束通过饱和放大器以放大探测信号能量,得到第i个经过放大的最优探测信号波束;我方雷达发射信号中的干扰信号Sjam通过调相器,得到干扰信号波束,干扰信号波束通过饱和放大器以放大干扰信号能量,得到经过放大的干扰信号波束;其中,i=1,2,…,Nt。
如图2b所示,当A=2Nt时:在第i个发射天线组中,我方雷达发射信号中的最优探测信号矩阵中第i个最优探测信号通过调相器,得到第i个最优探测信号波束,第i个最优探测信号波束通过饱和放大器以放大探测信号能量,得到第i个经过放大的最优探测信号波束;在Nt个干扰信号天线中,我方雷达发射信号中的干扰信号Sjam通过调相器,得到干扰信号波束,干扰信号波束通过饱和放大器以放大干扰信号能量,得到经过放大的干扰信号波束,所述经过放大的干扰信号波束包括Nt个经过放大的干扰信号;其中,i=1,2,…,Nt。
5.3如图2a所示,当A=Nt时:第i个经过放大的最优探测信号波束和经过放大的干扰信号波束通过功率合成网络,得到第i个合成后的雷达发射信号,然后利用发射波束置零的方法得到第i个发射天线的加权系数ωi;将第i个合成后的雷达发射信号乘以第i个发射天线的加权系数ωi后得到第i个最终的雷达发射信号,所述第i个最终的雷达发射信号经由第i个发射天线发射,以探测对方目标以及干扰对方目标中的对方雷达。
如图2b所示,当A=2Nt时:对第i个发射天线组利用发射波束置零的方法分别得到第i个探测信号天线的加权系数和第i个干扰信号天线的加权系数将第i个经过放大的最优探测信号波束乘以第i个探测信号天线的加权系数后得到第i个最终雷达探测信号,所述第i个最终雷达探测信号经由第i个发射天线组发射,以探测对方雷达;将第i个经过放大的干扰信号乘以第i个干扰信号天线的加权系数后得到第i个最终雷达干扰信号,所述第i个最终雷达干扰信号经由第i个发射天线组发射,以干扰对方目标中的对方雷达;将所述第i个最终雷达探测信号和所述第i个最终雷达干扰信号记为第i个最终的雷达发射信号。
5.4如果我方雷达发射第i个最终的雷达发射信号探测到对方目标后被对方目标反射,则将反射后得到的信号记为第g'个目标回波信号,雷达信号接收天线网络对第g'个目标回波信号进行接收,g'的初始值为1,并令g'的值加1;如果我方雷达发射第i个最终探测信号没有探测到对方目标,则忽略第i个最终的探测信号。
5.5在每一个接收天线中,对第g'个目标回波信号经过滤波处理,使得第g'个目标回波信号中干扰信号的能量降到最低,从而得到第g'个干扰低的雷达探测回波波束,将第g'个干扰低的雷达探测回波波束通过调相器后得到第g'个干扰低的雷达探测回波信号。
5.6令i的值分别取1至Nt,重复执行5.2至5.5,进而分别得到第1个干扰低的雷达探测回波信号至第Nr个干扰低的雷达探测回波信号,记为Nr个干扰低的雷达探测回波信号,g'=1,2,…,Nr,Nt>Nr。
步骤6,根据Nr个干扰低的雷达探测回波信号,从而得到对方目标位置,所述对方目标位置为一种射频掩护信号下的阵列雷达信号发射接收结果。
具体地,在雷达信号接收天线网络中,利用信号融合算法对Nr个干扰低的雷达探测回波信号进行融合积累,以增强检测性能,进而得到融合积累后的雷达探测回波信号;具体处理过程如图2a和图2b所示。
然后对融合积累后的雷达探测回波信号进行脉冲压缩处理,得到脉冲压缩处理后的雷达探测回波信号,其中脉冲压缩处理后的雷达探测回波信号主瓣位置,即为对方目标位置,所述对方目标位置为一种射频掩护信号下的阵列雷达信号发射接收结果。
通过仿真图对仿真的参数及效果本发明进行验证说明。
图3和图4中仿真的参数为:干扰信号为线性调频信号,其带宽为10MHz,时宽为12.8μs,采样率为20MHz,优化的雷达探测信号为相位编码信号,其阵元个数为1,其每一个阵元的码元长度为256;图3是调频信号与编码信号的互相关,其互相关旁瓣为-26.16dB;图4是编码信号的自相关,其自相关旁瓣为-26.16dB;可以看出此时调频信号与1个相位编码信号的互相关、编码信号的自相关的旁瓣较低,调频信号与编码信号的正交性较好。
图5和图6中仿真的参数为:调频信号的带宽为10MHz,时宽为12.8μs,采样率为20MHz,优化的相位编码信号的个数为3,其每一个的码元长度为256。图5是调频信号与编码信号的互相关,其互相关旁瓣为-24.74dB;图6是编码信号的自相关,其自相关旁瓣为-24.74dB;可以看出此时调频信号与3个相位编码信号的互相关、编码信号的自相关的旁瓣及编码信号之间的互相关的旁瓣比与调频信号与1个相位编码信号的情况下的旁瓣上升了1.42dB。
图7和图8中仿真的参数为:调频信号的带宽为10MHz,时宽为12.8μs,采样率为20MHz,优化的相位编码信号的个数为3,其每一个的码元长度为256。通过优化得到长度为256的3组编码信号后,改变调频信号的采样率为40MHz、80MHz,对编码信号进行相应采样,得到调频信号与编码信号的互相关。通过对图5、图7、图8进行对比,可以看出采样率越大,互相关旁瓣的波动就越大。
图9和图10中仿真的参数为:调频信号的带宽为10MHz,时宽为64μs,采样率为20MHz,优化的相位编码信号的个数为3,其每一个的码元长度为256。此时调频信号长度是编码信号长度的5倍,图9是调频信号与编码信号的互相关,其互相关旁瓣为-23.81dB,比原来升高了0.93dB左右;图10是编码信号的自相关,其自相关旁瓣为-23.81dB,比原来升高了0.93dB左右。可以看出当调频的长度与编码信号不等长时,虽然会增强其抗干扰性能,但会使其正交性变差。
图11和图12中仿真的参数为:调频信号的带宽为10MHz,时宽为12.8μs,采样率为20MHz,调频信号进行加窗,窗函数选择为海明窗,优化的相位编码信号的个数为3,其每一个的码元长度为256;图11是调频信号与编码信号的互相关,其互相关旁瓣为-25.58dB,比原来降低了0.84dB左右,但出现了类似于窗函数过渡带的下降区域;图12是编码信号的自相关,其互相关旁瓣为-25.58dB,比原来降低了0.84dB左右。可以看出加窗可以进一步降低其自相关和互相关旁瓣。
综上所述,仿真实验验证了本发明的正确性,有效性和可靠性。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围;这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。
Claims (6)
1.一种射频掩护信号下的阵列雷达信号发射接收方法,其特征在于,包括以下步骤:
步骤1,确定我方和对方,我方包括我方雷达,对方包括对方雷达和对方目标,对方雷达置于对方目标中;我方雷达发射信号,我方雷达发射的信号包括探测信号和干扰信号;
步骤2,确定我方雷达发射信号中的干扰信号;
步骤3,根据我方雷达发射信号中的探测信号和干扰信号,构建目标函数;
步骤4,对所述目标函数进行求解,得到我方雷达发射信号中的最优探测信号矩阵;
步骤5,设定雷达信号发射天线网络和雷达信号接收天线网络,然后将我方雷达发射信号中的最优探测信号矩阵和我方雷达发射信号中的干扰信号在所述雷达信号发射天线网络同时发射,经对方目标反射后得到干扰低的雷达探测回波信号;
步骤6,雷达信号接收天线网络对所述干扰低的雷达探测回波信号进行接收,并对接收到的干扰低的雷达探测回波信号进行融合积累后,得到融合积累后的雷达探测回波信号,对融合积累后的雷达探测回波信号进行脉冲压缩处理后,从而得到对方目标位置,所述对方目标位置为一种射频掩护信号下的阵列雷达信号发射接收结果;
在步骤2中,所述我方雷达发射信号中的干扰信号,其确定过程为:
首先将我方雷达发射信号中的干扰信号类型确定为线性调频信号,并得到我方雷达发射信号中的连续干扰信号sjam(t),其表达式为:
其中,rect()表示矩形窗函数,T表示我方雷达发射信号中的连续干扰信号时宽,exp()表示指数函数,j表示虚数单位,π表示圆周率,μ表示我方雷达发射信号中的连续干扰信号调频率,μ=B/T,B表示我方雷达发射信号中的连续干扰信号带宽,t表示时间变量;
然后对我方雷达发射信号中的连续干扰信号sjam(t)进行采样点数为Nt的采样后得到的离散信号,记为我方雷达发射信号中的干扰信号Sjam;Nt≥1;
步骤3的子步骤为:
3a)设定我方雷达发射信号中的探测信号包括Na个角频率,将所述Na个角频率分别进行归一化后得到Na个归一化角频率,记为我方雷达发射探测信号时的归一化角频率方向集合fa,fam表示我方雷达发射探测信号时的第m个归一化角频率,m=1,2,…,Na,Na表示我方雷达发射探测信号时的归一化角频率方向集合fa中包括的归一化角频率总个数;
根据我方雷达发射探测信号时的归一化角频率方向集合fa,得到我方雷达发射信号中探测信号的自相关峰值旁瓣电平PSLra:
其中,max表示取最大值操作,k=1,2,…,Ns-1,Ns表示我方雷达发射信号中每个探测信号的码元长度,∈表示属于,|·|表示取模值操作,at(fam)表示我方雷达发射信号中探测信号在fam方向形成的导向矢量,
at(fam)=[1,exp(j2πfam),…,exp(j2π(Nt-1)fam)]T,Nt表示我方雷达发射信号中的探测信号总个数,()H表示取共轭转置操作,()T表示取转置操作,S表示我方雷达发射信号中的探测信号,si表示我方雷达发射信号中的第i个探测信号,表示(Ns-k)×k维全零矩阵,0k×k表示k×k维全零矩阵,表示(Ns-k)×(Ns-k)维全1矩阵,表示k×(Ns-k)维全零矩阵;
3b)将我方雷达发射信号中的探测信号发射后经由对方目标反射回来、并由我方雷达接收到的信号,记为我方雷达接收的探测信号,所述我方雷达接收的探测信号包括Na个角频率,将所述Na个角频率分别归一化后得到我方雷达接收探测信号时的归一化角频率集合f′a,
然后计算得到我方雷达发射信号中的探测信号与我方雷达接收的探测信号之间的互相关距离旁瓣电平PSLrc:
其中,p=-(Ns-1),-(Ns-2),…,0,1,2,…,Ns-1,n=1,2,…,Nb,Nb表示我方雷达发射信号中的探测信号的归一化角频率组合集合中包括的归一化角频率组合总个数,at(fam)表示我方雷达发射信号中的探测信号在fam方向形成的导向矢量,at(f′am)表示我方雷达发射信号中的探测信号在f′am方向形成的导向矢量,ar(fam)表示我方雷达接收探测信号时在fam方向形成的导向矢量,ar(f′am)表示我方雷达接收探测信号时在f′am方向形成的导向矢量,
at(f′am)=[1,exp(j2πf′am),…,exp(j2π(Nt-1)f′am)]T,
ar(f′am)=[1,exp(j2πf′am),…,exp(j2π(Nr-1)f′am)]T,
3c)设定我方雷达发射信号中的干扰信号包括Na个角频率,将所述Na个角频率分别归一化后得到我方雷达发射干扰信号时的Na个归一化角频率fb,
然后计算得到我方雷达发射信号中的探测信号与干扰信号在相同方向的互相关旁瓣电平PSLjam1:
其中,PSLjam1表示我方雷达发射信号中的探测信号与干扰信号在相同方向的互相关旁瓣电平,p=-(Ns-1),-(Ns-2),…,0,1,2,…,Ns-1,Ns表示我方雷达发射信号中每个探测信号的码元长度,Sjam表示我方雷达发射信号中的干扰信号,S表示我方雷达发射信号中的探测信号,at(fbm)表示我方雷达发射信号中的探测信号在fbm方向形成的导向矢量,
at(fbm)=[1,exp(j2πfbm),…,exp(j2π(Nt-1)fbm)]T;
3d)将我方雷达发射信号中的干扰信号发射后经由对方目标反射回来、并由我方雷达接收到的信号,记为我方雷达接收的干扰信号,我方雷达接收的干扰信号包括Na个角频率,将所述Na个角频率分别归一化后得到我方雷达接收干扰信号时的归一化角频率集合f′b,
然后计算得到我方雷达发射信号中的探测信号与接收的干扰信号在不同方向的互相关旁瓣电平PSLjam2:
其中,at(f′bm)表示我方雷达发射信号中的探测信号在f′bm方向形成的导向矢量,ar(f′bm)表示我方雷达接收干扰信号时在f′bm方向形成的导向矢量,
at(f′bm)=[1,exp(j2πf′bm),…,exp(j2π(Nt-1)f′bm)]T,
ar(f′bm)=[1,exp(j2πf′bm),…,exp(j2π(Nr-1)f′bm)]T;
3e)构建目标函数:
其中,PSL表示我方雷达发射信号的峰值旁瓣电平值,min表示取最小值操作,P表示我方雷达发射信号中的探测信号相位矩阵,|| ||p表示取p范数操作,α、β、η分别表示正的加权系数。
2.如权利要求1所述的一种射频掩护信号下的阵列雷达信号发射接收方法,其特征在于,在步骤4中,所述我方雷达发射信号中的最优探测信号矩阵,其得到过程为:
4a)初始化:设置我方雷达发射信号中的探测信号总个数Nt以及我方雷达发射信号中每个探测信号的码元长度Ns;设置在求解时最为匹配的循环次数N;设置最高旁瓣标志为无穷大;令l表示第l次循环,l的初始值为1,l=1,2,…,N,N为大于1的正整数;
4b)确定第l次循环后我方雷达发射信号中的探测信号初始化相位矩阵Pl,即第l次循环后我方雷达发射信号中的探测信号初始化相位矩阵Pl中的每个元素均设置成一个随机相位值,每个相位值的取值范围均为0 2π;
4c)调用SQP算法优化第l次循环后我方雷达发射信号中的探测信号初始化相位矩阵Pl后,得到第k次优化后我方雷达发射信号中的探测信号相位矩阵进而计算得到第l次优化后我方雷达发射信号中的探测信号矩阵循环次数与优化次数相等且一一对应;
4d)根据第l次优化后我方雷达发射信号中的探测信号矩阵利用目标函数计算得到第l次循环后我方雷达发射信号中的探测信号峰值旁瓣电平值PSLl,比较第l次循环后我方雷达发射信号中的探测信号峰值旁瓣电平值PSLl与最高旁瓣标志的大小,如果第l次循环后我方雷达发射信号中的探测信号峰值旁瓣电平值PSLl小于最高旁瓣标志,则保存此次优化的结果,即将第k次优化后我方雷达发射信号中的探测信号相位矩阵作为第g个我方雷达发射信号中的探测信号优化相位矩阵Pg,将第l次优化后我方雷达发射信号中的探测信号矩阵作为第g个我方雷达发射信号中的优化探测信号矩阵g的初始值为1,并令g的值加1;然后令l的值加1,返回4b);
如果第l次循环后我方雷达发射信号中的探测信号峰值旁瓣电平值PSLl大于或等于最高旁瓣标志,则忽略此次优化的结果,然后令l的值加1,并返回4b);
直到N次循环结束后,此时得到g个我方雷达发射信号中的探测信号优化相位矩阵,[P1,P2,…,Pg],以及g个我方雷达发射信号中的优化探测信号矩阵然后根据目标函数分别计算第1个我方雷达发射信号中的探测信号峰值旁瓣电平值PSL1至第g个我方雷达发射信号中的探测信号峰值旁瓣电平值PSLg,记为我方雷达发射信号中探测信号的g个峰值旁瓣电平值,并将其中最小的峰值旁瓣电平值所对应的优化相位矩阵,记为我方雷达发射信号中的探测信号最优相位矩阵然后根据我方雷达发射信号中的探测信号最优相位矩阵计算得到我方雷达发射信号中的最优探测信号矩阵
3.如权利要求2所述的一种射频掩护信号下的阵列雷达信号发射接收方法,其特征在于,在步骤5中,Nr个干扰低的所述雷达探测回波信号,其得到过程为:
5.1确定我方雷达包括雷达信号发射天线网络和雷达信号接收天线网络,雷达信号发射天线网络包括A个发射天线,雷达信号接收天线网络是由Nr个相同的发射天线组成的天线阵列;
其中,A=Nt时,雷达信号发射天线网络包括Nt个发射天线;A=2Nt时,雷达信号发射天线网络包括2Nt个发射天线,所述2Nt个发射天线分为Nt个探测信号天线和Nt个干扰信号天线,进而得到Nt个发射天线组,每个发射天线组包括1个探测信号天线和1个干扰信号天线,且每个发射天线组分别对应我方雷达发射信号中的最优探测信号矩阵中1个最优探测信号;Nt>Nr;
5.2在第i个发射天线或在第i个发射天线组中,我方雷达发射信号中的最优探测信号矩阵中第i个最优探测信号通过调相处理后,得到第i个最优探测信号波束,第i个最优探测信号波束通过饱和放大处理后,得到第i个经过放大的最优探测信号波束;我方雷达发射信号中的干扰信号Sjam通过调相处理后,得到干扰信号波束,干扰信号波束通过饱和放大处理后,得到经过放大的干扰信号波束;其中,i=1,2,…,Nt;
5.3根据第i个经过放大的最优探测信号波束和经过放大的干扰信号波束,得到第i个最终的雷达发射信号;
5.4如果我方雷达发射第i个最终的探测信号探测到对方目标后被对方目标反射,则将反射后得到的信号记为第g'个目标回波信号,雷达信号接收天线网络对第g'个目标回波信号进行接收,g'的初始值为1,并令g'的值加1;如果我方雷达发射第i个最终的探测信号没有探测到对方目标,则忽略第i个最终的探测信号;
5.5对第g'个目标回波信号经过滤波处理后,得到第g'个干扰低的雷达探测回波波束,对第g'个干扰低的雷达探测回波波束通过调相处理后得到第g'个干扰低的雷达探测回波信号;
5.6令i的值分别取1至Nt,重复执行5.2至5.5,进而分别得到第1个干扰低的雷达探测回波信号至第Nr个干扰低的雷达探测回波信号,记为Nr个干扰低的雷达探测回波信号,g'=1,2,…,Nr,Nt>Nr。
4.如权利要求3所述的一种射频掩护信号下的阵列雷达信号发射接收方法,其特征在于,在5.3中,所述第i个最终的雷达发射信号,其得到过程为:
当A=Nt时:第i个经过放大的最优探测信号波束和经过放大的干扰信号波束通过功率合成网络,得到第i个合成后的雷达发射信号,然后利用发射波束置零的方法得到第i个发射天线的加权系数ωi;将第i个合成后的雷达发射信号乘以第i个发射天线的加权系数ωi后得到第i个最终的雷达发射信号,所述第i个最终的雷达发射信号经由第i个发射天线发射,以探测对方目标以及干扰对方目标中的对方雷达。
5.如权利要求3所述的一种射频掩护信号下的阵列雷达信号发射接收方法,其特征在于,在5.3中,所述第i个最终的雷达发射信号,其得到过程还包括:
当A=2Nt时:对第i个发射天线组利用发射波束置零的方法分别得到第i个探测信号天线的加权系数和第i个干扰信号天线的加权系数将第i个经过放大的最优探测信号波束乘以第i个探测信号天线的加权系数后得到第i个最终雷达探测信号,所述第i个最终雷达探测信号经由第i个发射天线组发射,以探测对方雷达;
6.如权利要求3所述的一种射频掩护信号下的阵列雷达信号发射接收方法,其特征在于,在步骤6中,所述对方目标位置,其得到过程为:
在雷达信号接收天线网络中,利用信号融合算法对Nr个干扰低的雷达探测回波信号进行融合积累,得到融合积累后的雷达探测回波信号;
然后对融合积累后的雷达探测回波信号进行脉冲压缩处理,得到脉冲压缩处理后的雷达探测回波信号,其中脉冲压缩处理后的雷达探测回波信号主瓣位置,即为对方目标位置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810650437.8A CN109061576B (zh) | 2018-06-22 | 2018-06-22 | 一种射频掩护信号下的阵列雷达信号发射接收方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810650437.8A CN109061576B (zh) | 2018-06-22 | 2018-06-22 | 一种射频掩护信号下的阵列雷达信号发射接收方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109061576A CN109061576A (zh) | 2018-12-21 |
CN109061576B true CN109061576B (zh) | 2022-11-04 |
Family
ID=64820758
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810650437.8A Active CN109061576B (zh) | 2018-06-22 | 2018-06-22 | 一种射频掩护信号下的阵列雷达信号发射接收方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109061576B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111781566A (zh) * | 2020-07-06 | 2020-10-16 | 上海志良电子科技有限公司 | 雷达抗有源干扰方法及设备 |
CN111880161B (zh) * | 2020-08-03 | 2021-06-01 | 湖南国天电子科技有限公司 | 一种便携式激光雷达测量方法及装置 |
CN112180331A (zh) * | 2020-09-29 | 2021-01-05 | 中国船舶重工集团公司第七二四研究所 | 一种自适应射频掩护脉冲频点策略调度方法 |
CN113311857A (zh) * | 2021-04-29 | 2021-08-27 | 重庆交通大学 | 一种基于无人机的环境感知与避障系统及方法 |
CN115639526A (zh) * | 2021-07-19 | 2023-01-24 | 湖南迈克森伟电子科技有限公司 | 用于探测雷达的无线电抗干扰方法及其抗干扰系统 |
CN115201807A (zh) * | 2022-07-08 | 2022-10-18 | 中山大学 | 一种高精度微波测距系统、装置与方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7403153B2 (en) * | 2004-12-15 | 2008-07-22 | Valeo Raytheon Systems, Inc. | System and method for reducing a radar interference signal |
CN104076334B (zh) * | 2014-07-08 | 2016-07-06 | 西安电子科技大学 | 一种设计mimo雷达波形和发射天线阵列的方法 |
CN106093877B (zh) * | 2016-07-19 | 2018-09-21 | 西安电子科技大学 | 正交宽主瓣相位编码信号和失配滤波器联合优化方法 |
CN106093876B (zh) * | 2016-07-19 | 2018-11-09 | 西安电子科技大学 | 分布式mimo雷达正交宽主瓣相位编码信号设计方法 |
CN107015210B (zh) * | 2017-03-16 | 2020-01-10 | 电子科技大学 | 一种基于相似性约束的雷达慢时间域恒模信号设计方法 |
-
2018
- 2018-06-22 CN CN201810650437.8A patent/CN109061576B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN109061576A (zh) | 2018-12-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109061576B (zh) | 一种射频掩护信号下的阵列雷达信号发射接收方法 | |
CN110520750B (zh) | 用于使用mimo雷达对物体进行无线检测的方法和系统 | |
CN104297734B (zh) | 基于频率分集阵列的mimo雷达的欺骗式干扰抑制方法 | |
CN103728614B (zh) | 基于机扫米波雷达的改进单脉冲测角方法 | |
CN106093931B (zh) | 基于数字阵列天线的通信雷达一体化收发方法 | |
Ma et al. | Accurate indoor ranging by broadband harmonic generation in passive NLTL backscatter tags | |
CN109901149A (zh) | 一种基于fda-mimo雷达的目标参数估计方法 | |
JP2008008887A (ja) | 信号到来角決定方法及び信号到来角決定装置、並びに信号の到来角を決定するためのシステム | |
CN104155647B (zh) | 基于机扫米波阵列雷达的估计目标方位角的方法 | |
CN109283497B (zh) | 双基地fda-mimo雷达距离欺骗干扰识别方法 | |
CN108845309A (zh) | 一种基于主动虚拟时反方法的水下目标doa估计方法 | |
CN105137410A (zh) | 基于ofdm的高分辨率雷达通信一体化的波形优化方法 | |
CN108132461B (zh) | 抑制调频连续波着陆雷达直流泄露的方法 | |
CN106066473A (zh) | 正交波形下mimo雷达目标的多波束联合测角和点迹融合方法 | |
CN105974390B (zh) | 基于多普勒信息的机扫米波雷达质量中心测角方法 | |
CN106772337A (zh) | 一种基于music算法的频率分集阵列两层介质目标定位方法 | |
CN109597041A (zh) | 基于相干fda的分段线性调频波形设计方法 | |
CN101471734B (zh) | 一种多发多收声定位网络系统及其定位方法 | |
CN113406573B (zh) | 基于多混频器的慢时间fda雷达信号处理方法、装置、介质及雷达系统 | |
CN113156380B (zh) | 基于收发失配处理的主瓣噪声压制干扰自适应抑制方法 | |
CN112881984B (zh) | 一种雷达信号抗干扰处理方法、装置及存储介质 | |
CN109188373A (zh) | 基于子阵阻塞矩阵预处理的抗主瓣干扰方法 | |
CN111505600B (zh) | 基于stpc的fda-mimo雷达信号处理方法、装置及介质 | |
CN109407098A (zh) | 一种处理正交脉冲串信号的并行处理算法的实现方法 | |
Fackelmeier et al. | Narrowband frequency scanning array antenna at 5.8 GHz for short range imaging |
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 |