一种卫星导航抗干扰天线数字信号正交分解方法

一种卫星导航抗干扰天线数字信号正交分解方法
【技术领域】
[0001] 本发明属于卫星导航抗干扰领域，是实现卫星导航抗干扰天线中频数字信号正交 分解的一种算法。
【背景技术】
[0002] 在卫星导航抗干扰天线的信号处理过程中，对AD采样后的数字实信号进行正交 分解是数字信号处理的第一步。数字信号正交分解得到的I、Q两路数字信号的幅度和相位 的平衡性是影响卫星导航抗干扰天线对干扰抑制能力的关键因素之一。因此，为了更好地 完成对卫星导航干扰信号的抑制，前提是设计一种可较好保持I、Q两路数字信号的幅度和 相位平衡性的数字信号正交分解方法。
[0003] 目前，对实信号进行数字正交分解的主要方法有：基于傅里叶变换的数字正交分 解方法、基于希尔伯特变换的数字正交分解方法、W特定的采样方式进行的数字正交分解 方法W及利用数字混频进行的数字正交分解方法。基于傅里叶变换的数字正交分解方法利 用两次傅里叶变换完成I、Q路的分解，可W得到较为理想的镜像抑制效果，但计算量较大。 基于希尔伯特变换的数字正交分解方法的计算量较小，但由于对输入信号采样点数为有限 字长，导致其分解后得到的I、Q路的不平衡性较大。W特定的采样方式进行的数字正交分 解方法工程实现较为方便，对窄带信号的分解效果较好，但对于宽带信号，其分解后得到的 I、Q路的不平衡性较大。利用数字混频进行的数字正交分解方法的采样率不受信号频率的 限制，对于窄带信号，其精度比W特定的采样方式进行的数字正交分解方法高，但对于宽带 信号其分解后得到的I、Q路的不平衡性依然较大。
[0004] 针对卫星导航信号的压制式干扰信号通常为覆盖整个卫星导航信号频带范围的 宽带信号。在卫星导航抗干扰天线处理中通常采用基于傅里叶变换的数字正交分解方法完 成数字信号的I、Q路分解。但是，由于卫星导航抗干扰天线包含多个通道，同时对多个通道 进行傅里叶变换及逆傅里叶变换的运算量非常大，需要耗费数字信号处理器中大量的乘法 器和存储器资源。

【发明内容】

[0005] 为了克服现有技术的不足，本发明提供一种卫星导航抗干扰天线数字信号正交分 解方法，采用偶对称FIR滤波器和奇对称FIR滤波器完成对模数转换后中频实数字信号的 正交分解，适用于对宽带信号的正交分解处理、计算量小且易于工程实现。
[0006] 本发明解决其技术问题所采用的技术方案包括W下步骤：
[0007] 第一步、计算FIR滤波器的归一化频率点，包括通带第一频率点
通带第二频率点
和阻带第二频率点
其中，f。是卫星导航抗干扰天线中频工作的中屯、频率，B是卫星导航 抗干扰天线中频工作的带宽，F,是模数转换器的采样频率，D是通带与阻带之间的过渡带 宽；
[0008] 第二步、设计FIR滤波器在归一化频率范围内的期望幅度响应向量
[0009] a= [a0,astopi,Spassi,apass2, 3別。。2,aj
[0010] 其中，a。表示在归一化频率点0上的期望幅度响应值；astwi表示在归一化频率点 fst。pl上的期望幅度响应值；apass康示在归一化频率点fpassl上的期望幅度响应值；apass2表 示在归一化频率点上的期望幅度响应值；表示在归一化频率点fstDp2上的期望幅 度响应值；曰1表示在归一化频率点1上的期望幅度响应值；
[0011] a。、3与1、apassi、3pass2、3加。2和a1的取值范围在0~1之间。
[0012] 第S步、根据归一化频率点和期望幅度响应向量，采用Parks-McClellan算法设 计K阶的滤波系数偶对称FIR滤波器和滤波系数奇对称滤波器，分别得到偶对称FIR滤波 器系数向量bw。。和奇对称FIR滤波器系数向量b"dd;将卫星导航抗干扰天线中频实数字信 号同时送入偶对称FIR滤波器和奇对称FIR滤波器，偶对称FIR滤波器输出与输入数字信 号同相的I路信号，奇对称FIR滤波器输出与输入数字信号正交的Q路信号。
[0013] 本发明的有益效果是：采用本发明方法实现的128阶偶对称FIR和奇对称FIR滤 波器对模数转换后中频实数字信号进行正交分解后对镜频的抑制效果与采用傅里叶变换 法进行8192点FFT运算得到的镜频抑制效果相当。但是，本发明方法需要的复乘和复加运 算共为254次，而8192点FFT运算则需要复乘和复加运算共为53248次。此外，傅里叶变 换数字信号正交分解方法采用块处理方式实现对每个数据块的正交分解，但数据块之间的 信号相位连续性无法保证，而本发明方法只要将偶对称FIR滤波器和奇对称FIR滤波器系 数确定后即可通过流水线形式实现对输入数字中频实信号的正交分解，可W保证输出正交 数字信号的相位连续性。因此，本发明方法比基于傅里叶变换的数字正交分解方法计算量 小，更适用于对卫星导航抗干扰天线信号处理中AD采样后的数字实信号进行正交分解。
【附图说明】
[0014] 图1是期望幅度响应向量与归一化频率点对应关系示意图；
[0015] 图2是卫星导航抗干扰天线数字信号正交分解过程示意图；
[0016] 图3是FIR滤波器期望幅度响应示意图；
[0017] 图4是偶对称FIR滤波器系数和奇对称FIR滤波器系数示意图，其中，（a)是偶对 称FIR滤波器系数，化）是奇对称FIR滤波器系数；
[0018] 图5是偶对称FIR滤波器和奇对称FIR滤波器幅度响应示意图，其中，（a)是偶对 称FIR滤波器幅度响应，化）是奇对称FIR滤波器幅度响应；
[0019] 图6是偶对称FIR滤波器和奇对称FIR滤波器相位响应。
【具体实施方式】
[0020] 下面结合附图和实施例对本发明进一步说明，本发明包括但不仅限于下述实施 例。
[0021] 本发明偶对称FIR滤波器和奇对称FIR滤波器的设计实现步骤如下：
[002引第一步；根据卫星导航抗干扰天线中频工作的中屯、频率f。、带宽B、模数转换器的 采样频率及通带与阻带之间的过渡带宽D(通常取2MHz~4MHz)计算FIR滤波器归 一化的通带第一频率点fpassl、通带第二频率点fpass2、阻带第一频率点fstDpl和阻带第二频率 点计算表达式如下：
[0027] 第二步：设计FIR滤波器在归一化频率范围内的期望幅度响应向量a,表达式如 下：
[0028] a= [a0, 3別。。1,Spassi,apass2, 3別。。2,aj 巧)
[0029] 其中，a。表示在归一化频率点0上的期望幅度响应值；a,tDpi表示在归一化频率点 fst。pl上的期望幅度响应值；apass康示在归一化频率点fpassl上的期望幅 度响应值；apass2表 示在归一化频率点上的期望幅度响应值；表示在归一化频率点fstDp2上的期望幅 度响应值；ai表示在归一化频率点1上的期望幅度响应值。a"、a日tDpi、a。。日日i、a。。日日2、屯tw2和a1 的取值范围在0~1之间。图1所示为期望幅度响应向量与归一化频率点的对应关系。
[0030] 第S步：根据第一步计算得到的归一化频率点和第二步设计得到的期望幅度响应 向量，采用Parks-McClellan算法设计阶数为K阶的滤波系数偶对称FIR滤波器和滤波系 数奇对称滤波器，分别得到偶对称FIR滤波器系数向量bw。。和奇对称FIR滤波器系数向量 b〇dd°
[0031] 将卫星导航抗干扰天线中频实数字信号同时送入滤波系数为bw。。和b"dd的FIR滤 波器，滤波系数为bw。。的FIR滤波器输出与输入数字信号同相的I路信号，滤波系数为b"dd 的FIR滤波器输出与输入数字信号正交的Q路信号。
[0032] 卫星导航抗干扰天线数字信号正交分解过程如图2所示。设接收天线阵列由N个 相同的天线单元构成。下变频处理将N单元天线阵列接收到的射频信号转变为中频模拟信 号。模数转换处理将中频模拟信号转换为中频实数数字信号后分别送入阶数等长的偶对称 FIR滤波器和奇对称FIR滤波器。偶对称FIR滤波器输出与输入实数数字信号同相的数字 信号，即I路信号；奇对称FIR滤波器输出与输入实数数字信号正交的数字信号，即Q路信 号。
[003引本发明适用于GPS、抓S和化0NASS卫星导航抗干扰天线中频数字信号的正交分解 处理。W抓SB3频点抗干扰天线中频数字信号的正交分解为例阐述本发明的具体实施方 式：
[0034] 步骤1 ;将抓2B3频点抗干扰天线中频中屯、频率片二15. 52MHz、采样率F,= 62MHz、工作带宽B= 20MHz和所设计FIR滤波器过渡带宽D=IMHz参数代入式（1)、式 似、式做和式…分别计算巧R滤波器归一化的通带第一频率点通带第二频率点 fpass2、阻带第一频率点fstopl和阻带第二频率点fSt。的：
[0039] 步骤2;设计FIR滤波器在归一化频率范围内的期望幅度响应向量a:
[0040]a= [a0,astopi, 过passl,SpassS, 3別。02,aj- [O, 0, 1, 1, 0,O] (10)
[0041] 得到FIR滤波器期望幅度响应如图3所示。
[0042] 步骤3;根据步骤1计算得到的归一化频率点和步骤2设计得到的期望幅度响应 向量a,采用Parks-McClellan算法设计阶数为16阶化=15)的偶对称FIR滤波器和奇对 称滤波器。具体实现采用Matl油软件Signalprocessing工具箱的fi巧m函数得到滤波系 数向量beven和bDdd;
[0043] beven=firpm(Xf,a) (11)
[0044] b〇dd=firpm(Xf,a,，hilbert，） (12)
[0045]其中，f= [0, 0. 1458, 0. 1781，0. 8232, 0. 8555, 1]为归一化频率向量，a= [0, 0, 1，1，0, 0]为与归一化频率向量对应的期望幅度响应向量。得到的偶对称FIR滤波器 系数和奇对称FIR滤波器系数分别如图4(a)和图4(b)所示。所设计的偶对称FIR滤波器 和奇对称FIR滤波器的幅度响应分别如图5 (a)和图5化）所示。图6为所设的偶对称FIR 滤波器和奇对称FIR滤波器的相位响应，其中实线为偶对称FIR滤波器的相位响应，虚线为 奇对称FIR滤波器的相位响应。
[0046] 将抓2B3频点卫星导航抗干扰天线中频采样信号同时送入由步骤3得到的偶对称 FIR滤波器和奇对称FIR滤波器进行滤波处理即可得到与输入信号同相的I路输出信号和 与输入信号正交的Q路输出信号。至此，完成对抓2B3频点卫星导航抗干扰天线中频数字 信号的正交分解处理。
【主权项】
1. 一种卫星导航抗干扰天线数字信号正交分解方法，其特征在于包括下述步骤： 第一步、计算FIR滤波器的归一化频率点，包括通带第一频率^通 带第二频率、阻带第一频率点和阻带第二频率点 .其中，&是卫星导航抗干扰天线中频工作的中心频率，B是卫星导航 抗干扰天线中频工作的带宽，Fs是模数转换器的采样频率，D是通带与阻带之间的过渡带 宽； 第二步、设计FIR滤波器在归一化频率范围内的期望幅度响应向量 ^ 一 [<3-〇> <3-stopl> ^passH ^pass2> ^stop2> 其中，Btl表示在归一化频率点O上的期望幅度响应值；a stopl表示在归一化频率点fs_ 上的期望幅度响应值；apassl表示在归一化频率点fpassll的期望幅度响应值；a pass2表示在归 一化频率点&_2上的期望幅度响应值；astop2表示在归一化频率点f st()p2上的期望幅度响应 值；S1表示在归一化频率点1上的期望幅度响应值； 第三步、根据归一化频率点和期望幅度响应向量，采用Parks-McClellan算法设计K阶 的滤波系数偶对称FIR滤波器和滤波系数奇对称滤波器，分别得到偶对称FIR滤波器系数 向量13_和奇对称FIR滤波器系数向量b 将卫星导航抗干扰天线中频实数字信号同时 送入偶对称FIR滤波器和奇对称FIR滤波器，偶对称FIR滤波器输出与输入数字信号同相 的I路信号，奇对称FIR滤波器输出与输入数字信号正交的Q路信号。2. 根据权利要求1所述的卫星导航抗干扰天线数字信号正交分解方法，其特征在于： aQ、a_、apassl、a pass2、a-和a i的取值范围在O~1之间。
【专利摘要】本发明提供了一种卫星导航抗干扰天线数字信号正交分解方法，根据归一化频率点和期望幅度响应向量，采用Parks-McClellan算法设计K阶的滤波系数偶对称FIR滤波器和滤波系数奇对称滤波器，将卫星导航抗干扰天线中频实数字信号同时送入偶对称FIR滤波器和奇对称FIR滤波器，偶对称FIR滤波器输出与输入数字信号同相的I路信号，奇对称FIR滤波器输出与输入数字信号正交的Q路信号。本发明方法比基于傅里叶变换的数字正交分解方法计算量小，更适用于对卫星导航抗干扰天线信号处理中AD采样后的数字实信号进行正交分解。
【IPC分类】G01S19/37
【公开号】CN104880721
【申请号】CN201510242516
【发明人】王晓宇, 谢斌斌, 金燕, 张骅
【申请人】中国电子科技集团公司第二十研究所
【公开日】2015年9月2日
【申请日】2015年5月13日