基于传播算子的卫星导航系统干扰源直接定位方法.pdf
《基于传播算子的卫星导航系统干扰源直接定位方法.pdf》由会员分享,可在线阅读,更多相关《基于传播算子的卫星导航系统干扰源直接定位方法.pdf(7页珍藏版)》请在咨信网上搜索。
1、第 21 卷 第 8 期2023 年 8 月Vol.21,No.8Aug.,2023太赫兹科学与电子信息学报Journal of Terahertz Science and Electronic Information Technology基于传播算子的卫星导航系统干扰源直接定位方法唐元春1,陈端云2,夏炳森1(1.国网福建省电力有限公司经济技术研究院,福建 福州 350000;2.国网福建省电力有限公司,福建 福州 350001)摘要:针对现有单星定位技术精确度不足问题,提出一种基于传播算子的卫星导航系统干扰源直接定位方法,解决现有技术无法快速有效地定位中高轨卫星干扰源的问题。对卫星接收信号
2、建模,计算每个时刻的接收信号协方差矩阵,通过对协方差矩阵分块,计算得到传播算子,进而得到导向矢量正交投影子空间的估计;然后建立损失函数,通过搜索得到干扰源位置估计结果。理论分析和仿真结果证明,该算法在复杂度大大降低的同时,定位精确度与传统基于多重信号分类(MUSIC)的直接定位算法相当,远优于现有的测向定位技术,具有重要的工程应用价值。关键词:导航卫星;传播算子;干扰源;直接定位中图分类号:TN96 文献标志码:Adoi:10.11805/TKYDA2021118Direct positioning method of interference source of satellite navi
3、gation system Direct positioning method of interference source of satellite navigation system based on propagation operatorbased on propagation operatorTANG Yuanchun1,CHEN Duanyun2,XIA Bingsen1(1.State Grid Fujian Economic Research Institue,Fuzhou Fujian 350000,China;2.State Grid Fujian Electric Pow
4、er Co.,Ltd,Fuzhou Fujian 350001,China)AbstractAbstract:To address the problem of insufficient precision of existing single-satellite positioning techniques,a direct positioning method of interference source for satellite navigation system based on propagation operator is proposed.Firstly,the receive
5、d signal of the satellite is modeled and the covariance matrix at each moment is calculated.After blocking the covariance matrix,the propagation operator is obtained to estimate the orthogonal projection subspace of the steering vector.Then the loss function can be established and the position of th
6、e interference source can be estimated by searching the loss function.Theoretical analysis and simulation results prove that the localization accuracy of the proposed method is comparable to that of the traditional Multiple Signal Classification(MUSIC)-based direct positioning algorithm,while the ca
7、lculation complexity is greatly reduced.Meanwhile,the localization accuracy is much higher than that of the existing bearing-only localization techniques,which is valuable for practical engineering applications.KeywordsKeywords:navigation satellite;propagation operator;interference source;direct pos
8、itioning卫星导航定位系统具有覆盖范围大,通信距离远,不受地理环境影响等优势。我国自主研制的北斗导航系统发展迅速,除了导航定位和授时外,还兼具短报文通信功能,已在交通、测绘、电信、气象等方面发挥着至关重要的作用,是我国综合国力的象征之一。由于卫星导航系统工作在一个开放的环境中,其导航定位等功能极易受到各种人为或无意的干扰。随着云计算、大数据、物联网、移动互联网、人工智能等新兴技术带动的新一轮科技革命和产业变革,无线电网络日益增多,卫星所处的电磁频谱环境日益恶化,复杂干扰信号的存在导致卫星导航质量下降,甚至会给国家带来巨大的经济损失。因此,快速有效地通过受扰卫星对干扰源进行定位是当务之急。
9、根据定位所需的卫星数量分类,卫星对地面干扰源的定位可分为单星定位体制和多星定位体制1-5。传统的双星和三星定位技术虽然定位精确度高,但需要多颗卫星协同工作,邻星的选择是一个难题。单星定位只需一文章编号:2095-4980(2023)08-0985-07收稿日期:2021-03-23;修回日期:2023-06-27太赫兹科学与电子信息学报第 21 卷颗卫星就能完成定位,不仅节省了有限的轨道资源,同时避免了邻星选择的困难。本文研究的主要对象为单星定位体制。现有的单星定位体制大都是测向交叉定位,整个定位过程分成参数估计和干扰源定位 2 个阶段。通过参数估计阶段对干扰源方向的估计,结合卫星与干扰源的空
10、间几何关系,进行干扰源位置解算6-7。但这种传统两步定位方法并没有充分考虑定位过程的整体性,每次参数估计步骤都是一个独立的过程,干扰源定位精确度完全依赖于参数估计的准确性,损失了整个定位过程中干扰源位置不变这一有效信息。由于低轨卫星的轨道高度相对较低,卫星可以截获功率相对较强的干扰源信号,因而可以获得较为准确的干扰源方位估计,卫星对干扰源的定位精确度容易达到需求。但对于卫星导航系统的中高轨卫星,由于卫星距离干扰源较远,一方面使测向带来的位置误差巨大;另一方面接收信号功率衰减严重,卫星接收到的干扰源信号较弱8-9,导致方位信息测量精确度不够,中间参数估计错误还会直接导致定位失败,不利于实际应用,
11、实现的技术难度非常大10。传统多重信号分类(MUSIC)等直接定位方法不需要估计定位参数过程,直接对原始信号进行采样处理,利用信号中的位置信息构建目标函数,通过穷尽搜索等算法实现定位11。本文针对现有的中高轨单星测向定位体制存在的问题,在干扰源位置不变的假设下,提出一种基于传播算子的卫星导航系统直接定位方法,系统性地考虑整个定位过程,利用传播算子及正交投影子空间的性质构造定位损失函数,通过直接搜索损失函数得到干扰源位置估计。算法避免了中间参数的估计步骤,且由于直接从原始数据层抽取位置信息,避免了信息的二次损耗,相比于传统测向定位算法,定位精确度得到了大幅提升。本文所提方法避免了传统基于 MUS
12、IC 等直接定位算法的特征值分解步骤,大大降低了计算复杂度。文中也给出了算法的计算复杂度分析,仿真结果表明该算法具有较高的估计精确度,计算效率高,是一种适用于中高轨卫星的快速有效的干扰源定位方法,同时也可拓展到无线电监测、雷达、勘探等领域。1系统模型系统模型如图 1 所示,建立地固坐标系OeX4Y4Z4,以地心为原点,OeX4由地心向外指向格林尼治子午线,OeZ4垂直于地球赤道平面与地球自转轴重合指向北极,OeY4与OeX4和OeZ4构成右手直角坐标系。在地理坐标系中,采用经度、纬度和高程来描述空间的位置,表示为H。地球半径为R,根据地理坐标系与地固坐标系之间的转换关系,可以计算出地固坐标系下
13、的空间位置为:x=(R+H)coscosy=(R+H)sincosz=(R+H)sin(1)干扰源地理位置矢量为p k=kkHkT,根据卫星星历得到卫星地固位置矢量为ul=xlylzlT,卫星接收机上集成了一阵元数为M的 L 阵,该 L 阵由 2 个阵元数相同的均匀线阵组成,所有阵元均被视为无差别的点,且对于各个方向的接收增益相等。传统的测向定位算法通常需要先从卫星截获的信号中估计出干扰源信号的方位角和俯仰角,利用卫星与干扰源之间的相对运动,通过空间几何关系才能实现定位。由于利用了方向信息,传统测向定位对卫星的姿态测量和控制精确度要求较高。本文采用的定位原理不需要通过中间参数的估计,通过构造损
14、失函数就能直接从原始数据层实现干扰源位置估计,大大避免了中间参数估计带来的误差。2定位算法2.1 算法描述为了简化模型,假设卫星接收的干扰源信号为窄带信号,且一段时间内对卫星存在连续干扰,则卫星在第l个位置的快拍t的接收信号为:xl(t)=Alsl(t)+nl(t)(2)式中:t=12J,J为快拍数;sl(t)=bl1s1(t)bl2s2(t)blKsK(t)T为阵列天线的参考阵元接收到的信号矢量,blk OeX4Z4Y4 u1ulpkFig.1 System model图1 系统模型986第 8 期唐元春等:基于传播算子的卫星导航系统干扰源直接定位方法为卫星接收信号的传输衰减系数;sk(t)
15、为干扰源k的发射信号,k=12K,K为干扰源个数;nl(t)(l=12L)为相互独立的零均值加性高斯白噪声,方差为2,且信号与噪声之间互不相关;Al=al(p1)al(p2)al(pK)为阵列流型矩阵,al(pk)=ejkTl(pk)d1ejkTl(pk)d2ejkTl(pk)dMT,其中,pk为干扰源地固位置矢量,dm表示阵列的第m个阵元相对于参考阵元的位置矢量,kl(pk)为卫星在位置ul处的波数向量,表示为:kl(pk)=2pk-ulpk-ul(3)式中:为波长;表示向量的 2 范数。由于实际信号采样长度有限,接收信号协方差矩阵由采样协方差矩阵代替,可计算为:Rl=1Jt=1Txl(t)
16、xHl(t)=AlRssAHl+2IM(4)式中:Rss=Esl(t)sHl(t)为信源信号协方差矩阵;IM为M M维单位矩阵。考虑到阵列流型矩阵Al为M K维列满秩矩阵,将Al分块为:Al=A1;A2(5)式中:A1 CK K为满秩矩阵;A2 C(M-K)K。显然,矩阵A2可以通过矩阵A1的线性变换得到:PA1=A2(6)式中 P 为(M-K)K维传播算子。将式(5)代入式(4),可得Rl=A1RssAHlA2RssAHl+2IM(7)考虑无噪声环境下的协方差矩阵Rl=Rl-2IM,令Rl=Ga;Gb,其中Ga为矩阵Rl的前K行,Gb为矩阵Rl的后(M-K)行,利用式(6)的关系,则有PGa
17、=Gb(8)由于实际环境中不可避免地存在噪声,因此,根据最小二乘理论,传播算子矩阵的估计为:P=minGHaPH-GHb2(9)式中:Ga和Gb分别由矩阵Rl的前K行和后(M-K)行构成,由式(9)可解得传播算子矩阵的估计值为:P=(GaGHa)-1GaGHbH(10)构造矩阵Pc=IK;P(11)式中IK为K K维单位矩阵,显然有PcA1=Al(12)根据式(12)可得矩阵Pc和矩阵Al由相同的列空间组成,即spanPc=spanAl。因此,可得导向矢量的正交投影子空间的估计为:Ql=Pc(PHcPc)-1PHc(13)利用导向矢量正交投影子空间的性质,可以构造出损失函数:f-1(p)=l=
18、1LaHl(p)(IM-Ql)al(p)(14)将搜索区域内的位置坐标 p 分别代入式(14),对上述损失函数进行搜索,找到式(14)的最大值,该最大值对应的位置坐标即为干扰源的位置估计结果。987太赫兹科学与电子信息学报第 21 卷2.2 搜索区域的选取由于地表范围很大,如果对整个地球表面进行一次全搜索,则定位效率很低。因此考虑利用卫星的对地覆盖区域来缩小定位搜索区域,从而加快卫星对干扰源的定位速度,这也是成功定位干扰源的保障。对于位于地面的静止干扰源(即干扰源在地理坐标系下的高度为 0),其可能存在的地固空间是无限的,而全球经纬度范围是有限的,因而考虑将干扰源可能存在的区域根据经纬度范围网
19、格化,如图 2所示。根据图 2 所示的卫星覆盖区域示意图,卫星位置点为S,干扰源可能存在的位置点为E,显然卫星覆盖区域内的位置点到卫星与地心连线的投影大于卫星与地球的切点到卫星与地心连线的投影,即向量OeE到向量OeS的投影应大于dmin。从图 2 所示的几何关系中,可以计算出dmin=R2(R+H)(15)同时通过前面介绍的系统模型可知地心Oe为坐标系原点,显然可得OeE=piOeS=ul(16)式中pi为干扰源可能存在的位置矢量,由离散化的经纬度网格点p i转换的地固坐标系下的坐标,则卫星在时刻l的对地覆盖区域可表示为:areal=pi|pieiR2(R+H)(17)式中:piei表示向量
20、pi与向量ei的点积,ei=ulul为卫星与地心连线方向上的单位向量。取L个时刻卫星覆盖区域的交集,可得卫星覆盖区域的重叠区域为:Area=pi|piarea1area2areaL2.3 算法步骤综上所述,本文所提算法的流程可归纳如下:步骤 1:选取卫星在L个不同时刻接收的干扰源信号,分别计算每个时刻的接收信号协方差矩阵;步骤 2:对协方差矩阵分块,计算得到传播算子,进而得到导向矢量正交投影子空间的估计;步骤 3:根据卫星的轨道信息计算出被干扰时卫星覆盖区域的重叠区域;步骤 4:利用导向矢量正交投影子空间的性质建立损失函数,将L个时刻的导向矢量正交投影子空间融合,通过搜索重叠区域得到干扰源位置
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 传播 算子 卫星 导航系统 干扰 直接 定位 方法
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。