CN108333586A - 基于双接收机垂直交叉扫描的无源雷达测向方法及系统 - Google Patents
基于双接收机垂直交叉扫描的无源雷达测向方法及系统 Download PDFInfo
- Publication number
- CN108333586A CN108333586A CN201810135996.5A CN201810135996A CN108333586A CN 108333586 A CN108333586 A CN 108333586A CN 201810135996 A CN201810135996 A CN 201810135996A CN 108333586 A CN108333586 A CN 108333586A
- Authority
- CN
- China
- Prior art keywords
- receiver
- movement locus
- signal
- matrix
- finding
- 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
Links
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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开一种双接收机垂直交叉扫描的无源雷达测向方法及系统,所述方法包括:获取第一接收机的第一回波信号;获取第二接收机的第二回波信号,根据第一接收机的运动轨迹函数构建第一信号识别矩阵;根据第一回波信号和第一信号识别矩阵确定第一测向预处理向量;根据第二接收机的运动轨迹函数构建第二信号识别矩阵;根据第二回波信号和第二信号识别矩阵确定第二测向预处理向量;根据第一测向预处理向量和第二测向预处理向量确定待测目标的方向。本发明通过双接收机椭圆垂直交叉扫描能够减小系统占地面积,对旁瓣进行有效抑制,并且优化了系统分辨率,提升了系统的测向性能。
Description
技术领域
本发明涉及无源雷达测向领域,特别是涉及一种基于双接收机垂直交叉扫描的无源雷达测向方法及系统。
背景技术
目前的固定目标雷达测向方法主要分主动雷达测向和无源雷达测向两种。主动测向雷达包括实孔径雷达或者阵列天线测向,其方位分辨率均为,其中为信号波长,而为雷达孔径或者阵列天线长度。为了提升测向精度,主动测向雷达系统一般采用大孔径天线或者长阵列天线,所以系统造价高昂,且占地面积较大。另外由于主动进行电磁照射,在战场上其生存能力堪忧。无源雷达测向方法方法中,以一种基于民用窄带外辐射源信号、接收机围着其转动中心做匀速圆周运动的雷达测向方法为代表,该方法的方位分辨率为,其中为接收机转动半径。该测向系统的结构简单,系统分辨能力较好。但是该系统测向结果的最大旁瓣高达-7.9dB,会出现对多目标测向时强旁瓣可能淹没弱主瓣而导致目标丢失的情况。
发明内容
本发明的目的是提供一种基于双接收机垂直交叉扫描的无源雷达测向方法及系统,提升了目标测向的分辨性能。
为实现上述目的,本发明提供了如下方案:
一种基于双接收机垂直交叉扫描的无源雷达测向方法,所述无源雷达测向方法包括:
获取第一接收机的第一回波信号;
获取第二接收机的第二回波信号,所述第二接收机和所述第一接收机围绕同一中心作椭圆运动,所述第一接收机的运动轨迹所在平面与所述第二接收机的运动轨迹所在平面互相垂直,所述第一接收机的运动轨迹的长轴与所述第二接收机的运动轨迹的短轴的方向相同;
根据所述第一接收机的运动轨迹函数构建第一信号识别矩阵;
根据所述第一回波信号和所述第一信号识别矩阵确定第一测向预处理向量;
根据所述第二接收机的运动轨迹函数构建第二信号识别矩阵;
根据所述第二回波信号和所述第二信号识别矩阵确定第二测向预处理向量;
根据所述第一测向预处理向量和所述第二测向预处理向量确定待测目标的方向。
可选的,所述根据所述第一接收机的运动轨迹函数构建第一信号识别矩阵,具体包括:
根据所述第一接收机的运动轨迹确定待测目标相对所述第一接收机的第一斜距历程;
根据所述第一斜距历程得到第一离散回波信号矩阵;
根据所述第一离散回波信号矩阵构造第一信号识别矩阵。
可选的,所述根据所述第二接收机的运动轨迹函数构建第二信号识别矩阵,具体包括:
根据所述第二接收机的运动轨迹确定待测目标相对所述第二接收机的第二斜距历程;
根据所述第二斜距历程得到第二离散回波信号矩阵;
根据所述第二离散回波信号矩阵构造第二信号识别矩阵。
可选的,所述第一信号识别矩阵的形式为:
其中,G1(m1,n)为第一信号识别矩阵,m1=1,2,…,M1,M1为第一接收机在成像过程中的采样次数,n=1,2,…,N,N为遍历次数,a1为所述第一接收机运动轨迹的长半轴,b1为所述第一接收机运动轨迹的短半轴,λ是系统所采用外辐射源信号波长,θ1(m1)为第一接收机在椭圆转动轨迹上的角度位置,Δα为遍历步长,N=2π/Δα,exp{j……}为复数的指数形式。
可选的,所述第二信号识别矩阵的形式为:
其中,G2(m2,n)为第二信号识别矩阵,m2=1,2,…,M2,M2为第二接收机在成像过程中的采样次数,a2为所述第二接收机运动轨迹的长半轴,b2为所述第二接收机运动轨迹的短半轴,λ是系统所采用外辐射源信号波长,θ2(m2)为第二接收机在椭圆转动轨迹上的角度位置。
可选的,所述根据所述第一离散回波信号矩阵和所述第一信号识别矩阵确定第一测向预处理向量,具体包括:
根据公式计算第一测向预处理向量;
其中,F1(n)为第一测向预处理向量,g1(m1)为第一离散回波信号矩阵,G1(m1,n)为第一信号识别矩阵,m1=1,2,…,M1,M1为第一接收机在成像过程中的采样次数,n=1,2,…,N,N为遍历次数。
可选的,所述根据所述第二离散回波信号矩阵和所述第二信号识别矩阵确定第二测向预处理向量,具体包括:
根据公式计算第二测向预处理向量;
其中,F2(n)为第二测向预处理向量,g2(m2)为第二离散回波信号矩阵,G2(m2,n)为第二信号识别矩阵,m2=1,2,…,M2,M2为第二接收机在成像过程中的采样次数。
可选的,所述根据所述第一测向预处理向量和所述第二测向预处理向量确定待测目标的方向,具体包括:
根据公式F(n)=F1(n)*F2(n)计算待测目标的位置函数;其中,F(n)为待测目标的位置函数,F1(n)为第一测向预处理向量,F2(n)为第二测向预处理向量,*为Schur积,n=1,2,…,N,N为遍历次数。
确定所述位置函数的峰值;
根据所述位置函数的峰值对应的n值,确定待测目标的方位角为nΔα,其中,Δα为遍历步长。
可选的,a1=b2,a2=b1,b1=0.63a1,其中,a1为所述第一接收机运动轨迹的长半轴,a2为所述第二接收机运动轨迹的长半轴,b1为所述第一接收机运动轨迹的短半轴,b2为所述第二接收机运动轨迹的短半轴。
一种基于双接收机垂直交叉扫描的无源雷达测向系统,所述无源雷达测向系统包括:
第一回波信号获取模块,用于获取第一接收机的第一回波信号;
第二回波信号获取模块,用于获取第二接收机的第二回波信号,所述第二接收机和所述第一接收机围绕同一中心作椭圆运动,所述第一接收机的运动轨迹所在平面与所述第二接收机的运动轨迹所在平面互相垂直,所述第一接收机的运动轨迹的长轴与所述第二接收机的运动轨迹的短轴的方向相同;
第一信号识别矩阵构建模块,用于根据所述第一接收机的运动轨迹函数构建第一信号识别矩阵;
第一测向预处理向量确定模块,用于根据所述第一回波信号和所述第一信号识别矩阵确定第一测向预处理向量;
第二信号识别矩阵构建模块,用于根据所述第二接收机的运动轨迹函数构建第二信号识别矩阵;
第二测向预处理向量确定模块,用于根据所述第二回波信号和所述第二信号识别矩阵确定第二测向预处理向量;
待测目标的方向确定模块,用于根据所述第一测向预处理向量和所述第二测向预处理向量确定待测目标的方向。
根据本发明提供的具体实施例,本发明公开了以下技术效果:
本发明首次提出一种双接收机椭圆垂直交叉扫描的无源雷达测向方法,通过双接收机椭圆扫描能够减小系统占地面积,对旁瓣进行有效抑制,并且优化了系统分辨率,提升了系统的测向性能。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明基于双接收机垂直交叉扫描的无源雷达测向方法的流程图;
图2为本发明测向系统模型示意图;
图3为各方位3dB分辨率;
图4为本发明对15个随机目标的测向结果;
图5为圆形扫描对15个随机目标的测向结果;
图6为不同频率外辐射源对测向结果的影响;
图7为a取不同值时的测向结果;
图8为本发明基于双接收机垂直交叉扫描的无源雷达测向系统的结构图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的目的是提供一种基于双接收机垂直交叉扫描的无源雷达测向方法及系统,提升了目标测向的分辨性能。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1为本发明基于双接收机垂直交叉扫描的无源雷达测向方法的流程图,如图1所示,所述无源雷达测向方法包括:
步骤101:获取第一接收机的第一回波信号。
步骤102:获取第二接收机的第二回波信号,所述第二接收机和所述第一接收机围绕同一中心作椭圆运动,所述第一接收机的运动轨迹所在平面与所述第二接收机的运动轨迹所在平面互相垂直,所述第一接收机的运动轨迹的长轴与所述第二接收机的运动轨迹的短轴的方向相同;a1=b2,a2=b1,b1=0.63a1,其中,a1为所述第一接收机运动轨迹的长半轴,a2为所述第二接收机运动轨迹的长半轴,b1为所述第一接收机运动轨迹的短半轴,b2为所述第二接收机运动轨迹的短半轴。
步骤103:根据所述第一接收机的运动轨迹函数构建第一信号识别矩阵;具体包括:
步骤1031:根据所述第一接收机的运动轨迹确定待测目标相对所述第一接收机的第一斜距历程。
步骤1032:根据所述第一斜距历程得到第一离散回波信号矩阵。
步骤1033:根据所述第一离散回波信号矩阵构造第一信号识别矩阵。所述第一信号识别矩阵的形式为:
其中,G1(m1,n)为第一信号识别矩阵,m1=1,2,…,M1,M1为第一接收机在成像过程中的采样次数,n=1,2,…,N,N为遍历次数,a1为所述第一接收机运动轨迹的长半轴,b1为所述第一接收机运动轨迹的短半轴,λ是系统所采用外辐射源信号波长,θ1(m1)为第一接收机在椭圆转动轨迹上的角度位置,Δα为遍历步长,N=2π/Δα,exp{j……}为复数的指数形式。
步骤104:根据所述第一回波信号和所述第一信号识别矩阵确定第一测向预处理向量;具体包括:
根据公式计算第一测向预处理向量;
其中,F1(n)为第一测向预处理向量,g1(m1)为第一离散回波信号矩阵,G1(m1,n)为第一信号识别矩阵,m1=1,2,…,M1,M1为第一接收机在成像过程中的采样次数,n=1,2,…,N,N为遍历次数。
步骤105:根据所述第二接收机的运动轨迹函数构建第二信号识别矩阵;具体包括:
步骤1051:根据所述第二接收机的运动轨迹确定待测目标相对所述第二接收机的第一斜距历程。
步骤1052:根据所述第二斜距历程得到第二离散回波信号矩阵。
步骤1053:根据所述第二离散回波信号矩阵构造第二信号识别矩阵。所述第二信号识别矩阵的形式为:
其中,G2(m2,n)为第二信号识别矩阵,m2=1,2,…,M2,M2为第二接收机在成像过程中的采样次数,a2为所述第二接收机运动轨迹的长半轴,b2为所述第二接收机运动轨迹的短半轴,λ是系统所采用外辐射源信号波长,θ2(m2)为第二接收机在椭圆转动轨迹上的角度位置。
步骤106:根据所述第二回波信号和所述第二信号识别矩阵确定第二测向预处理向量,具体包括:
根据公式计算第二测向预处理向量;
其中,F2(n)为第二测向预处理向量,g2(m2)为第二离散回波信号矩阵,G2(m2,n)为第二信号识别矩阵,m2=1,2,…,M2,M2为第二接收机在成像过程中的采样次数。
步骤107:根据所述第一测向预处理向量和所述第二测向预处理向量确定待测目标的方向,具体包括:
步骤1071:根据公式F(n)=F1(n)*F2(n)计算待测目标的位置函数;其中,F(n)为待测目标的位置函数,F1(n)为第一测向预处理向量,F2(n)为第二测向预处理向量,*为Schur积,n=1,2,…,N,N为遍历次数。
步骤1072:确定所述位置函数的峰值;
步骤1073:根据所述位置函数的峰值对应的n值,确定待测目标的方位角为nΔα,其中,Δα为遍历步长。
本发明首次提出一种双接收机椭圆垂直交叉扫描的无源雷达测向方法,通过双椭圆扫描模式,减小了所需占地面积,极大降低了系统旁瓣,并降低了系统分辨率,使得系统测向性能有了显著提升。
以下是本发明的理论推理过程:
步骤1:建立测向系统模型。
图2为本发明测向系统模型示意图。本测向系统为无源测向系统,在接收端采用了双接收机模式,两个接收机均绕着同一转动中心转动,其转动轨迹均为椭圆。其中接收机1的长半轴为a,短半轴为b,且有b=0.63a。以接收机1的长轴方向为X轴,以短轴方向为Y轴,以转动中心为坐标原点建立直角坐标系。接收机2的长半轴与短半轴与接收机1相同,但是接收机2的椭圆轨迹与接收机1的椭圆轨迹垂直交叉,也即将接收机1的运动轨迹整体逆时针转动90°,即可得到接收机2的运动轨迹,所以接收机2的长轴在上述直角坐标系的Y轴上,而短轴在X轴上。在t时刻有两接收机的椭圆角度参数为θi(t)=ωit+θi0,i=1,2,其中ωi为转动角速度,θi0为初始角度参数,下标i=1时表示接收机1,i=2时表示接收机2,下面也采用此标记方法。在t时刻有两接收机的位置为(xir(t),yir(t)),则有:
和
其中,a1=a,b1=b,a2=b,b2=a。
步骤2:推导监测区域内固定目标相对于两个接收机的斜距历程:
Ri(t)=RT+R0-aicosαcos(ωit+θi0)-bisinαsin(ωit+θi0);
其中,RT是目标与发射机之间的距离,R0是目标到坐标原点之间的距离,α是在监测区域内所选任一固定目标的方位角。
斜距历程推导过程为如下:
假设发射机位于(xt,yt),在观测区域存在一固定目标,其直角坐标和极坐标分别为(x0,y0)和(R0,α),则有:
上式中的近似处理基于目标到转动中心的距离远大于接收机转动尺寸而展开。
步骤3:推导解调后回波信号的离散采样形式:
其中,σ是该固定目标散射系数,λ是系统所采用外辐射源信号波长,θi(mi)=ωimiΔti+θi0,Δti是接收机采样的时间步长,接收机在成像过程中的采样次数Mi,
所用民用信号一般为窄带信号,假设其为:f(t)=exp{j2πft}。
则目标回波信号为:
其中,C指光速。所以解调并离散采样后,其回波信号为:
其中,σ是该固定目标散射系数,λ是系统所采用外辐射源信号波长,θi(mi)=ωimiΔti+θi0,Δti是接收机采样的时间步长,接收机在成像过程中的采样次数Mi,
步骤4:根据回波信号形式构造信号识别矩阵:
其中,Δα为遍历步长,共遍历N次,N=2π/Δα。为了便于后续处理,两个接收机的信号识别矩阵中的Δα和N要一致。
步骤5:求测向预处理向量:
其中,J0(·)为零阶贝塞尔函数。
Fi(n)的推导过程为:
其中, J0(·)为零阶贝塞尔函数。
在3dB波束内,变化很小,所以有因此所以在3dB波束内有:
由上式可知,当nΔα=α,也即n=α/Δα时,Fi(n)出现峰值Mi|σ|,可以根据该峰值确定目标方向。由贝塞尔函数的性质可知,其分辨率为0.36λ/Bi(α)rad。由于Bi(α)是α的函数,所以两个接收机对目标的分辨率与目标所在方位有关。当目标位于0.5πrad或者1.5πrad(也即目标位于Y轴正向或者负向上)时,接收机1处理结果的分辨性能最好,分辨率为0.36λ/arad,而接收机2处理结果的分辨性能最差,分辨率为0.36λ/brad。当目标位于πrad或者2πrad(也即目标位于X轴正向或者负向上)时,接收机2处理结果的分辨性能最好,分辨率为0.36λ/arad,而接收机1处理结果的分辨性能最差,分辨率为0.36λ/brad。
通过求解F1(n)和F2(n)的Schur积得到F(n)。由于Fi(n)为零阶贝塞尔函数,根据贝塞尔函数性质可得,Fi(n)的第一零点和第一旁瓣与主瓣之间角度差分别为0.383λ/Bi(α)rad和0.61λ/Bi(α)rad。因此只要满足a2=0.63a1,在0.5πrad或1.5πrad附近就能够使接收机2处理结果的第一零点位置与接收机1第一旁瓣的位置相同,则通过求F1(n)和F2(n)的Schur积可以使F(n)在此处出现零点。同理,满足b1=0.63b2,则在πrad或2πrad附近就能够使接收机1处理结果的第一零点位置与接收机2第一旁瓣的位置相同,F1(n)和F2(n)的Schur积F(n)在此处出现零点。所以在建立系统模型时系统参数有a1=b2=a,b1=a2=b=0.63a。分析表明,F(n)的分辨率在整个遍历区间(0rad,2πrad]上基本为常数0.309λ/arad,波动很小,图3是根据步骤5和步骤6作出的各方位3dB分辨率,可以看出,虽然目标位于不同方位角时其分辨率会有波动,但是其波动非常小,基本可以忽略不计。当目标方位角位于0.25πrad、0.75πrad、1.25πrad、1.75πrad时,其最大旁瓣为-15.8dB,而目标方位角位于0.5πrad、πrad、1.5πrad、2πrad时,其最大旁瓣为-19.43dB,其余位置F(n)的最大旁瓣介于-15.8dB和-19.43dB之间。
本发明的效果通过以下仿真实验进一步说明。
仿真条件
外辐射源坐标为(20000m,15000m),转动过程中接收机1采样2000次,椭圆角度参数的采样步长为0.001πrad,接收机2采样1000次,椭圆角度参数的采样步长为0.002πrad。接收机1和接收机2的角度参数初始值分别为0.125πrad和0.1πrad,不失一般性,仿真中目标散射系数均为1。
仿真内容
仿真1:在上述仿真条件基础上,接收机转动参数a=20m,采用频率为300MHz的民用信号,监测区域存在有15个随机目标,分别用本发明所提方法和圆形扫描(接收机转动半径为20m)进行测向,测向结果如图4和图5所示,其中“*”表示目标的真实方位角。从仿真结果可以看出,本发明的测向结果明显要比圆形扫描测向结果的旁瓣要低,而且其分辨率也更好,所以本发明的测向结果更为干净,其分辨效果更优。
仿真2:接收机转动参数a=20m,外辐射源频率分别为20MHz、40MHz、90MHz时,对位于(5000m,πrad)固定目标测向,其仿真结果如图6所示。从图6可以看出,系统的方位分辨性能会随着外辐射源信号频率的提高而提升。
仿真3:当外辐射源频率为50MHz时,接收机转动参数a分别取10m、20m、30m时对位于(5000m,πrad)固定目标测向,其仿真结果如图7所示。从图7可以看出,分辨性能会随着长半轴的增长而提升。
图8为本发明基于双接收机垂直交叉扫描的无源雷达测向系统的结构图。如图8所示,所述无源雷达测向系统包括:
第一回波信号获取模块801,用于获取第一接收机的第一回波信号;
第二回波信号获取模块802,用于获取第二接收机的第二回波信号,所述第二接收机和所述第一接收机围绕同一中心作椭圆运动,所述第一接收机的运动轨迹所在平面与所述第二接收机的运动轨迹所在平面互相垂直,所述第一接收机的运动轨迹的长轴与所述第二接收机的运动轨迹的短轴的方向相同;
第一信号识别矩阵构建模块803,用于根据所述第一接收机的运动轨迹函数构建第一信号识别矩阵;
第一测向预处理向量确定模块804,用于根据所述第一回波信号和所述第一信号识别矩阵确定第一测向预处理向量;
第二信号识别矩阵构建模块805,用于根据所述第二接收机的运动轨迹函数构建第二信号识别矩阵;
第二测向预处理向量确定模块806,用于根据所述第二回波信号和所述第二信号识别矩阵确定第二测向预处理向量;
待测目标的方向确定模块807,用于根据所述第一测向预处理向量和所述第二测向预处理向量确定待测目标的方向。
对于实施例公开的系统而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。
Claims (10)
1.一种基于双接收机垂直交叉扫描的无源雷达测向方法,其特征在于,所述无源雷达测向方法包括:
获取第一接收机的第一回波信号;
获取第二接收机的第二回波信号,所述第二接收机和所述第一接收机围绕同一中心作椭圆运动,所述第一接收机的运动轨迹所在平面与所述第二接收机的运动轨迹所在平面互相垂直,所述第一接收机的运动轨迹的长轴与所述第二接收机的运动轨迹的短轴的方向相同;
根据所述第一接收机的运动轨迹函数构建第一信号识别矩阵;
根据所述第一回波信号和所述第一信号识别矩阵确定第一测向预处理向量;
根据所述第二接收机的运动轨迹函数构建第二信号识别矩阵;
根据所述第二回波信号和所述第二信号识别矩阵确定第二测向预处理向量;
根据所述第一测向预处理向量和所述第二测向预处理向量确定待测目标的方向。
2.根据权利要求1所述的无源雷达测向方法,其特征在于,所述根据所述第一接收机的运动轨迹函数构建第一信号识别矩阵,具体包括:
根据所述第一接收机的运动轨迹确定待测目标相对所述第一接收机的第一斜距历程;
根据所述第一斜距历程得到第一离散回波信号矩阵;
根据所述第一离散回波信号矩阵构造第一信号识别矩阵。
3.根据权利要求1所述的无源雷达测向方法,其特征在于,所述根据所述第二接收机的运动轨迹函数构建第二信号识别矩阵,具体包括:
根据所述第二接收机的运动轨迹确定待测目标相对所述第二接收机的第二斜距历程;
根据所述第二斜距历程得到第二离散回波信号矩阵;
根据所述第二离散回波信号矩阵构造第二信号识别矩阵。
4.根据权利要求2所述的无源雷达测向方法,其特征在于,所述第一信号识别矩阵的形式为:
其中,G1(m1,n)为第一信号识别矩阵,m1=1,2,…,M1,M1为第一接收机在成像过程中的采样次数,n=1,2,…,N,N为遍历次数,a1为所述第一接收机运动轨迹的长半轴,b1为所述第一接收机运动轨迹的短半轴,λ是系统所采用外辐射源信号波长,θ1(m1)为第一接收机在椭圆转动轨迹上的角度位置,Δα为遍历步长,N=2π/Δα,exp{j……}为复数的指数形式。
5.根据权利要求3所述的无源雷达测向方法,其特征在于,所述第二信号识别矩阵的形式为:
其中,G2(m2,n)为第二信号识别矩阵,m2=1,2,…,M2,M2为第二接收机在成像过程中的采样次数,a2为所述第二接收机运动轨迹的长半轴,b2为所述第二接收机运动轨迹的短半轴,λ是系统所采用外辐射源信号波长,θ2(m2)为第二接收机在椭圆转动轨迹上的角度位置。
6.根据权利要求2所述的无源雷达测向方法,其特征在于,所述根据所述第一离散回波信号矩阵和所述第一信号识别矩阵确定第一测向预处理向量,具体包括:
根据公式计算第一测向预处理向量;
其中,F1(n)为第一测向预处理向量,g1(m1)为第一离散回波信号矩阵,G1(m1,n)为第一信号识别矩阵,m1=1,2,…,M1,M1为第一接收机在成像过程中的采样次数,n=1,2,…,N,N为遍历次数。
7.根据权利要求3所述的无源雷达测向方法,其特征在于,所述根据所述第二离散回波信号矩阵和所述第二信号识别矩阵确定第二测向预处理向量,具体包括:
根据公式计算第二测向预处理向量;
其中,F2(n)为第二测向预处理向量,g2(m2)为第二离散回波信号矩阵,G2(m2,n)为第二信号识别矩阵,m2=1,2,…,M2,M2为第二接收机在成像过程中的采样次数。
8.根据权利要求1所述的无源雷达测向方法,其特征在于,所述根据所述第一测向预处理向量和所述第二测向预处理向量确定待测目标的方向,具体包括:
根据公式F(n)=F1(n)*F2(n)计算待测目标的位置函数;
其中,F(n)为待测目标的位置函数,F1(n)为第一测向预处理向量,F2(n)为第二测向预处理向量,*为Schur积,n=1,2,…,N,N为遍历次数。
确定所述位置函数的峰值;
根据所述位置函数的峰值对应的n值,确定待测目标的方位角为nΔα,其中,Δα为遍历步长。
9.根据权利要求1所述的无源雷达测向方法,其特征在于,a1=b2,a2=b1,b1=0.63a1,其中,a1为所述第一接收机运动轨迹的长半轴,a2为所述第二接收机运动轨迹的长半轴,b1为所述第一接收机运动轨迹的短半轴,b2为所述第二接收机运动轨迹的短半轴。
10.一种基于双接收机垂直交叉扫描的无源雷达测向系统,其特征在于,所述无源雷达测向系统包括:
第一回波信号获取模块,用于获取第一接收机的第一回波信号;
第二回波信号获取模块,用于获取第二接收机的第二回波信号,所述第二接收机和所述第一接收机围绕同一中心作椭圆运动,所述第一接收机的运动轨迹所在平面与所述第二接收机的运动轨迹所在平面互相垂直,所述第一接收机的运动轨迹的长轴与所述第二接收机的运动轨迹的短轴的方向相同;
第一信号识别矩阵构建模块,用于根据所述第一接收机的运动轨迹函数构建第一信号识别矩阵;
第一测向预处理向量确定模块,用于根据所述第一回波信号和所述第一信号识别矩阵确定第一测向预处理向量;
第二信号识别矩阵构建模块,用于根据所述第二接收机的运动轨迹函数构建第二信号识别矩阵;
第二测向预处理向量确定模块,用于根据所述第二回波信号和所述第二信号识别矩阵确定第二测向预处理向量;
待测目标的方向确定模块,用于根据所述第一测向预处理向量和所述第二测向预处理向量确定待测目标的方向。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810135996.5A CN108333586B (zh) | 2018-02-09 | 2018-02-09 | 基于双接收机垂直交叉扫描的无源雷达测向方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810135996.5A CN108333586B (zh) | 2018-02-09 | 2018-02-09 | 基于双接收机垂直交叉扫描的无源雷达测向方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108333586A true CN108333586A (zh) | 2018-07-27 |
CN108333586B CN108333586B (zh) | 2020-03-27 |
Family
ID=62929207
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810135996.5A Expired - Fee Related CN108333586B (zh) | 2018-02-09 | 2018-02-09 | 基于双接收机垂直交叉扫描的无源雷达测向方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108333586B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111693959A (zh) * | 2020-06-09 | 2020-09-22 | 周口师范学院 | 一种双圆周复合扫描被动雷达测向方法 |
CN111693957A (zh) * | 2020-06-09 | 2020-09-22 | 周口师范学院 | 一种圆周与简谐复合扫描被动雷达测向方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103499974A (zh) * | 2013-09-27 | 2014-01-08 | 中国人民解放军空军工程大学 | 一种双机协同被动雷达探测路径规划方法 |
US20140210661A1 (en) * | 2011-04-07 | 2014-07-31 | Mrinal Singh Balaji | Systems and methods for calibrating dual polarization radar systems |
CN105182283A (zh) * | 2015-08-17 | 2015-12-23 | 周口师范学院 | 一种无源雷达固定目标时域测向方法 |
US9482752B1 (en) * | 2013-03-14 | 2016-11-01 | Enterprise Electronics Corporation | Process for measuring circular depolarization ratios in a weather radar |
CN106249217A (zh) * | 2016-08-11 | 2016-12-21 | 周口师范学院 | 一种单发多收无源雷达随机目标的定位方法及装置 |
CN106291450A (zh) * | 2016-08-11 | 2017-01-04 | 周口师范学院 | 一种单发双收无源雷达的测向方法及装置 |
CN106353750A (zh) * | 2016-08-11 | 2017-01-25 | 周口师范学院 | 一种多发单收无源雷达的测向方法及装置 |
CN106483514A (zh) * | 2016-09-23 | 2017-03-08 | 电子科技大学 | 一种基于eemd和支持向量机的飞机运动模式识别方法 |
-
2018
- 2018-02-09 CN CN201810135996.5A patent/CN108333586B/zh not_active Expired - Fee Related
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140210661A1 (en) * | 2011-04-07 | 2014-07-31 | Mrinal Singh Balaji | Systems and methods for calibrating dual polarization radar systems |
US9482752B1 (en) * | 2013-03-14 | 2016-11-01 | Enterprise Electronics Corporation | Process for measuring circular depolarization ratios in a weather radar |
CN103499974A (zh) * | 2013-09-27 | 2014-01-08 | 中国人民解放军空军工程大学 | 一种双机协同被动雷达探测路径规划方法 |
CN103499974B (zh) * | 2013-09-27 | 2016-08-17 | 中国人民解放军空军工程大学 | 一种双机协同被动雷达探测路径规划方法 |
CN105182283A (zh) * | 2015-08-17 | 2015-12-23 | 周口师范学院 | 一种无源雷达固定目标时域测向方法 |
CN106249217A (zh) * | 2016-08-11 | 2016-12-21 | 周口师范学院 | 一种单发多收无源雷达随机目标的定位方法及装置 |
CN106291450A (zh) * | 2016-08-11 | 2017-01-04 | 周口师范学院 | 一种单发双收无源雷达的测向方法及装置 |
CN106353750A (zh) * | 2016-08-11 | 2017-01-25 | 周口师范学院 | 一种多发单收无源雷达的测向方法及装置 |
CN106483514A (zh) * | 2016-09-23 | 2017-03-08 | 电子科技大学 | 一种基于eemd和支持向量机的飞机运动模式识别方法 |
Non-Patent Citations (2)
Title |
---|
HONGWEI LI等: ""Passive coherent radar tracking algorithm based on particle filter and multiple TDOA measurements"", 《2009 2ND INTERNATIONAL CONGRESS ON IMAGE AND SIGNAL PROCESSING》 * |
王宗全等: ""一种二维无源交叉定位方法"", 《雷达科学与技术》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111693959A (zh) * | 2020-06-09 | 2020-09-22 | 周口师范学院 | 一种双圆周复合扫描被动雷达测向方法 |
CN111693957A (zh) * | 2020-06-09 | 2020-09-22 | 周口师范学院 | 一种圆周与简谐复合扫描被动雷达测向方法 |
CN111693959B (zh) * | 2020-06-09 | 2024-04-19 | 周口师范学院 | 一种双圆周复合扫描被动雷达测向方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108333586B (zh) | 2020-03-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105388465B (zh) | 基于海浪谱模型的海杂波模拟方法 | |
CN104777453B (zh) | 舰船线谱噪声源定位的波束域时频分析方法 | |
CN104035095B (zh) | 基于空时最优处理器的低空风切变风速估计方法 | |
CN108363046A (zh) | 基于双接收机椭圆扫描的无源雷达测向方法及系统 | |
CN103235297B (zh) | 基于宽带雷达观测的空间章动目标参数估计方法 | |
CN109031219A (zh) | 基于相位测距的宽带雷达弹道目标微动几何参数估计方法 | |
CN102997988B (zh) | 大型潜标矢量水听器低频声学指向性的水池测试方法 | |
CN104133213B (zh) | 一种结合rm算法与bp算法的柱面近场三维rcs成像方法 | |
CN105738891B (zh) | 一种机载数字阵列雷达对弱机动目标角度进行跟踪的方法 | |
CN108333586A (zh) | 基于双接收机垂直交叉扫描的无源雷达测向方法及系统 | |
CN108169511A (zh) | 三维空间来风的风速测量系统及方法 | |
CN109738894A (zh) | 一种雷达的成像方法 | |
CN106371083A (zh) | 雷达监测点选取方法及装置以及回波信号处理方法及装置 | |
WO2016142071A1 (de) | Ultraschallwindmesser | |
CN107942314A (zh) | 基于lasso特征提取的多普勒穿墙雷达定位方法 | |
CN108919174A (zh) | 非规则天线阵列结构的短波无线电测向系统和方法 | |
CN106249217B (zh) | 一种单发多收无源雷达随机目标的定位方法及装置 | |
CN107390170A (zh) | 基于定向电磁波及旋转角度进行三维空间定位的方法 | |
CN103064084A (zh) | 基于距离频域的解模糊方法 | |
CN104597446B (zh) | 一种星载合成孔径雷达地距分辨率表征与参数设计方法 | |
CN108562872A (zh) | 一种超短基线水下声学定位校准时异常值探测的方法 | |
CN108363045B (zh) | 基于单接收机复合扫描的无源雷达测向方法及系统 | |
CN108344983A (zh) | 一种无源雷达测向方法及系统 | |
CN108303686A (zh) | 一种外辐射源雷达测向方法及系统 | |
Piniés et al. | Fast localization of avalanche victims using sum of gaussians |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200327 Termination date: 20210209 |