发明名称 一种基于改进Rake结构的频域抗干扰方法
摘要 本发明属于卫星导航抗干扰领域,具体涉及一种基于改进Rake结构的频域抗干扰方法,主要包括以下步骤:(S1)输入中频采样信号,计算频域抗干扰后信号的相关函数,获取码时延粗估计;(S2)获取通道0的直达信号跟踪结果;(S3)获取除通道0外其余通道的多径信号跟踪结果;(S4)计算去除多径后的直达信号。本发明的频域抗干扰方法在抑制干扰的同时不可避免会导致信号质量下降。尤其在面对宽带干扰时,本发明能够在一定程度上减小抗干扰损耗,提升频域抗干扰后信号的信噪比,减小导航数据解调的误码率,提升接收机的通信性能。
申请公布号 CN106443719A 申请公布日期 2017.02.22
申请号 CN201610813745.9 申请日期 2016.09.09
申请人 中国人民解放军国防科学技术大学 发明人 王飞雪;李建;聂俊伟;耿正霖;鲁祖坤;陈飞强;伍微;彭竞;王勇
分类号 G01S19/21(2010.01)I 主分类号 G01S19/21(2010.01)I
代理机构 北京中济纬天专利代理有限公司 11429 代理人 陈立新
主权项 一种基于改进Rake结构的频域抗干扰方法,其特征在于,包括以下步骤:(S1)输入中频采样信号,计算频域抗干扰后信号的相关函数,获取码时延粗估计,R(τ)表示相关函数,<maths num="0001"><math><![CDATA[<mrow><mi>R</mi><mrow><mo>(</mo><mi>&tau;</mi><mo>)</mo></mrow><mo>=</mo><mo>|</mo><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><mi>x</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mi>c</mi><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mi>&tau;</mi><mo>)</mo></mrow><msup><mi>e</mi><mrow><mo>-</mo><mi>j</mi><mrow><mo>(</mo><mn>2</mn><msub><mi>&pi;f</mi><mi>c</mi></msub><mi>n</mi><mo>/</mo><msub><mi>f</mi><mi>s</mi></msub><mo>+</mo><mi>&theta;</mi><mo>)</mo></mrow></mrow></msup><mo>|</mo></mrow>]]></math><img file="FDA0001112255120000011.GIF" wi="799" he="153" /></maths>其中,x(n)为中频输入采样信号,n为采样点数,c(n‑τ)为接收机本地参考伪码,其中τ是指伪码时延,f<sub>c</sub>为导航信号载频,θ为导航信号初相,f<sub>s</sub>为采样频率,N为积累点数,π为圆周率,j为虚数单位;通过搜索R(τ)的极大值获取码时延粗估计<img file="FDA0001112255120000012.GIF" wi="166" he="67" />m表示通道号,取值为m=0,1,2,...,M,m取0时表示对应通道是直达信号,m取1至M时表示对应通道是多径信号;(S2)获取通道0的直达信号跟踪结果将各个跟踪通道跟踪结果初始值<img file="FDA0001112255120000013.GIF" wi="125" he="71" />全部置零,令通道号m=0,入口信号为<img file="FDA0001112255120000014.GIF" wi="457" he="63" />初始时码跟踪延为<img file="FDA0001112255120000015.GIF" wi="138" he="58" />其中<img file="FDA0001112255120000016.GIF" wi="112" he="61" />指除通道0外所有通道跟踪结果之和,即<maths num="0002"><math><![CDATA[<mrow><msub><mover><mi>s</mi><mo>&OverBar;</mo></mover><mn>0</mn></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mover><mi>s</mi><mo>^</mo></mover><mi>k</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001112255120000017.GIF" wi="333" he="120" /></maths><img file="FDA0001112255120000018.GIF" wi="124" he="72" />为第k个跟踪通道的跟踪结果;启动通道的载波跟踪和码跟踪,码跟踪环和载波跟踪环锁定后,码跟踪环提供的码时延精确估计为<img file="FDA0001112255120000019.GIF" wi="39" he="51" />和载波跟踪环提供的载波相位精确估计为<img file="FDA00011122551200000110.GIF" wi="65" he="63" />利用下式计算信号幅度估计:<maths num="0003"><math><![CDATA[<mrow><msub><mover><mi>a</mi><mo>^</mo></mover><mn>0</mn></msub><mo>=</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>n</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><msub><mover><mi>s</mi><mo>^</mo></mover><mrow><mi>i</mi><mi>n</mi><mo>,</mo><mn>0</mn></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mi>c</mi><mrow><mo>(</mo><mi>n</mi><mo>-</mo><msub><mover><mi>&tau;</mi><mo>^</mo></mover><mn>0</mn></msub><mo>)</mo></mrow><msup><mi>e</mi><mrow><mo>-</mo><mi>j</mi><mrow><mo>(</mo><mn>2</mn><msub><mi>&pi;f</mi><mi>c</mi></msub><mi>n</mi><mo>/</mo><msub><mi>f</mi><mi>s</mi></msub><mo>+</mo><msub><mover><mi>&theta;</mi><mo>^</mo></mover><mn>0</mn></msub><mo>)</mo></mrow></mrow></msup></mrow>]]></math><img file="FDA00011122551200000111.GIF" wi="760" he="121" /></maths>获得通道0的直达信号跟踪结果:<maths num="0004"><math><![CDATA[<mrow><msub><mover><mi>s</mi><mo>^</mo></mover><mn>0</mn></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><msub><mover><mi>a</mi><mo>^</mo></mover><mn>0</mn></msub><mi>c</mi><mrow><mo>(</mo><mi>n</mi><mo>-</mo><msub><mover><mi>&tau;</mi><mo>^</mo></mover><mn>0</mn></msub><mo>)</mo></mrow><msup><mi>e</mi><mrow><mo>-</mo><mi>j</mi><mrow><mo>(</mo><mn>2</mn><msub><mi>&pi;f</mi><mi>c</mi></msub><mi>n</mi><mo>/</mo><msub><mi>f</mi><mi>s</mi></msub><mo>+</mo><msub><mover><mi>&theta;</mi><mo>^</mo></mover><mn>0</mn></msub><mo>)</mo></mrow></mrow></msup><mo>;</mo></mrow>]]></math><img file="FDA0001112255120000021.GIF" wi="815" he="87" /></maths>(S3)获取除通道0外其余通道的多径信号跟踪结果通道号分别取m=1,2,...,M,入口信号为<img file="FDA0001112255120000022.GIF" wi="534" he="76" />初始时码跟踪延为<img file="FDA0001112255120000023.GIF" wi="168" he="67" />其中<img file="FDA0001112255120000024.GIF" wi="112" he="62" />指除通道m外所有通道跟踪结果之和,<img file="FDA0001112255120000025.GIF" wi="435" he="145" />启动通道的载波跟踪和码跟踪,码跟踪环和载波跟踪环锁定后,码环提供的码时延精确估计为<img file="FDA0001112255120000026.GIF" wi="55" he="58" />和载波跟踪环提供的载波相位精确估计为<img file="FDA0001112255120000027.GIF" wi="79" he="71" />利用下式计算信号幅度估计:<maths num="0005"><math><![CDATA[<mrow><msub><mover><mi>a</mi><mo>^</mo></mover><mi>m</mi></msub><mo>=</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>n</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><msub><mover><mi>s</mi><mo>^</mo></mover><mrow><mi>i</mi><mi>n</mi><mo>,</mo><mi>m</mi></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mi>c</mi><mrow><mo>(</mo><mi>n</mi><mo>-</mo><msub><mover><mi>&tau;</mi><mo>^</mo></mover><mi>m</mi></msub><mo>)</mo></mrow><msup><mi>e</mi><mrow><mo>-</mo><mi>j</mi><mrow><mo>(</mo><mn>2</mn><msub><mi>&pi;f</mi><mi>c</mi></msub><mi>n</mi><mo>/</mo><msub><mi>f</mi><mi>s</mi></msub><mo>+</mo><msub><mover><mi>&theta;</mi><mo>^</mo></mover><mi>m</mi></msub><mo>)</mo></mrow></mrow></msup></mrow>]]></math><img file="FDA0001112255120000028.GIF" wi="928" he="146" /></maths>获得通道m的多径信号跟踪结果:<maths num="0006"><math><![CDATA[<mrow><msub><mover><mi>s</mi><mo>^</mo></mover><mi>m</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><msub><mover><mi>a</mi><mo>^</mo></mover><mi>m</mi></msub><mi>c</mi><mrow><mo>(</mo><mi>n</mi><mo>-</mo><msub><mover><mi>&tau;</mi><mo>^</mo></mover><mi>m</mi></msub><mo>)</mo></mrow><msup><mi>e</mi><mrow><mo>-</mo><mi>j</mi><mrow><mo>(</mo><mn>2</mn><msub><mi>&pi;f</mi><mi>c</mi></msub><mi>n</mi><mo>/</mo><msub><mi>f</mi><mi>s</mi></msub><mo>+</mo><msub><mover><mi>&theta;</mi><mo>^</mo></mover><mi>m</mi></msub><mo>)</mo></mrow></mrow></msup><mo>,</mo><mi>m</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mo>,</mo><mi>M</mi><mo>;</mo></mrow>]]></math><img file="FDA0001112255120000029.GIF" wi="1119" he="95" /></maths>(S4)计算去除多径后的直达信号,根据步骤(S3)得到的M个多径信号跟踪结果,更新<img file="FDA00011122551200000210.GIF" wi="129" he="72" /><maths num="0007"><math><![CDATA[<mrow><msub><mover><mi>s</mi><mo>&OverBar;</mo></mover><mn>0</mn></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mover><mi>s</mi><mo>^</mo></mover><mi>k</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA00011122551200000211.GIF" wi="358" he="138" /></maths>去除多径后的直达信号估计<img file="FDA00011122551200000212.GIF" wi="101" he="53" />为:<maths num="0008"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><msubsup><mover><mi>s</mi><mo>^</mo></mover><mn>0</mn><mo>&prime;</mo></msubsup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mi>x</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>-</mo><msub><mover><mi>s</mi><mo>&OverBar;</mo></mover><mn>0</mn></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><mi>x</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mover><mi>s</mi><mo>^</mo></mover><mi>k</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow></mtd></mtr></mtable><mo>.</mo></mrow>]]></math><img file="FDA00011122551200000213.GIF" wi="623" he="239" /></maths>
地址 410073 湖南省长沙市开福区德雅路109号