发明名称 一种快速收敛的卫星导航天线阵抗干扰方法
摘要 本发明公开了一种快速收敛的卫星导航天线阵抗干扰方法,其特征在于:包括下述步骤:(1)对阵列天线接收的射频信号分别进行模拟下变频、A/D采样数字化和数字正交下变频,生成数字基带信号;对于阵元数目为N的天线阵,总共可以得到N路数字基带信号,分别为x<sub>1</sub>(n),x<sub>2</sub>(n),…,x<sub>N</sub>(n);(2)选取第N路基带数字信号x<sub>N</sub>(n)作为参考信号,并对其他N‑1路数字基带信号进行滤波处理;本发明在干扰数目或干扰功率突然减少的情况下,提高算法的收敛速度,减小抗干扰处理对卫星信号的损耗,从而提高接收机捕获跟踪和定位解算性能。
申请公布号 CN106443721A 申请公布日期 2017.02.22
申请号 CN201610992497.9 申请日期 2016.11.11
申请人 深圳市天弓导航科技有限公司 发明人 罗雷波
分类号 G01S19/21(2010.01)I 主分类号 G01S19/21(2010.01)I
代理机构 北京华仲龙腾专利代理事务所(普通合伙) 11548 代理人 黄玉珏
主权项 一种快速收敛的卫星导航天线阵抗干扰方法,其特征在于:包括下述步骤:(1)对阵列天线接收的射频信号分别进行模拟下变频、A/D采样数字化和数字正交下变频,生成数字基带信号;对于阵元数目为N的天线阵,总共可以得到N路数字基带信号,分别为x<sub>1</sub>(n),x<sub>2</sub>(n),…,x<sub>N</sub>(n);(2)选取第N路基带数字信号x<sub>N</sub>(n)作为参考信号,并对其他N‑1路数字基带信号进行滤波处理,滤波过程可表示为:<maths num="0001"><math><![CDATA[<mrow><msub><mi>y</mi><mi>a</mi></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><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><msubsup><mi>w</mi><mi>k</mi><mo>*</mo></msubsup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><msub><mi>x</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001150359490000011.GIF" wi="446" he="134" /></maths>其中,y<sub>a</sub>(n)为n时刻得到的滤波输出信号,w<sub>k</sub>(n)为n时刻阵元k对应的阵列权值,k=1,2,…,N,阵列权值的初始值设定为0,(·)<sup>*</sup>表示取共轭;(3)求n时刻的阵列输出信号,阵列输出信号按下式计算:y(n)=x<sub>N</sub>(n)‑y<sub>a</sub>(n)(4)估计参考信号当前的功率:<maths num="0002"><math><![CDATA[<mrow><mi>p</mi><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mi>L</mi></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>n</mi><mo>=</mo><mi>m</mi><mi>L</mi><mo>+</mo><mn>1</mn></mrow><mrow><mi>m</mi><mi>L</mi><mo>+</mo><mi>L</mi></mrow></munderover><msubsup><mi>x</mi><mi>N</mi><mo>*</mo></msubsup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><msub><mi>x</mi><mi>N</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001150359490000012.GIF" wi="517" he="127" /></maths>式中,p(m)为参考信号当前的功率,L为估计功率所用的数据长度;(5)对干扰功率进行门限检测,看下式是否成立:<maths num="0003"><math><![CDATA[<mrow><mi>z</mi><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>p</mi><mrow><mo>(</mo><mi>m</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow><mrow><mi>p</mi><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow></mrow></mfrac><mo>&gt;</mo><mi>T</mi></mrow>]]></math><img file="FDA0001150359490000013.GIF" wi="414" he="126" /></maths>其中,z(m)为检测量,p(m‑1)为参考信号前一时刻的功率,可按照步骤(4)中的公式进行计算,T为检测门限;(6)更新阵列权值:当检测量小于等于门限时,按下式更新阵列权值:w<sub>k</sub>(n+1)=w<sub>k</sub>(n)‑μx<sub>k</sub>(n)y<sup>*</sup>(n) k=1,2,...,N‑1.其中,μ为步长因子,步长因子的选取是信号处理领域的公知常识;而当检测量大于门限时,按下式对阵列权值进行复位处理:w<sub>k</sub>(n+1)=0 k=1,2,...,N‑1.(7)令n=n+1,重复步骤(2)~(6)进行迭代,得到n+1时刻的阵列输出信号。
地址 518000 广东省深圳市宝安区西乡街道共和工业路107号华丰互联网+创意园A座(华丰科技园)五楼520-521
您可能感兴趣的专利