发明名称 北斗卫星用户设备接收抗干扰方法
摘要 本发明提供一种北斗卫星用户设备接收抗干扰方法,以解决北斗导航系统实际使用中的抗干扰问题,涉及卫星导航及智能天线技术领域,结合时空自适应滤波和卡尔曼滤波的抗干扰算法,通过计算通道的幅相误差校正权值,对数字中频信号数据进行校正,协方差矩阵估计,计算自适应权值和更新自适应滤波器权系数,可以在不知道干扰和信号方向的情况下在北斗卫星用户设备接收天线的方向图上对干扰进行精确的零陷,抑制各类干扰,并保持信号的完整性;本发明计算量小、鲁棒性好、抗干扰能力强、对信号不损伤。
申请公布号 CN101887129A 申请公布日期 2010.11.17
申请号 CN201010154321.9 申请日期 2010.04.23
申请人 西安航空电子科技有限公司 发明人 陈雪峰
分类号 G01S19/21(2010.01)I;H04B7/08(2006.01)I 主分类号 G01S19/21(2010.01)I
代理机构 西安西交通盛知识产权代理有限责任公司 61217 代理人 罗永娟
主权项 1.北斗卫星用户设备接收抗干扰方法,其特征在于包括以下基本步骤:第一步:计算通道的幅相误差校正权值;在接收天线阵列的法线方向通过发射天线辐射单点频的校准信号s(w),第k通道的传递函数表示为s<sub>k</sub>(w)=α<sub>k</sub>(w)exp[jβ<sub>k</sub>(w)]*s(w),以任一通道的数据为参考信号s<sub>1</sub>(w),那么k通道和参考通道之间的差异为校正权值,用复数除法得到校正权值η<sub>k</sub>,<maths num="0001"><![CDATA[<math><mrow><msub><mi>&eta;</mi><mi>k</mi></msub><mo>=</mo><mfrac><mrow><msub><mi>s</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>w</mi><mo>)</mo></mrow></mrow><mrow><msub><mi>s</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>w</mi><mo>)</mo></mrow></mrow></mfrac><mo>=</mo><mfrac><mrow><msub><mi>&alpha;</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>w</mi><mo>)</mo></mrow></mrow><mrow><msub><mi>&alpha;</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>w</mi><mo>)</mo></mrow></mrow></mfrac><mi>exp</mi><mo>[</mo><mi>j</mi><mrow><mo>(</mo><msub><mi>&beta;</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>w</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>&beta;</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>w</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>]</mo></mrow></math>]]></maths>第二步:将阵列天线接收到的北斗卫星信号进行低噪声放大和低通滤波,然后下变频为模拟中频信号,再经数模转换为数字中频信号,供后续数字域的信号处理;第三步:对数字中频信号数据进行校正;对第二步中得到每个通道数字中频信号数据用校正权值进行补偿,以消除不同通道间的幅相误差影响,第k通道校正后的输出信号y′<sub>k</sub>=y<sub>k</sub>·η<sub>k</sub>,第四步:协方差矩阵估计;抗干扰算法中采用n个阵元、p个抽头,则在时刻k接收的空时采样快拍向量为x(k)=[(x<sup>0</sup>(n))<sup>T</sup>,(x<sup>1</sup>(n))<sup>T</sup>,...,(x<sup>p</sup>(n))<sup>T</sup>]<sup>T</sup>,考虑收集M个空时快拍,则接收的空时(n×p)×M维采样矩阵为X=[x(1),x(2),...,x(M)],则空时协方差矩阵按如下计算:R=XX<sup>H</sup>/M第五步:计算自适应权值w<sub>k</sub>;获得协方差矩阵逆后,通过如下公式计算空时自适应权值,<maths num="0002"><![CDATA[<math><mrow><msub><mi>w</mi><mi>k</mi></msub><mo>=</mo><msubsup><mi>R</mi><mi>k</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><msub><mi>a</mi><mi>st</mi></msub><mo>/</mo><mrow><mo>(</mo><msubsup><mi>a</mi><mi>st</mi><mi>H</mi></msubsup><msubsup><mi>R</mi><mi>k</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><msub><mi>a</mi><mi>st</mi></msub><mo>)</mo></mrow></mrow></math>]]></maths>其中<img file="FDA0000020850670000021.GIF" wi="260" he="63" />a<sub>s</sub>=[1 1 1 1 ]<sup>T</sup>,a<sub>t</sub>=[1 0 0 0]<sup>T</sup>,R<sub>k</sub>为k时刻的协方差矩阵,H表示共轭操作。
地址 710075 陕西省西安市高新六路28号