CN104199053B - 一种基于卫星信号到达角约束的稳健波束形成方法 - Google Patents
一种基于卫星信号到达角约束的稳健波束形成方法 Download PDFInfo
- Publication number
- CN104199053B CN104199053B CN201410486111.8A CN201410486111A CN104199053B CN 104199053 B CN104199053 B CN 104199053B CN 201410486111 A CN201410486111 A CN 201410486111A CN 104199053 B CN104199053 B CN 104199053B
- Authority
- CN
- China
- Prior art keywords
- signal
- satellite
- angle
- array antenna
- convex optimization
- 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
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)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
本发明涉及一种基于卫星信号到达角约束的稳健波束形成方法,其特征在于:步骤1:建立阵列天线接收信号模型,计算卫星信号和干扰信号的导向矢量;步骤2:根据卫星信号到达角的范围,确定阵列天线幅值响应约束条件;步骤3:将步骤2中的非凸优化约束条件转化为凸优化约束条件;步骤4:根据噪声功率来确定阵列天线所接收信号协方差矩阵的对角加载因子γ;步骤5:建立基于角度约束的波束形成代价函数,采用凸优化工具求解计算阵列天线最优权值,获取卫星信号方向上的波束和干扰信号方向的零陷。
Description
技术领域
本发明涉及一种基于卫星信号到达角约束的稳健波束形成方法。
背景技术
由于GPS信号是在众所周知的频率上发射的,其调制特征广为人知,信噪比又比较低,因而很容易进行干扰或欺骗。GPS信号依靠其扩频体制能够为其获得43dB的处理增益,然而随着人为干扰技术的提升,卫星信号只依靠其扩频体制进行抗干扰已经不能满足用户需求。
目前,由于基于阵列天线的波束形成技术能够在空域对干扰进行很好的抑制,因此,该技术得到了许多学者的关注。天线阵在接收到信号后,通过对天线阵的每个阵元上引入可以调整的权系数,来实现对天线阵主波束方向的灵活控制,达到对整个空域信号的接收,这个天线阵元上的权系数的调整也即实现了“空域滤波”功能。
然而,基于阵列天线的波束形成技术在GPS抗干扰领域仍然存在一些问题,如载体与卫星之间的相对移动,载体上的卫星接收机所形成的波束不能够准确的对准空间卫星位置。从而,阵列天线可能对卫星信号造成失真接收,对基带信号处理部分的捕获和跟踪造成严重影响。因此,一些稳健的波束形成算法被提出来解决此类问题。H.Cox等人提出了一种对角加载的方法(Robustadaptivebeamforming,IEEETransactionsonAcoustic,Speech,andSignalProcessing,35,October,1987,1365–1376.),通过对接收信号的协方差矩阵进行修正来提高波束形成算法的稳健性。L.Chang和C.C.Yeh等人提出了一种子空间投影的方法(PerformanceofDMIandeigenspace-basedbeamformers,IEEETransactionsonAntennasandPropagation,40,November,1992,1336–1347.)。这种方法将接收到的信号投影到与干扰信号相正交的子空间来抑制干扰,从而增强波束形成的稳健性。当信噪比比较低的时候,此类方法容易产生信号子空间交换,得到的与干扰信号相正交的子空间不够准确。现有稳健的方法均是在以上两类方法的基础之上发展而来的,但是以上方法主要考虑的是对干扰信号进行抑制,没有对期望信号进行约束,因此,期望信号容易失真接收。
发明内容
本发明目的在于提供一种基于卫星信号到达角约束的稳健波束形成方法,能够实现卫星信号的无失真接收。
实现本发明目的技术方案:
一种基于卫星信号到达角约束的稳健波束形成方法,其特征在于:
步骤1:建立阵列天线接收信号模型,计算卫星信号和干扰信号的导向矢量;
步骤2:根据卫星信号到达角的范围,确定阵列天线幅值响应约束条件;
步骤3:将步骤2中的非凸优化约束条件转化为凸优化约束条件;
步骤4:根据噪声功率来确定阵列天线所接收信号协方差矩阵的对角加载因子γ;:
步骤5:建立基于角度约束的波束形成代价函数,采用凸优化工具求解计算阵列天线最优权值,获取卫星信号方向上的波束和干扰信号方向的零陷。
步骤2中所说的幅值响应约束条件为,
L(θ)≤|wHa(θ)|≤U(θ),θ∈Θ
式中,|·|和()H分别为绝对值和共轭转置操作,θ为期望卫星信号到达角,Θ为θ的一个取值范围,w为天线阵列权值,L(θ)和U(θ)分别为幅值响应|wHa(θ)|的下限值和上限值。
步骤4中,对角加载因子γ=10δ2,δ2为噪声功率。
步骤5中,基于角度约束的波束形成代价函数表示为
L(θ)≤wHa(θ)≤U(θ),θ∈Θ
式中,为天线阵列接收信号的协方差矩阵,I为与同维数的单位矩阵,γ为对角加载因子。
本发明具有的有益效果:
本发明根据卫星信号到达角的范围,来确定阵列幅值响应区间,将此阵列幅值响应区间作为一个约束条件,从而保证了卫星信号的无失真接收。本发明将非凸优化形式的阵列幅值响应约束条件转化为凸优化的形式,使用CVX工具箱进行求解,有效简化算法,避免了采用其他复杂的优化算法进行求解。本发明直接根据噪声功率来确定阵列天线所接收信号协方差矩阵加载因子,避免了采用H.Cox等人所提出的复杂对角加载算法,进一步简化了算法。
附图说明
图1是本发明基于卫星信号到达角约束的稳健波束形成方法流程图;
图2是基于最小方差准则与角度约束的阵列天线示意图;
图3是波束宽度=10°时,阵列输出的信干噪比随着输入信噪比的变化曲线图。
具体实施方式
如图1所示,本发明基于卫星信号到达角约束的稳健波束形成方法包括以下步骤:
步骤1:建立阵列天线接收信号模型,计算卫星信号和干扰信号的导向矢量;
1)阵列天线接收卫星信号、干扰信号和空间白噪声,建立阵列天线接收信号模型
假设有M个天线阵元,相邻阵元的间隔为λ/2,λ是GPS信号的波长。假设期望信号和干扰信号的到达角度分别为θ0和θk(k=1、2,…P),则阵列天线接收信号可以表示为
式(1)中,s0(t)代表卫星直达信号,sk(t)为第k个干扰信号,n(t)代表相互独立的零均值高斯白噪声,噪声功率表示为δ2,n(t)可以表示为n(t)=[n0(t),n1(t),…,np(t)]T。
2)由卫星信号和干扰信号的达到角,计算卫星信号和干扰信号的导向矢量;
在式(1)中,a(θ0)为卫星信号的导向矢量,a(θk)为第k个干扰信号的导向矢量。当θ0和θk(k=1、2,…P)为确定值时,a(θ0)和a(θk)可以分别表示为 和
3)在公式(1)中,GPS信号s0(t)为由伪随机码(C/A码)、L1载波(载波频率为1575.42MHz得余弦波)和导航数据电文组成,可以表示为
s0(t)=AC(t)D(t)cos(ωt+φ0)(2)
式(2)中,A为信号幅值,C(t)表示C/A码,D(t)表示导航电文数据,ω为L1载波角频率,φ0为L1频段初始载波相位。
步骤2:根据卫星信号到达角的范围,确定阵列天线幅值响应约束条件;
1)当卫星信号的到达角精确已知时;
当估计的期望卫星信号到达角准确时,基于最小方差准则的波束形成代价函数可以表示为
(3)
s.t.aH(θ0)w=1
式中,w为阵列天线权值,x(n)为离散化的x(t),N为采样点个数,(·)H为共轭转置,s.t.代表约束条件。
2)当卫星信号的到达角存在偏差时;
然而,由于载体与卫星之间的相对移动,假设的期望信号到达角度不再是确定的θ0,而是一个区间范围Θ,即θ0=θ∈Θ。图2给出了基于最小方差准则与角度约束的阵列天线示意图。
为了保证卫星信号在此空间角度都能够被无失真的接收,本发明中采用下式中的幅值响应约束
L(θ)≤|wHa(θ)|≤U(θ),θ∈Θ(4)
式中,|·|为绝对值操作,L(θ)和U(θ)分别为幅值响应|wHa(θ)|的下限值和上限值。因此,(3)中的代价函数可以表示为
(5)
L(θ)≤|wHa(θ)|≤U(θ),θ∈Θ
步骤3:将步骤2中的非凸优化约束条件转化为凸优化约束条件;
将步骤2中的非凸优化约束条件转化为凸优化形式,采用CVX工具箱进行求解;
在式(5)中,L(θ)≤|wHa(θ)|≤U(θ)为非凸优化约束,需要将其转化为凸优化约束,进而通过CVX工具箱进行计算。
以下凸优化变换是在天线阵元满足几何对称的条件下进行计算的。本发明考虑均匀线性阵列,现假设M为奇数阵元,以中间阵元为参考阵元。则卫星信号的导向矢量可以表示为
各阵元上面的权系数具有下面的关系
式(7)中,()*为共轭操作。从而,阵列幅值响应|wHa(θ)|可以表示为
式中,Re{·}为取实数部分。从式(8)中可以看出|wHa(θ)|=wHa(θ)。因此,代价函数(5)可以变换为
(9)
L(θ)≤wHa(θ)≤U(θ),θ∈Θ
式(9)中的约束条件为凸优化形式,可以采用CVX工具箱直接进行计算。
步骤4:根据噪声功率来确定阵列天线所接收信号协方差矩阵的对角加载因子γ;
由于卫星信号到达角的不确定,在进行(3)式运算时,协方差矩阵因a(θ0)的误差而不准确。先假定准确的协方差矩阵为因a(θ0)而导致的协方差矩阵误差为Δ,则
在(10)中,||Δ||≤γ,||·||为Frobenius范数操作,γ为较小的常数。所以(9)中的代价函数在最坏的情况下表示为
L(θ)≤wHa(θ)≤U(θ),θ∈Θ(11)
||Δ||≤γ
为了求解式(11),需要先求解
(12)
||Δ||≤γ
式(12)可以表示为
(13)
||Δ||≤γ
根据拉格朗日乘子法,Δ可以求解为这里γ取10δ2,δ2为噪声功率。所以式(11)中,
步骤5:建立基于角度约束的波束形成代价函数,采用凸优化工具求解计算阵列天线最优权值,获取卫星信号方向上的波束和干扰信号方向的零陷。
根据步骤2所确定的阵列幅值响应约束条件、步骤3所进行的凸优化变换和步骤4所计算的对角加载因子,建立基于卫星信号到达角约束的代价函数,并采用凸优化工具求解计算阵列最优权值,获取卫星信号方向上的波束和干扰信号方向的零陷。
1)采用凸优化工具箱计算最优阵列权值w;
根据步骤2、步骤3和步骤4的描述,将式(13)和(14)代入(11)式中,即可得到本发明所述方法的代价函数
(15)
L(θ)≤wHa(θ)≤U(θ),θ∈Θ
在数学应用方面,式(15)为凸优化形式的最小值问题。采用CVX最优算法求解式(15)中的阵列最优权值w。
2)根据最优化阵列权值,获取卫星信号方向上的波束和干扰信号方向的零陷。
在获得阵列最优权值之后,可以得到阵列的输出表达是为
y(t)=wHx(t)(16)
阵列权值w与卫星信号的导向矢量a(θ),θ∈Θ相乘即可获取卫星信号方向的波束F(θ),
F(θ)=wHa(θ),θ∈Θ(17)
式(17)表明,阵列天线能够在卫星信号的空间角度Θ范围内形成波束,从而保持卫星信号无失真接收。
同理,若干扰信号的到达角范围为Ω,阵列权值w与干扰信号的导向矢量a(θ),θ∈Ω相乘即可获取干扰信号方向的零陷G(θ),
G(θ)=wHa(θ),θ∈Ω(18)
式(18)表明,阵列天线能够在干扰信号的空间角度Ω范围内形成零陷,从而抑制干扰信号的接收。
为证明本发明的有效性,图3给出了当Θ=[-50,50],即波束宽度为10°时,本发明所提方法的阵列输出信干噪比(SINR)随着输入信噪比(SNR)的变化曲线,图3证明了本发明方法具有更好的输出SINR。
Claims (4)
1.一种基于卫星信号到达角约束的稳健波束形成方法,其特征在于:
步骤1:建立阵列天线接收信号模型,计算卫星信号和干扰信号的导向矢量;
步骤2:根据卫星信号到达角的范围,确定阵列天线幅值响应的非凸优化约束条件;
步骤3:将步骤2中的非凸优化约束条件转化为凸优化约束条件;
步骤4:根据噪声功率来确定阵列天线所接收信号协方差矩阵的对角加载因子γ;
步骤5:建立基于角度约束的波束形成代价函数,采用凸优化工具求解计算阵列天线最优权值,获取卫星信号方向上的波束和干扰信号方向的零陷。
2.根据权利要求1所述的基于卫星信号到达角约束的稳健波束形成方法,其特征在于:
步骤2中所述的阵列天线幅值响应的非凸优化约束条件为,
L(θ)≤|wHa(θ)|≤U(θ),θ∈Θ
式中,|·|和()H分别为绝对值和共轭转置操作,θ为期望卫星信号到达角,Θ为θ的一个取值范围,w为天线阵列权值,L(θ)和U(θ)分别为幅值响应|wHa(θ)|的下限值和上限值,a(θ)为卫星信号的导向矢量。
3.根据权利要求2所述的基于卫星信号到达角约束的稳健波束形成方法,
其特征在于:步骤4中,对角加载因子γ=10δ2,δ2为噪声功率。
4.根据权利要求3所述的基于卫星信号到达角约束的稳健波束形成方法,其特征在于:
步骤5中,基于角度约束的波束形成代价函数表示为
L(θ)≤wHa(θ)≤U(θ),θ∈Θ
式中,为天线阵列接收信号的协方差矩阵,I为与同维数的单位矩阵,γ为对角加载因子。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410486111.8A CN104199053B (zh) | 2014-09-22 | 2014-09-22 | 一种基于卫星信号到达角约束的稳健波束形成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410486111.8A CN104199053B (zh) | 2014-09-22 | 2014-09-22 | 一种基于卫星信号到达角约束的稳健波束形成方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104199053A CN104199053A (zh) | 2014-12-10 |
CN104199053B true CN104199053B (zh) | 2016-06-29 |
Family
ID=52084365
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410486111.8A Active CN104199053B (zh) | 2014-09-22 | 2014-09-22 | 一种基于卫星信号到达角约束的稳健波束形成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104199053B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104536018B (zh) * | 2015-01-06 | 2017-02-22 | 中国人民解放军国防科学技术大学 | 一种使用阵列天线抗干扰技术的gnss多星联合捕获方法 |
CN105655727B (zh) * | 2015-09-28 | 2018-07-24 | 中国电子科技集团公司第二十九研究所 | 一种基于凸优化的增益损失恒定波束的形成方法及装置 |
CN107306244B (zh) * | 2016-04-21 | 2021-08-31 | 中兴通讯股份有限公司 | 一种对角加载方法和装置 |
CN108663668B (zh) * | 2018-05-18 | 2022-03-22 | 西安电子科技大学 | 基于iaa的干扰加噪声协方差矩阵重构稳健波束形成方法 |
CN110895327B (zh) * | 2019-11-08 | 2022-10-14 | 电子科技大学 | 基于直接凸优化建模的鲁棒性自适应波束形成方法 |
CN112379334B (zh) * | 2020-11-02 | 2023-03-14 | 上海无线电设备研究所 | 一种自适应波束形成方法及装置 |
CN112543047B (zh) * | 2020-11-04 | 2022-02-22 | 西安交通大学 | 一种多波束卫星干扰抑制方法、存储介质及计算设备 |
CN113050030B (zh) * | 2021-03-26 | 2021-11-05 | 江苏蓝策电子科技有限公司 | 基于到达角度测距的定位方法及装置 |
CN113050127B (zh) * | 2021-04-13 | 2024-01-02 | 广州海格通信集团股份有限公司 | 信号处理方法、装置、计算机设备和存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6195037B1 (en) * | 1999-06-01 | 2001-02-27 | Motorola, Inc. | Method and apparatus for increased system capacity using antenna beamforming |
EP1717969A2 (en) * | 2005-04-27 | 2006-11-02 | Delphi Technologies, Inc. | Method and system for digital beam forming in communications systems |
CN103293517A (zh) * | 2013-05-13 | 2013-09-11 | 西安电子科技大学 | 基于脊参数估计的对角加载稳健自适应雷达波束形成方法 |
CN109344624A (zh) * | 2018-10-26 | 2019-02-15 | 深信服科技股份有限公司 | 基于云端协作的渗透测试方法、平台、设备及存储介质 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103944624B (zh) * | 2014-03-25 | 2017-07-18 | 电子科技大学 | 基于迭代算法的稳健波束形成方法 |
-
2014
- 2014-09-22 CN CN201410486111.8A patent/CN104199053B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6195037B1 (en) * | 1999-06-01 | 2001-02-27 | Motorola, Inc. | Method and apparatus for increased system capacity using antenna beamforming |
EP1717969A2 (en) * | 2005-04-27 | 2006-11-02 | Delphi Technologies, Inc. | Method and system for digital beam forming in communications systems |
CN103293517A (zh) * | 2013-05-13 | 2013-09-11 | 西安电子科技大学 | 基于脊参数估计的对角加载稳健自适应雷达波束形成方法 |
CN109344624A (zh) * | 2018-10-26 | 2019-02-15 | 深信服科技股份有限公司 | 基于云端协作的渗透测试方法、平台、设备及存储介质 |
Non-Patent Citations (2)
Title |
---|
一种新的稳健Capon波束形成算法;徐定杰等;《宇航学报》;20121031;第33卷(第10期);第1492-1497页 * |
基于范数优化的对角加载稳健自适应波束形成;史英春等;《电路与***学报》;20130228;第18卷(第1期);第218-223页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104199053A (zh) | 2014-12-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104199053B (zh) | 一种基于卫星信号到达角约束的稳健波束形成方法 | |
CN104199052B (zh) | 一种基于范数约束的波束旁瓣抑制方法 | |
Li et al. | A robust anti-jamming navigation receiver with antenna array and GPS/SINS | |
CN106569181A (zh) | 基于协方差矩阵重构稳健Capon波束形成的算法 | |
CN103245956B (zh) | 一种基于稳健波束形成算法的gps抗多径方法 | |
CN102353970A (zh) | 一种高抗干扰性能gps/sins组合导航***及实现方法 | |
CN103513250B (zh) | 一种基于鲁棒自适应波束形成原理的模基定位方法及*** | |
CN105629206B (zh) | 导向矢量失配下的机载雷达稳健空时波束形成方法及*** | |
CN109143193B (zh) | 一种基于二阶峰的高频地波雷达浪场反演方法 | |
CN104345300B (zh) | 杂波空时谱线性补偿的机载非正侧视阵雷达stap方法 | |
CN103954950A (zh) | 一种基于样本协方差矩阵稀疏性的波达方向估计方法 | |
CN103353596A (zh) | 基于压缩感知的波束空间域米波雷达测高方法 | |
CN104408278A (zh) | 一种基于干扰噪声协方差矩阵估计的稳健波束形成方法 | |
CN109799495A (zh) | 一种用于高保真阵列处理的宽带时延估计方法 | |
CN103267966A (zh) | 一种单个二维压差式矢量水听器方位估计的方法 | |
CN102393525A (zh) | 子空间投影的导航干扰抑制与信号增强方法 | |
Zheng et al. | Joint towed array shape and direction of arrivals estimation using sparse Bayesian learning during maneuvering | |
CN104049262A (zh) | 一种基于矢量跟踪环路的波束形成抗干扰方法 | |
CN103698753A (zh) | 一种小型阵列的无源通道校正方法 | |
CN102749067A (zh) | 自适应多级计算的精确反演海面高程遥测装置 | |
CN106443594A (zh) | 一种基于稀疏约束的雷达天线阵列稳健波束形成方法 | |
JP6362262B2 (ja) | 角度推定装置および角度推定方法 | |
CN103902830A (zh) | 一种圆环形阵列稳健旁瓣控制超指向性波束形成方法 | |
RU2515179C1 (ru) | Способ определения направления на гидроакустический маяк-ответчик в условиях многолучевого распространения навигационного сигнала | |
CN109405954A (zh) | 一种基于矢量水听器阵列的uuv移动水声通信技术 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |