CN108152836A - 一种基于fda天线抗三星频差定位体制的欺骗方法 - Google Patents
一种基于fda天线抗三星频差定位体制的欺骗方法 Download PDFInfo
- Publication number
- CN108152836A CN108152836A CN201810003504.7A CN201810003504A CN108152836A CN 108152836 A CN108152836 A CN 108152836A CN 201810003504 A CN201810003504 A CN 201810003504A CN 108152836 A CN108152836 A CN 108152836A
- Authority
- CN
- China
- Prior art keywords
- fda
- frequency
- samsung
- satellite
- frequency difference
- 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.)
- Pending
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
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/21—Interference related issues ; Issues related to cross-correlation, spoofing or other methods of denial of service
- G01S19/215—Interference related issues ; Issues related to cross-correlation, spoofing or other methods of denial of service issues related to spoofing
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明属于目标定位技术领域,具体的说是一种基于FDA天线抗三星频差定位体制的欺骗方法。本发明在相控阵天线的基础上,在相邻阵元对发射信号附加了一个远小于工作载频的频率增量,使得卫星侦察接收机无法准确测量目标辐射源与侦察接收机之间的多普勒频率,从而导致卫星侦察系统采用基于FDOA的定位方法时,无法获得目标辐射源的真实位置,即利用FDA有效解决了上述现有技术中存在的问题。
Description
技术领域
本发明属于目标定位技术领域,具体的说是一种利用FDA(FDA,FrequencyDiverse Array频率分集阵列)天线辐射的信号进行抗三星频差定位体制的欺骗方法。
背景技术
相控阵天线通常每个阵元发射的是同一个信号,通过在每个阵元的输出端接入移相器来实现波束方向控制,调整移相器的相移量便可实现波束的空域扫描。
相控阵雷达因其灵活的波束扫描特性被广泛的应用到军事和民用领域。然而,因其每个阵元发射的信号频率相同,当相控阵辐射源暴露在由三颗侦察卫星组成的侦察系统侦收区域时,侦察系统利用现有的测频技术能够准确获取相控阵辐射源的工作频率。这就意味着多个卫星侦察接收机将得到目标辐射源的多普勒频差信息,结合侦察卫星自身位置,利用FDOA(Frequency Difference of Arrival,频率到达差)定位方法就能计算目标辐射源的位置。暴露目标辐射源的真实位置是非常危险的,尤其在战争时期,这将同时危及周边重要军事设施的安全。如何隐蔽目标辐射源的真实位置,针对卫星侦察系统实现定位欺骗显得尤为重要。
发明内容
本发明的目的,就是针对上述问题,提出一种利用FDA天线抗三星频差测量的定位欺骗方法。本发明所采用的频率分集阵列天线与相控阵天线不同之处在于:FDA是相控阵天线的扩展,其在相邻阵元上对发射信号附加了一个远小于工作载频的频率增量,即每个阵元辐射信号的频率不同。因频率增量较小,且不同阵元发射的信号频率不同,考虑到信号在频域分析中不可避免的产生栅栏效应和频率泄露现象。这一改变使得传统测频方法的测频误差增大,从而使得卫星侦察接收机难以有效获取目标辐射源发射信号的频率真实值,也就无法在多个卫星侦察接收机之间形成准确的多普勒频差信息。那么,当卫星侦察系统采用基于FDOA的三星频差定位方法时,也就无法完成对目标辐射源的准确定位,即利用FDA天线辐射的信号可以实现抗三星频差测量的定位欺骗。
本发明所采用的技术方案为:
一种基于FDA天线抗三星频差定位体制的欺骗方法,其特征在于,包括以下步骤:
a、对卫星侦察接收机接收的FDA信号s(θ,r,t)进行中频采样后,再进行快速傅里叶变换,其中,θ为波束指向角,r为目标辐射源第一个阵元到卫星侦察接收机的距离,t为目标辐射源发射FDA信号的时间;
b、对经步骤a处理后的信号进行CFAR(Constant False Alarm Ratio,衡虚警)检测,记录超过CFAR门限Th所对应的频谱峰值;其中,Th由系统参数设定;
c、对经步骤b获得的所有频点进行平均,得到粗的频率测量值;
d、以粗的频率测量值为中心,在频率范围内,利用Chirp-Z(一种频率估计方法)变换得到细化后的频率测量值;其中,由系统参数设定;
e、重复步骤a-d,直至获得第i部卫星侦察接收机的频率测量值fi,i=2,3;以第一部卫星接收机为基准,得到第i部卫星接收机相对于第一部卫星接收机的到达频率差fi1,i=2,3;
f、对经步骤e获得的到达频率差信息建立频差定位方程,采用最大似然估计法进行求解,从而得到目标辐射源的位置坐标。
具体的,所述步骤a中接收的FDA信号s(θ,r,t)具有如下表达式:
其中,M为阵元个数,f0为FDA天线第一个阵元辐射的频率,λ0为第一个阵元发射信号波长,d为阵元间距,△f为相邻阵元的载频增量,c为光速。
具体的,所述步骤b中的CFAR检测方法为OS-CFAR(有序统计量恒虚警检测),其检测门限Th由系统参数设定。
具体的,所述步骤d中Chirp-Z变换方法为:
其中,A0表示起始取样点的半径长度,θ0表示起始取样点的相角,W0表示螺旋线的伸展率,表示相邻两点间的等分角;A0,θ0,W0,由Chirp-Z变换的参数设定。
具体的,所述步骤f中频差方程具有如下形式:
其中,pi,vi,(i=1,2,3)分别为第i个卫星的位置矢量和速度矢量,q为目标的位置矢量
具体的,所述步骤f中最大似然估计的方法为:
其中,Ω为卫星覆盖区域,为待估计的频差向量,f(q)=【f21(q),f31(q)]T。
本发明的有益效果为,与相控阵目标辐射源相比,由三颗卫星组成的侦察系统不能通过多普勒频差对使用FDA天线的目标辐射源进行准确定位。也就是说,本发明对三星频差定位系统具有定位欺骗性,特别适合一些特定场合的应用。
附图说明
图1为实施例中FDA天线发射阵列示意图;
图2为实施例中接收的FDA信号的频谱图;
图3为实施例中Chirp-Z变换后的频谱图;
图4为实施例中三星频差定位示意图;
图5为实施例中基于FDOA方法的目标定位结果。
具体实施方式
下面结合附图和具体设施方式对本发明进行详细的描述,以便本领域的技术人员更好地理解本发明。
实施例
本实施例是以由三颗侦察卫星组成的侦察系统对使用FDA天线的目标辐射源进行测频定位为例。FDA天线发射阵列示意图如图1所示,卫星侦察接收机接收的FDA信号为:
其中,θ为波束指向角,r为目标辐射源到卫星侦察接收机第一个阵元的距离,t为目标辐射源发射FDA信号的时间,M为阵元个数,f0为FDA天线第一个阵元辐射的频率,λ为第一个阵元发射信号波长,d为阵元间距,△f为相邻阵元的载频增量,c为光速。
根据上述信号模型,使用Matlab(一种计算机算法语言)进行仿真验证,具体的仿真参数如下:
系统参数设置:FDA天线第一个阵元的辐射源频率f0=8GHz,频率增量△f=6KHz,阵元间隔d=λ/2,阵元个数M=8,中频频率fIF=30MHz,采样率fs=100MHz,OS-CFAR的检测门限Th=0.5。
侦察接收机参数设置:搭载侦察接收机的三颗卫星平台位置坐标依次为(3909.0,4499.4,3452.7),(3892.3,4527.8,3454.5),(3868.3,4545.9,3457.7),单位:km;三颗卫星的速度坐标依次为(-1.667,-3.553,6.517),(-1.706,-3.510,6.520),(-1.752,-3.476,6.529),单位:km/s。
目标参数设置:目标位置二维坐标为(3362.1,4006.7),单位:km。目标与卫星侦察接收机位于同一坐标系下。
具体定位过程如下:
A.对卫星侦察接收机接收的FDA信号进行中频采样,再进行快速傅里叶变换。快速傅里叶变换的点数为16384,其变换结果如图2所示,可以看出:由于FDA信号本身的特殊性,接收信号的频谱包含多个峰值。
B.采用OS-CFAR方法对快速傅里叶变换的结果进行检测,将门限Th与信号频谱的幅度比较,记录大于门限Th所对应的频谱峰值。其中,Th=0.5。
C.对所有大于检测门限的频谱峰值对应的信号频率进行平均,得到粗的频率测量值fco。
D.以粗的频率测量值fco为中心,在频率范围内,利用Chirp-Z变换得到细化后的频率测量值f1。Chirp-Z变换公式为:
其中,N=2048,A0=1,θ0=-1.27rad,W0=1,经Chirp-Z变换后频谱如图3所示。
E.重复步骤A~D,获得第i部卫星侦察接收机的频率测量值fi,i=2,3;根据下式得到第i部卫星接收机相对于第一部卫星接收机的到达频率差
fi1=fi-f1 i=2,3
F.利用到达频率差信息建立频差定位方程如下:
其中,p1=(3909.0,4499.4,3452.7)T p2=(3892.3,4527.8,3454.5)T,
p3=(3868.3,4545.9,3457.7)T,单位:km;f21=1085Hz,f31=2075Hz;
v1=(-1.667,-3.553,6.517)T,v2=(-1.706,-3.510,6.520)T,v3=(-1.752,-3.476,6.529)T,单位km/s。
接下来,利用最大似然法对目标位置进行估计:
为获得准确的目标位置坐标,在一定区域内进行搜索。当上式出现最小值时,即得到估计的目标位置坐标。
为验证本发明的性能,采用蒙特卡洛仿真的方法对10dB信噪比下的定位性能进行统计。设蒙特卡洛仿真次数为400次,忽略目标高度误差带来的影响。采用三星频差定位体制(如图4所示),得到目标的定位结果如图5所示。其中,“〇”表示目标真实位置,“☆”表示FDA辐射源目标的定位结果,“*”表示相控阵辐射源目标的定位结果,“+”表示三颗卫星的位置。可以看出,利用三星频差定位方法估计得到的相控阵辐射源目标位置与目标真实位置偏差很小,而估计的FDA辐射源目标位置与目标真实位置偏差很大。
表1列出了三颗卫星侦察接收机分别接收相控阵辐射源信号与FDA辐射信号条件下,采用FDOA方法获得的目标定位结果以及定位结果的RMSE(Root Mean Square Error)。可以看出:利用频差定位方法估计得到的FDA辐射源目标定位结果的RMSE比相控阵辐射源目标定位结果的RMSE大7倍以上,这说明本发明对三星频差定位功能的卫星侦察系统具有良好的定位欺骗性。
表1三颗卫星侦察接收机分别接收相控阵辐射源信号与FDA辐射信号条件下,采用FDOA
方法获得的目标定位结果以及定位结果的RMSE
目标位置 | (3362.1,4006.7) |
相控阵辐射源的定位结果(km) | (3367.2,4012.6) |
FDA辐射源的定位结果(km) | (3396.0,4038.5) |
相控阵辐射源定位的RMSE(km) | 6.0 |
FDA辐射源定位定位的RMSE(km) | 46.5 |
Claims (6)
1.一种基于FDA天线抗三星频差定位体制的欺骗方法,其特征在于,包括以下步骤:
a、对卫星侦察接收机接收的FDA信号s(θ,r,t)进行中频采样后,再进行快速傅里叶变换,其中,θ为波束指向角,r为目标辐射源第一个阵元到卫星侦察接收机的距离,t为目标辐射源发射FDA信号的时间;
b、对经步骤a处理后的信号进行CFAR检测,记录超过CFAR门限Th所对应的频谱峰值;其中,Th由系统参数设定;
c、对经步骤b获得的所有频点进行平均,得到粗的频率测量值;
d、以粗的频率测量值为中心,在频率范围内,利用Chirp-Z变换得到细化后的频率测量值;其中,由系统参数设定;
e、重复步骤a-d,直至获得第i部卫星侦察接收机的频率测量值fi,i=2,3;以第一部卫星接收机为基准,得到第i部卫星接收机相对于第一部卫星接收机的到达频率差fi1,i=2,3;
f、对经步骤e获得的到达频率差信息建立频差定位方程,采用最大似然估计法进行求解,从而得到目标辐射源的位置坐标。
2.根据权利要求1所述的一种基于FDA天线抗三星频差定位体制的欺骗方法,其特征在于,所述步骤a中接收的FDA信号s(θ,r,t)具有如下表达式:
其中,M为阵元个数,f0为FDA天线第一个阵元辐射的频率,λ0为第一个阵元发射信号波长,d为阵元间距,△f为相邻阵元的载频增量,c为光速。
3.根据权利要求2所述的一种基于FDA天线抗三星频差定位体制的欺骗方法,其特征在于,所述步骤b中的CFAR检测方法为有序统计量恒虚警检测,其检测门限Th由系统参数设定。
4.根据权利要求3所述的一种基于FDA天线抗三星频差定位体制的欺骗方法,其特征在于,所述步骤d中Chirp-Z变换方法为:
其中,A0表示起始取样点的半径长度,θ0表示起始取样点的相角,W0表示螺旋线的伸展率,表示相邻两点间的等分角;A0,θ0,W0,由Chirp-Z变换的参数设定。
5.根据权利要求4所述的一种基于FDA天线抗三星频差定位体制的欺骗方法,其特征在于,所述步骤f中频差方程具有如下形式:
其中,pi,vi,(i=1,2,3)分别为第i个卫星的位置矢量和速度矢量,q为目标的位置矢量。
6.根据权利要求5所述的一种基于FDA天线抗三星频差定位体制的欺骗方法,其特征在于,所述步骤f中最大似然估计的方法为:
其中,Ω为卫星覆盖区域,为待估计的频差向量,f(q)=[f21(q),f31(q)]T。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810003504.7A CN108152836A (zh) | 2018-01-03 | 2018-01-03 | 一种基于fda天线抗三星频差定位体制的欺骗方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810003504.7A CN108152836A (zh) | 2018-01-03 | 2018-01-03 | 一种基于fda天线抗三星频差定位体制的欺骗方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108152836A true CN108152836A (zh) | 2018-06-12 |
Family
ID=62460518
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810003504.7A Pending CN108152836A (zh) | 2018-01-03 | 2018-01-03 | 一种基于fda天线抗三星频差定位体制的欺骗方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108152836A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109031188A (zh) * | 2018-06-14 | 2018-12-18 | 中国人民解放军战略支援部队信息工程大学 | 一种基于蒙特卡罗的窄带辐射源频差估计方法及装置 |
CN109581450A (zh) * | 2018-12-18 | 2019-04-05 | 电子科技大学 | 基于频控阵天线的双站定位欺骗方法 |
CN111273223A (zh) * | 2020-01-20 | 2020-06-12 | 西安交通大学 | 一种基于双天线时延无源定位对抗方法 |
FR3105833A1 (fr) * | 2019-12-31 | 2021-07-02 | Thales | Procede anti-geolocalisation d'un premier terminal de telecommunication par satellite, dispositif et programme d'ordinateur associes |
CN117078863A (zh) * | 2023-08-30 | 2023-11-17 | 电子科技大学 | 一种异构阵列天线的快速波束赋形方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103645485A (zh) * | 2013-10-28 | 2014-03-19 | 中国科学院国家授时中心 | 一种基于双星时差频差无源定位的伪距差分方法 |
CN104849738A (zh) * | 2015-04-28 | 2015-08-19 | 中国电子科技集团公司第三十六研究所 | 一种卫星定位系统和定位方法 |
CN105044689A (zh) * | 2015-04-14 | 2015-11-11 | 电子科技大学 | 一种基于频控阵的射频隐身方法及设备 |
CN106353744A (zh) * | 2016-10-28 | 2017-01-25 | 中国人民解放军信息工程大学 | 基于双基地fda‑mimo雷达的多参数联合估计方法 |
CN107037410A (zh) * | 2017-04-17 | 2017-08-11 | 电子科技大学 | 一种对雷达进行干扰的方法、装置及频控阵干扰机 |
-
2018
- 2018-01-03 CN CN201810003504.7A patent/CN108152836A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103645485A (zh) * | 2013-10-28 | 2014-03-19 | 中国科学院国家授时中心 | 一种基于双星时差频差无源定位的伪距差分方法 |
CN105044689A (zh) * | 2015-04-14 | 2015-11-11 | 电子科技大学 | 一种基于频控阵的射频隐身方法及设备 |
CN104849738A (zh) * | 2015-04-28 | 2015-08-19 | 中国电子科技集团公司第三十六研究所 | 一种卫星定位系统和定位方法 |
CN106353744A (zh) * | 2016-10-28 | 2017-01-25 | 中国人民解放军信息工程大学 | 基于双基地fda‑mimo雷达的多参数联合估计方法 |
CN107037410A (zh) * | 2017-04-17 | 2017-08-11 | 电子科技大学 | 一种对雷达进行干扰的方法、装置及频控阵干扰机 |
Non-Patent Citations (2)
Title |
---|
徐海源 等: "三星多普勒频差无源定位方法及定位精度分析", 《宇航学报》 * |
王文钦 等: "频控阵雷达概念、原理与应用", 《电子与信息学报》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109031188A (zh) * | 2018-06-14 | 2018-12-18 | 中国人民解放军战略支援部队信息工程大学 | 一种基于蒙特卡罗的窄带辐射源频差估计方法及装置 |
CN109031188B (zh) * | 2018-06-14 | 2021-01-29 | 中国人民解放军战略支援部队信息工程大学 | 一种基于蒙特卡罗的窄带辐射源频差估计方法及装置 |
CN109581450A (zh) * | 2018-12-18 | 2019-04-05 | 电子科技大学 | 基于频控阵天线的双站定位欺骗方法 |
FR3105833A1 (fr) * | 2019-12-31 | 2021-07-02 | Thales | Procede anti-geolocalisation d'un premier terminal de telecommunication par satellite, dispositif et programme d'ordinateur associes |
CN111273223A (zh) * | 2020-01-20 | 2020-06-12 | 西安交通大学 | 一种基于双天线时延无源定位对抗方法 |
CN117078863A (zh) * | 2023-08-30 | 2023-11-17 | 电子科技大学 | 一种异构阵列天线的快速波束赋形方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108152836A (zh) | 一种基于fda天线抗三星频差定位体制的欺骗方法 | |
US9341706B2 (en) | Radar system | |
CN104678368B (zh) | 一维相扫三坐标雷达空域自适应干扰抑制方法 | |
CN104515971B (zh) | 宽带多目标机载单站无源定位方法 | |
CN110850400B (zh) | 基于干涉仪测向的lfmcw雷达多目标检测方法 | |
CN106054165B (zh) | 一种实现分布式无源雷达目标检测的方法 | |
CN105467369B (zh) | 一种目标回波仿真方法和装置 | |
CN102819008A (zh) | 一种基于非线性最小二乘的非合作雷达辐射源定位方法 | |
CN105738887B (zh) | 基于多普勒通道划分的机载雷达杂波功率谱的优化方法 | |
CN104730503A (zh) | 确定高分辨率sar参考目标rcs对定标影响的方法及补偿方法 | |
CN114002673A (zh) | 星载无源sar的非合作信号感知系统与多维参数估计方法 | |
CN105891817A (zh) | 一种无直达波条件下分布式无源雷达目标检测的方法 | |
RU2557808C1 (ru) | Способ определения наклонной дальности до движущейся цели пассивным моностатическим пеленгатором | |
CN112666543A (zh) | 一种稀疏阵列tdm-mimo雷达及其校正方法 | |
CN109521418B (zh) | 基于干涉场的地基雷达测角方法 | |
CN107907860A (zh) | 一种利用频控阵技术抗时差测量的定位欺骗方法 | |
CN112147593B (zh) | 一种高速密集爆炸物破片目标的四维参数估计方法 | |
CN112859028A (zh) | 一种外场地物散射时域回波采集及谱分析方法及系统 | |
CN104020465A (zh) | 基于八单元小孔径圆阵天线的外辐射源雷达测角方法 | |
KR102229191B1 (ko) | 방향 탐지 장치를 위한 보정 시스템 | |
US10845475B2 (en) | Method of measuring azimuth of radar target | |
CN115877346B (zh) | 一种基于二维相控阵雷达的无人机载脱靶量矢量检测方法 | |
Yu et al. | Single RFI localization based on conjugate cross-correlation of dual-channel SAR signals | |
Liu et al. | A clutter cancellation method for airborne passive radar | |
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 | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20180612 |
|
WD01 | Invention patent application deemed withdrawn after publication |