CN109633585B - 分布式机会阵雷达非合作目标动态回波的高精度计算方法 - Google Patents
分布式机会阵雷达非合作目标动态回波的高精度计算方法 Download PDFInfo
- Publication number
- CN109633585B CN109633585B CN201811471434.4A CN201811471434A CN109633585B CN 109633585 B CN109633585 B CN 109633585B CN 201811471434 A CN201811471434 A CN 201811471434A CN 109633585 B CN109633585 B CN 109633585B
- Authority
- CN
- China
- Prior art keywords
- scattering
- range gate
- radar
- equivalent
- target
- 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/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
本发明提供一种分布式机会阵雷达非合作目标动态回波的高精度计算方法,首先根据飞行器的航迹和雷达视线在飞行器坐标系上的时变姿态角,得到特定航迹下的非合作目标动态RCS序列,以该非合作目标动态RCS序列为基础,将回波信号分成若干个小平面,利用目标小平面单元模型,kirchhoff近似模型法和逆傅里叶计算得到每个距离门的等效散射中心的时域散射特性,将所述的时域散射特性与雷达发射信号的多距离门载频信号做卷积运算,获得每个等效散射中心的时域回波信号;对上述时域回波信号进行多普勒频率调制,延迟叠加之后得到宽带雷达目标的动态回波信号。本发明能够准确描述分布式机会阵雷达非合作目标的动态回波信号。
Description
技术领域
本发明涉及雷达信号处理领域,特别涉及一种分布式机会阵雷达非合作目标动态回波的高精度计算方法。
背景技术
分布式机会阵雷达系统的目标电磁散射特性和目标回波模型是分布式机会阵雷达系统的核心技术之一,分布式机会阵雷达由于雷达平台和目标之间的高速运动和姿态变化等使雷达目标回波信号随时间、空间和视线变化等因素具有瞬时动态变化大、不确定性高等特点,准确描述动态回波特性难度大。目前利用复杂目标的静态RCS(雷达散射截面积)计算是当前比较成熟的方法,但静态目标的RCS不能够精准的反映雷达动态目标的电磁散射。
发明内容
发明目的:为了能够准确的描述动态回波特性,本发明提供一种分布式机会阵雷达非合作目标动态回波的高精度计算方法。
技术方案:本发明提供一种分布式机会阵雷达非合作目标动态回波的高精度计算方法,具体包括如下步骤:
步骤1:根据飞行器的几何模型,建立飞行器全空域静态RCS序列数据库;
步骤2:根据飞行器的飞行航迹和在飞行过程中的姿态扰动影响,计算飞行器在雷达坐标系中的航迹和姿态变化,得到雷达视线在飞行器坐标系上的时变姿态角,根据时变姿态角在飞行器全空域静态RCS序列数据库中对应的RCS值,得到特定航迹下的非合作目标动态RCS序列;
步骤3:基于上述非合作目标动态RCS序列,根据雷达收、发平台位置和波束参数,确定收发波束的重叠区域,从而提取各散射点对应的RCS值和各散射点的距离,并根据各散射点的距离,得到散射点的回波相位和各散射点所在的距离门,叠加同一距离门内所有散射点的散射特性,得到该距离门的等效散射中心的频域散射特性,从而得该距离门的等效散射中心的时域散射特性;
步骤4:根据步骤3所述的距离门的等效散射中心的时域散射特性与雷达发射信号的多距离门载频信号,获得每个等效散射中心的时域回波信号;
步骤5:对每个等效散射中心的时域回波信号进行多普勒频率调制,延迟叠加之后得到宽带雷达目标的动态回波信号。
进一步的,所述步骤1采用GRECO(图形电磁)算法,根据飞行器与雷达之间的距离、散射场的强度和入射场的强度,利用GRECO算法得到飞行器全空域静态RCS序列数据库。
进一步的,所述步骤3中距离门的等效散射中心的时域散射特性的具体计算步骤为:
步骤3.1:将回波信号分成若干个小平面单元,根据复杂目标小平面单元模型,得到各小平面单元的中心坐标;
步骤3.2:根据kirchhoff近似模型,得到各小平面单元的后向散射系数;
步骤3.3:在距离向上根据天线、各小平面单元的几何关系、雷达收、发平台的位置、双基角,确定双基角平分线方向上的距离门所属的小平面单元;
步骤3.4:根据双基角平分线方向上的距离门所属的小平面单元的后向散射系数和叠加的同一距离门内所有散射点的散射特性,得到该距离门的等效散射中心的频域散射特性;
步骤3.5:对上述距离门的等效散射中心的频域散射特性进行逆傅里叶变换计算,得到该距离门的等效散射中心的时域散射特性。
进一步的,所述距离门的等效散射中心的时域散射特性与雷达发射信号的多距离门载频信号做卷积运算,得到每个等效散射中心对应的目标回波信号sbR(t);
其中hbn(t)为第n个距离门的等效散射中心的时域散射特性,sc(t)为雷达发射信号的多距离门载频信号:
其中an为第n个距离门的等效散射中心对应的雷达天线增益和空间衰减总和;sT(t)为雷达线性调频信号,B为雷达发射信号的带宽;n=1,2,3,…,N,N为回波信号中距离门的个数;
根据sT(t)=rect(t/TP)exp[j2πf0+jπμt2],其中TP为雷达发射信号的带宽的脉宽;μ为调频斜率μ=B/Tp,f0为雷达发射信号的起始频率,得到:
其中Hbn(f)为第n个距离门的等效散射中心的散射频域特性,PT为雷达发射功率;和分别为第n个距离门的等效散射中心在发射平台与目标视线方向接收平台与目标视线方向的发射、接收天线增益;θTn、为第n个距离门的等效散射中心在发射平台的仰俯角和方位角;θRn、为第n个距离门的等效散射中心在接收平台的仰俯角和方位角;RTn(tTn)和RRn(tRn)分别为发射、接收平台到第n个距离门的散射中心的距离,tTn与tRn分别为到达发射、接受平台的时间延迟;c为光速;j为复数。
进一步的,所述多普勒频率调制的具体方法为:设第n个距离门的等效散射中心的多普勒频率为fdn,则:
srn(t)=In(t)·cos[2πfdnt]-Qn(t)·sin[2πfdnt]
其中,In(t)和Qn(t)分别为第n个距离门的等效散射中心的回波I分量和Q分量,Srn(t)为fdn调制后的第n个距离门的等效散射中心回波信号。
有益效果:本发明是在静态目标RCS计算的基础上,通过全空间的静态目标RCS 数据序列来构建动态目标的RCS数据序列。在目标动态RCS的基础上,进一步计算目标的回波信号,本发明能够准确描述分布式机会阵雷达非合作目标的动态回波信号,是实现雷达系统资源管理(调度)、目标检测和目标跟踪的重要理论和技术基础。
附图说明
图1为本发明的目标动态RCS仿真流程图;
图2为本发明的飞行器的目标几何建模;
图3为本发明的计算等效散射中心的时域散射特性的方法。
具体实施方式
构成本发明的一部分的附图用来提供对本发明的进一步理解,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。
如图1-3所示,本发明提供一种分布式机会阵雷达非合作目标动态回波的高精度计算方法,具体包括如下步骤:
步骤1:根据飞行器的几何模型,建立飞行器全空域静态RCS序列数据库;
飞机外形重建应该遵守原有几何设计准则,以保证飞机的理论外形,像有些曲面之间是切矢,曲率连续,不能简单用点连续来代替,这样构造出来的曲面就会产生明显的分界线,造成不必要的误差。同时,构造曲面时应严格按所绘制的轮廓线去定义划分区域,不能简单的以整张曲面去替换所在区域的全部曲面,这样会造成曲面过渡问题,使几何图形失真。
飞机外形的各个部件的特点都不相同,像机身是由纵向轮廓线,横向截面定义构成,进气道在纵向上为按一定规律的横向变截面的曲面组成,在构造曲面时,应抓住部件的具体特点,进行分块构建。为确保分块建模精确形成,必须先绘制纵向控制线和截面线构成曲面的轮廓线,再定义曲面。
利用图形电磁(GRECO)计算法,得到飞行器目标RCS:
其中R为飞机与雷达之间的距离;ES为散射场;Ei为入射场。
步骤2:根据飞行器的飞行航迹和在飞行过程中的姿态扰动影响,计算飞行器在雷达坐标系中的航迹和姿态变化,得到雷达视线在飞行器坐标系上的时变姿态角,根据时变姿态角在飞行器全空域静态RCS序列数据库中对应的RCS值,得到特定航迹下的非合作目标动态RCS序列;
步骤3:基于上述非合作目标动态RCS序列,将每个距离门的RCS值分解成二维或三维的等效散射中心,根据雷达收、发平台位置和波束参数,确定收发波束的重叠区域,从而提取各散射点对应的RCS值和各散射点的距离,根据RCS值获得散射点的电磁散射特性,并根据各散射点的距离,得到散射点的回波相位和各散射点所在的距离门,叠加同一距离门内所有散射点的散射特性,得到该距离门的等效散射中心的频域散射特性,从而得到该距离门的等效散射中心的时域散射特性;
步骤4:根据步骤3所述的距离门的等效散射中心的时域散射特性与雷达发射信号的多距离门载频信号,获得每个等效散射中心对应的目标时域回波信号;
步骤5:对每个等效散射中心对应的目标时域回波信号进行多普勒频率调制,延迟叠加之后得到宽带雷达目标的动态回波信号。
所述步骤3中距离门的等效散射中心的时域散射特性的具体计算步骤为:
步骤3.1:将回波信号分成若干个小平面单元,根据复杂目标小平面单元模型,得到各小平面单元的中心坐标;
步骤3.2:根据基尔霍夫(kirchhoff)近似模型,得到各小平面单元的后向散射系数;
步骤3.3:在距离向上根据天线、各小平面单元的几何关系、雷达收、发平台的位置、双基角,确定双基角平分线方向上的距离门所属的小平面单元,所述双基角等效为单机地雷达的仰俯角、方位角;
根据双基角平分线方向上的距离门所属的小平面单元的后向散射系数和叠加的同一距离门内所有散射点的散射特性,得到该距离门的等效散射中心的频域散射特性;
步骤3.5:对上述距离门的等效散射中心的频域散射特性进行逆傅里叶变换计算,得到该距离门的等效散射中心的时域散射特性;
hbn(t)=FT-1[Hbn(f)]
其中FT-1为逆傅里叶变换,hbn(t)为第n个距离门的等效散射中心的时域散射特性,其中Hbn(f)为第n个距离门的等效散射中心的频域散射特性。
步骤4中,所述距离门的等效散射中心的时域散射特性与雷达发射信号的多距离门载频信号做卷积运算,得到每个等效散射中心对应的目标回波信号sbR(t);
sc(t)为雷达发射信号的多距离门载频信号:
其中an为第n个距离门的等效散射中心对应的雷达天线增益和空间衰减总和;sT(t)为雷达线性调频信号,B为雷达发射信号的带宽;N为回波信号中距离门的个数;
其中Hbn(f)为第n个距离门的等效散射中心的散射频域特性,PT为雷达发射功率;和分别为第n个距离门的等效散射中心在发射平台与目标视线方向接收平台与目标视线方向的发射、接收天线增益;θTn、为第n个距离门的等效散射中心在发射平台的仰俯角和方位;θRn、为第n个距离门的等效散射中心在接受平台的仰俯角和方位,tTn与tRn分别为到达发射、接受平台的时间延迟;c为光速;j为复数,FT-1为逆傅里叶变换。
所述多普勒频率调制的具体方法为:设第n个距离门的等效散射中心的多普勒频率为fdn,则:
srn(t)=In(t)·cos[2πfdnt]-Qn(t)·sin[2πfdnt]
其中,In(t)和Qn(t)分别为第n个距离门的等效散射中心的回波I分量和Q分量,Srn(t)为fdn调制后的第n个距离门的等效散射中心回波信号。
另外需要说明的是,在上述具体实施方式中所描述的各个具体技术特征,在不矛盾的情况下,可以通过任何合适的方式进行组合。为了避免不必要的重复,本发明对各种可能的组合方式不再另行说明。
Claims (3)
1.分布式机会阵雷达非合作目标动态回波的高精度计算方法,其特征在于,具体包括如下步骤:
步骤1:根据飞行器的几何模型,建立飞行器全空域静态RCS序列数据库;
步骤2:根据飞行器的飞行航迹和在飞行过程中的姿态扰动影响,计算飞行器在雷达坐标系中的航迹和姿态变化,得到雷达视线在飞行器坐标系上的时变姿态角,根据时变姿态角在飞行器全空域静态RCS序列数据库中对应的RCS值,得到特定航迹下的非合作目标动态RCS序列;
步骤3:基于上述非合作目标动态RCS序列,将每个距离门的RCS值分解成二维或三维的等效散射中心,根据雷达收、发平台位置和波束参数,确定收发波束的重叠区域,从而提取各散射点对应的RCS值和各散射点的距离,根据RCS值获得散射点的电磁散射特性,并根据各散射点的距离,得到散射点的回波相位和各散射点所在的距离门,叠加同一距离门内所有散射点的散射特性,得到该距离门的等效散射中心的频域散射特性,从而得到该距离门的等效散射中心的时域散射特性;
步骤4:根据步骤3所述的距离门的等效散射中心的时域散射特性与雷达发射信号的多距离门载频信号,获得每个等效散射中心的目标时域回波信号;
步骤5:对每个等效散射中心的目标时域回波信号进行多普勒频率调制,延迟叠加之后得到宽带雷达目标的动态回波信号;
所述步骤3中距离门的等效散射中心的时域散射特性的具体计算步骤为:
步骤3.1:将回波信号分成若干个小平面单元,根据复杂目标小平面单元模型,得到各小平面单元的中心坐标;
步骤3.2:根据kirchhoff近似模型,得到各小平面单元的后向散射系数;
步骤3.3:在距离向上根据天线、各小平面单元的几何关系、雷达收、发平台的位置、双基角,确定双基角平分线方向上的距离门所属的小平面单元;
步骤3.4:根据双基角平分线方向上的距离门所属的小平面单元的后向散射系数和叠加的同一距离门内所有散射点的散射特性,得到该距离门的等效散射中心的频域散射特性;
步骤3.5:对上述距离门的等效散射中心的频域散射特性进行逆傅里叶变换计算,得到该距离门的等效散射中心的时域散射特性;
所述距离门的等效散射中心的时域散射特性与雷达发射信号的多距离门载频信号做卷积运算,得到每个等效散射中心对应的目标回波信号sbR(t);
其中hbn(t)为第n个距离门的等效散射中心的时域散射特性,sc(t)为雷达发射信号的多距离门载频信号:
其中an为第n个距离门的等效散射中心对应的雷达天线增益和空间衰减总和;sT(t)为雷达线性调频信号,B为雷达发射信号的带宽;n=1,2,3,…,N,N为回波信号中距离门的个数;
根据sT(t)=rect(t/TP)exp[j2πf0+jπμt2],其中TP为雷达发射信号的带宽的脉宽;μ为调频斜率μ=B/Tp,f0为雷达发射信号的起始频率,得到:
2.根据权利要求1所述的分布式机会阵雷达非合作目标动态回波的高精度计算方法,其特征在于,所述步骤1采用GRECO算法,根据飞行器与雷达之间的距离、散射场的强度和入射场的强度,利用GRECO算法得到飞行器全空域静态RCS序列数据库;
3.根据权利要求1所述的分布式机会阵雷达非合作目标动态回波的高精度计算方法,其特征在于,所述多普勒频率调制的具体方法为:设第n个距离门的等效散射中心的多普勒频率为fdn,则:
srn(t)=In(t)·cos[2πfdnt]-Qn(t)·sin[2πfdnt]
其中,In(t)和Qn(t)分别为第n个距离门的等效散射中心的回波I分量和Q分量,Srn(t)为fdn调制后的第n个距离门的等效散射中心回波信号。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811471434.4A CN109633585B (zh) | 2018-12-04 | 2018-12-04 | 分布式机会阵雷达非合作目标动态回波的高精度计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811471434.4A CN109633585B (zh) | 2018-12-04 | 2018-12-04 | 分布式机会阵雷达非合作目标动态回波的高精度计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109633585A CN109633585A (zh) | 2019-04-16 |
CN109633585B true CN109633585B (zh) | 2023-01-13 |
Family
ID=66070852
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811471434.4A Active CN109633585B (zh) | 2018-12-04 | 2018-12-04 | 分布式机会阵雷达非合作目标动态回波的高精度计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109633585B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114485676B (zh) * | 2022-02-18 | 2023-07-21 | 四川大学 | 一种分布式飞行雷达平台的航迹规划方法 |
CN117572423B (zh) * | 2024-01-19 | 2024-04-16 | 逻安科技(天津)有限责任公司 | 多普勒孔径同收发机会阵雷达 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5281972A (en) * | 1992-09-24 | 1994-01-25 | Hughes Aircraft Company | Beam summing apparatus for RCS measurements of large targets |
CN103376447A (zh) * | 2012-04-24 | 2013-10-30 | 中国人民解放军海军航空工程学院 | 一种实现非合作双多基地雷达目标三维定位的方法 |
CN104865562A (zh) * | 2015-06-12 | 2015-08-26 | 西安电子科技大学 | 基于混合模型的雷达非合作目标的识别方法 |
CN107656264A (zh) * | 2017-08-02 | 2018-02-02 | 南京航空航天大学 | 杂波环境下的机会阵雷达多目标跟踪的功率资源管理方法 |
CN108693510A (zh) * | 2018-05-18 | 2018-10-23 | 西安电子科技大学 | 基于gpu的知识辅助机载机会阵地杂波快速仿真方法 |
-
2018
- 2018-12-04 CN CN201811471434.4A patent/CN109633585B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5281972A (en) * | 1992-09-24 | 1994-01-25 | Hughes Aircraft Company | Beam summing apparatus for RCS measurements of large targets |
CN103376447A (zh) * | 2012-04-24 | 2013-10-30 | 中国人民解放军海军航空工程学院 | 一种实现非合作双多基地雷达目标三维定位的方法 |
CN104865562A (zh) * | 2015-06-12 | 2015-08-26 | 西安电子科技大学 | 基于混合模型的雷达非合作目标的识别方法 |
CN107656264A (zh) * | 2017-08-02 | 2018-02-02 | 南京航空航天大学 | 杂波环境下的机会阵雷达多目标跟踪的功率资源管理方法 |
CN108693510A (zh) * | 2018-05-18 | 2018-10-23 | 西安电子科技大学 | 基于gpu的知识辅助机载机会阵地杂波快速仿真方法 |
Non-Patent Citations (1)
Title |
---|
Time resource management of OAR based on fuzzy logic priority for multiple target tracking;Han qignhua等;《JournalofSystemsEngineeringandElectronics》;20180831;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN109633585A (zh) | 2019-04-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104007440B (zh) | 一种加速分解后向投影聚束合成孔径雷达成像方法 | |
CN101526614B (zh) | 基于子孔径和等效散射体的sar回波快速模拟方法 | |
CN103176170B (zh) | 一种基于gpu并行计算的sar回波模拟方法 | |
CN107976660B (zh) | 弹载多通道雷达超低空目标分析与多径回波建模方法 | |
CN103207387B (zh) | 一种机载相控阵pd雷达杂波的快速模拟方法 | |
CN109613532B (zh) | 一种机载雷达实时多普勒波束锐化超分辨成像方法 | |
CN104914415A (zh) | 基于目标距离像模板匹配的单脉冲雷达相干干扰方法 | |
CN107831480A (zh) | 弹载雷达和差通道稳健自适应杂波抑制方法 | |
CN106291489B (zh) | 适用于多种发射信号波形的合成孔径雷达回波仿真方法 | |
CN107942295B (zh) | 一种前视阵列sar系统的稀疏天线 | |
CN106054187B (zh) | 基于斜距模型下的大斜视sar曲线轨迹波数域成像方法 | |
CN109633585B (zh) | 分布式机会阵雷达非合作目标动态回波的高精度计算方法 | |
CN112859074B (zh) | 多频带多视角isar融合成像方法 | |
CN108387884A (zh) | 基于知识辅助稀疏渐进最小方差的机载雷达杂波抑制方法 | |
CN105738887A (zh) | 基于多普勒通道划分的机载雷达杂波功率谱的优化方法 | |
CN107271997A (zh) | 机载多通道cssar地面运动目标运动参数估计方法 | |
CN113567978B (zh) | 一种多基分布式雷达协同成像方法 | |
CN113376625B (zh) | 目标物体的偏离角度获得方法、装置、电子设备及存储介质 | |
Sun et al. | Polar format algorithm for spotlight bistatic SAR with arbitrary geometry configuration | |
CN111708015B (zh) | 一种多径效应下的低空目标跟踪滤波方法 | |
CN103675777B (zh) | 基于拟合法的机载雷达杂波模拟方法及装置 | |
CN104035078A (zh) | 一种基于阵元阶数递推的降维空时自适应权值计算方法 | |
CN103869297A (zh) | 高超声速平台载雷达的杂波生成系统及其方法 | |
CN115480245A (zh) | 一种基于二维超分辨算法的机载雷达前视三维成像技术 | |
CN109212525A (zh) | 一种基于keystone变换的合成孔径雷达逆投影成像方法 |
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 |