发明名称 一种基于子空间投影的导航接收机欺骗式干扰消除方法
摘要 一种基于子空间投影的导航接收机欺骗式干扰消除方法,涉及导航接收机抗干扰技术。本发明的目的是利用单天线即可有效的对欺骗式干扰进行抑制,同时不降低可用卫星的数量。确定未被欺骗的M<sub>1</sub>个授权卫星信号以及M<sub>3</sub>个欺骗式干扰信号;对已经完成接收的欺骗式干扰信号的载波频率信息、扩频码信息和导航电文信息进行复现;利用观测到的欺骗式干扰信号的参数信息对欺骗式干扰信号矩阵进行重建;利用得到的欺骗式干扰信号矩阵构建欺骗式干扰信号子空间的投影矩阵;计算欺骗式干扰子空间的补空间;将接收的中频信号在补空间上进行投影,即完成干扰抑制运算。本发明能确保导航接收机在授权卫星信号与欺骗式干扰信号并存时能够正常对授权信号进行接收。
申请公布号 CN104459723A 申请公布日期 2015.03.25
申请号 CN201410834899.7 申请日期 2014.12.29
申请人 哈尔滨工业大学 发明人 韩帅;陈雷;孟维晓
分类号 G01S19/21(2010.01)I 主分类号 G01S19/21(2010.01)I
代理机构 哈尔滨市松花江专利商标事务所 23109 代理人 杨立超
主权项 一种基于子空间投影的导航接收机欺骗式干扰消除方法,定义导航接收机接收到的信号经过下变频后得到的中频信号为r<sub>IF</sub>,包含M<sub>1</sub>个未被欺骗的授权卫星信号、M<sub>2</sub>个被欺骗的授权卫星信号和M<sub>3</sub>个欺骗式干扰信号,其中有M<sub>2</sub>≤M<sub>3</sub>,n<sub>0</sub>为高斯白噪声;<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>r</mi><mi>IF</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>M</mi><mn>1</mn></msub></munderover><msubsup><mi>r</mi><mrow><mi>real</mi><mn>1</mn></mrow><mi>j</mi></msubsup><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>M</mi><mn>2</mn></msub></munderover><msubsup><mi>r</mi><mrow><mi>real</mi><mn>2</mn></mrow><mi>j</mi></msubsup><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>M</mi><mn>3</mn></msub></munderover><msubsup><mi>r</mi><mi>spoof</mi><mi>i</mi></msubsup><mo>+</mo><msub><mi>n</mi><mn>0</mn></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000646297160000011.GIF" wi="838" he="158" /></maths>M<sub>3</sub>个欺骗式干扰信号对M<sub>1</sub>个未被欺骗的授权卫星信号的解算结果产生影响,将这M<sub>3</sub>个欺骗式干扰信号消除,其特征在于,欺骗式干扰信号的消除过程为:步骤一、导航接收机开机对信号进行接收,同时对欺骗式干扰信号进行检测,确定未被欺骗的M<sub>1</sub>个授权卫星信号;确定M<sub>3</sub>个欺骗式干扰信号;步骤二、对已经完成接收的欺骗式干扰信号的载波频率信息<img file="FDA0000646297160000012.GIF" wi="364" he="85" />扩频码信息<img file="FDA0000646297160000013.GIF" wi="340" he="85" />和导航电文信息<img file="FDA0000646297160000014.GIF" wi="332" he="86" />进行复现;而这三种参数可在欺骗式干扰信号的观测通道的载波跟踪环和码跟踪环得到;步骤三、利用观测到的欺骗式干扰信号的参数信息对欺骗式干扰信号矩阵进行重建,假设每次运算的数据长度为K,则欺骗式干扰信号矩阵按下面的公式进行重建:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>S</mi><mo>=</mo><msub><mrow><mo>[</mo><msub><mi>S</mi><mn>1</mn></msub><mo>,</mo><msub><mi>S</mi><mn>2</mn></msub><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msub><mi>S</mi><msub><mi>M</mi><mn>3</mn></msub></msub><mo>]</mo></mrow><mrow><mi>K</mi><mo>&times;</mo><msub><mi>M</mi><mn>3</mn></msub></mrow></msub></mrow>]]></math><img file="FDA0000646297160000015.GIF" wi="476" he="77" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>S</mi><mi>j</mi></msub><mo>=</mo><msub><mfenced open='[' close=']'><mtable><mtr><mtd><msubsup><mi>d</mi><mi>spoof</mi><mi>j</mi></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mn>1</mn></msub><mo>)</mo></mrow><msubsup><mi>C</mi><mi>spoof</mi><mi>j</mi></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mn>1</mn></msub><mo>)</mo></mrow><msup><mi>e</mi><mrow><mi>j</mi><mrow><mo>(</mo><mn>2</mn><mi>&pi;</mi><msubsup><mi>f</mi><mi>spoof</mi><mi>j</mi></msubsup><msub><mi>t</mi><mn>1</mn></msub><mo>)</mo></mrow></mrow></msup></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><msubsup><mi>d</mi><mi>spoof</mi><mi>j</mi></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mi>K</mi></msub><mo>)</mo></mrow><msubsup><mi>C</mi><mi>spoof</mi><mi>j</mi></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mi>K</mi></msub><mo>)</mo></mrow><msup><mi>e</mi><mrow><mi>j</mi><mrow><mo>(</mo><mn>2</mn><mi>&pi;</mi><msubsup><mi>f</mi><mi>spoof</mi><mi>j</mi></msubsup><msub><mi>t</mi><mi>K</mi></msub><mo>)</mo></mrow></mrow></msup></mtd></mtr></mtable></mfenced><mrow><mi>K</mi><mo>&times;</mo><mn>1</mn></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000646297160000016.GIF" wi="1249" he="287" /></maths>式中参数含义如步骤二所示,t<sub>i</sub>代表第i个采样时刻;步骤四、利用得到的欺骗式干扰信号矩阵按照公式(3)构建欺骗式干扰信号子空间的投影矩阵:H<sub>spoof</sub>=S(S<sup>H</sup>S)<sup>‑1</sup>S<sup>H</sup>    (3)式中,S<sup>H</sup>是S的共轭转置矩阵;步骤五、计算欺骗式干扰子空间的补空间:H<sub>c</sub>=I‑H<sub>spoof</sub>    (4)步骤六、将接收的中频信号在公式(4)得到的补空间上进行投影,即完成干扰抑制运算:r=(I‑H<sub>spoof</sub>)·r<sub>IF</sub>    (5)式中,I是单位矩阵;步骤七、M<sub>3</sub>个欺骗式干扰信号被消除,剩下的M<sub>2</sub>个被欺骗的授权卫星信号变为正常的授权卫星信号;得到的信号r近似为包含M<sub>1</sub>+M<sub>2</sub>个授权卫星信号和噪声的混合信号,即可对其进行捕获、跟踪、解调和解算。
地址 150001 黑龙江省哈尔滨市南岗区西大直街92号