CN111239777A - 一种基于位置指纹的卫星集群分级定位方法 - Google Patents

一种基于位置指纹的卫星集群分级定位方法 Download PDF

Info

Publication number
CN111239777A
CN111239777A CN202010014588.1A CN202010014588A CN111239777A CN 111239777 A CN111239777 A CN 111239777A CN 202010014588 A CN202010014588 A CN 202010014588A CN 111239777 A CN111239777 A CN 111239777A
Authority
CN
China
Prior art keywords
satellite
star
coordinate system
time
positioning
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
Application number
CN202010014588.1A
Other languages
English (en)
Other versions
CN111239777B (zh
Inventor
马琳
黄鹏飞
谭学治
王孝
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harbin Institute of Technology
Original Assignee
Harbin Institute of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Harbin Institute of Technology filed Critical Harbin Institute of Technology
Priority to CN202010014588.1A priority Critical patent/CN111239777B/zh
Publication of CN111239777A publication Critical patent/CN111239777A/zh
Application granted granted Critical
Publication of CN111239777B publication Critical patent/CN111239777B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13Receivers
    • G01S19/24Acquisition or tracking or demodulation of signals transmitted by the system
    • G01S19/25Acquisition or tracking or demodulation of signals transmitted by the system involving aiding data received from a cooperating element, e.g. assisted GPS
    • G01S19/258Acquisition or tracking or demodulation of signals transmitted by the system involving aiding data received from a cooperating element, e.g. assisted GPS relating to the satellite constellation, e.g. almanac, ephemeris data, lists of satellites in view
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/16Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02DCLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00Reducing energy consumption in communication networks
    • Y02D30/70Reducing energy consumption in communication networks in wireless communication networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Computing Systems (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

本发明提出一种基于位置指纹的卫星集群分级定位方法,所述方法首先以较大的参考点间隔计算稀疏Radio Map,以使得其占用较小的存储空间,实现对成员星的粗定位;再以粗定位结果为中心,确定精定位区域,并对精定位区域进行精细的参考点划分,根据当前时刻计算精定位区域的Radio Map,对成员星进行精定位,以得到成员星最终的定位结果。该方法能够以增加在线计算时间来减小存储空间,并减小时间偏移对定位精度的影响。所述方法能够减小时间偏移对定位精度的影响同时极大地减小存储空间,能够更好地运用于实际工程中。

Description

一种基于位置指纹的卫星集群分级定位方法
技术领域
本发明属于信号处理技术领域,特别是涉及一种基于位置指纹的卫星集群分级定位方法。
背景技术
卫星集群是指由多个卫星组成的一个分布式卫星***。其中,选取一颗或若干颗卫星作为基准星,其余卫星为成员星。基准星沿着预设的轨道飞行,成员星跟随基准星进行伴飞。作为一种分布式卫星***,在卫星集群中,卫星之间保持着一个相对稳定的分布,卫星之间通过星间链路,实现信息交互,从而使得卫星之间相互协作,能够共同完成一个或多个任务。相比于大卫星的设计周期长、成本高等缺点,小卫星具有集成度高、设计周期短、成本低等优点,因此小卫星有着不可替代的优势。
传统的卫星集群自主定位主要针对两个卫星或小规模卫星集群。随着航天任务的任务需求的不断发展,卫星集群成员数量不断增多,体积进一步减小,功能趋向专一化,因此需要研究一种快速、低成本的,并具有集群自主性的定位导航技术。基于位置指纹的定位方法是一种较为成熟的定位方法,指纹定位***具有成本低、结构简单的特点,因此将指纹定位方法扩展至卫星集群定位中,能够有效减小研制成本,并对卫星集群的协同控制以及卫星技术的发展具有重要意义。
发明内容
本发明目的是为了解决现有技术中的问题,提出了一种基于位置指纹的卫星集群分级定位方法。
本发明是通过以下技术方案实现的,本发明提出一种基于位置指纹的卫星集群分级定位方法,所述方法分为三个阶段:构建坐标系阶段、离线阶段和在线阶段;在构建坐标系阶段建立LVLH坐标系和基准星本体坐标系;在离线阶段,由地面计算出稀疏Radio Map,并储存在成员星中;在线阶段分为三个子阶段:成员星粗定位阶段、指纹库在线计算阶段和成员星精定位阶段;在成员星粗定位阶段,成员星根据当前时刻和稀疏Radio Map选取最近时刻的RSS矩阵,并结合实时获取到的RSS向量进行粗定位;在指纹库在线计算阶段,成员星根据粗定位结果计算稠密Radio Map;在成员星精定位阶段,成员星根据实时获取到的RSS向量和计算出的稠密Radio Map进行精定位,得到成员星的最终估计位置。
进一步地,所述构建坐标系阶段具体为:
在卫星集群中基准星的数量为4,其余卫星为成员星,基准星成对进行绕飞,绕飞的中心称为“中心星”,中心星为一虚拟的点,实际上并不是卫星;成员星则在基准星附近进行伴飞;
LVLH坐标系原点OL位于卫星集群中的“中心星”,坐标轴之间满足右手定则;OLXL轴方向为中心星运动的切线方向,OLYL轴由中心星指向地球质心,OLZL轴垂直于XLYL平面;
基准星本体坐标系原点OB位于基准星的质心,坐标轴方向与基准星姿态有关;OBXB轴为天线方向图0度方向,OBYB轴为天线方向图90度方向,OBZB轴垂直于XBYB平面,OBXB轴和OBYB轴满足右手定则;并且需要对4个基准星都建立对应的基准星本体坐标系。
进一步地,所述离线阶段具体步骤为:
步骤1、在LVLH坐标系下确定边长为l0的粗定位区域,并将其均匀划分为n个边长为l1的网格,将每个边长为l1的网格的中心位置作为粗定位区域的参考点,得到参考点位置向量L:
Figure BDA0002358390680000021
其中,(xj,yj,zj),j=1,2,3…n为粗定位区域中第j个参考点在LVLH坐标系下的坐标;
步骤2、将基准星的运动周期T以间隔ΔT划分为m个时间段,每个时间段的中心时刻为:
T1,T2,…,Tm (2)
根据基准星运动方程获取Ti,i=1,2,…m时刻4个基准星位置向量:
Figure BDA0002358390680000022
步骤3、通过预先设置的参数获取第k个基准星,k=1,2,3,4,在Ti时刻的发射功率
Figure BDA0002358390680000031
第k个基准星在LVLH坐标系下的坐标为(xk(Ti),yk(Ti),zk(Ti)),第k个基准星本体坐标系相对于LVLH坐标系旋转角度为
Figure BDA0002358390680000032
Figure BDA0002358390680000033
θk(Ti),ψk(Ti)分别为三个坐标轴相应的旋转角度;
计算在Ti时刻粗定位区域内第j个参考点相对于第k个基准星的位置:
Figure BDA0002358390680000034
步骤4、计算在Ti时刻LVLH坐标系变换到第k个基准星本体坐标系的旋转矩阵:
Figure BDA0002358390680000035
其中,
Figure BDA0002358390680000036
Figure BDA0002358390680000037
Figure BDA0002358390680000038
步骤5、计算第j个参考点在Ti时刻第k个基准星本体坐标系下的直角坐标:
Figure BDA0002358390680000039
为了使用天线方向图计算传输损耗,需要将直角坐标变换为球坐标,从直角坐标到球坐标的变换公式为:
Figure BDA00023583906800000310
于是,第j个参考点在Ti时刻第k个基准星本体坐标系下的球坐标表示为:
Figure BDA0002358390680000041
其中,rk(Ti)是第j个参考点在Ti时刻与第k个基准星之间的距离;θG,k(Ti)为第j个参考点在Ti时刻相对于第k个基准星的天顶角;
Figure BDA0002358390680000042
为第j个参考点在Ti时刻相对于第k个基准星的方位角;
步骤6、根据球坐标计算在Ti时刻天线增益:
Figure BDA0002358390680000043
其中,Gmax为发射天线最大增益,
Figure BDA0002358390680000044
为天线方向性函数,fmax是天线方向性函数最大值;
计算基准星发射的信号到达成员星的传输损耗:
Figure BDA0002358390680000045
其中,f是天线发射频率,GR是成员星接收天线增益;
于是,Ti时刻在第j个参考点处接收到第k个基准星的接收信号强度RSS为:
Figure BDA0002358390680000046
步骤7、重复步骤1至步骤6可得到Ti时刻在所有参考点处接收到4个基准星的RSS,构成如下Ti时刻RSS矩阵:
Figure BDA0002358390680000047
其中,
Figure BDA0002358390680000048
为Ti时刻在第j个参考点处接收到4个基准星的RSS构成的向量;
步骤8、重复步骤1至步骤7可得到T1,T2,…,Tm时刻的RSS矩阵:
R1,R2,…,Rm (16)
于是得到稀疏Radio Map:
M=[L R1 R2 … Rm] (17)
稀疏Radio Map由地面提前计算好并存储至每一个成员星中,成员星进行定位时调用。
进一步地,所述成员星粗定位阶段具体为:
成员星获取当前时刻接收到的4个基准星的RSS:
Figure BDA0002358390680000051
根据卫星时钟获取当前时刻t,并在时刻T1,T2,…,Tm中选择与当前时刻t最接近的时刻T*,根据式(17)得到T*时刻的稀疏Radio Map:
Figure BDA0002358390680000052
将成员星在线获得的RSS与T*时刻的稀疏Radio Map中的RSS进行匹配,根据kNN算法获得成员星的粗定位结果为:
Figure BDA0002358390680000053
进一步地,所述指纹库在线计算阶段和成员星精定位阶段具体为:
根据成员星粗定位结果确定成员星精定位区域:以成员星粗定位结果为中心,确定边长为l2的正方体区域作为成员星精定位区域;
将精定位区域按照均匀划分为w个边长为l3的网格,将每个边长为l3的网格的中心位置作为精定位区域的参考点;
得到精定位区域参考点位置向量:
Figure BDA0002358390680000054
根据基准星运行参数得到当前时刻4个基准星的发射功率、在LVLH坐标系下的坐标以及基准星本体坐标系相对于LVLH坐标系旋转角度,从而计算精定位区域在该当前时刻的RSS矩阵:
Figure BDA0002358390680000061
因此,稠密Radio Map表示为:
Figure BDA0002358390680000062
将成员星在线获得的RSS向量与当前时刻的稠密Radio Map中的每个参考点的RSS进行匹配,根据最近邻算法获得成员星的精定位结果为:
Figure BDA0002358390680000063
本发明提出的基于位置指纹的卫星集群分级定位方法首先以较大的参考点间隔计算稀疏Radio Map,以使得其占用较小的存储空间,实现对成员星的粗定位;再以粗定位结果为中心,确定精定位区域,并对精定位区域进行精细的参考点划分,根据当前时刻计算精定位区域的Radio Map,对成员星进行精定位,以得到成员星最终的定位结果。该方法能够以增加在线计算时间来减小存储空间,并减小时间偏移对定位精度的影响。实验结果表明,当精定位区域边长分别为20m、30m、40m、50m时,定位误差小于5m的概率分别超过60%、75%、80%、85%;并且该方法能够减小时间偏移对定位精度的影响同时极大地减小存储空间,能够更好地运用于实际工程中。
附图说明
图1为卫星集群示意图;
图2为坐标系示意图;
图3位粗定位区域网格化和参考点选取示意图;
图4为基准星本体坐标系下直角坐标和球坐标示意图;
图5为精定位区域选取示意图;
图6为精定位区域网格化和参考点选取示意图;
图7为基于位置指纹的卫星集群分级定位方法流程图;
图8为时间偏移量示意图;
图9为定位误差CDF曲线示意图;
图10为时间偏移对定位误差的影响示意图。
具体实施方式
下面将结合本发明实施例中的附图对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
结合图7,本发明提出一种基于位置指纹的卫星集群分级定位方法,所述方法分为三个阶段:构建坐标系阶段、离线阶段和在线阶段;在构建坐标系阶段建立LVLH坐标系和基准星本体坐标系;在离线阶段,由地面计算出稀疏Radio Map,并储存在成员星中;在线阶段分为三个子阶段:成员星粗定位阶段、指纹库在线计算阶段和成员星精定位阶段;在成员星粗定位阶段,成员星根据当前时刻和稀疏Radio Map选取最近时刻的RSS矩阵,并结合实时获取到的RSS向量进行粗定位;在指纹库在线计算阶段,成员星根据粗定位结果计算稠密Radio Map;在成员星精定位阶段,成员星根据实时获取到的RSS向量和计算出的稠密RadioMap进行精定位,得到成员星的最终估计位置。
所述构建坐标系阶段具体为:
在卫星集群中,基准星负责整个卫星集群的协调和控制,成员星负责数据的采集和转发。其中,基准星体积、质量大,一般能够通过星历或者安装高精度定位设备确定准确的轨道位置;而成员星体积、质量小得多,且载荷有限,因此所能承载的定位设备有限。本发明所述方法在卫星集群中基准星的数量为4,其余卫星为成员星,基准星成对进行绕飞,绕飞的中心称为“中心星”,中心星为一虚拟的点,实际上并不是卫星;成员星则在基准星附近进行伴飞,如图1所示;
结合图2,LVLH坐标系原点OL位于卫星集群中的“中心星”,坐标轴之间满足右手定则;OLXL轴方向为中心星运动的切线方向,OLYL轴由中心星指向地球质心,OLZL轴垂直于XLYL平面;
基准星本体坐标系原点OB位于基准星的质心,坐标轴方向与基准星姿态有关;OBXB轴为天线方向图0度方向,OBYB轴为天线方向图90度方向,OBZB轴垂直于XBYB平面,OBXB轴和OBYB轴满足右手定则;并且需要对4个基准星都建立对应的基准星本体坐标系。卫星集群自主定位算法就是要确定成员星在卫星集群中的相对位置,因此本发明中基准星和成员星的位置均在LVLH坐标系下表示。
所述离线阶段具体步骤为:
步骤1、在LVLH坐标系下确定边长为l0的粗定位区域,并将其均匀划分为n个边长为l1的网格,如图3所示,将每个边长为l1的网格的中心位置作为粗定位区域的参考点,得到参考点位置向量L:
Figure BDA0002358390680000081
其中,(xj,yj,zj),j=1,2,3…n为粗定位区域中第j个参考点在LVLH坐标系下的坐标;
步骤2、将基准星的运动周期T以间隔ΔT划分为m个时间段,每个时间段的中心时刻为:
T1,T2,…,Tm (2)
根据基准星运动方程获取Ti,i=1,2,…m时刻4个基准星位置向量:
Figure BDA0002358390680000082
步骤3、通过预先设置的参数获取第k个基准星,k=1,2,3,4,在Ti时刻的发射功率
Figure BDA0002358390680000083
第k个基准星在LVLH坐标系下的坐标为(xk(Ti),yk(Ti),zk(Ti)),第k个基准星本体坐标系相对于LVLH坐标系旋转角度为
Figure BDA0002358390680000084
Figure BDA0002358390680000085
θk(Ti),ψk(Ti)分别为三个坐标轴相应的旋转角度;
计算在Ti时刻粗定位区域内第j个参考点相对于第k个基准星的位置:
Figure BDA0002358390680000086
步骤4、计算在Ti时刻LVLH坐标系变换到第k个基准星本体坐标系的旋转矩阵:
Figure BDA0002358390680000087
其中,
Figure BDA0002358390680000088
Figure BDA0002358390680000091
Figure BDA0002358390680000092
步骤5、计算第j个参考点在Ti时刻第k个基准星本体坐标系下的直角坐标:
Figure BDA0002358390680000093
为了使用天线方向图计算传输损耗,需要将直角坐标变换为球坐标,如图4所示,从直角坐标到球坐标的变换公式为:
Figure BDA0002358390680000094
于是,第j个参考点在Ti时刻第k个基准星本体坐标系下的球坐标表示为:
Figure BDA0002358390680000095
其中,rk(Ti)是第j个参考点在Ti时刻与第k个基准星之间的距离;θG,k(Ti)为第j个参考点在Ti时刻相对于第k个基准星的天顶角;
Figure BDA0002358390680000096
为第j个参考点在Ti时刻相对于第k个基准星的方位角;
步骤6、根据球坐标计算在Ti时刻天线增益:
Figure BDA0002358390680000097
其中,Gmax为发射天线最大增益,
Figure BDA0002358390680000098
为天线方向性函数,fmax是天线方向性函数最大值;
计算基准星发射的信号到达成员星的传输损耗:
Figure BDA0002358390680000106
其中,f是天线发射频率,GR是成员星接收天线增益;
于是,Ti时刻在第j个参考点处接收到第k个基准星的接收信号强度RSS为:
Figure BDA0002358390680000101
步骤7、重复步骤1至步骤6可得到Ti时刻在所有参考点处接收到4个基准星的RSS,构成如下Ti时刻RSS矩阵:
Figure BDA0002358390680000102
其中,
Figure BDA0002358390680000103
为Ti时刻在第j个参考点处接收到4个基准星的RSS构成的向量;
步骤8、重复步骤1至步骤7可得到T1,T2,…,Tm时刻的RSS矩阵:
R1,R2,…,Rm (16)
于是得到稀疏Radio Map:
M=[L R1 R2 … Rm] (17)
稀疏Radio Map由地面提前计算好并存储至每一个成员星中,成员星进行定位时调用。
所述成员星粗定位阶段具体为:
成员星获取当前时刻接收到的4个基准星的RSS:
Figure BDA0002358390680000104
根据卫星时钟获取当前时刻t,并在时刻T1,T2,…,Tm中选择与当前时刻t最接近的时刻T*,根据式(17)得到T*时刻的稀疏Radio Map:
Figure BDA0002358390680000105
将成员星在线获得的RSS与T*时刻的稀疏Radio Map中的RSS进行匹配,根据kNN算法获得成员星的粗定位结果为:
Figure BDA0002358390680000111
所述指纹库在线计算阶段和成员星精定位阶段具体为:
根据成员星粗定位结果确定成员星精定位区域:以成员星粗定位结果为中心,确定边长为l2的正方体区域作为成员星精定位区域;如图5所示;
将精定位区域按照均匀划分为w个边长为l3的网格,将每个边长为l3的网格的中心位置作为精定位区域的参考点;如图6所示;
得到精定位区域参考点位置向量:
Figure BDA0002358390680000112
根据基准星运行参数得到当前时刻4个基准星的发射功率、在LVLH坐标系下的坐标以及基准星本体坐标系相对于LVLH坐标系旋转角度,从而计算精定位区域在该当前时刻的RSS矩阵:
Figure BDA0002358390680000113
因此,稠密Radio Map表示为:
Figure BDA0002358390680000114
将成员星在线获得的RSS向量与当前时刻的稠密Radio Map中的每个参考点的RSS进行匹配,根据最近邻算法获得成员星的精定位结果为:
Figure BDA0002358390680000115
仿真实验
在成员星自主定位过程中,成员星所处的时刻和稀疏Radio Map构建的离散时刻会存在一定的时间偏移量,如图8所示。
图8中,三角形表示稀疏Radio Map建立的时刻,菱形表示当前定位时刻,ΔT表示稀疏Radio Map建立的时间间隔。而成员星会选择一个和当前定位时刻最接近的稀疏RadioMap时刻进行定位。实际定位过程中,因为ΔT始终大于0,因此当前成员星自主定位时刻和稀疏Radio Map中的时刻存在一个偏差Δt。时间偏移越小,则定位时刻越接近稀疏RadioMap的建立时刻,定位时刻基准星的位置就越接近稀疏Radio Map建立时基准星的位置。
令l0=1000m,l1=20m,l3=2.5m,时间段的个数m=186,当时间偏移量为0时,本发明所述方法和传统方法的定位误差累计分布函数(CDF)曲线和存储空间随精定位区域边长l2的变化情况如图9所示。
由图9可知,当时间偏移量为0时,传统方法的定位精度最高,本发明所述方法的定位精度低于传统方法,且精度随着精定位区域范围的增大而提高,但是由表1可知,本发明所述方法能够大大减小存储空间,并且随着精定位区域边的变化存储空间基本不变。这是因为本发明所述方法有两次建立Radio Map的过程。稀疏Radio Map的建立和传统方式一样,只不过此时参考点的间隔很大,无法得到精确的定位结果,只能得到一个大致的范围,但是稀疏Radio Map的存储空间较小;而稠密Radio Map的参考点间距很小,但其范围局限在粗定位确定的定位范围之中,该范围远远小于初始的定位范围,因此稠密Radio Map也只会占用很小的存储空间。从而本发明提出的定位方法整个的存储空间将会很小,可以运用于实际当中。
表1本发明所述方法和传统方法的存储空间
Figure BDA0002358390680000121
固定精定位区域边长l2=30m,ΔT=30s,时间偏移量Δt从-15s至15s,统计在置信概率为67%时的定位误差,并和传统方法进行对比,如图10所示。
从图10可以看出,传统方法受时间偏移的影响更大,时间偏移越大,其定位误差越大。这是因为受时间偏移的影响,定位时刻基准星的位置和其Radio Map构建时刻基准星所处位置不同,从而产生定位误差。时间偏移越大,基准星的真实位置和Radio Map构建时刻基准星的位置相差越大,从而造成定位误差增大。而本发明所述方法虽然一定程度上也会受到时间偏移的影响,但从图中可以看出,本发明所述方法受时间偏移的影响要小得多,尤其是当时间偏移较大的时候。这是因为:稀疏Radio Map虽然是以一定时间间隔建立的,粗定位结果会受到时间偏移的影响,但是稠密Radio Map则是根据当前定位时刻建立的,因此在一定程度上能够减轻由时间偏移带来的影响。
以上对本发明所提出的一种基于位置指纹的卫星集群分级定位方法,进行了详细介绍,本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发明的限制。

Claims (5)

1.一种基于位置指纹的卫星集群分级定位方法,其特征在于:所述方法分为三个阶段:构建坐标系阶段、离线阶段和在线阶段;在构建坐标系阶段建立LVLH坐标系和基准星本体坐标系;在离线阶段,由地面计算出稀疏Radio Map,并储存在成员星中;在线阶段分为三个子阶段:成员星粗定位阶段、指纹库在线计算阶段和成员星精定位阶段;在成员星粗定位阶段,成员星根据当前时刻和稀疏Radio Map选取最近时刻的RSS矩阵,并结合实时获取到的RSS向量进行粗定位;在指纹库在线计算阶段,成员星根据粗定位结果计算稠密Radio Map;在成员星精定位阶段,成员星根据实时获取到的RSS向量和计算出的稠密Radio Map进行精定位,得到成员星的最终估计位置。
2.根据权利要求1所述的方法,其特征在于:所述构建坐标系阶段具体为:
在卫星集群中基准星的数量为4,其余卫星为成员星,基准星成对进行绕飞,绕飞的中心称为“中心星”,中心星为一虚拟的点,实际上并不是卫星;成员星则在基准星附近进行伴飞;
LVLH坐标系原点OL位于卫星集群中的“中心星”,坐标轴之间满足右手定则;OLXL轴方向为中心星运动的切线方向,OLYL轴由中心星指向地球质心,OLZL轴垂直于XLYL平面;
基准星本体坐标系原点OB位于基准星的质心,坐标轴方向与基准星姿态有关;OBXB轴为天线方向图0度方向,OBYB轴为天线方向图90度方向,OBZB轴垂直于XBYB平面,OBXB轴和OBYB轴满足右手定则;并且需要对4个基准星都建立对应的基准星本体坐标系。
3.根据权利要求2所述的方法,其特征在于:所述离线阶段具体步骤为:
步骤1、在LVLH坐标系下确定边长为l0的粗定位区域,并将其均匀划分为n个边长为l1的网格,将每个边长为l1的网格的中心位置作为粗定位区域的参考点,得到参考点位置向量L:
Figure FDA0002358390670000011
其中,(xj,yj,zj),j=1,2,3…n为粗定位区域中第j个参考点在LVLH坐标系下的坐标;
步骤2、将基准星的运动周期T以间隔ΔT划分为m个时间段,每个时间段的中心时刻为:
T1,T2,…,Tm (2)
根据基准星运动方程获取Ti,i=1,2,…m时刻4个基准星位置向量:
Figure FDA0002358390670000021
步骤3、通过预先设置的参数获取第k个基准星,k=1,2,3,4,在Ti时刻的发射功率
Figure FDA0002358390670000029
第k个基准星在LVLH坐标系下的坐标为(xk(Ti),yk(Ti),zk(Ti)),第k个基准星本体坐标系相对于LVLH坐标系旋转角度为(
Figure FDA0002358390670000027
θk(Ti),ψk(Ti));
Figure FDA0002358390670000028
θk(Ti),ψk(Ti)分别为三个坐标轴相应的旋转角度;
计算在Ti时刻粗定位区域内第j个参考点相对于第k个基准星的位置:
Figure FDA0002358390670000022
步骤4、计算在Ti时刻LVLH坐标系变换到第k个基准星本体坐标系的旋转矩阵:
Figure FDA0002358390670000023
其中,
Figure FDA0002358390670000024
Figure FDA0002358390670000025
Figure FDA0002358390670000026
步骤5、计算第j个参考点在Ti时刻第k个基准星本体坐标系下的直角坐标:
Figure FDA0002358390670000031
为了使用天线方向图计算传输损耗,需要将直角坐标变换为球坐标,从直角坐标到球坐标的变换公式为:
Figure FDA0002358390670000032
于是,第j个参考点在Ti时刻第k个基准星本体坐标系下的球坐标表示为:
Figure FDA0002358390670000033
其中,rk(Ti)是第j个参考点在Ti时刻与第k个基准星之间的距离;θG,k(Ti)为第j个参考点在Ti时刻相对于第k个基准星的天顶角;
Figure FDA0002358390670000037
为第j个参考点在Ti时刻相对于第k个基准星的方位角;
步骤6、根据球坐标计算在Ti时刻天线增益:
Figure FDA0002358390670000034
其中,Gmax为发射天线最大增益,f(θG,k(Ti),
Figure FDA0002358390670000038
)为天线方向性函数,fmax是天线方向性函数最大值;
计算基准星发射的信号到达成员星的传输损耗:
Figure FDA0002358390670000035
其中,f是天线发射频率,GR是成员星接收天线增益;
于是,Ti时刻在第j个参考点处接收到第k个基准星的接收信号强度RSS为:
Figure FDA0002358390670000036
步骤7、重复步骤1至步骤6可得到Ti时刻在所有参考点处接收到4个基准星的RSS,构成如下Ti时刻RSS矩阵:
Figure FDA0002358390670000041
其中,
Figure FDA0002358390670000045
为Ti时刻在第j个参考点处接收到4个基准星的RSS构成的向量;
步骤8、重复步骤1至步骤7可得到T1,T2,…,Tm时刻的RSS矩阵:
R1,R2,…,Rm (16)
于是得到稀疏Radio Map:
M=[L R1 R2 … Rm] (17)
稀疏Radio Map由地面提前计算好并存储至每一个成员星中,成员星进行定位时调用。
4.根据权利要求3所述的方法,其特征在于:所述成员星粗定位阶段具体为:
成员星获取当前时刻接收到的4个基准星的RSS:
Figure FDA0002358390670000042
根据卫星时钟获取当前时刻t,并在时刻T1,T2,…,Tm中选择与当前时刻t最接近的时刻T*,根据式(17)得到T*时刻的稀疏Radio Map:
Figure FDA0002358390670000043
将成员星在线获得的RSS与T*时刻的稀疏Radio Map中的RSS进行匹配,根据kNN算法获得成员星的粗定位结果为:
Figure FDA0002358390670000044
5.根据权利要求4所述的方法,其特征在于:所述指纹库在线计算阶段和成员星精定位阶段具体为:
根据成员星粗定位结果确定成员星精定位区域:以成员星粗定位结果为中心,确定边长为l2的正方体区域作为成员星精定位区域;
将精定位区域按照均匀划分为w个边长为l3的网格,将每个边长为l3的网格的中心位置作为精定位区域的参考点;
得到精定位区域参考点位置向量:
Figure FDA0002358390670000051
根据基准星运行参数得到当前时刻4个基准星的发射功率、在LVLH坐标系下的坐标以及基准星本体坐标系相对于LVLH坐标系旋转角度,从而计算精定位区域在该当前时刻的RSS矩阵:
Figure FDA0002358390670000052
因此,稠密Radio Map表示为:
Figure FDA0002358390670000053
将成员星在线获得的RSS向量与当前时刻的稠密Radio Map中的每个参考点的RSS进行匹配,根据最近邻算法获得成员星的精定位结果为:
Figure FDA0002358390670000054
CN202010014588.1A 2020-01-07 2020-01-07 一种基于位置指纹的卫星集群分级定位方法 Active CN111239777B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010014588.1A CN111239777B (zh) 2020-01-07 2020-01-07 一种基于位置指纹的卫星集群分级定位方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010014588.1A CN111239777B (zh) 2020-01-07 2020-01-07 一种基于位置指纹的卫星集群分级定位方法

Publications (2)

Publication Number Publication Date
CN111239777A true CN111239777A (zh) 2020-06-05
CN111239777B CN111239777B (zh) 2023-07-25

Family

ID=70879860

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010014588.1A Active CN111239777B (zh) 2020-01-07 2020-01-07 一种基于位置指纹的卫星集群分级定位方法

Country Status (1)

Country Link
CN (1) CN111239777B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111836188A (zh) * 2020-06-17 2020-10-27 华中科技大学 一种基于Wi-Fi RSS的在线协同定位及***

Citations (35)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA714200A (en) * 1965-07-20 J. Herzl Peter Coarse-fine positioning system
EP0401973A2 (en) * 1989-05-05 1990-12-12 International Business Machines Corporation Reliable track crossing detection in data storage apparatus
JPH05504100A (ja) * 1990-04-06 1993-07-01 インターナシヨナル・ビジネス・マシーンズ・コーポレーシヨン 超高速電気力学式x、y、θ位置決めステージ
US5745869A (en) * 1995-09-28 1998-04-28 Lockheed Missiles & Space Company, Inc. Techniques for optimizing an autonomous star tracker
EP0843245A2 (de) * 1996-11-15 1998-05-20 Oerlikon Contraves AG Verfahren und Anordnung zur Positionshaltung eines geostationären Satellitenschwarms unter Benutzung einer optischen Satellitenverbindung
CN1202053A (zh) * 1997-04-25 1998-12-16 洛克希德马丁公司 利用gps卫星***快速而精确的蜂窝式电话机地理定位
TW200510751A (en) * 2003-08-13 2005-03-16 Arbitron Inc Universal display exposure monitor using personal locator service
CN1689353A (zh) * 2002-09-30 2005-10-26 摩托罗拉公司 用于在定位***中协助定位的***和方法
EP1626288A1 (fr) * 2004-08-13 2006-02-15 Frédéric Henry Système de géolocalisation hybride
CN101025439A (zh) * 1999-02-01 2007-08-29 施耐普特拉克股份有限公司 用于卫星定位***(sps)信号测量值处理的方法和装置
JP2007263637A (ja) * 2006-03-28 2007-10-11 Mitsubishi Electric Corp 測位装置及び測位方法及びプログラム
US20090303115A1 (en) * 2008-06-06 2009-12-10 Skyhook Wireless, Inc. Methods and systems for stationary user detection in a hybrid positioning system
TW201009300A (en) * 2008-08-21 2010-03-01 Mstar Semiconductor Inc Identification tag navigation system and method of the same
CN101672653A (zh) * 2008-09-12 2010-03-17 晨星软件研发(深圳)有限公司 识别标签导航***、导航装置及导航方法
JP2010112834A (ja) * 2008-11-06 2010-05-20 Tokyo Univ Of Science 変位計測システム、衛星信号受信装置および衛星測位方法
CN102138079A (zh) * 2007-03-30 2011-07-27 SiRF技术公司 改善的卫星时钟预测
CN102204396A (zh) * 2008-09-10 2011-09-28 科姆拉布斯公司 广域定位***
WO2012011624A1 (ko) * 2010-07-21 2012-01-26 (주)브이아이소프트 위치 추정에 사용자 이동성을 활용하는 와이파이 라디오 맵 기반 실내 네비게이션 시스템 및 방법
CN102679945A (zh) * 2012-06-05 2012-09-19 哈尔滨工业大学 基于三点反射合作的卫星指向与姿态测量方法与装置
CN102741701A (zh) * 2009-12-04 2012-10-17 诺基亚公司 使用压缩指纹档案的设备上定位的方法和设备
CN103238085A (zh) * 2010-11-04 2013-08-07 李雨燮 利用卫星信号发生装置的室内位置测量***及其方法
CN103323864A (zh) * 2011-10-06 2013-09-25 波音公司 来自卫星***的精确的绝对时间传递
EP2685753A1 (en) * 2012-07-11 2014-01-15 Forsvarets Forskningsinstitutt Large-scale peer-to-peer discovery mechanism for frequency allocation
KR20140099875A (ko) * 2011-11-02 2014-08-13 샤이 론넨 위치지문맵의 생성 및 이용
KR20150035745A (ko) * 2012-06-26 2015-04-07 더 거버닝 카운실 오브 더 유니버시티 오브 토론토 라디오 맵의 동적 생성을 위한 시스템, 방법 그리고 컴퓨터 프로그램
CN105044740A (zh) * 2008-04-15 2015-11-11 高通股份有限公司 通过混合sps轨道数据进行定位的方法和装置
WO2017021218A1 (en) * 2015-07-31 2017-02-09 Avanti Communications Group Plc Satellite operations support system
GB201704216D0 (en) * 2017-03-16 2017-05-03 Ranplan Wireless Network Design Ltd WIFI multi-band fingerprint-based indoor positioning
EP3230766A1 (fr) * 2014-12-12 2017-10-18 Thales Procédé et système de validation de géolocalisation par satellite
CN107728169A (zh) * 2008-11-04 2018-02-23 知维科技有限公司 对于被援助性定位***的改进
CN109459016A (zh) * 2018-11-15 2019-03-12 上海航天控制技术研究所 一种基于位置指纹的微纳卫星集群相对定位方法
US20190265365A1 (en) * 2018-02-28 2019-08-29 Robert Bosch Gmbh Method for Determining a Data Profile for the Satellite-Based Determination of a Position of a Vehicle
US10422882B1 (en) * 2018-11-01 2019-09-24 Beihang University ARAIM subset selection method and system based on BeiDou constellation
WO2019208592A1 (ja) * 2018-04-27 2019-10-31 国立大学法人静岡大学 衛星測位システム
CN110536235A (zh) * 2019-08-30 2019-12-03 哈尔滨工程大学 一种基于位置指纹筛选的室内定位方法

Patent Citations (37)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA714200A (en) * 1965-07-20 J. Herzl Peter Coarse-fine positioning system
EP0401973A2 (en) * 1989-05-05 1990-12-12 International Business Machines Corporation Reliable track crossing detection in data storage apparatus
JPH05504100A (ja) * 1990-04-06 1993-07-01 インターナシヨナル・ビジネス・マシーンズ・コーポレーシヨン 超高速電気力学式x、y、θ位置決めステージ
US5745869A (en) * 1995-09-28 1998-04-28 Lockheed Missiles & Space Company, Inc. Techniques for optimizing an autonomous star tracker
EP0843245A2 (de) * 1996-11-15 1998-05-20 Oerlikon Contraves AG Verfahren und Anordnung zur Positionshaltung eines geostationären Satellitenschwarms unter Benutzung einer optischen Satellitenverbindung
CN1202053A (zh) * 1997-04-25 1998-12-16 洛克希德马丁公司 利用gps卫星***快速而精确的蜂窝式电话机地理定位
CN101025439A (zh) * 1999-02-01 2007-08-29 施耐普特拉克股份有限公司 用于卫星定位***(sps)信号测量值处理的方法和装置
CN1689353A (zh) * 2002-09-30 2005-10-26 摩托罗拉公司 用于在定位***中协助定位的***和方法
TW200510751A (en) * 2003-08-13 2005-03-16 Arbitron Inc Universal display exposure monitor using personal locator service
EP1626288A1 (fr) * 2004-08-13 2006-02-15 Frédéric Henry Système de géolocalisation hybride
JP2007263637A (ja) * 2006-03-28 2007-10-11 Mitsubishi Electric Corp 測位装置及び測位方法及びプログラム
CN102138079A (zh) * 2007-03-30 2011-07-27 SiRF技术公司 改善的卫星时钟预测
CN105044740A (zh) * 2008-04-15 2015-11-11 高通股份有限公司 通过混合sps轨道数据进行定位的方法和装置
US20090303115A1 (en) * 2008-06-06 2009-12-10 Skyhook Wireless, Inc. Methods and systems for stationary user detection in a hybrid positioning system
TW201009300A (en) * 2008-08-21 2010-03-01 Mstar Semiconductor Inc Identification tag navigation system and method of the same
CN104122563A (zh) * 2008-09-10 2014-10-29 耐克斯特纳威公司 广域定位***
CN102204396A (zh) * 2008-09-10 2011-09-28 科姆拉布斯公司 广域定位***
CN104122564A (zh) * 2008-09-10 2014-10-29 耐克斯特纳威公司 广域定位***
CN101672653A (zh) * 2008-09-12 2010-03-17 晨星软件研发(深圳)有限公司 识别标签导航***、导航装置及导航方法
CN107728169A (zh) * 2008-11-04 2018-02-23 知维科技有限公司 对于被援助性定位***的改进
JP2010112834A (ja) * 2008-11-06 2010-05-20 Tokyo Univ Of Science 変位計測システム、衛星信号受信装置および衛星測位方法
CN102741701A (zh) * 2009-12-04 2012-10-17 诺基亚公司 使用压缩指纹档案的设备上定位的方法和设备
WO2012011624A1 (ko) * 2010-07-21 2012-01-26 (주)브이아이소프트 위치 추정에 사용자 이동성을 활용하는 와이파이 라디오 맵 기반 실내 네비게이션 시스템 및 방법
CN103238085A (zh) * 2010-11-04 2013-08-07 李雨燮 利用卫星信号发生装置的室内位置测量***及其方法
CN103323864A (zh) * 2011-10-06 2013-09-25 波音公司 来自卫星***的精确的绝对时间传递
KR20140099875A (ko) * 2011-11-02 2014-08-13 샤이 론넨 위치지문맵의 생성 및 이용
CN102679945A (zh) * 2012-06-05 2012-09-19 哈尔滨工业大学 基于三点反射合作的卫星指向与姿态测量方法与装置
KR20150035745A (ko) * 2012-06-26 2015-04-07 더 거버닝 카운실 오브 더 유니버시티 오브 토론토 라디오 맵의 동적 생성을 위한 시스템, 방법 그리고 컴퓨터 프로그램
EP2685753A1 (en) * 2012-07-11 2014-01-15 Forsvarets Forskningsinstitutt Large-scale peer-to-peer discovery mechanism for frequency allocation
EP3230766A1 (fr) * 2014-12-12 2017-10-18 Thales Procédé et système de validation de géolocalisation par satellite
WO2017021218A1 (en) * 2015-07-31 2017-02-09 Avanti Communications Group Plc Satellite operations support system
GB201704216D0 (en) * 2017-03-16 2017-05-03 Ranplan Wireless Network Design Ltd WIFI multi-band fingerprint-based indoor positioning
US20190265365A1 (en) * 2018-02-28 2019-08-29 Robert Bosch Gmbh Method for Determining a Data Profile for the Satellite-Based Determination of a Position of a Vehicle
WO2019208592A1 (ja) * 2018-04-27 2019-10-31 国立大学法人静岡大学 衛星測位システム
US10422882B1 (en) * 2018-11-01 2019-09-24 Beihang University ARAIM subset selection method and system based on BeiDou constellation
CN109459016A (zh) * 2018-11-15 2019-03-12 上海航天控制技术研究所 一种基于位置指纹的微纳卫星集群相对定位方法
CN110536235A (zh) * 2019-08-30 2019-12-03 哈尔滨工程大学 一种基于位置指纹筛选的室内定位方法

Non-Patent Citations (11)

* Cited by examiner, † Cited by third party
Title
KANARIS, L: ""Improving RSS Fingerprint-based Localization Using Directional Antennas"", 《2014 8TH EUROPEAN CONFERENCE ON ANTENNAS AND PROPAGATION (EUCAP)》 *
MA, LIN: ""3GPE: An Energy Efficient Probabilistic Fingerprint-Assisted Localization in Indoor Wi-Fi Areas"", 《2012 IEEE INTERNATIONAL CONFERENCE ON COMMUNICATIONS (ICC)》 *
MA, LIN: ""A SOCP-Based Automatic Visual Fingerprinting Method for Indoor Localization System"", 《IEEE ACCESS》 *
MA, LIN: ""Green wireless local area network received signal strength dimensionality reduction and indoor localization based on fingerprint algorithm"", 《INTERNATIONAL JOURNAL OF COMMUNICATION SYSTEMS》 *
SHIT, R.C: ""Probabilistic RSS fingerprinting for localization in smart platforms"", 《2018 INTERNATIONAL CONFERENCE ON INFORMATION TECHNOLOGY (ICIT)》 *
孙程君等: "倒GPS技术在小卫星群外测中的应用", 《高技术通讯》 *
张小红: ""低轨导航增强GNSS发展综述"", 《测绘学报》 *
李晨光等: "编队飞行卫星群相对轨道测量研究", 《北京航空航天大学学报》 *
闻新: ""卫星集群***的应用现状与发展动态"", 《人民论坛学术前沿》 *
马琳: ""基于图像关键帧的Visual-Depth Map建立方法"", 《哈尔滨工业大学学报》 *
马琳: ""高性能卫星导航接收机信号处理技术"", 《科技成果》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111836188A (zh) * 2020-06-17 2020-10-27 华中科技大学 一种基于Wi-Fi RSS的在线协同定位及***
CN111836188B (zh) * 2020-06-17 2021-05-18 华中科技大学 一种基于Wi-Fi RSS的在线协同定位及***

Also Published As

Publication number Publication date
CN111239777B (zh) 2023-07-25

Similar Documents

Publication Publication Date Title
CN108508918B (zh) 一种静轨遥感卫星数传天线高精度实时对地指向控制方法
US8775078B1 (en) Vehicle navigation using cellular networks
Wang et al. A GNSS/INS integrated navigation algorithm based on Kalman filter
CN109539903B (zh) 一种固体运载火箭椭圆转移轨道迭代制导控制方法
CN108897338B (zh) 基于pio的圆轨道航天器编队重构型防碰撞路径规划方法
CN109375172B (zh) 一种相控阵雷达解耦方法
CN103915673A (zh) 船载a-e-c三轴卫星通信天线波束指向跟踪控制方法
CN107300700B (zh) 敏捷合成孔径雷达卫星聚束模式姿态机动需求计算方法
CN111273312A (zh) 一种智能车辆定位与回环检测方法
CN110940310A (zh) 弹箭载中继测控终端相控阵天线波束指向角度计算方法
CN110456304A (zh) 机载测向定位方法
CN111239777B (zh) 一种基于位置指纹的卫星集群分级定位方法
CN109459016B (zh) 一种基于位置指纹的微纳卫星集群相对定位方法
CN112130590A (zh) 一种基于瞬时惯性系下速度补偿的星载天线对地指向确定方法
CN102023301B (zh) 一种中轨道卫星搜救***的选星方法
CN113359167A (zh) 一种通过惯性测量参数将gps与激光雷达融合定位的方法
CN103048553A (zh) 一种利用面积加权精确计算导航天线相位中心的方法
CN112213754A (zh) 一种利用无人机进行地面目标定位跟踪的方法
CN115632699A (zh) 一种目标天线的指向控制***、方法、设备及存储介质
Corazzini et al. Onboard pseudolite augmentation system for relative navigation
CN113252039B (zh) 面向地形辅助导航的粒子群快速匹配方法
CN109541587A (zh) 一种天基预警雷达功率孔径积优化设计方法
Wallace et al. Cooperative relative UAV attitude estimation using DoA and RF polarization
Yang et al. A Layout Method of Space-Based Pseudolite System Based on GDOP Geometry
DeVon et al. Innovation-based fusion of multiple satellite positioning systems for minimizing uncertainty

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