CN104155663B - 一种应用于卫星导航系统的反欺骗攻击的方法及系统 - Google Patents

Abstract

本发明公开一种应用于卫星导航系统的反欺骗攻击的方法及系统，所述方法包括：S1.对卫星信号进行捕获并处理，获取所述卫星信号的码多普勒信息和载波多普勒信息；S2.对所述码多普勒信息和载波多普勒信息进行一致性检测，得到一致性比对结果；S3.根据所述一致性比对结果，判断所述卫星信号是否为欺骗信号，若是，则发出报警信息，若否，则将当前码相位和载波频率信息送入预设的跟踪环路。本发明的方法及系统不仅能判断接收到的信号中是否存在欺骗信号，而且还能判断出哪个是欺骗信号、哪个是真实信号，即可达到第二层次，具有欺骗信号抑制功能，具有较强的反欺骗能力。

Description

一种应用于卫星导航系统的反欺骗攻击的方法及系统

技术领域

本发明涉及卫星导航技术领域，具体涉及一种应用于卫星导航系统的反欺骗攻击的方法。

背景技术

卫星导航系统，包括GPS、GLONASS、Galileo、北斗等，由于能提供精确的时间和位置信息，已被广泛运用于交通、通信、电力、航空、航天等领域。然而，由于卫星播发的导航信号功率微弱、无线信道十分复杂以及信号结构公开等因素的影响，卫星导航系统极易受到干扰和欺骗攻击。欺骗攻击能在不被察觉的情况下使目标接收设备输出错误的位置、时间信息，因此欺骗攻击对卫星导航系统的安全使用构成了巨大威胁。

反欺骗攻击的方法可以分为两个层次。第一次层次是欺骗检测，即判断所接收到的信号中是否存在欺骗信号；第二层次是欺骗抑制，即识别出所接收信号中哪个是欺骗信号、哪个是真实信号。

现有的反欺骗攻击方法，如“多峰检测”与“信号质量监测”，一般仅能达到反欺骗攻击第一层次，还无法达到第二层次。

发明内容

本发明所要解决的技术问题是现有的反欺骗攻击的方法仅能判断接收到的信号中是否存在欺骗信号，无法识别出哪个是欺骗信号、哪个是真实信号。

为此目的，第一方面，本发明提供一种应用于卫星导航系统的反欺骗攻击的方法，该方法包括：

S1.对卫星信号进行捕获并处理，获取所述卫星信号的码多普勒信息和载波多普勒信息；

S2.对所述码多普勒信息和载波多普勒信息进行一致性检测，得到一致性比对结果；

S3.根据所述一致性比对结果，判断所述卫星信号是否为欺骗信号，若是，则发出报警信息，若否，则将当前码相位和载波频率信息送入预设的跟踪环路。

可选地，所述步骤S1包括：

S11.对卫星信号进行捕获并处理，并将处理后得到的相关值与预设捕获门限进行对比，记录所有大于所述预设捕获门限的相关值所对应的码相位参数、载波多普勒参数以及当前时刻；

S12.每隔预设时间Δt，重复执行步骤S11，得到集合Q，所述Q＝{(t_l(n),τ_l(n)，f_l(n))|1≤l≤L_n，1≤n≤N}；

其中，所述t_l(n)为第n次捕获对应的时刻，所述τ_l(n)为码相位参数，所述f_l(n)为载波多普勒参数；l为大于所述预设捕获门限的相关值数目下标，表示第l个大于所述预设捕获门限的相关值；L_n为第n次捕获中大于所述预设捕获门限的相关值数目，N为预设捕获次数；

S13.根据所述集合Q，通过霍夫变换HT，得到所述卫星信号的码多普勒信息和载波多普勒信息。

可选地，在步骤S13中，所述卫星信号的码多普勒信息和载波多普勒信息的获取步骤包括：

根据所述集合Q，得到子集合Q_τ，所述Q_τ＝{(t_l(n),τ_l(n))|1≤l≤L_n，1≤n≤N}；

根据预设码多普勒最大值d_max、预设码多普勒最小值d_min以及预设码多普勒分辨率Δd，得到初始码多普勒集合D，所述D为：

D＝{d_m|d_m＝d_min+m×Δd，0≤m<M}；

其中所述为不小于的最小整数；

根据所述子集合Q_τ，得到初始码相位集B＝{b(n,l,m)|1≤n≤N，1≤l≤L_n，0≤m<M}，所述b(n,l,m)为第n个初始码相位子集：b(n,l,m)＝-t_l(n)×d_m+τ_l(n)；

根据所述初始码相位集B的最大值b_max和最小值b_min以及预设初始码相位分辨率Δb，得到离散的初始码相位集：

$\hat{B}=\left\{{\hat{b}}_j\right|{\hat{b}}_j=b_{\min }+j\&times;\mathrm{\&Delta;b},0\&le;j<J\};$

其中 为不小于的最小整数；

根据所述初始码多普勒集合D、所述初始码相位集及离散的初始码相位集，得到卫星信号的码多普勒信息和载波多普勒信息。

第二方面，本发明还提出一种应用于卫星导航系统的反欺骗攻击的系统，所述系统包括：

信息获取模块，用于对卫星信号进行捕获并处理，获取所述卫星信号的码多普勒信息和载波多普勒信息；

一致性检测模块，用于对所述码多普勒信息和载波多普勒信息进行一致性检测，得到一致性比对结果；

判断模块，用于根据所述一致性比对结果，判断所述卫星信号是否为欺骗信号，若是，则发出报警信息，若否，则将当前码相位和载波频率信息送入预设的跟踪环路。

可选地，所述信息获取模块包括：

信号处理单元，用于对卫星信号进行捕获并处理，并将处理后得到的相关值与预设捕获门限进行对比，记录所有大于所述预设捕获门限的相关值所对应的码相位参数、载波多普勒参数以及当前时刻；

每隔预设时间Δt，重复执行步骤S11，得到集合Q，所述Q＝{(t_l(n),τ_l(n)，f_l(n))|1≤l≤L_n，1≤n≤N}；

其中，所述t_l(n)为第n次捕获对应的时刻，所述τ_l(n)为码相位参数，所述f_l(n)为载波多普勒参数；l为大于所述预设捕获门限的相关值数目下标，表示第l个大于所述预设捕获门限的相关值；L_n为第n次捕获中大于所述预设捕获门限的相关值数目，N为预设捕获次数；

信息获取单元，用于根据所述集合Q，通过霍夫变换HT，得到所述卫星信号的码多普勒信息和载波多普勒信息。

可选地，在所述信息获取单元中，所述卫星信号的码多普勒信息和载波多普勒信息的获取步骤包括：

根据所述集合Q，得到子集合Q_τ，所述Q_τ＝{(t_l(n),τ_l(n))|1≤l≤L_n，1≤n≤N}；

根据预设码多普勒最大值d_max、预设码多普勒最小值d_min以及预设码多普勒分辨率Δd，得到初始码多普勒集合D，所述D为：

D＝{d_m|d_m＝d_min+m×Δd，0≤m<M}；

其中所述为不小于的最小整数；

根据所述子集合Q_τ，得到初始码相位集B＝{b(n,l,m)|1≤n≤N，1≤l≤L_n，0≤m<M}，所述b(n,l,m)为第n个初始码相位子集：b(n,l,m)＝-t_l(n)×d_m+τ_l(n)；

根据所述初始码相位集B的最大值b_max和最小值b_min以及预设初始码相位分辨率Δb，得到离散的初始码相位集：

$\hat{B}=\left\{{\hat{b}}_j\right|{\hat{b}}_j=b_{\min }+j\&times;\mathrm{\&Delta;b},0\&le;j<J\};$

其中 为不小于的最小整数；

根据所述初始码多普勒集合D、所述初始码相位集及离散的初始码相位集，得到卫星信号的码多普勒信息和载波多普勒信息。

相比于现有技术，本发明的应用于卫星导航系统的反欺骗攻击的方法及系统通过对接收到的卫星信号进行多次捕获以获取信号码多普勒和载波多普勒信息，然后通过检测码多普勒和载波多普勒的一致性来判断所接收到的信号是欺骗信号还是真实信号。若为欺骗信号则发出报警信息，若为真实信号则将当前码相位和载波频率信息送入跟踪环路以进行正常信号跟踪，可见，本发明的方法及系统不仅能判断接收到的信号中是否存在欺骗信号，而且还能判断出哪个是欺骗信号、哪个是真实信号(即可达到第二层次)，具有欺骗信号抑制功能，具有较强的反欺骗能力。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1示出了实施例一的应用于卫星导航系统的反欺骗攻击的方法流程图；

图2示出了实施例一中对卫星信号进行捕获并处理的流程图；

图3示出了真实信号信噪比为-21dB，不同欺信比以及不同欺骗信号个数时的欺骗信号检测率示意图；

图4示出了在信噪比为-21dB，欺信比为0dB，不同捕获次数N以及不同欺骗信号个数下的欺骗信号检测率示意图；

图5示出了实施例二的应用于卫星导航系统的反欺骗攻击的系统结构图。

具体实施方式

为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

实施例一：

本实施例公开一种应用于卫星导航系统的反欺骗攻击的方法，如图1所示，该方法包括：

S1.对卫星信号进行捕获并处理，获取所述卫星信号的码多普勒信息和载波多普勒信息；

S2.对所述码多普勒信息和载波多普勒信息进行一致性检测，得到一致性比对结果；

S3.根据所述一致性比对结果，判断所述卫星信号是否为欺骗信号，若是，则发出报警信息，若否，则将当前码相位和载波频率信息送入预设的跟踪环路。

所述步骤S1包括：

S11.对卫星信号进行捕获并处理，并将处理后得到的相关值与预设捕获门限进行对比，记录所有大于所述预设捕获门限的相关值所对应的码相位参数、载波多普勒参数以及当前时刻；

S12.每隔预设时间Δt，重复执行步骤S11，得到集合Q，所述Q＝{(t_l(n),τ_l(n)，f_l(n))|1≤l≤L_n，1≤n≤N}；

其中，所述t_l(n)为第n次捕获对应的时刻，所述τ_l(n)为码相位参数，所述f_l(n)为载波多普勒参数；l为大于所述预设捕获门限的相关值数目下标，表示第l个大于所述预设捕获门限的相关值；L_n为第n次捕获中大于所述预设捕获门限的相关值数目，N为预设捕获次数；

S13.根据所述集合Q，通过霍夫变换HT，得到所述卫星信号的码多普勒信息和载波多普勒信息。

具体地，如图2所示，对经接收机前端处理后的中频信号进行A/D采样并进行正常的伪码捕获(伪码捕获可采用现有的伪码捕获模块来实现)，将所得相关值(所述相关值是伪码捕获时本地信号与接收信号进行相关所产生的数据，在信号处理领域都称之为相关值)与预设捕获门限进行对比，记录所有大于所述预设捕获门限的相关值对应的码相位和载波多普勒参数。经过预设时间Δt间隔时间后，再次进行正常的伪码捕获，并记录当前所有大于预设捕获门限的相关值对应的码相位和载波多普勒参数。每过预设时间Δt间隔再次进行正常的伪码捕获，如此反复进行N次，并记录下每次大于预设捕获门限的相关值对应的码相位和载波多普勒参数。这样将得到一个捕获时刻、码相位和载波多普勒参数集合Q。

在步骤S13中，所述卫星信号的码多普勒信息和载波多普勒信息的获取步骤包括：

根据所述集合Q，得到子集合Q_τ，所述Q_τ＝{(t_l(n),τ_l(n))|1≤l≤L_n，1≤n≤N}；

根据预设码多普勒最大值d_max、预设码多普勒最小值d_min以及预设码多普勒分辨率Δd，得到初始码多普勒集合D，所述D为：

D＝{d_m|d_m＝d_min+m×Δd，0≤m<M}；

其中所述为不小于的最小整数；

根据所述子集合Q_τ，得到初始码相位集B＝{b(n,l,m)|1≤n≤N，1≤l≤L_n，0≤m<M}，所述b(n,l,m)为第n个初始码相位子集：b(n,l,m)＝-t_l(n)×d_m+τ_l(n)；

根据所述初始码相位集B的最大值b_max和最小值b_min以及预设初始码相位分辨率Δb，得到离散的初始码相位集：

$\hat{B}=\left\{{\hat{b}}_j\right|{\hat{b}}_j=b_{\min }+j\&times;\mathrm{\&Delta;b},0\&le;j<J\};$

其中 为不小于的最小整数；

根据所述初始码多普勒集合D、所述初始码相位集及离散的初始码相位集，得到卫星信号的码多普勒信息和载波多普勒信息。

具体地，构造一个M行J列的二维累加矩阵A＝[a_mj]_M×J，并将该矩阵初始化为0，对于初始码相位集B中的每一个元素b(n,l,m)，求得离散初始码相位集中与其最接近的元素对应的下标j，并将二维累加矩阵A中的元素a_mj加一。按照以上方式遍历初始码相位集B中的每一个元素后，完成对二维累加矩阵A的更新。

按照列优先的原则搜索二维累加矩阵A，求得其最大的L个值，表示所有捕获中大于预设捕获门限的信号数目的最大值，L_n为第n次捕获中大于所述预设捕获门限的相关值数目，并记录对应的行下标和列下标，然后按如下方法求得对应的码多普勒和初始码相位：若二维累加矩阵A中第m行j列元素a_mj为最大值，其行下标为m则对应的码多普勒记为码多普勒集合D中元素d_m，其列下标为j则对应的初始码相位为离散的初始码相位集中元素这样可求得L个信号的码多普勒参数和初始码相位参数，并将其分别记为码多普勒矢量初始码相位矢量

$\stackrel{\&RightArrow;}{r}=[r_1,r_2,...,r_L];$

${\stackrel{\&RightArrow;}{\mathrm{\&tau;}}}_{\mathrm{init}}=[{\mathrm{\&tau;}}_1,{\mathrm{\&tau;}}_2,...,{\mathrm{\&tau;}}_L].$

其中，码多普勒矢量的L个元素表示L个信号的码多普勒。利用求得的码多普勒矢量和初始码相位矢量在捕获起始时刻和码相位构成的二维空间中拟合出多条直线，并按照同一条直线上对应的码相位和载波多普勒表示同一信号的原则，将码相位和载波多普勒参数集合Q划分为L个子集。每一个子集的元素表示同一信号不同捕获时刻对应的码相位和载波多普勒。对这L个子集中的第l个子集，对其N个载波多普勒f_l(n)(1≤n≤N)求解均值得到该信号对应的平均载波多普勒

${\stackrel{\&OverBar;}{f}}_l=\frac{1}{N}\underset{n=1}{\overset{N}{\mathrm{\&Sigma;}}}{f}_l\left(n\right),1\&le;l\&le;L$

这样，可以求得L个信号各自的平均载波多普勒。所求得的L个平均载波多普勒构成了载波多普勒矢量：

$\stackrel{\&RightArrow;}{f}=[{\stackrel{\&OverBar;}{f}}_1,{\stackrel{\&OverBar;}{f}}_2,...,{\stackrel{\&OverBar;}{f}}_L].$

载波多普勒矢量的L个元素表示L个信号的载波多普勒。

至此，完成对所有信号的码多普勒和载波多普勒的计算。

在获得所有信号的码多普勒和载波多普勒以后，按如下方法进行欺骗信号检测与抑制，并完成正常伪码捕获工作：

首先，根据接收信号标称的载波频率f_carrier和伪码速率f_code计算比例系数δ

$\mathrm{\&delta;}=\frac{f_{\mathrm{carrier}}}{f_{\mathrm{code}}};$

其次，利用比例系数δ对码多普勒矢量中的L个元素分别构造检测统计量，构造公式如下：

$T_l=r_l*\mathrm{\&delta;}-{\stackrel{\&OverBar;}{f}}_l,1\&le;l\&le;L$

再次，假定预设的欺骗信号检测判决门限为V_JD，如果T_l≥V_JD，则将第l个信号判为欺骗信号并发出告警，进而实现对欺骗信号的检测与抑制；如果T_l＜V_JD，则将第l个信号判为真实信号，并将该信号对应的伪码相位和载波频率参数送给预设的跟踪环路，进而完成正常的伪码捕获工作。

在具体应用中，如图2所示的对卫星信号进行捕获并处理的流程如下：

步骤1：初始化捕获次数n＝1。

步骤2：对经接收终端前端处理后的中频信号进行A/D采样，将所得信号与本地载波和伪码进行相关，得到同相支路相关积分值I和正交支路相关积分值Q。

步骤3：对同相支路相关积分值I和正交支路相关积分值Q进行非相干累加，得到非相干累加值I²+Q²。将所得非相干累加值与预设的捕获门限行比较，如果超过门限，则记录下当前时刻t₁(1)、码相位参数τ₁(1)和载波多普勒参数f₁(1))，并记为(t₁(1)，τ₁(1)，f₁(1))，然后则转入步骤4；否则，调整当前码相位和载波频率，然后转入步骤2继续搜索下一个码相位和载波频率格。

步骤4：判断是否搜索完所有的码相位和载波频率格：如果已完成所有搜素，则转入步骤5；否则，调整当前码相位和载波频率，然后转入步骤2继续搜索后续码相位和载波频率格，并记录下其它大于预设捕获门限的相关值对应的当前时刻t_l(1)、码相位参数τ_l(1)和载波多普勒参数f_l(1)，直至搜索完所有的码相位与载波频率格。

步骤5：经过Δt间隔时间后，将捕获次数n加1，并通过重复步骤2至步骤4再次进行正常的伪码捕获，并记录下所有超过捕获门限的相关值对应的当前时刻t_l(n)、码相位参数τ_l(n)和载波多普勒参数f_l(n)。

步骤6：每经过Δt间隔时间后，重复步骤5，直至捕获次数n＝N。至此，得到一个由捕获时刻、码相位和载波多普勒参数构成的集合Q，Q＝{(t_l(n),τ_l(n)，f_l(n))|1≤l≤L_n，1≤n≤N}；

其中，所述t_l(n)为第n次捕获对应的时刻，所述τ_l(n)为码相位参数，所述f_l(n)为载波多普勒参数；l为大于所述预设捕获门限的相关值数目下标，表示第l个大于所述预设捕获门限的相关值；L_n为第n次捕获中大于所述预设捕获门限的相关值数目，N为预设捕获次数。利用参数表示所有捕获中大于预设捕获门限的信号数目的最大值。

步骤7：计算所有信号的码多普勒和载波多普勒。

步骤8：进行欺骗信号检测与抑制，并完成正常伪码捕获工作。

本实施例的应用于卫星导航系统的反欺骗攻击方法与信号质量监测方法相比，首先是捕获阶段，所述捕获阶段是卫星导航系统基带信号处理的第一个阶段，相比于在跟踪阶段实现的信号质量监测方法，本实施例的方法可以更早地发现信号异常——当存在欺骗信号时，可以更早地发出欺骗告警。同时，本实施例的方法可利用卫星导航接收终端已有的捕获模块实现，无需增加额外的硬件成本。

本实施例的方法与多峰检测方法相比，本实施例的方法是根据接收信号的码多普勒和载波多普勒的一致性来检测、抑制欺骗信号。对于真实信号，码多普勒和载波多普勒都是由于接收机和卫星间的相对运动产生的，满足一致性关系；而对于欺骗信号，其码相位和载波多普勒一致性一般难以满足。因此，该方法不但能检测出是否存在欺骗信号，而且能识别出哪个是欺骗信号、哪个是真实信号，可及时发出告警，具有欺骗抑制能力。

图3给出了真实信号信噪比为-21dB，不同欺信比(欺信比表示欺骗信号与真实信号功率比)以及不同欺骗信号个数时的欺骗信号检测率(即正确地检测出欺骗信号的概率)。可以看出，当欺骗信号和真实信号功率差在一定范围内时，该方法能保持较高的欺骗信号检测率。同时，随着欺骗信号个数的增多，要获得相同的欺骗信号检测率所要求的欺骗信号和真实信号间的功率差将进一步减小。

图4给出了在信噪比为-21dB，欺信比为0dB，不同捕获次数N以及不同欺骗信号个数下的欺骗信号检测率。可以看出，随着捕获次数的增加，欺骗信号识别正确率呈上升趋势。即当欺骗信号增多时，要达到相同的欺骗信号检测率所需要的捕获次数也将增加。要使欺骗信号检测率达到0.9以上，当存在一个欺骗信号时所需捕获次数应大于14次，当存在4个欺骗信号时所需捕获次数应大于43次。

实施例二：

本实施例公开一种应用于卫星导航系统的反欺骗攻击的系统，如图5所示，所述系统可包括：

信息获取模块，用于对卫星信号进行捕获并处理，获取所述卫星信号的码多普勒信息和载波多普勒信息；

一致性检测模块，用于对所述码多普勒信息和载波多普勒信息进行一致性检测，得到一致性比对结果；

判断模块，用于根据所述一致性比对结果，判断所述卫星信号是否为欺骗信号，若是，则发出报警信息，若否，则将当前码相位和载波频率信息送入预设的跟踪环路。

具体应用中，所述信息获取模块可包括图5中未示出的以下单元：

信号处理单元，用于对卫星信号进行捕获并处理，并将处理后得到的相关值与预设捕获门限进行对比，记录所有大于所述预设捕获门限的相关值所对应的码相位参数、载波多普勒参数以及当前时刻；

每隔预设时间Δt，重复执行步骤S11，得到集合Q，所述Q＝{(t_l(n),τ_l(n)，f_l(n))|1≤l≤L_n，1≤n≤N}；

其中，所述t_l(n)为第n次捕获对应的时刻，所述τ_l(n)为码相位参数，所述f_l(n)为载波多普勒参数；l为大于所述预设捕获门限的相关值数目下标，表示第l个大于所述预设捕获门限的相关值；L_n为第n次捕获中大于所述预设捕获门限的相关值数目，N为预设捕获次数；

信息获取单元，用于根据所述集合Q，通过霍夫变换HT，得到所述卫星信号的码多普勒信息和载波多普勒信息。

在所述信息获取单元中，所述卫星信号的码多普勒信息和载波多普勒信息的获取步骤包括：

根据所述集合Q，得到子集合Q_τ，所述Q_τ＝{(t_l(n),τ_l(n))|1≤l≤L_n，1≤n≤N}；

根据预设码多普勒最大值d_max、预设码多普勒最小值d_min以及预设码多普勒分辨率Δd，得到初始码多普勒集合D，所述D为：

D＝{d_m|d_m＝d_min+m×Δd，0≤m<M}；

其中所述为不小于的最小整数；

根据所述子集合Q_τ，得到初始码相位集B＝{b(n,l,m)|1≤n≤N，1≤l≤L_n，0≤m<M}，所述b(n,l,m)为第n个初始码相位子集：b(n,l,m)＝-t_l(n)×d_m+τ_l(n)；

根据所述初始码相位集B的最大值b_max和最小值b_min以及预设初始码相位分辨率Δb，得到离散的初始码相位集：

$\hat{B}=\left\{{\hat{b}}_j\right|{\hat{b}}_j=b_{\min }+j\&times;\mathrm{\&Delta;b},0\&le;j<J\};$

其中 为不小于的最小整数；

根据所述初始码多普勒集合D、所述初始码相位集及离散的初始码相位集，得到卫星信号的码多普勒信息和载波多普勒信息。

与现有技术相比，本实施例的应用于卫星导航系统的反欺骗攻击系统可以更早地发现信号异常——当存在欺骗信号时，可以更早地发出欺骗告警。同时，本实施例的系统可利用卫星导航接收终端已有的捕获模块实现，无需增加额外的硬件成本。

本实施例的系统根据接收信号的码多普勒和载波多普勒的一致性来检测、抑制欺骗信号。对于真实信号，码多普勒和载波多普勒都是由于接收机和卫星间的相对运动产生的，满足一致性关系；而对于欺骗信号，其码相位和载波多普勒一致性一般难以满足。因此，该系统不但能检测出是否存在欺骗信号，而且能识别出哪个是欺骗信号、哪个是真实信号，可及时发出告警，具有欺骗抑制能力。

虽然结合附图描述了本发明的实施方式，但是本领域技术人员可以在不脱离本发明的精神和范围的情况下做出各种修改和变型，这样的修改和变型均落入由所附权利要求所限定的范围之内。

Claims (6)

1.一种应用于卫星导航系统的反欺骗攻击的方法，其特征在于，该方法包括：

S1.对卫星信号进行捕获并处理，获取所述卫星信号的码多普勒信息和载波多普勒信息；

S2.对所述码多普勒信息和载波多普勒信息进行一致性检测，得到一致性比对结果；

在获得所有信号的码多普勒和载波多普勒以后，按如下方法进行欺骗信号检测与抑制，并完成正常伪码捕获工作：

首先，根据接收信号标称的载波频率f_carrier和伪码速率f_code计算比例系数δ：

$\mathrm{\&delta;}=\frac{f_{carrier}}{f_{code}};$

其次，利用比例系数δ对码多普勒矢量中的L个元素分别构造检测统计量，构造公式如下：

$T_l=r_l*\mathrm{\&delta;}-{\stackrel{\&OverBar;}{f}}_l,1\&le;l\&le;L$

再次，假定预设的欺骗信号检测判决门限为V_JD，如果T_l≥V_JD，则将第l个信号判为欺骗信号并发出告警，进而实现对欺骗信号的检测与抑制；如果T_l＜V_JD，则将第l个信号判为真实信号，并将该信号对应的伪码相位和载波频率参数送给预设的跟踪环路，进而完成正常的伪码捕获工作；

其中，为第l个信号对应的平均载波多普勒：

${\stackrel{\&OverBar;}{f}}_l=\frac{1}{N}\underset{n=1}{\overset{N}{\&Sigma;}}{f}_l\left(n\right),1\&le;l\&le;L$

N为预设捕获次数，f_l(n)(1≤n≤N)为第l个信号对应的第n次捕获所得载波多普勒；

S3.根据所述一致性比对结果，判断所述卫星信号是否为欺骗信号，若是，则发出报警信息，若否，则将当前码相位和载波频率信息送入预设的跟踪环路。

2.根据权利要求1所述的方法，其特征还在于，所述步骤S1包括：

S11.对卫星信号进行捕获并处理，并将处理后得到的相关值与预设捕获门限进行对比，记录所有大于所述预设捕获门限的相关值所对应的码相位参数、载波多普勒参数以及当前时刻；

S12.每隔预设时间Δt，重复执行步骤S11，得到集合Q，所述Q＝{(t_l(n),τ_l(n)，f_l(n))|1≤l≤L_n，1≤n≤N}；

其中，所述t_l(n)为第n次捕获对应的时刻，所述τ_l(n)为码相位参数，所述f_l(n)为载波多普勒参数；l为大于所述预设捕获门限的相关值数目下标，表示第l个大于所述预设捕获门限的相关值；L_n为第n次捕获中大于所述预设捕获门限的相关值数目，N为预设捕获次数；

S13.根据所述集合Q，通过霍夫变换HT，得到所述卫星信号的码多普勒信息和载波多普勒信息。

3.根据权利要求2所述的方法，其特征在于，在步骤S13中，所述卫星信号的码多普勒信息和载波多普勒信息的获取步骤包括：

根据所述集合Q，得到子集合Q_τ，所述Q_τ＝{(t_l(n),τ_l(n))|1≤l≤L_n，1≤n≤N}；

根据预设码多普勒最大值d_max、预设码多普勒最小值d_min以及预设码多普勒分辨率Δd，得到初始码多普勒集合D，所述D为：

D＝{d_m|d_m＝d_min+m×Δd，0≤m<M}；

其中所述为不小于的最小整数；

根据所述子集合Q_τ，得到初始码相位集B＝{b(n,l,m)|1≤n≤N，1≤l≤L_n，0≤m<M}，所述b(n,l,m)为第n个初始码相位子集：b(n,l,m)＝-t_l(n)×d_m+τ_l(n)；

根据所述初始码相位集B的最大值b_max和最小值b_min以及预设初始码相位分辨率Δb，得到离散的初始码相位集：

$\hat{B}=\left\{{\hat{b}}_j\right|{\hat{b}}_j=b_{min}+j\&times;\mathrm{\&Delta;}b,0\&le;j<J\};$

其中为不小于的最小整数；

根据所述初始码多普勒集合D、所述初始码相位集及离散的初始码相位集，得到卫星信号的码多普勒信息和载波多普勒信息。

4.一种应用于卫星导航系统的反欺骗攻击的系统，其特征在于，所述系统包括：

信息获取模块，用于对卫星信号进行捕获并处理，获取所述卫星信号的码多普勒信息和载波多普勒信息；

一致性检测模块，用于对所述码多普勒信息和载波多普勒信息进行一致性检测，得到一致性比对结果；

在获得所有信号的码多普勒和载波多普勒以后，按如下方法进行欺骗信号检测与抑制，并完成正常伪码捕获工作：

首先，根据接收信号标称的载波频率f_carrier和伪码速率f_code计算比例系数δ：

$\mathrm{\&delta;}=\frac{f_{carrier}}{f_{code}};$

其次，利用比例系数δ对码多普勒矢量中的L个元素分别构造检测统计量，构造公式如下：

$T_l=r_l*\mathrm{\&delta;}-{\stackrel{\&OverBar;}{f}}_l,1\&le;l\&le;L$

再次，假定预设的欺骗信号检测判决门限为V_JD，如果T_l≥V_JD，则将第l个信号判为欺骗信号并发出告警，进而实现对欺骗信号的检测与抑制；如果T_l＜V_JD，则将第l个信号判为真实信号，并将该信号对应的伪码相位和载波频率参数送给预设的跟踪环路，进而完成正常的伪码捕获工作；

其中，为第l个信号对应的平均载波多普勒：

${\stackrel{\&OverBar;}{f}}_l=\frac{1}{N}\underset{n=1}{\overset{N}{\&Sigma;}}{f}_l\left(n\right),1\&le;l\&le;L$

N为预设捕获次数，f_l(n)(1≤n≤N)为第l个信号对应的第n次捕获所得载波多普勒；

判断模块，用于根据所述一致性比对结果，判断所述卫星信号是否为欺骗信号，若是，则发出报警信息，若否，则将当前码相位和载波频率信息送入预设的跟踪环路。

5.根据权利要求4所述的系统，其特征在于，所述信息获取模块包括：

信号处理单元，用于对卫星信号进行捕获并处理，并将处理后得到的相关值与预设捕获门限进行对比，记录所有大于所述预设捕获门限的相关值所对应的码相位参数、载波多普勒参数以及当前时刻；

每隔预设时间Δt，重复执行步骤S11，得到集合Q，所述Q＝{(t_l(n),τ_l(n)，f_l(n))|1≤l≤L_n，1≤n≤N}；

其中，所述t_l(n)为第n次捕获对应的时刻，所述τ_l(n)为码相位参数，所述f_l(n)为载波多普勒参数；l为大于所述预设捕获门限的相关值数目下标，表示第l个大于所述预设捕获门限的相关值；L_n为第n次捕获中大于所述预设捕获门限的相关值数目，N为预设捕获次数；

信息获取单元，用于根据所述集合Q，通过霍夫变换HT，得到所述卫星信号的码多普勒信息和载波多普勒信息。

6.根据权利要求5所述的系统，其特征在于，在所述信息获取单元中，所述卫星信号的码多普勒信息和载波多普勒信息的获取步骤包括：

根据所述集合Q，得到子集合Q_τ，所述Q_τ＝{(t_l(n),τ_l(n))|1≤l≤L_n，1≤n≤N}；

根据预设码多普勒最大值d_max、预设码多普勒最小值d_min以及预设码多普勒分辨率Δd，得到初始码多普勒集合D，所述D为：

D＝{d_m|d_m＝d_min+m×Δd，0≤m<M}；

其中所述为不小于的最小整数；

根据所述子集合Q_τ，得到初始码相位集B＝{b(n,l,m)|1≤n≤N，1≤l≤L_n，0≤m<M}，所述b(n,l,m)为第n个初始码相位子集：b(n,l,m)＝-t_l(n)×d_m+τ_l(n)；

根据所述初始码相位集B的最大值b_max和最小值b_min以及预设初始码相位分辨率Δb，得到离散的初始码相位集：

$\hat{B}=\left\{{\hat{b}}_j\right|{\hat{b}}_j=b_{\min }+j\&times;\mathrm{\&Delta;}b,0\&le;j<J\};$

其中为不小于的最小整数；

根据所述初始码多普勒集合D、所述初始码相位集及离散的初始码相位集，得到卫星信号的码多普勒信息和载波多普勒信息。