1.基于改进Ratio的多径环境下的欺骗干扰检测方法,其特征在于,包括如下步骤:步骤1.通过差值形式构建DIQR检测量,即计算超前支路与滞后支路的相关值之和以及与即时码支路相关值之间的差值,并将差值结果作为DIQR检测量;
步骤2.设定第一检测门限,第一检测门限包括第一检测门限的下限和第一检测门限的上限;对步骤1得到的DIQR检测量采取加窗处理措施,得到多个第一检测窗;
将每个第一检测窗时间内的DIQR检测量的数量,作为第一预设数量;
统计每个第一检测窗时间内DIQR检测量低于第一检测门限的下限以及高于第一检测门限的上限的数量之和,并将二者的数量之和定义为第一数量;
若当前第一检测窗时间内的第一数量与第一预设数量比值超过第一预设阈值,则判断该第一检测窗内存在干扰;否则,判断该第一检测窗不存在干扰;
步骤3.针对步骤1得到的DIQR检测量设定一个固定长度的滑动窗口,提取DIQR检测量的子集,计算子集中数值的平方平均值与子集均值平方的差值,得到当前子集的方差;
滑动窗口以预设的固定步长向前移动,选取新的数据子集,并重复上述方差计算过程,在整个DIQR检测量上依次进行,最终生成一系列方差值,即DIQR_MV检测量;
步骤4.设定第二检测门限,第二检测门限包括第二检测门限的下限和第二检测门限的上限;对步骤3得到的整个DIQR_MV检测量采取加窗处理措施,得到多个第二检测窗;
其中第二检测窗与第一检测窗在时间序列上是一一对应的;
针对步骤2中存在干扰的第一检测窗,选取与存在干扰的第一检测窗在时间序列上一一对应的第二检测窗,进一步利用DIQR_MV检测量进行干扰判断,过程如下:将对应第二检测窗时间内的DIQR_MV检测量的数量,作为第二预设数量;
统计当前第二检测窗时间内DIQR_MV检测量低于第二检测门限的下限以及高于第二检测门限的上限的数量之和,并将二者的数量之和定义为第二数量;
若当前第二检测窗时间内第二数量与第二预设数量比值超过第二预设阈值,则判断该第二检测窗内存在欺骗干扰;否则,判断该第二检测窗内存在多径干扰。
2.根据权利要求1所述的基于改进Ratio的多径环境下的欺骗干扰检测方法,其特征在于,所述步骤1中,DIQR检测量的计算公式为:DIQR(k)=E(k)+L(k)‑αP(k) (1)式中,E(k)、L(k)和P(k)分别为在第k个相干积分内相关器支路的超前、滞后和即时输出,α代表互相关函数在其峰值处的斜率,用于调节DIQR检测量中即时码输出的影响程度。
3.根据权利要求2所述的基于改进Ratio的多径环境下的欺骗干扰检测方法,其特征在于,所述步骤1中,E(k)、L(k)和P(k)的计算过程如公式(2)至(4)所示;
式中Ie(k)、Qe(k)分别为在第k个相干积分内相关器同相支路、正交支路的超前输出;
式中Il(k)、Ql(k)分别为在第k个相干积分内相关器同相支路、正交支路的滞后输出;
式中Ip(k)、Qp(k)分别为在第k个相干积分内相关器同相支路、正交支路的即时输出。
4.根据权利要求2所述的基于改进Ratio的多径环境下的欺骗干扰检测方法,其特征在于,所述步骤2中,第一检测门限值的确定方法如下:定义Pf1表示第一虚警概率,Pf1的计算公式如公式(5)所示;
其中T2为第一检测门限的下限,T1为第一检测门限的上限,num(·)为满足括号内条件的DIQR检测量个数,M为无干扰的检测样本个数;
第一检测门限值根据DIQR检测量的分布特性和预设的第一虚警概率Pf1计算得出。
5.根据权利要求1所述的基于改进Ratio的多径环境下的欺骗干扰检测方法,其特征在于,所述步骤3中,DIQR_MV检测量的求解过程如下:定义第n个滑动窗口的移动方差DIQR_MV为DIQR_MV(n),其表达式为:式中,DIQR(i)为第i个DIQR检测量,w为移动窗口长度,m为滑动区间大小,定义N表示数据处理过程中滑动窗口将覆盖的总次数,则n=1,2,...,N。
6.根据权利要求5所述的基于改进Ratio的多径环境下的欺骗干扰检测方法,其特征在于,所述步骤4中,第二检测门限值的确定方法如下:建立检验统计量MV,其表达式如下:
其中 表示方差,定义第二虚警概率Pf2,Pf2通过以下公式(9)计算得到;
其中 f(MV)为MV的概率密度函数;T4为第二检测门限的下限,T3为第二检测门限的上限,F(Tm,h;k)为第二检测门限的上限的累积分布函数,F(Tm,l;k)为第二检测门限的下限的累积分布函数;
第二检测门限值根据DIQR_MV检测量的分布特性和预设第二虚警概率Pf2计算得出。