一种GPS干扰检测方法技术

本发明专利技术公开了一种GPS干扰检测方法，包括步骤：一、射频信号x(t)接收；二、求取射频信号x(t)的协方差矩阵Rxx；三、对协方差矩阵Rxx进行酉变换，得到酉变换后的实数矩阵Rl；四、接收机首先求取实数矩阵Rl的逆矩阵然后求取矩阵的m次方，得到五、根据MUSIC算法构造空间谱函数；六、计算干扰波达方向角本发明专利技术方法步骤简单，运算量小，能够快速、准确地得到GPS干扰波达方向角，适用范围广，实用性强，便于推广使用。

【技术实现步骤摘要】
【专利摘要】本专利技术公开了一种GPS干扰检测方法，包括步骤：一、射频信号x(t)接收；二、求取射频信号x(t)的协方差矩阵Rxx；三、对协方差矩阵Rxx进行酉变换，得到酉变换后的实数矩阵Rl；四、接收机首先求取实数矩阵Rl的逆矩阵然后求取矩阵的m次方，得到五、根据MUSIC算法构造空间谱函数；六、计算干扰波达方向角本专利技术方法步骤简单，运算量小，能够快速、准确地得到GPS干扰波达方向角，适用范围广，实用性强，便于推广使用。【专利说明】—种GPS干扰检测方法
本专利技术属于GPS抗干扰
，具体涉及一种GPS干扰检测方法。
技术介绍
在GPS抗干扰技术中，检测干扰波达方向(DOA)是空时自适应抗干扰信号处理的一个重要内容，近几年来受到了广泛关注。学者提出了很多高分辨率子空间类算法，其中多重信号分类(MUSIC)算法是最重要的一种。这个算法利用子空间的正交特性能获得很好的分辨性能。但这种算法包含大量的复数运算，同时它还要进行特征分解，因此运算量比较大，硬件实现复杂。为了降低运算量，曾经有人提出了均匀线阵模式的酉变换实值化方法，但这种方法还是要进行特征分解；还有人提出了一种无特征分解的快速子空间DOA方法，但这种方法还是要进行大量的复数运算。
技术实现思路
本专利技术所要解决的技术问题在于针对上述现有技术中的不足，提供一种GPS干扰检测方法，其方法步骤简单，运算量小，能够快速、准确地得到GPS干扰波达方向角，适用范围广，实用性强，便于推广使用。为解决上述技术问题，本专利技术采用的技术方案是:一种GPS干扰检测方法，其特征在于该方法包括以下步骤:步骤一、射频信号接收:接收机通过天线阵接收射频信号x(t)；步骤二、求取射频信号的协方差矩阵，具体过程如下:步骤201、接收机对射频信号x(t)进行下变频，变为中频信号；步骤202、接收机对中频信号进行采样，得到采样时刻η对应的采样值矩阵χ (η)；步骤203、接收机根据公式:【权利要求】1.一种GPS干扰检测方法，其特征在于该方法包括以下步骤: 步骤一、射频信号接收:接收机通过天线阵接收射频信号x(t)； 步骤二、求取射频信号的协方差矩阵，具体过程如下: 步骤201、接收机对射频信号x(t)进行下变频，变为中频信号； 步骤202、接收机对中频信号进行采样，得到采样时刻η对应的采样值矩阵X (η)； 步骤203、接收机根据公式: 2.按照权利要求1所述的一种GPS干扰检测方法，其特征在于:步骤一中所述天线阵为均匀线阵。3.按照权利要求1所述的一种GPS干扰检测方法，其特征在于:步骤四中所述m为3~8的正整数。4.按照权利要求1所述的一种GPS干扰检测方法，其特征在于:步骤五中所述Θ的取值范围为-90°~90°。【文档编号】G01S19/21GK103926598SQ201410182106【公开日】2014年7月16日 申请日期:2014年5月2日 优先权日:2014年5月2日 【专利技术者】葛大江, 李文元, 梁茂国, 周光彬, 张亚妮 申请人:中国人民解放军西安通信学院本文档来自技高网...

【技术保护点】
一种GPS干扰检测方法，其特征在于该方法包括以下步骤：步骤一、射频信号接收：接收机通过天线阵接收射频信号x(t)；步骤二、求取射频信号的协方差矩阵，具体过程如下：步骤201、接收机对射频信号x(t)进行下变频，变为中频信号；步骤202、接收机对中频信号进行采样，得到采样时刻n对应的采样值矩阵x(n)；步骤203、接收机根据公式:Rxx≈1snΣn=1snx(n)xH(n)---(1)]]>计算采样值矩阵x(n)的协方差矩阵Rxx，公式(1)中，sn为快拍数，xH(n)为x(n)的共轭转置矩阵；步骤三、接收机根据公式：Rl=QHRxxQ+(Q*)HRxx*Q*---(2)]]>对协方差矩阵Rxx进行酉变换，得到酉变换后的实数矩阵Rl；公式(2)中，Q为维数为ML×ML的酉变换矩阵，QH为Q的共轭转置矩阵，Q*为Q的共轭矩阵，(Q*)H为Q*的共轭转置矩阵，为Rxx的共轭矩阵，M为所述天线阵的阵元数，L为所述天线阵的抽头数；当ML为偶数时Q=12IjIF-jF---(3)]]>公式(3)中，I为维数为[ML/2]×[ML/2]的单位矩阵，F为维数为[ML/2]×[ML/2]的反对角置换矩阵；当ML为奇数时Q=12I0jI0‾20‾F0-jF---(4)]]>公式(4)中，I为维数为[(ML‑1)/2]×[(ML‑1)/2]的单位矩阵，F为维数为[(ML‑1)/2]×[(ML‑1)/2]的反对角置换矩阵，为维数为1×[(ML‑1)/2]的全零行向量；步骤四、接收机首先求取实数矩阵Rl的逆矩阵然后求取矩阵的m次方，得到矩阵其中，m为正整数；步骤五、接收机根据MUSIC算法构造空间谱函数：Pl=1a~H(θ)Rl-ma~(θ)---(5)]]>公式(5)中，θ为射频信号x(t)的入射角，为的共轭转置矩阵，a(θ)为θ的空时导向矢量且为Kronecker积的符号，as(θ)为θ的空域导向矢量且as(θ)＝[1e‑jφ… e‑j(M‑1)φ]T，φ为所述天线阵的空间相位且φ＝2πdsinθ/λ，d为所述天线阵的阵元间的距离，λ为射频信号x(t)的波长；at(θ)为θ的时域导向矢量且为所述天线阵的时间延时且ft为多普勒频率且vr为所述天线阵的径向速度分量，T′为时间抽头的时延；步骤六、接收机根据公式：θ~l=argminθ1Pl---(6)]]>计算得到干扰波达方向角...

【技术特征摘要】

【专利技术属性】
技术研发人员：葛大江，李文元，梁茂国，周光彬，张亚妮，
申请(专利权)人：中国人民解放军西安通信学院，
类型：发明
国别省市：陕西;61