CN110501727A - 一种基于空频自适应滤波的卫星导航抗干扰方法 - Google Patents
一种基于空频自适应滤波的卫星导航抗干扰方法 Download PDFInfo
- Publication number
- CN110501727A CN110501727A CN201910748274.1A CN201910748274A CN110501727A CN 110501727 A CN110501727 A CN 110501727A CN 201910748274 A CN201910748274 A CN 201910748274A CN 110501727 A CN110501727 A CN 110501727A
- Authority
- CN
- China
- Prior art keywords
- data
- frequency point
- matrix
- point data
- frequency
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/21—Interference related issues ; Issues related to cross-correlation, spoofing or other methods of denial of service
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/35—Constructional details or hardware or software details of the signal processing chain
- G01S19/37—Hardware or software details of the signal processing chain
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Signal Processing (AREA)
- Noise Elimination (AREA)
Abstract
本发明属于卫星导航抗干扰领域,公开了一种空频自适应处理方法。该方法包括接收M个阵元的中频数据;在每个阵元的中频数据中,每N个点形成一个数据段,对每个数据段进行加窗处理得到加窗数据段;对加窗数据段做25%的重叠处理得到重叠数据段;对每个重叠数据段进行N点傅立叶变换得到对应的频点数据;利用频点数据,计算与频点数据中的每个数据段的协方差矩阵;对协方差矩阵进行自适应对角加载得到加载矩阵;对加载矩阵进行求逆得出逆矩阵;利用逆矩阵求得频点数据的最优权值;利用最优权值对频点数据进行加权处理得到加权数据;对加权数据进行逆傅立叶变换并输出处理结果。
Description
技术领域
本发明属于卫星导航抗干扰技术领域,具体涉及一种空频自适应的信号处理方法。
背景技术
卫星信号到达地面时只有-130dBm,极易受到干扰。在复杂的电磁环境下,卫星导航的脆弱性成为制约其持续稳定定位的主要因素。卫星导航抗干扰技术,利用方向图可控天线阵(CRPA),实时自适应地生成空域滤波权矢量,并做加权处理,使抗干扰天线在空间来向未知的干扰方向自适应地形成零陷,或在卫星信号入射方向形成增益,获得良好的信噪比。这对于战时复杂电磁环境下惯性/卫星组合导航的生存能力,有至关重要的意义。因为干扰个数、带宽、时变性各不同,对抗干扰阵列天线的自由度、各通道的一致性及抗干扰算法的运算速度都有严苛要求,因此抗干扰算法及其高效实现成为卫星导航专业中非常复杂和必要的研究方向。
发明内容
发明目的:
针对卫星导航易受干扰的问题,提出一种空域自适应滤波的数字信号处理方法,实现卫星导航抗干扰的功能。
技术方案:
第一方面,提出了一种基于空频自适应滤波的卫星导航抗干扰方法,其特征在于,包括:
接收M个阵元的中频数据;
在每个阵元的中频数据中,每N个点形成一个数据段,对每个数据段进行汉明窗加窗处理得到加窗数据段;
将每个加窗数据段的后25%×N点数据与后一加窗数据段的前25%×N点数据相加并用加和之后的数据替换所述当前加窗数据段的后25%×N点数据得到重叠数据段;
对每个重叠数据段进行N点快速傅立叶变换得到对应的频点数据;
将所述频点数据缓存在双口RAM中;
利用所述频点数据,计算与所述频点数据中的每个数据段中的前13个频点以及后13个频点对应的协方差矩阵;
对所述协方差矩阵进行自适应对角加载得到加载矩阵;
对所述加载矩阵进行求逆得出逆矩阵;
利用逆矩阵,根据线性约束最小方差准则求得所述频点数据的最优权值;
从所述双口RAM中读取所述频点数据;
利用所述最优权值对所述频点数据进行加权处理得到加权数据;
对加权数据进行逆快速傅立叶变换并输出处理结果,其中,所述M和N均为自然数。
可选地,对所述加载矩阵进行求逆是利用FPGA对所述加载矩阵进行求逆。
可选地,所述协方差矩阵是二阶的。
可选地,对所述协方差矩阵进行自适应对角加载是指对每个协方差矩阵,自适应地加上第一行第一列数据除以8192的数值。
可选地,其中根据线性约束最小方差准则求得所述频点数据的最优权值包括对所述频点数据中的每个数据段中的第14至第N-13个频点的权值取零。
可选地,所利用的最优权值是与其进行加权处理的频点数据计算得到。
可选地,所述N为64。
第二方面,提出了一种基于空频自适应滤波的卫星导航抗干扰装置,其特征在于,包括:
接收模块,用于接收M个阵元的中频数据;
加窗模块,用于在每个阵元的中频数据中,每N个点形成一个数据段,对每个数据段进行汉明窗加窗处理得到加窗数据段;
叠加模块,用于将每个加窗数据段的后16点数据与后一加窗数据段的前16点数据相加并用加和之后的数据替换所述当前加窗数据段的后16点数据得到重叠数据段;
变换模块,用于对每个重叠数据段进行N点快速傅立叶变换得到对应的频点数据;
存储模块,用于将所述频点数据缓存在双口RAM中;
协方差矩阵计算模块,用于利用所述频点数据,计算与所述频点数据中的每个数据段中的前13个频点以及后13个频点对应的协方差矩阵;
对角加载模块,用于对所述协方差矩阵进行自适应对角加载得到加载矩阵;
求逆模块,用于对所述加载矩阵进行求逆得出逆矩阵;
最优权值计算模块,用于利用逆矩阵,根据线性约束最小方差准则求得所述频点数据的最优权值;
读取模块,用于从所述双口RAM中读取所述频点数据;
加权处理模块,用于利用所述最优权值对所述频点数据进行加权处理得到加权数据;
逆变换模块,用于对加权数据进行逆快速傅立叶变换并输出处理结果
有益效果:
相较于传统的时域L抽头的空时自适应处理STAP,本发明将LM×LM维矩阵运算降至2M×2M维矩阵,减小了运算量,提高了运算实时性,其中M是天线阵元数量。同时,对部分频点的弃用,运算实时性进一步提高。并且,因为二阶FFT、时域加窗和重叠的引入,缓解了现有SFAP零限杂乱的问题。第三,自适应对角加载的引用提高了抗干扰处理的数值稳健性,抗干扰性能得到提高。
附图说明
图1是本发明方法的流程图。
具体实施方式
下面结合附图对本发明进行详细描述
本发明采用二阶空频自适应处理(Space Frequency Adaptive Processing,SFAP)方法,通过FFT(快速傅里叶变换)将输入信号划分为N个频点,并在每个频点上采用矩阵求逆的方法计算各个阵元的最优权值。缓存FFT后的数据,与该数据计算出的最优权值做加权处理,再进行IFFT(快速傅里叶逆变换)后达到抑制干扰的目的。
假设阵列有M个阵元,每个频点的操作在2M×2M维矩阵中进行,包括协方差矩阵的统计平均计算、矩阵求逆和权值归一化三个数字信号处理步骤。本发明选取N=64,为减小运算量,提高抗干扰实时性,只计算前13个频点和后13个频点的最优权值,其余频点的最优权值取零。
同时,为减少频谱泄漏带来的影响,除采用二阶FFT以外,本发明采用时域加窗和25%叠加的方法,使宽带干扰达到理想中的均匀分割,缓解了SFAP比STAP(SpaceFrequency Adaptive Processing)零限较杂乱的问题。
而且,当干扰较大时,各阵元间的自适应矩阵病态度高,数值稳健性差。本发明采用对协方差矩阵进行自适应对角加载的方法,降低矩阵求逆的误差,提高抗干扰性能。
具体步骤如下:
1)对M个阵元的中频采样处理,选择合适的N点,本设计N=64。每N个点形成一个数据段,对每个数据段进行Hamming窗加窗处理得到加窗数据段;
2)相邻的每个数据段做25%重叠。即把每个数据段的最后25%×N点数据和下一段数据的前25%×N点数据重叠相加,并用加和之后的数据替换当前加窗数据段的后25%×N点数据;
3)对个数据段做N点FFT得到对应的频点数据;
其中m=1,2,……M,k=1,2,……N.
xm(n)为对应的时域中,第M个阵元后的第n个采样值,ωn是时域窗函数的第n个系数,而Xm(k)为频域中第m个阵元后的第k个采样值。
4)将频点数据缓存在双口RAM中;
5)计算协方差矩阵;
M个阵元组成的第k个频点的矢量为:
X(k)=(x1(k),x2(k),……xM(k))T,增加为二阶矢量得到:
X’(k)=(x1(k),x1(k+1),x2(k),x2(k+1),……xM(k),xM(k+1))T,
协方差矩阵为R(k)=E{X’(k)X’h(k)},是一个2M×2M维矩阵。期望的实现形式是对170个协方差矩阵的统计平均。为简化运算,只计算k=[1…13]和k=[52…64]区间内的协方差矩阵。
6)取出每个频点协方差矩阵的第一行第一列数据,将其除以8192,将结果与该协方差矩阵对角线上的元素相加;
7)对协方差矩阵进行FPGA硬件矩阵求逆,并根据线性约束最小方差准则(LCMV),求得k=[1…13]和k=[52…64]频点的最优权值;
其中k=1,2,……13,52,53,……,64。
u(k)是信号从期望方向入射时,在第k个频点形成的导向矢量,最优权值h(k)是一个2M维矢量。
8)取步骤4中RAM中的数据,与7步骤计算得到的权值做加权处理。注意时序控制,确保利用的最优权值是与其进行加权处理的频点数据计算得到。
其中k=1,2,……13,52,53,……,64。
Claims (8)
1.一种基于空频自适应滤波的卫星导航抗干扰方法,其特征在于,包括:
接收M个阵元的中频数据;
在每个阵元的中频数据中,每N个点形成一个数据段,对每个数据段进行汉明窗加窗处理得到加窗数据段;
将每个加窗数据段的后25%×N点数据与后一加窗数据段的前25%×N点数据相加,并用加和之后的数据替换当前加窗数据段的后25%×N点数据得到重叠数据段;
对每个重叠数据段进行N点快速傅立叶变换得到对应的频点数据;
将所述频点数据缓存在双口RAM中;
利用所述频点数据,计算与所述频点数据中的每个数据段中的前13个频点以及后13个频点对应的协方差矩阵;
对所述协方差矩阵进行自适应对角加载得到加载矩阵;
对所述加载矩阵进行求逆得出逆矩阵;
利用逆矩阵,根据线性约束最小方差准则求得所述频点数据的最优权值;
从所述双口RAM中读取所述频点数据;
利用所述最优权值对所述频点数据进行加权处理得到加权数据;
对加权数据进行逆快速傅立叶变换并输出处理结果,其中,所述M和N均为自然数。
2.根据权利要求1所述的方法,其特征在于,对所述加载矩阵进行求逆是利用FPGA对所述加载矩阵进行求逆。
3.根据权利要求1或2所述的方法,其特征在于,所述协方差矩阵是二阶的。
4.根据权利要求1或2中任一项所述的方法,其特征在于,对所述协方差矩阵进行自适应对角加载是指对每个协方差矩阵,自适应地加上第一行第一列数据除以8192的数值。
5.根据权利要求1或2中任一项所述的方法,其特征在于,根据线性约束最小方差准则求得所述频点数据的最优权值,包括:
对所述频点数据中的每个数据段中的第14至第N-13个频点的权值取零。
6.根据权利要求1或2中任一项所述的方法,其特征在于,所利用的最优权值是与其进行加权处理的频点数据计算得到。
7.根据权利要求1或2中任一项所述的方法,其特征在于,所述N为64。
8.一种基于空频自适应滤波的卫星导航抗干扰装置,其特征在于,包括:
接收模块,用于接收M个阵元的中频数据;
加窗模块,用于在每个阵元的中频数据中,每N个点形成一个数据段,对每个数据段进行汉明窗加窗处理得到加窗数据段;
叠加模块,用于将每个加窗数据段的后16点数据与后一加窗数据段的前16点数据相加并用加和之后的数据替换当前加窗数据段的后16点数据得到重叠数据段;
变换模块,用于对每个重叠数据段进行N点快速傅立叶变换得到对应的频点数据;
存储模块,用于将所述频点数据缓存在双口RAM中;
协方差矩阵计算模块,用于利用所述频点数据,计算与所述频点数据中的每个数据段中的前13个频点以及后13个频点对应的协方差矩阵;
对角加载模块,用于对所述协方差矩阵进行自适应对角加载得到加载矩阵;
求逆模块,用于对所述加载矩阵进行求逆得出逆矩阵;
最优权值计算模块,用于利用逆矩阵,根据线性约束最小方差准则求得所述频点数据的最优权值;
读取模块,用于从所述双口RAM中读取所述频点数据;
加权处理模块,用于利用所述最优权值对所述频点数据进行加权处理得到加权数据;
逆变换模块,用于对加权数据进行逆快速傅立叶变换并输出处理结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910748274.1A CN110501727B (zh) | 2019-08-13 | 2019-08-13 | 一种基于空频自适应滤波的卫星导航抗干扰方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910748274.1A CN110501727B (zh) | 2019-08-13 | 2019-08-13 | 一种基于空频自适应滤波的卫星导航抗干扰方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110501727A true CN110501727A (zh) | 2019-11-26 |
CN110501727B CN110501727B (zh) | 2023-10-20 |
Family
ID=68587336
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910748274.1A Active CN110501727B (zh) | 2019-08-13 | 2019-08-13 | 一种基于空频自适应滤波的卫星导航抗干扰方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110501727B (zh) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104280745A (zh) * | 2014-10-28 | 2015-01-14 | 北京北斗星通导航技术股份有限公司 | 一种改进的空频信号处理方法及处理装置 |
CN104698472A (zh) * | 2015-03-18 | 2015-06-10 | 北京北斗星通导航技术股份有限公司 | 改进宽带干扰抑制处理方法及处理装置 |
CN104716982A (zh) * | 2015-03-25 | 2015-06-17 | 北京北斗星通导航技术股份有限公司 | 一种扩频系统稳健抗干扰处理方法和装置 |
CN104898132A (zh) * | 2015-05-14 | 2015-09-09 | 电子科技大学 | 一种阈值处理和空频自适应算法结合的导航抗干扰算法 |
CN105005035A (zh) * | 2015-06-25 | 2015-10-28 | 西安电子科技大学 | 基于二维滑窗稳健空时自适应处理的目标检测方法 |
WO2016119388A1 (zh) * | 2015-01-30 | 2016-08-04 | 华为技术有限公司 | 一种基于语音信号构造聚焦协方差矩阵的方法及装置 |
CN107907891A (zh) * | 2017-11-16 | 2018-04-13 | 北京卫星信息工程研究所 | 卫星导航应用系统中空频抗干扰方法及应用系统 |
US10243761B1 (en) * | 2018-03-29 | 2019-03-26 | Rockwell Collins, Inc. | Scalable space frequency adaptive processing (SFAP) |
-
2019
- 2019-08-13 CN CN201910748274.1A patent/CN110501727B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104280745A (zh) * | 2014-10-28 | 2015-01-14 | 北京北斗星通导航技术股份有限公司 | 一种改进的空频信号处理方法及处理装置 |
WO2016119388A1 (zh) * | 2015-01-30 | 2016-08-04 | 华为技术有限公司 | 一种基于语音信号构造聚焦协方差矩阵的方法及装置 |
CN104698472A (zh) * | 2015-03-18 | 2015-06-10 | 北京北斗星通导航技术股份有限公司 | 改进宽带干扰抑制处理方法及处理装置 |
CN104716982A (zh) * | 2015-03-25 | 2015-06-17 | 北京北斗星通导航技术股份有限公司 | 一种扩频系统稳健抗干扰处理方法和装置 |
CN104898132A (zh) * | 2015-05-14 | 2015-09-09 | 电子科技大学 | 一种阈值处理和空频自适应算法结合的导航抗干扰算法 |
CN105005035A (zh) * | 2015-06-25 | 2015-10-28 | 西安电子科技大学 | 基于二维滑窗稳健空时自适应处理的目标检测方法 |
CN107907891A (zh) * | 2017-11-16 | 2018-04-13 | 北京卫星信息工程研究所 | 卫星导航应用系统中空频抗干扰方法及应用系统 |
US10243761B1 (en) * | 2018-03-29 | 2019-03-26 | Rockwell Collins, Inc. | Scalable space frequency adaptive processing (SFAP) |
Non-Patent Citations (4)
Title |
---|
NIU QIANGJUN等: "Wide-gap Frequency-Hopping Sequences Generated By Improved Dual Sub-band Method", 《TELECOMMUNICATION ENGINEERING》, vol. 59, no. 1 * |
张提升;郭文飞;郑建生;: "GNSS接收机中一种新的射频干扰抑制级联方法", 宇航学报, no. 07 * |
郭文飞: "抗干扰GPS接收系统关键技术研究与实现", 《中国博士学位论文全文数据库》 * |
郭文飞: "抗干扰GPS接收系统关键技术研究与实现", 《中国博士学位论文全文数据库》, 30 June 2014 (2014-06-30), pages 008 - 1 * |
Also Published As
Publication number | Publication date |
---|---|
CN110501727B (zh) | 2023-10-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102608588B (zh) | 基于子带分解的宽带子阵自适应波束形成方法 | |
CN101369014B (zh) | 应用于多输入多输出雷达的双边约束自适应波束形成方法 | |
CN108462521A (zh) | 自适应阵列天线的抗干扰实现方法 | |
Xiao et al. | MVDR algorithm based on estimated diagonal loading for beamforming | |
CN107966677B (zh) | 一种基于空间稀疏约束的圆阵模态域方位估计方法 | |
CN1752771A (zh) | 单通道相控阵接收信号重构及空间信号处理方法 | |
CN113156383B (zh) | 基于fpga的雷达信号抗干扰方法、系统及数据处理设备 | |
CN111880198B (zh) | 基于交替极化敏感阵列的空时极化抗干扰方法 | |
CN109633692A (zh) | Gnss导航卫星信号抗干扰处理方法 | |
CN111693971A (zh) | 一种用于弱目标检测的宽波束干扰抑制方法 | |
CN109188373B (zh) | 基于子阵阻塞矩阵预处理的抗主瓣干扰方法 | |
CN108828586B (zh) | 一种基于波束域的双基地mimo雷达测角优化方法 | |
CN106443672B (zh) | 一种方位多通道sar信号自适应重构方法 | |
CN110501727A (zh) | 一种基于空频自适应滤波的卫星导航抗干扰方法 | |
CN104346532A (zh) | 一种mimo雷达降维自适应波束形成方法 | |
CN114563760B (zh) | 一种基于sca阵型的二阶超波束形成方法、设备及介质 | |
CN109597034A (zh) | 一种基于欧几里得距离的空时自适应处理方法 | |
CN111965617B (zh) | 基于gpu的时分mimo雷达信号处理方法 | |
CN110208830B (zh) | 一种基于空时二维稀疏阵列的导航抗干扰方法 | |
CN115085784A (zh) | 一种分级加权稀疏子阵天线阵列稳健波束合成方法 | |
CN110146847B (zh) | 一种基于伪框架理论的强干扰信号抑制方法 | |
CN107135026A (zh) | 未知互耦存在时基于矩阵重构的稳健波束形成方法 | |
Li et al. | Adaptive pulse compression of MIMO radar based on GSC | |
CN102006105A (zh) | 深空接收天线组阵相关加权方法及系统 | |
CN107515409B (zh) | 抗脉冲干扰的调零天线波束合成方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |