CN106338741B - 一种基于迭代思想的空时导航抗干扰方法 - Google Patents
一种基于迭代思想的空时导航抗干扰方法 Download PDFInfo
- Publication number
- CN106338741B CN106338741B CN201610727439.3A CN201610727439A CN106338741B CN 106338741 B CN106338741 B CN 106338741B CN 201610727439 A CN201610727439 A CN 201610727439A CN 106338741 B CN106338741 B CN 106338741B
- Authority
- CN
- China
- Prior art keywords
- time
- iteration
- interference
- navigation
- vector
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 47
- 239000011159 matrix material Substances 0.000 claims abstract description 16
- 238000012545 processing Methods 0.000 claims description 11
- 238000004891 communication Methods 0.000 abstract description 3
- 238000012937 correction Methods 0.000 abstract description 3
- 238000004422 calculation algorithm Methods 0.000 description 8
- 238000004088 simulation Methods 0.000 description 4
- 241000712899 Lymphocytic choriomeningitis mammarenavirus Species 0.000 description 3
- 230000003044 adaptive effect Effects 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 238000010276 construction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
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
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radio Transmission System (AREA)
Abstract
本发明涉及属于卫星导航通信领域,提供一种基于迭代思想的空时导航抗干扰方法,用于克服现有空时导航抗干扰方法在干扰数目和时间抽头增加的情况下实时性急剧变差的问题;本发明首先根据输入信号计算自相关矩阵RX,然后计算计算误差向量e,最后根据迭代公式:w(k+1)=Cwk+s+μe求得最优导航权向量。本发明在最优抗干扰权向量求解过程中,在高斯‑赛德尔迭代的过程中加入修正量μe,大大提升迭代速度,从而有效提高空时导航抗干扰的实时性,保证实际应用需求。
Description
技术领域
本发明涉及属于卫星导航通信领域,特别涉及卫星导航***中抗干扰方法的优化,具体为一种基于迭代思想的空时导航抗干扰方法。
背景技术
全球卫星导航***利用导航卫星发射导航卫星信号,地面、海洋和空间的用户通过接收机接收卫星导航信息从而享受连续、实时、高精度的定位、导航和授时服务。近几十年,卫星导航***发展迅速,在许多领域得到广泛应用,同时产生了显著地社会效益和经济效益。建设和发展导航***能够在一定程度上促进国家的政治、经济和文化建设,所以各国都投入大量人力、物力和财力对卫星导航***进行研究。
通常导航卫星均被设计成发射功率仅有几毫瓦的弱信号卫星,这对于卫星降低造价和延长使用寿命是很有效的。但由于卫星信号弱,所以非常容易受到干扰。除了恶意的干扰外,一些频率较高的商业电视台、航空卫星通信和机动卫星***终端都可能削弱导航卫星信号,而且自然界所发生的一些现象也会引起信号干扰。
一旦卫星信号被干扰就可能中断其使用,使定位误差增大甚至完全无法实现导航功能。目前,社会生活对导航技术的依赖越来越大,导航接收机抗干扰性能的要求也越来越高。专门针对导航接收机的抗干扰技术的研究就具有了现实意义。
目前针对导航抗干扰采取的主要方法是空时自适应信号处理(STAP,Space TimeAdapt ive Processing)的办法,传统的STAP算法的结构如图1所示:
天线阵元个数为M,时间延迟单元数为P,时延间隔为Δ(一般要求Δ<1/B,B是信号处理带宽),假设第m个阵元上的输入信号为xm(n),则阵元m的第p个延迟抽头数接收到的信号为:
xm(n-(p-1)Δ)(m=1,2,...,M;p=1,2,...,P) (1)
若第m个阵元的第p个抽头的权值为wmp,则滤波器的输出为:
空时二维滤波的输出信号可表示为:
y=wHx (3)
其中,
w=[w11,...,w1P,...,wM1,...,wMP]T (4)
x=[x1(n),...,x1(n-(P-1)Δ),...,xM(Δ),...xM(n-(P-1)Δ)]T (5)
将上述权向量w中第i个元素表示为wi,则w=[w1,w2,…,wMP]T;
由此可以看出,纯空域天线阵列在时域上扩展后便构成了空时二维联合处理器;对于每个天线阵元来说,增加的延时抽头数可构成时域滤波器,可对不同频率的干扰源形成分辨能力;而对于不同天线阵元的相同的延迟节点来说,其天线阵元则构成了纯空域滤波器,从而对不同来向的空间干扰源形成分辨力;而对整个空时二维滤波器来说,在增加滤波器自由度的同时,又使天线阵元个数保持不变,在不以牺牲硬件成本为代价的同时,又能比单纯的空域滤波滤除更多的干扰;
对于M阵元的均匀线阵,若每个阵元接收到的信号的快拍数为N,而时域滤波器的延迟单元数为P,构成的空时二维滤波器的滤波器阶数为M×P,则在进行算法仿真的时候,空时二维滤波器的快拍数则为N-P+1,空时二维接收到的数据的表示见式(6);其中的每一列则依次作为第n时刻的空时二维滤波器的输入信号:
由Frost提出的线性约束方差准则(LCMV)的基本思想就是使滤波器的输出功率(即信号方差)最小,即
其中,Pout为滤波器输出功率,E(X)表示X的期望(均值);
若直接求其最小输出功率,但是没有约束条件,则会在w=0时取得,此时的输出信号y(n)也变为0,没有了实际工程意义;因此需加上一种约束条件,即
wHCX=f (8)
其中,CX是约束矩阵,f为相应的约束响应向量;
构造实值代价函数:
L(w)=wHRXw-λ(fH-CX Hw) (9)
其中,λ为拉格朗日因子、是个常数,RX=E{X*XH}为输入信号数据的自相关矩阵;
令求得LCMV的最优权向量为
wLCMV=RX -1CX(CX HRX -1CX)-1fH (10)
最小方差无失真响应(MVDR)准则是在有用信号的约束矢量s已知的情况下,使其约束响应f等于固定值1;因此它是LCMV准则的一个特例,可描述为:
求得MVDR最优权向量为:
可以写为:
w=RX -1s (13)
RXw=s
由式(13)可知,最优抗干扰权值的求取等同于解一个MP维的线性方程组;在解线性方程组的数值方法中,包括直接求逆法与迭代法;直接求逆法运算过程复杂,且计算时容易产生不可控状态;在迭代计算方法中,简单迭代法利于并行实现,具体过程为:
由式(13)基于线性约束最小方差准则可得线性迭代方法的权值求解过程所用的公式为:
w=(I-RX)w+s=Cw+s (14)
其中,C=I-RX,为迭代算法的系数矩阵;I为单位阵;M为天线阵元数目,P为时间单元数目,s为MP×1维的约束向量;简单迭代法可表示为:
式中,cij为矩阵C的第i行、第j列元素,为w=[w1,w2,…,wMP]T中第i个权值的第k次迭代结果,si为s中的第i个元素的值;
还有一种高斯—赛德尔迭代法可以用于解式(13),其算法分析如下:
解方程RXw=s中,将系数矩阵RX中主对角元***,rij为RX中第i行、第j列元素,设RX=D-L-U;
其中:D=diag(r11,r22,…,rMPMP) (16)
高斯-赛德尔迭代法的矩阵形式为:
也可以写为:
其中
高斯—赛德尔迭代法可表示为:
在迭代计算方法中,简单迭代法利于并行实现,但是迭代过程可能会发散;而高斯—赛德尔迭代法对于输入信号自相关矩阵(对称正定)而言一定收敛,其缺点在于并行实现能力较弱;因此,提出一种将简单迭代法与高斯—赛德尔迭代法相结合的线性化迭代方法,该算法可实现以上两种方法的优势互补,而简单迭代法与高斯—赛德尔迭代法相结合的线性化迭代方法具体计算公式为:
其中,A为同时进行更新的权值数目(A>1),有MP/A=Z,Z为将空时处理维数MP分成的段数;其迭代算法方案如图2所示。
上述方法在空时二维处理维数(空时处理维数等于天线阵元数目M与时间单元数目P的乘积)增大时,虽可获得较好抗干扰效果,并增加可抗干扰的数目;但是其运算量急剧增加,大运算量意味着算法需要很长的运算时间和很高的***复杂度,这将影响自适应算法的实时性,即导致空时导航抗干扰方法的实时性急剧变差;针对这个问题,我们提出了一种新的直接解线性方程组的方法。
发明内容
本发明的目的在于针对现有空时导航抗干扰方法在干扰数目和时间抽头增加的情况下实时性急剧变差的问题,提供一种基于迭代思想的空时导航抗干扰方法,该方法在求解最优抗干扰权向量过程中,在高斯-赛德尔迭代的过程中加入合适的修正量以加快迭代速度。
为实现上述目的,本发明采取的技术方案为:
一种基于迭代思想的空时导航抗干扰方法,包括以下步骤:
步骤a.将每个阵元的输入信号进行处理,求得自相关矩阵RX;
步骤b.设定初始权值,令权向量w中所有元素为0:w=0;
步骤c.计算误差向量e:e=s-RXw,其中s为约束矢量;
步骤d.按照以下迭代公式进行迭代:
其中,μ为步长因子,满足λmax为自相关矩阵RX的最大特征值;ei为误差向量e的第i个元素。
本发明的有益效果在于:
本发明提供一种基于迭代思想的空时导航抗干扰方法,该方法在最优抗干扰权向量的迭代求解过程中,引入修正量μe,从而大大降低了迭代次数,从而有效提高空时导航抗干扰的实时性,保证实际应用需求。
附图说明
图1为现有的空时自适应算法结构框图。
图2为现有空时权矢量w的迭代过程示意图。
图3为本发明与现有空时权矢量w的迭代速度的比较图。
具体实施方式
以下结合实施例和附图对本发明的上述内容再作进一步的详细说明,但不应将此理解为本发明上述主题的范围仅限于以下的实例;在不脱离本发明上述技术思想情况下,根据本领域普通技术知识和惯用手段做出的各种替换或变更,均应包括在本发明的范围内。本发明说明书中未作详细描述的内容属于本领域专业技术人员的公知技术。
本实施例提供基于迭代思想的空时导航抗干扰方法,包括以下步骤:
a.将每个阵元的输入信号进行处理,求得自相关矩阵RX;
RX=E{X*XH}
b.设定初始权值,令权向量w中所有元素为0;
c.计算出误差向量e:e=s-RXw,其中s为约束矢量;
d.在迭代中加入修正量μe,w(k+1)=Cwk+s+μ[s-RXwk]=Cwk+s+μe;其中,μ为步
长因子,满足λmax为自相关矩阵RX的最大特征值;
其中,ei为误差向量e的第i个元素。
在该方案中,将时空处理维数分成若干段,每一段含有A个数据(A是可选的);在每一次权值更新时,只计算本次数据相乘时的乘法结果和误差,而乘法结果的累加值在下一次权值更新时进行计算,这是因为浮点乘法与浮点加法在硬件中执行一次所需要的时间基本相同;这样在每次计算权值时,实际上是使用了前两次及所有以前的权值更新结果,而前一次的权值更新结果尚未计算出来;这样做既充分利用硬件资源,又最大程度地提高了时空抗干扰处理的实时性。
采用MATLAB对本实施例进行仿真测试,并设置现有空时权矢量w的迭代方法作为对比,本实施例设定的仿真参数为:
自相关矩阵RX和约束矢量s分别为:
步长因子μ=0.03,A=Z=2,其仿真结果如图3所示,从图中可以看到,本发明迭代过程中误差达到足够小所需迭代次数要远小于现有的迭代方法的,从而可知,本发明基于迭代思想的空时导航抗干扰方法在在干扰数目和时间抽头增加的情况下具备更好的实时性。
以上所述,仅为本发明的具体实施方式,本说明书中所公开的任一特征,除非特别叙述,均可被其他等效或具有类似目的的替代特征加以替换;所公开的所有特征、或所有方法或过程中的步骤,除了互相排斥的特征和/或步骤以外,均可以任何方式组合。
Claims (1)
1.一种基于迭代思想的空时导航抗干扰方法,包括以下步骤:
步骤a.将每个阵元的输入信号进行处理,求得自相关矩阵RX;
步骤b.设定初始权值,令权向量w中所有元素为0:w=0;
步骤c.计算误差向量e:e=s-RXw,其中s为约束矢量;
步骤d.按照以下迭代公式进行迭代:
其中,μ为步长因子,满足λmax为自相关矩阵RX的最大特征值;ei为误差向量e的第i个元素;M为天线阵元数目,P为时间单元数目,A为同时进行更新的权值数目,Z为将空时处理维数MP分成的段数;cij为系数矩阵C的第i行、第j列元素,为w中第j个权值的第k次迭代结果,si为约束矢量s中的第i个元素的值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610727439.3A CN106338741B (zh) | 2016-08-25 | 2016-08-25 | 一种基于迭代思想的空时导航抗干扰方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610727439.3A CN106338741B (zh) | 2016-08-25 | 2016-08-25 | 一种基于迭代思想的空时导航抗干扰方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106338741A CN106338741A (zh) | 2017-01-18 |
CN106338741B true CN106338741B (zh) | 2019-10-01 |
Family
ID=57825313
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610727439.3A Active CN106338741B (zh) | 2016-08-25 | 2016-08-25 | 一种基于迭代思想的空时导航抗干扰方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106338741B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100965100B1 (ko) * | 2008-01-15 | 2010-06-25 | 중앙대학교 산학협력단 | 적응적 빔포밍 장치 및 그 방법 |
CN101937088A (zh) * | 2010-07-20 | 2011-01-05 | 中国电子科技集团公司第五十四研究所 | 一种卫星导航接收设备空时二维抗干扰实现方法 |
CN103944624A (zh) * | 2014-03-25 | 2014-07-23 | 电子科技大学 | 基于迭代算法的稳健波束形成方法 |
CN104808219A (zh) * | 2015-04-24 | 2015-07-29 | 太原理工大学 | 一种新型的空时联合抗干扰方法 |
CN104914451A (zh) * | 2015-05-05 | 2015-09-16 | 中国人民解放军国防科学技术大学 | 一种块Toeplitz矩阵低复杂度求逆的空时抗干扰方法 |
-
2016
- 2016-08-25 CN CN201610727439.3A patent/CN106338741B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100965100B1 (ko) * | 2008-01-15 | 2010-06-25 | 중앙대학교 산학협력단 | 적응적 빔포밍 장치 및 그 방법 |
CN101937088A (zh) * | 2010-07-20 | 2011-01-05 | 中国电子科技集团公司第五十四研究所 | 一种卫星导航接收设备空时二维抗干扰实现方法 |
CN103944624A (zh) * | 2014-03-25 | 2014-07-23 | 电子科技大学 | 基于迭代算法的稳健波束形成方法 |
CN104808219A (zh) * | 2015-04-24 | 2015-07-29 | 太原理工大学 | 一种新型的空时联合抗干扰方法 |
CN104914451A (zh) * | 2015-05-05 | 2015-09-16 | 中国人民解放军国防科学技术大学 | 一种块Toeplitz矩阵低复杂度求逆的空时抗干扰方法 |
Non-Patent Citations (3)
Title |
---|
导航终端空时二维联合抗干扰实现方法;邓志鑫;《无线电工程》;20121231;第42卷(第5期);第30-33页 * |
近似最佳加速因子改进高斯-赛德尔法潮流计算;易水平等;《江西电力》;20081231;第5-7页 * |
通信***中空时二维抗干扰方法研究;杨亮;《综合电子信息技术》;20131231;第39卷(第6期);第74-77页 * |
Also Published As
Publication number | Publication date |
---|---|
CN106338741A (zh) | 2017-01-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wen et al. | Angle estimation and mutual coupling self-calibration for ULA-based bistatic MIMO radar | |
CN103941220B (zh) | 一种基于稀疏重构的网格外目标波达方向估计方法 | |
CN103245956B (zh) | 一种基于稳健波束形成算法的gps抗多径方法 | |
CN108562866B (zh) | 基于矩阵填充的双基地mimo雷达角度估算方法 | |
CN111337893A (zh) | 一种基于实值稀疏贝叶斯学习的离格doa估计方法 | |
CN108717174B (zh) | 基于信息论的预测快速协方差交互融合无源协同定位方法 | |
Zhang et al. | Robust beamforming for coherent signals based on the spatial-smoothing technique | |
CN114338301B (zh) | 一种基于压缩感知的ris辅助毫米波***的信道估计方法 | |
CN108111208A (zh) | 一种灵活码本配置的毫米波***角度域信道估计方法 | |
CN108872932A (zh) | 基于神经网络的超视距目标直接定位结果纠偏方法 | |
CN114172547A (zh) | 基于智能反射面的无线携能通信混合预编码设计方法 | |
Rahman et al. | Ising model formulation of outlier rejection, with application in wifi based positioning | |
CN111337873A (zh) | 一种基于稀疏阵的doa估计方法 | |
CN104539340A (zh) | 一种基于稀疏表示和协方差拟合的稳健波达角估计方法 | |
CN107255796A (zh) | 一种非均匀噪声下窄带近场信号源定位方法 | |
CN111580042B (zh) | 一种基于相位优化的深度学习测向方法 | |
CN110954860B (zh) | 一种doa和极化参数估计方法 | |
CN113671439A (zh) | 基于非均匀智能超表面阵列的无人机集群测向***及方法 | |
CN103378888A (zh) | 下行波束赋形方法及设备 | |
CN104716982A (zh) | 一种扩频***稳健抗干扰处理方法和装置 | |
CN105068090B (zh) | 一种利用单快拍数据进行gnss天线阵干扰抑制的方法 | |
Yang et al. | Underdetermined DOA estimation for moving array with reduced mutual coupling in unknown nonuniform noise environment | |
CN104698472A (zh) | 改进宽带干扰抑制处理方法及处理装置 | |
CN106338741B (zh) | 一种基于迭代思想的空时导航抗干扰方法 | |
CN104777492A (zh) | 卫星导航***抗干扰处理方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |