法律状态公告日
法律状态信息
法律状态
2019-01-08
授权
授权
2016-12-07
实质审查的生效 IPC(主分类):G01S7/295 申请日:20160531
实质审查的生效
2016-11-09
公开
公开
技术领域
本发明涉及雷达下视尾后探测技术领域,具体涉及一种雷达下视尾后探测下主副瓣地海杂波主分布计算方法。
背景技术
地海杂波是雷达下视尾后探测和跟踪目标性能的最大制约因素。下视尾后探测空中目标,目标回波信号落入副瓣杂波区,雷达发射中重复频率波形用于探测目标。相对高重复频率波形,采用中重复频率波形能够减小副瓣杂波距离混叠程度,并进一步通过距离维切割降低副瓣杂波功率,使雷达具备下视尾后探测能力。但由于副瓣杂波在雷达处理的距离-多普勒域上并非均匀分布,而是存在强弱分布,当目标回波信号落入强副瓣杂波区域时,一方面大大降低了检测信杂比,影响雷达探测距离;另一方面极大增加了雷达错误锁定或者跟踪副瓣杂波的概率,导致丢失目标,这种情况亟需改进。
发明内容
本发明的目的在于提供一种雷达下视尾后探测下主副瓣地海杂波主分布计算方法,能够在线计算获取地海杂波在距离-多普勒域上的主分布准确位置,解决强副瓣杂波严重影响雷达下视尾后探测和跟踪目标性能的问题,为雷达下视尾后探测目标时使目标回波信号避开强副瓣杂波,避免落入强副瓣杂波区域提供必要参考条件,有效提升雷达下视尾后探测威力和跟踪性能。
为了达到上述目的,本发明通过以下技术方案实现:一种雷达下视尾后探测下主副瓣地海杂波主分布计算方法,其特点是,包含以下步骤:
S1、根据雷达下视尾后探测波束各方向增益强弱,确定主副瓣地海杂波主分布对应区域;
S2、在主副瓣地海杂波主分布对应区域内选取相应的杂波点,并计算所有杂波点的位置坐标;
S3、根据雷达的位置坐标、雷达的速度矢量参数及所有杂波点的位置坐标,计算所有杂波点与雷达的相对距离和相对多普勒频率;
S4、结合雷达发射信号的最大不模糊距离和雷达发射信号的脉冲重复频率参数,计算所有杂波点的视在距离和视在多普勒频率;
S5、量化所有杂波点的视在距离和视在多普勒频率,并将所有杂波点的量化距离和量化多普勒频率映射至距离-多普勒域,得到主副瓣地海杂波主分布的准确位置。
所述的步骤S1包含:
S1.1、分别建立基准坐标系和平台坐标系;
S1.2、根据实际雷达天线方向图强弱分布特征,以平台坐标系下的指向角,设置主副瓣地海杂波远距边界点和近距边界点与雷达天线相对指向,并根据雷达在基准坐标系下的位置和姿态分别计算远距边界点和近距边界点在基准坐标系下的坐标;
S1.3、在基准坐标系下,由远距边界点和近距边界点的连线构成主副瓣地海杂波主分布对应区域的计算区间。
所述的步骤S1.1中基准坐标系的原点设定在地面/海面上,X轴为北向,Y轴为天向,Z轴为东向,雷达所处位置定义为PR,在基准坐标系下的坐标为(xR,yR,zR),按照右手定则确定雷达姿态,雷达姿态定义为(Ψ,θ,γ),Ψ表示偏航角,θ表示俯仰角,γ表示滚转角。
所述的步骤S1.1中的平台坐标系是将基准坐标系原点平移至雷达所处位置,再将其依次按照Ψ、θ、γ姿态角进行三次转动,得到平台坐标系PRX″Y″Z″,指向角(α,β)表示雷达天线波束、地面/海面位置点与雷达在平台坐标系中的相对位置关系,α表示方位角,β表示俯仰角。
所述的步骤S1.2中远距边界点在基准坐标系下的坐标计算方法包含:
在平台坐标系下,由远距边界点指向角计算得到指向单位矢量;
通过欧拉坐标转换,将平台坐标系下的指向单位矢量转换到基准坐标系中,计算得到基准坐标系下指向单位矢量;
通过基准坐标系下指向单位矢量、平台坐标计算得到远距边界点在基准坐标系下的坐标。
所述的步骤S1.2中近距边界点在基准坐标系下的坐标计算方法包含:
在平台坐标系下,由近距边界点指向角计算得到指向单位矢量;
通过欧拉坐标转换,将平台坐标系下的指向单位矢量转换到基准坐标系中,计算得到基准坐标系下指向单位矢量;
通过基准坐标系下指向单位矢量、平台坐标计算得到近距边界点在基准坐标系下的坐标。
所述的步骤S2包含:
在基准坐标系下,根据远距边界点的坐标和近距边界点的坐标计算远距边界点和近距边界点的相对距离;
在基准坐标系下,通过远距边界点的坐标和近距边界点的坐标计算远距边界点到近距边界点连线与Z轴的夹角;
设定远距边界点和近距边界点连线上取点个数,依次计算各选取点相对于近距边界点的距离及各选取点在基准坐标系中的(X,Z)坐标;
依次根据各选取点在基准坐标系中的(X,Z)坐标,从高程数据中查询各选取点的Y坐标,以得到所有杂波点在基准坐标系下的位置坐标。
所述的步骤S3中所有杂波点与雷达的相对距离的计算方法为:
在基准坐标系下,利用杂波点的位置坐标和雷达的位置坐标,根据两点间距离计算公式可得所有杂波点与雷达的相对距离。
所述的步骤S3中所有杂波点与雷达的相对多普勒频率的计算方法为:
在基准坐标系下,利用杂波点的位置坐标和雷达的位置坐标,计算得到指向矢量;
利用指向矢量和雷达速度矢量计算得到各杂波点和雷达连线与速度方向的夹角;
利用平台合成速度值计算得打杂波点与雷达的相对多普勒频率。
所述的远距边界点和近距边界点连线上取点个数N,需满足以下条件:
其中,ΔR表示相邻取点直接的距离间隔,RFN表示远距边界点和近距边界点直接的相对距离,N表示取点个数,ΔRR表示雷达距离分辨率。
本发明一种雷达下视尾后探测下主副瓣地海杂波主分布计算方法与现有技术相比具有以下优点:通过确定主副瓣杂波主分布与地面/海面区域的对应关系简化杂波区域选择,缩小计算区域,计算区域内所有杂波单元与雷达的相对距离和相对多普勒频率,并映射至距离-多普勒域,得到杂波主分布的准确位置,能够在雷达下视尾后探测目标时,快速计算主副瓣杂波回波信号在距离-多普勒域主分布位置,为雷达利用目标信号与强副瓣杂波信号相对距离和频率关系判断或者避免目标信号落入或者穿越强副瓣杂波区提供核心参数和重要依据;本发明计算实时性强,能够满足机载、弹载雷达对处理实时性的要求。
附图说明
图1为本发明一种雷达下视尾后探测下主副瓣地海杂波主分布计算方法的流程图;
图2为雷达下视尾后探测的几何关系框图;
图3为本发明实施例中主副瓣地海杂波主分布计算结果示意图。
具体实施方式
以下结合附图,通过详细说明一个较佳的具体实施例,对本发明做进一步阐述。
如图1所示,一种雷达下视尾后探测下主副瓣地海杂波主分布计算方法,包含以下步骤:
S1、根据雷达下视尾后探测波束各方向增益强弱,确定主副瓣地海杂波主分布对应区域。
步骤S1具体包含:
S1.1、分别建立基准坐标系和平台坐标系,如图2所示;
其中,基准坐标系的原点设定在地面/海面上,X轴为北向,Y轴为天向,Z轴为东向,雷达所处位置定义为PR,在基准坐标系下的坐标为(xR,yR,zR),按照右手定则确定雷达姿态,雷达姿态定义为(Ψ,θ,γ),Ψ表示偏航角,θ表示俯仰角,γ表示滚转角;平台坐标系是将基准坐标系原点平移至雷达所处位置,再将其依次按照Ψ、θ、γ姿态角进行三次转动,得到平台坐标系PRX″Y″Z″,指向角(α,β)表示雷达天线波束、地面/海面位置点与雷达在平台坐标系中的相对位置关系,α表示方位角,β表示俯仰角;
S1.2、根据实际雷达天线方向图强弱分布特征,以平台坐标系下的指向角,设置主副瓣地海杂波远距边界点P1和近距边界点PN与雷达天线相对指向,远距边界点P1和近距边界点PN的指向角分别为(αFar,βFar)和(αNear,βNear),并根据雷达在基准坐标系下的位置和姿态分别计算远距边界点P1和近距边界点PN在基准坐标系下的坐标;
远距边界点P1在基准坐标系下的坐标计算方法包含:
在平台坐标系下,由远距边界点指向角(αFar,βFar)计算得到指向单位矢量(xP1,yP1,zP1);
xP1=cos(βFar)×cos(αFar)
yP1=cos(βFar)×sin(αFar)
zP1=sin(-βFar)
通过欧拉坐标转换,将平台坐标系下的指向单位矢量(xP1,yP1,zP1)转换到基准坐标系中,计算得到基准坐标系下指向单位矢量(xU1,yU1,zU1);
通过基准坐标系下指向单位矢量、平台坐标计算得到远距边界点P1在基准坐标系下的坐标(x1,y1,z1);
同理,近距边界点PN在基准坐标系下的坐标计算方法包含:
在平台坐标系下,由近距边界点指向角计算得到指向单位矢量;
通过欧拉坐标转换,将平台坐标系下的指向单位矢量转换到基准坐标系中,计算得到基准坐标系下指向单位矢量;
通过基准坐标系下指向单位矢量、平台坐标计算得到近距边界点PN在基准坐标系下的坐标(xN,yN,zN);
S1.3、在基准坐标系下,由远距边界点P1和近距边界点PN的连线构成主副瓣地海杂波主分布对应区域的计算区间。
S2、在主副瓣地海杂波主分布对应区域内选取相应的杂波点,并计算所有杂波点的位置坐标。
步骤S2具体包含:
在基准坐标系下,根据远距边界点的坐标(x1,y1,z1)和近距边界点的坐标(xN,yN,zN)计算远距边界点P1和近距边界点PN的相对距离RFN;
在基准坐标系下,通过远距边界点的坐标(x1,y1,z1)和近距边界点的坐标(xN,yN,zN)计算远距边界点P1到近距边界点PN连线与Z轴的夹角ε;
设定远距边界点P1和近距边界点PN连线上取点个数N,依次计算各选取点相对于近距边界点PN的距离RnN及各选取点在基准坐标系中的(X,Z)坐标(xn,zn);
其中,取点个数N需满足条件
其中ΔRR为雷达距离分辨率
RnN=(N-n+1)×ΔR
xn=RnN×sin(ε)+xN
zn=RnN×cos(ε)+zN
依次根据各选取点在基准坐标系中的(X,Z)坐标(xn,zn),从高程数据中查询各选取点的Y坐标yn,其中yn=HGn,HGn为选取点所在地面/海面高度,以得到所有杂波点在基准坐标系下的位置坐标(xn,yn,zn),n由1到N取值。
S3、根据雷达的位置坐标、雷达的速度矢量参数及所有杂波点的位置坐标,计算所有杂波点与雷达的相对距离和相对多普勒频率。
步骤S3具体包含:
在基准坐标系下,利用杂波点的位置坐标(xn,yn,zn)和雷达的位置坐标(xR,yR,zR),根据两点间距离计算公式可得所有杂波点与雷达的相对距离Rn;
在基准坐标系下,利用杂波点的位置坐标(xn,yn,zn)和雷达的位置坐标(xR,yR,zR),计算得到指向矢量
利用指向矢量和雷达速度矢量计算得到各杂波点和雷达连线与速度方向的夹角φn;
利用雷达合成速度值计算得到杂波点与雷达的相对多普勒频率Fn;
其中,λ为雷达发射信号波长,雷达速度矢量由惯性装置获得,机载雷达或者弹载雷达即对应的平台的速度矢量。
S4、结合雷达发射信号的最大不模糊距离和雷达发射信号的脉冲重复频率参数,计算所有杂波点的视在距离和视在多普勒频率。
步骤S4中视在距离ΔRn和视在多普勒频率ΔFn的计算公式分别为:
ΔRn=mod(Rn,Rmax)
其中,Rmax为发射信号最大不模糊距离;
ΔFn=mod(Fn,PRF)
其中,PRF为发射信号的脉冲重复频率。
S5、量化所有杂波点的视在距离和视在多普勒频率,并将所有杂波点的量化距离和量化多普勒频率映射至距离-多普勒域(ΔQRn,ΔQFn),(n=1,2…N),得到主副瓣地海杂波主分布的准确位置。
步骤S5中,量化所有杂波点的视在距离ΔRn和视在多普勒频率ΔFn,得到量化距离ΔQRn和量化多普勒频率ΔQFn,其中
其中,ΔRR为雷达距离分辨率;
其中,ΔFR雷达频率分辨率。
如图3所示,给出了雷达坐标(61310,14958,52736),雷达姿态(-38,-22,0),速度矢量(709,-371,564),远距边界点指向(-0.5,-10),近距边界点指向(0,-80),雷达发射信号的脉冲重复频率41kHz条件下,一个主副瓣地海杂波主分布计算结果图。
尽管本发明的内容已经通过上述优选实施例作了详细介绍,但应当认识到上述的描述不应被认为是对本发明的限制。在本领域技术人员阅读了上述内容后,对于本发明的多种修改和替代都将是显而易见的。因此,本发明的保护范围应由所附的权利要求来限定。
机译: 一种包括主单元和从单元的照明系统,其中,在从单元接管主单元的状态下,主单元可以切换到睡眠模式。
机译: 一种包括主单元和从单元的照明系统,其中,在从单元接管主单元的状态下,主单元可以切换到睡眠模式。
机译: 使用12V直流电(DC)电动机的电动风扇,包括主风扇主体和太阳能电池,其中所述主风扇主体包含外壳,风扇叶片和三相无刷直流电流(BLDC)电动机,以及在所说的BLDC电动机驱动风扇叶片周向旋转的情况下,其特征在于所说的BLDC电动机通过导电线直接连接到太阳能电池,而所说的BLDC电动机在其转子和定子磁芯以及其芯线所在的位置同样大小和规格的牙齿,甚至分布在