欺骗式干扰信号数据仿真方案.docx

1GNSS欺骗式干扰系统组成

GNSS欺骗式干扰系统由GNSS实时星历/时间接收系统、目标运动信息探测系统、欺骗式GNSS信号源组成，设备的连接关系如图1所示。

系统工作时，GNSS实时星历/时间接收系统接收GNSS真实信号，解调出导航卫星广播星历，计算出当前时间；目标运动信息探测系统探测欺骗目标的实时位置速度信息；欺骗式GNSS信号源根据前两个系统的输入信息，采用仿真模型计算出被欺骗目标上的GNSS接收机天线口面处的伪距、载波相位、信号功

率、导航电文等GNSS信号数据，然后进行信号调制、上变频等处理，生成GNSS欺骗RF信号，以强功率向被欺骗目标播发；被欺骗目标接收欺骗信号，在欺骗信号的引导下，导向错误位置。

GNSS信号接收天线

(()

探测信号接收天线

欺骗信号发射天线

目标运动信息探测系统

星历/时间注入接收机

目标位置速度

实时星历/时间欺骗式干扰信号源

GNSS欺骗系统

GNSS欺骗信号

图1、GNSS欺骗式干扰系统示意图

2欺骗GNSS信号参数仿真模型

欺骗GNSS信号参数仿真的数学模型及其信息流程图如下：

欺骗位置、速度仿

欺骗位置、速度仿真模型

GNSs实时星历

卫星轨道仿真模型

大气参数仿真模型

对流层延迟仿真模型

电离层延退仿真模型

卫星钟差模型

相对论效应模型

时空丽统模型

导航电文

欺骗GNSS信号参数数学仿真

真距仿真模型

日标位置遵度

实时时间

图2欺骗GNSS信号参数仿真模型信息流

2.1时空系统仿真模型

(1)UTC与导航系统时间转换模型

考虑模拟器的本地时间(年月日时分秒)和模拟器所在的时区，将本地时间应转化为格林尼治零时区的时间，再将格林尼治零时区的时间转化为导航系统时。以GPS系统时GPST为例，UTC转导GPST的模型如下：

如果(y-1980)%4=0或m≤2,则D?-1

tsec=(d?%7)×86400+(h,-Tinezone)×3600+mm×60+s,+13

如果tsec≤0,则WN-1

tsec=tse+604800,Z_Count=INT(sec/6)+1

其中y、m、d、h、m、s为目前年、月、日、时、分、秒

tse为系统时的周内秒数%为模运算

d。(12)={0,31,59,90,120,151,181,212,243,273,304,334}Timezom为模拟器模拟的时区

(2)坐标系统

数学仿真中，将涉及的坐标系统主要包括：BDS参考坐标系统CGS2000、GPS参考坐标系统WGS84。坐标系之间采用7参数转换模型：

式中，[xYz]si、[xYz]s?为不同坐标系下的坐标；k为坐标系比例尺系数；α、β、γ分别为两个坐标系的夹角；[△X△Y△Z]]为两个坐标系坐标原点偏差。

2.2导航卫星轨道仿真模型

卫星轨道仿真需要计算导航卫星在地固坐标系下的位置和速度。卫星轨道仿真模型采用基于广播星历轨道根数的运动学仿真。

A.卫星位置仿真模型算法

卫星轨道位置仿真模型算法如下：

采用WGS-84坐标系的地球引力常数：GM=H=3986005×1014(m3/s2)(1)计算卫星运行的平

采用WGS-84坐标系的地球引力常数：GM=H=3986005×1014(m3/s2)

为地球引力常数

(2)计算t时刻卫星的平近点角

M(t)=Mo+n.(t-to);

M?为星历设置中给出的参考时刻t。的平近点角；M的单位为弧度。n的单位为弧度/秒。

(3)计算偏近点角

E(t)=M(t)+e×sinE(t)

观测历元t的平近点角M已经得出，e由星历中给出，其中E的单位为弧度。

(4)计算真近点角

真近点角的单位为弧度。

(5)计算升交距角(未经改正的)Φ(t)=θ(t)+w

w为轨道近地点角距，由星历设定。

(6)计算卫星向径r(t)=A(1-e·cosE(t))

(7)计算摄动改正项

δu(t)=Cc·cos(2u(t))+C·sin(2u(t))δr(t)=C·cos(2u(t))+C,·sin(2u(t))

δi(t)=Cic·cos(2u(t))+C·sin(2u(t))

(8)进行摄动改正

u(t)=u(t)+δu(t)r(t)=r(t)+δr(t)

i(t)=i?+i·(t-toe)+δi(t)

(9)计算卫星在轨道平面坐标系中的位置

x(t)=r(t)·cosu(t)y(t)=r(t)·sinu(t)

(10)计算升交点经度

Ω(t)=Ω?+(Ω-w)·(t-toe)-weto

(11)计算卫星在地固坐标