基于和差波束测角的角度跟踪仿真研究
2023-01-09
来源:易榕旅网
第1期 2017年2月 雷达科学与技术 Radar Scfence and Technology Vo1.15 No.1 February 2017 基于和差波束测角的角度跟踪仿真研究 马越。李朝海 (电子科技大学电子工程学院,四川成都611731) 摘 要:飞机、导弹等高速运动载体与卫星系统实时通信过程中,载体的姿态变化,会影响通信质量。 为使载体上的相控阵天线始终对准卫星方向,保持载体与卫星正常通信,提出了一种角度跟踪系统设计方 案,对系统关键部分的功能进行说明;针对均匀对称的矩形阵列,在单脉冲测角理论和卡尔曼滤波算法的基 础上进行研究和改进,提出了基于和差波束测角和自适应a 滤波的角度跟踪算法,实现了对来波的角度测 量与跟踪。通过仿真,对系统测角精度与跟踪精度进行说明,结合算法原理分析影响系统跟踪精度的因素。 关键词:相控阵天线;和差波束;a—J9滤波;角度跟踪 中图分类号:TN953 .5;TN958.92 文献标志码:A 文章编号:1672—2337(2017)O1—0068—05 Simulation Research of Angle Tracking Based on Angle Measurement of Sum and Difference Beams MAYLie,LIChaohal (School of Electronic Engineering,University of Electronic Science and Technology of China,Chengdu 611731,China) Abstract:In the process of real-time communications between high-speed motion platforms such as air— craft,missiles and a repeater satellite,the postures of platforms change all the time,which affects the quali— ty of communication seriously.For keeping communications with the satellite,the phased array antenna on them must aim at the direction of satellite all the time.The design of angle tracking system is given,and the key functions of the system are introduced.For the rectangular array with symmetrical and uniform charac— teristies,by studying and improving the traditional theory of monopulse angle measurement and Kalman fil— ter,an angle tracking algorithm based on sum—difference angle measurement and adaptive a filtering is pro— posed.Through simulation,the angle measurement and tracking precision is given,and the factors affecting the system tracking precision is analysed. Key words:phased array antenna;sum and difference beams;口印filtering;angle tracking 见的角度估计方法主要有直接比相法、子阵相关 O 引言 运动载体卫星通信系统,即“动中通”,具有机 动灵活、覆盖范围广、抗干扰能力突出等特点,在 法 。 以及空间谱估计方法等。本文研究的和差波 束测角方法,测角精度更高,并通过匹配滤波提高 系统抗干扰能力。角度跟踪的经典方法是卡尔曼 滤波 以及扩展卡尔曼滤波等。但其涉及增益矩 阵的计算,运算量大,工程实现比较复杂。a 滤 波l5 能够克服这一问题,工程实现简单。本文在a— 军事、民用领域上扮演着不可替代的角色。其系统 设计的关键技术在于稳定准确的卫星跟踪能力。 高精度的跟踪系统,能有效地隔离天线平台的运 动,自身姿态的旋转以及平台周围位置环境对卫 滤波算法的基础上,提出自适应a 滤波算法,根 据方差的变化实时调整滤波器增益,提高算法的 稳健性。 星通信的影响,使天线始终对准卫星,保证平台在 高速运动中能够与卫星保持连续正常通信口]。 卫星跟踪系统按功能主要分为天线波束控 制、资源调度、信号搜索捕获、角度估计与跟踪等。 角度估计与角度跟踪是整个系统的关键部分。常 l系统总体设计方案 运动载体角度跟踪系统由波控单元、收发阵 收稿日期:2016—07—14;修回日期:2016—08—23 2017年第1期 马越:基于和差波束测角的角度跟踪仿真研究 69 列、射频通路以及基带信号处理单元等部分组成。 将接收阵列天线划分成4个子阵,每个子阵包含 8×8个均匀分布的间距为半波长的阵元。各子阵 输出的微弱信号经射频通路作放大以及模拟下变 频处理,得到4路中频信号,通过SMA接口送入 数字板。数字板基带信号处理单元通过4路ADC 对各子阵输出信号进行采样,采样频率为 100 MHz,得到基带数字信号。信号经过FIFO,数 字下变频(DDC),通道幅相校正等预处理模块后, 结合载体平台波控单元提供的系统初始引导信 图2搜索/捕获实现逻辑框图 息,完成目标搜索/捕获,利用和差波束测角对卫 系统波位顺序设定了两种状态_6],如图3所 星进行角度估计和跟踪滤波,把角度预测值传给 示。系统初始化获取引导信息,将引导信息中的指 运动载体波控单元,实时调整天线波束指向,从而 向作为初始波束指向的中心波位,按照状态l的顺 保证阵列天线一直对准卫星来波方向。系统整体 序调整波束指向进行搜索。若遍历状态l中5个 架构如图1所示。 波位均未搜索/捕获信号,说明此时系统初始引导 信息精度较差,则将状态1中的中心波位(波位1) -予接N收cl子阵D1阵列—}= 波控单兀 按照状态2的顺序进行调整,在状态2中9个中心 牺 — 一 波位依次按照状态1方式形成5个搜索波位,进行 曼 搜索/捕获。若遍历状态2波位仍未捕获信号,则 系统重新初始化,更新引导信息。 在角度跟踪系统工作的开始阶段,载体波控 单元提供的初始引导信息跟实际卫星方向存在一 状态1波位排列方式 状态2波位排列方式 定偏差,若直接进入角度跟踪状态,系统角度测量 图3波位编排两种状态 和跟踪精度也会受到极大的影响。因此,需要先对 来波方向进行搜索/捕获,成功捕获信号后转入跟 系统跟踪状态主要包括和差波束角度测量和 踪状态。 自适应a 滤波两部分。跟踪状态硬件实现逻辑 搜索/捕获实现逻辑框图如图2所示。载体波 框图如图4所示。 控单元根据引导信息,给出初始天线波束指向。并 根据系统搜索捕获状态下给出的波位调整信息, 实时地对波束指向作出调整。对每一个波位,系 统对四子阵输出信号A/D采样后,经预处理模块 后对子阵间信号作互相关运算,若超过检测门限, 则对其和波束进行解扩、解调等,确认是期望的协 作通信信号,则捕获成功,系统转入角度跟踪状 态;否则,捕获失败,反馈信息至波控单元,按照设 图4跟踪状态实现逻辑框图 定的波位编排方案调整波束指向,继续对来波方 各子阵接收信号经数字板上A/D采样、数字 向进行搜索/捕获。 下变频、低通滤波与抽取等,得到基带数字信号 。 70 雷达科学与技术 第15卷第1期 经过和差波束形成与匹配滤波,在提高信噪比的 参考子阵的相移分别为 n sin ≯2一-__si27tDy 同时,保留子阵信号间的相位差信息。完成和差波 束法角度测量,并将测量结果( , )实时地送入 DSP中进行跟踪滤波,得到下一时刻来波方向的 预测值( 。,0。),反馈到波控单元,用于调整阵列天 (3) 声。一 一声z+声。 in (4) (5) 线波束实时地对准卫星。 式中,D =Md /z,D,=Nd /z分别为 方向和 Y方向的子阵间距,M和N分别为z方向和Y方向 的阵元个数。 同理,记波束指向为( 。,0。),则第2,3,4子阵 对应的相移量分别为 2nDy 2。一_2 和差波束测角算法 针对上节所述的系统角度跟踪状态中角度测 量模块,本文提出一种基于和差波束测角的角度 测量方法,主要包括和差波束形成、相位差计算和 角度测量三部分。 2.1 阵列天线模型 sin 。sin 。 sin 。c。 。 (6) (7) (8) 2rtD ̄声3。::=_4o一 20+ 3o 因此,移相之后,各个子阵信号之间的相位差 分别为 △庐z一声z十≯ o一 均匀矩形面阵的几何结构以及阵列天线角度 坐标定义如图5所示。相控阵雷达接收阵列位于 xoy平面上,设阵元数为M×N,平行于z轴和Y 轴方向均匀分布的阵元且间距分别为d 和d 。记 [sin0 sin ̄--sin0 0 Sin 。](9) △ 3一声3+ 30一 入射信号在xoy平面上的投影与z轴夹角为方位 角 ,取值范围为一9O。~90。;记入射信号与z轴夹 角为俯仰角 ,取值范围为0。~9O。。 ————一1[一 sin0 c。 —icos — sin0 n0 coscocoscp p] ol(1)0) △ 一△ 2+△ 3 (11) 设子阵1的接收信号为 向 Y1( )一 1(咒)+"O1( ) (12) 式中, ( )为子阵1输出噪声。则第2,3,4子阵 接收信号为 fy 2(咒)= l( )eJ ̄ ̄2+7.3 2( ) Y 3( )= 1( )eJa ̄3+ 3( ) (13) 【 4( )一 1( )ej△ 4+ 4( ) 式中, ( ), 。( ), ( )分别为第2,3,4子阵 输出噪声。 2.2相位差求解 图5 均匀矩形面阵的几何结构 对各子阵输出信号进行和差波束形成得到: 信号的方向向量为 r:一[sin0 coscp sin0 sincp cos0] 分别是 P。=Eo,D ,O] P。一[D ,0,O] P 一ED ,D ,O] 和波束输出信号为 (1) y ( )一Ey ( )+ ( )]+[ 。( )+ ( )]+ ( )= 】以第1子阵为参考子阵,第2,3,4子阵的坐标 ( )(1+ej )(1+ej 3)+ :( ) (14) 俯仰差波束输出信号为 E△( )一Ey (咒)+ 。( )]一[ 2( )+ ( )]+ E△( )一 ( )(1+el 。)(1一 z)+ E△( ) (15) (2) 方位差波束输出信号为 Y (咒)一Ey1( )+ ( )]一[ 。( ) ( )]+ A△( )一 (16) ( )(1+ 郇z)(1一e) s)+ A ( ) 若人射信号先到达第2,3,4子阵,根据相控阵 空域滤波相关理论,可以得到第2,3,4子阵相对于 2017年第1期 式中, :( ), ( ), 马越:基于和差波束测角的角度跟踪仿真研究 ( )分别为和波束、俯仰 71 差波束、方位差波束的输出噪声。 通过匹配滤波,提高和差波束信号信噪比的 同时不改变原来的各子阵信号之间的相位关系, 因此有效地提高角度测量的精度。匹配滤波后和 3 跟踪滤波 系统的跟踪过程是通过将和差波束测角得到 的目标角度值( ,0)送给跟踪滤波器,通过跟踪滤 波束通道输出信号为 s : 1∑Y:(7z) ( )一 俯仰差波束通道输出信号为 一 ∑Y E△(n)z (72)一 EE△(1十e s)(1一eJa ̄2)+ (18) 方位差波束通道输出信号为 SAA一 ∑ ( )z (n)一 式中:Ⅳ为信号采样点快拍数;E ,E ,E 分别 为和波束信号、俯仰差波束信号和方位差波束信 号的能量; ,73 , 分别为和波束信号、俯仰 差波束信号和方位差波束信号的输出噪声幅度。 一 S=(,2) 1 + ej3 , s=( ) === 1+e 2 ㈣, 化简,提取出相位差△ ,△ 。: 2一谍1---'A△ 舔sz(n)(21) △≯。一 A (22 LX ̄。= A … 理论上,累积快拍数N足够大,角度测量的精 度比较理想。 2.3 角度计算 将2.1节的式(9)和式(10)进行化简可得 sinO sine:sinO 0sine o一△ 2 /2=D (23) sinO cos 一sinO ocos 0一△ 3 /27rD (24) 将2.2节得到的相位差△ 。,△ 。代入式(23) 和式(24)中,求解出角度值( ,0)如下:  ̄=ararctan(an sl —:-— ————_n(]0 cos6 _—————二——— =—÷n一A击 /Z7c /A ) (25) 0一±arcsin((sinO o sine o一△ 2A/2nD ) + (sin0 ocos ̄0~△ 3 /27cD )。)… (26) 波得到下一时刻来波方向的预测值( 。,0。),作为 下一时刻角度测量的输入参数,完成一次闭环的 自适应a一口跟踪滤波。 假设观测信息为目标角度测量值( ,0),通常 用方向余弦作为状态向量。目标的方位角和俯仰 角信息与方向余弦的对应关系如下: Jc]f-sinOc唧] . LY I2fJ l—j s inc0sosi0n9 J l (27) 式中,z ,Y 和 分别为目标在惯性坐标系下3 个方向的方向余弦。本文以目标在3个方向的 方向余弦及其变化的速度为状态向量 [z x Y 乏 ] ,其中的z ,Y 和z 均可 由( ,0)求得。主 , 和乏 可由下式求得: f. z ( )~ ( 一1) fX c一——— -—一 Yc一—————_ ———一一 (28)L 1. 2 ( )一 ( 一1) l 一——— -一 式中,T为跟踪滤波采样周期。 以一个通道的方向余弦估计为例进行说明, 其余两通道类似。 假设目标运动状态方程为 X( )一F( , 一1) ( 一1)+ r( , 一1)W(7z一1) (29) 式中,X( )一[z ( )主 (n)] 为目标z方向 的方向余弦、方向余弦变化速度构成的状态向 量,F( ,挖一1): T1]为状态转移矩阵, r( ,,z一1)一f-T2T/2]为状态噪声输入矩阵, ( ) 为该模型的过程噪声,服从零均值、方差为Q(n)的 高斯分布。 观测方程为 z( ):C( ) ( )+V( ) (3O) 式中,c( )一[1 O]为观测矩阵,v(n)为测量噪 声且服从零均值、方差为R(n)的高斯分布,它与 过程噪声w( )相互独立。 基于以上运动方程和观测方程,一次自适应 72 雷达科学与技术 第l5卷第1期 a一口滤波的过程如下: 步骤1:状态一步预测 ( l -1)一F( , 一1) ( 一1 l --1)(31) 式中, ( l n一1)表示用观测量 (1), (2),…, ( 一1)对信号 ( )进行的最小均方误差估计, 下同。 步骤2:由观测信号计算新息过程 ( )一z( )一三( I 一1)=== z( )一C( ) ( I 一1) (32) 步骤3:计算卡尔曼增益K(n),具体计算过程 由下文详细介绍。 步骤4:状态估计 ( 『,z)一 ( l 一1)+K(n)口( ) (33) 步骤5:重复步骤1~4,进行递推滤波计算。 步骤3中,K )一l T}由其中各个通道的 增益根据相应的机动指数计算获得,定义目标机 动指数为 (34) 式中, 和 :分别为机动加速度和测量误差的 方差。 a( )和 ( )与r的关系为 a(n)-- c√H )(3 ( )一2(2一 (n))一4 ̄/ 二 根据抽样统计知识以及目标机动的判定时间 区间,可得到机动加速度的方差 2 的计算公式如 下,其中V(n)是新息值,N可取3到5: 2( )一 [∑ 。( — )] (36) 在建立起目标运动轨迹模型的基础上,对目 标信号方向进行测量与跟踪,设子阵接收信号的 信噪比为一8 dB,和差波束匹配滤波模块的累积快 拍数为8 192,系统对其进行角度测量与跟踪的结 果如图6和图7所示。 图中观测误差表示角度测量结果与真实轨迹 的来波方向的偏差,跟踪误差表示系统跟踪滤波 过程得到的下一时刻角度预测值与真实值之间的 偏差。仿真结果表明,用和波束、俯仰差波束、方位 差波束信号相比的方法进行角度测量,结合自适 应a一 滤波进行角度跟踪,能有效地提高系统角度 测量与跟踪精度。 自 1无 一 三一 、 椭 援 椒 - x:0.422 二 0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 观测时间/s 图6方位角测量与跟踪误差 蕞耋I I 一。)、猁 晕 们 % ¨ ∞ 0 X 0102 Y:02637 c● i'X;0 178 ,揍 0 02 04 0.6 0 8 1.0 1.2 I_4 I 6 1.8 2.0 观测时间/s 图7俯仰角测量与跟踪误差 4结束语 本文针对高速运动载体卫星实时通信的背景 下,提出角度测量与跟踪系统的设计方案。结合阵 列天线子阵划分的特点,在传统单脉冲测角理论 的基础上,提出一种新的和差波束测角方法。并针 对载体的运动特性,研究和发展了自适应 滤波 算法,保证了系统角度测量跟踪的实时性和稳健 性。综上,本文搭建的角度测量与跟踪的闭环系 统,是实现高速运动载体卫星跟踪技术的关键,对 工程实践具有一定的指导意义。 参考文献: E1]赵来定.一种车载卫星通信跟踪技术EJ].南京邮电 大学学报(自然科学版),2014,34(5):80—86. [2]姜义成,喻春曦.一种改善单脉冲雷达测角精度的新 方法[J].雷达科学与技术,2009,7(5):380—383. (下转第80页) 80 雷达科学与技术 表3雷达抗干扰方法效能评估 第15卷第1期 型分析的雷达在复杂电磁环境下的抗干扰性能定 量评估方法。该方法将试验数据映射到高维量测 趋势[J].火力与指挥控制,2015,40(4):1—5. [7]聂红霞.雷达抗干扰效能评估方法探讨[J].舰船电 子工程,2Ol3,33(4):75-77. 特征空间中,通过统计理论分析量测值与真实目 标特性的差异,提出了MEIF、MSIF和JIF三种指 [8]朱宁,赵肖肖,唐庆华.多元正态性检验三种方法的比 较及sAS程序设计[J].苏州大学学报(自然科学版), 2012,28(3):20—25. 标,既可以描述干扰对于雷达的影响,也可以反映 雷达的抗干扰能力,并且对于单一干扰、复合干扰 不敏感,对复杂电磁干扰环境具有良好的适用性。 通过仿真,可以看出该方法理论可行,对雷达的抗 干扰能力建设具有指导性意义。 参考文献: [1]李浩,邱超凡,赵小亮.雷达抗干扰效能的多层次模糊 评估方法[J].雷达科学与技术,2012,10(2):143— 149,155. [9]梁小筠.我国正在制定“正态性检验”的新标准[J].应 用概率统计,2002,18(3):269—276. r1O]BOX G E P,COX D R.An Analysis of Transforma一 tions[J].Journal of the Royal Statistical Society, 1964,26(2):211-252. [11] 李长江,邓文平,曹元元,等.基于Box-Cox变换与 Johnson变换非正态过程能力分析[J].齐齐哈尔大 学学报(自然科学版),2015,3i(1):66—70. 作者简介 陈文东 男,1987年出生于山西永 济,在读博士研究生,主要研究方向为 [2]徐沙,张洁.一种基于层次分析法的雷达抗干扰效能 评估方法[J].航天电子对抗,2016,32(3):13—15. [3]刘松涛,姜宁,刘振兴.舰载电子对抗干扰效果在线评 估方法[J].电子信息对抗技术,2016,31(4):63—67. [43周振宇,贺志毅,刘建欣.末制导雷达抗干扰性能评估系 统工程实现EJ].现代防御技术,2011,39(2):138—141. [5]刘江波,舒敬环.复杂电磁环境下的雷达抗干扰效能评 雷达总体设计与雷达电子反对抗 技术。E-mail;chen_eecs@sina.corn 估指标体系构建[J].信息化研究,2015,41(5):1-5. [6]李亚南,韩壮志.雷达抗干扰性能评价指标体系及发展 汤 斌男,1977年出生,高级工程师,主要研究方向 为雷达总体设计与雷达电子反对抗技术。 (上接第72页) 作者简介 马越男,1991年生,山东烟台人, [3]杨雪亚,刘张林.子阵级和差波束形成及测角方法研究 口].中国电子科学研究院学报,2015,10(1):82—86. [4]姜雯献,李朝海.相控阵跟踪系统测角与角度跟踪算 法[J].雷达科学与技术,2015,13(2):190—194. 硕士研究生,主要研究方向为雷达信 号处理。 E-mail:woshimyue@163.corn [5]张海成,杨江平,王晗中.大型跟踪测量雷达的卫星标定 方法研究口].雷达科学与技术,2014,12(5):470—472. [6]LI Yang,LYu Huihui,PENG Sun,et a1.Study on Search Performance of Subarray Multi—Channel Phased 李朝海男,1972年生,四川I南充人, Array Radar Based on Multiple Received Beams[C] } IET Internationa1 Radar Conference 2013,Xi’an:IET, 2013:卜6. 电子科技大学电子工程学院教授级高 级工程师,主要研究方向为高速实时 信号处理、雷达系统及信号处理。 [7]吴越,严济鸿,何子述.基于FPGA的多通道高速数据 采集系统[J].雷达科学与技术,2012,10(6):671—676.