CN111239777A - 一种基于位置指纹的卫星集群分级定位方法 - Google Patents
一种基于位置指纹的卫星集群分级定位方法 Download PDFInfo
- 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
Links
Images
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/24—Acquisition or tracking or demodulation of signals transmitted by the system
- G01S19/25—Acquisition or tracking or demodulation of signals transmitted by the system involving aiding data received from a cooperating element, e.g. assisted GPS
- G01S19/258—Acquisition 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE 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/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing 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:
其中,(xj,yj,zj),j=1,2,3…n为粗定位区域中第j个参考点在LVLH坐标系下的坐标;
步骤2、将基准星的运动周期T以间隔ΔT划分为m个时间段,每个时间段的中心时刻为:
T1,T2,…,Tm (2)
根据基准星运动方程获取Ti,i=1,2,…m时刻4个基准星位置向量:
步骤3、通过预先设置的参数获取第k个基准星,k=1,2,3,4,在Ti时刻的发射功率第k个基准星在LVLH坐标系下的坐标为(xk(Ti),yk(Ti),zk(Ti)),第k个基准星本体坐标系相对于LVLH坐标系旋转角度为 θk(Ti),ψk(Ti)分别为三个坐标轴相应的旋转角度;
计算在Ti时刻粗定位区域内第j个参考点相对于第k个基准星的位置:
步骤4、计算在Ti时刻LVLH坐标系变换到第k个基准星本体坐标系的旋转矩阵:
其中,
步骤5、计算第j个参考点在Ti时刻第k个基准星本体坐标系下的直角坐标:
为了使用天线方向图计算传输损耗,需要将直角坐标变换为球坐标,从直角坐标到球坐标的变换公式为:
于是,第j个参考点在Ti时刻第k个基准星本体坐标系下的球坐标表示为:
步骤6、根据球坐标计算在Ti时刻天线增益:
计算基准星发射的信号到达成员星的传输损耗:
其中,f是天线发射频率,GR是成员星接收天线增益;
于是,Ti时刻在第j个参考点处接收到第k个基准星的接收信号强度RSS为:
步骤7、重复步骤1至步骤6可得到Ti时刻在所有参考点处接收到4个基准星的RSS,构成如下Ti时刻RSS矩阵:
步骤8、重复步骤1至步骤7可得到T1,T2,…,Tm时刻的RSS矩阵:
R1,R2,…,Rm (16)
于是得到稀疏Radio Map:
M=[L R1 R2 … Rm] (17)
稀疏Radio Map由地面提前计算好并存储至每一个成员星中,成员星进行定位时调用。
进一步地,所述成员星粗定位阶段具体为:
成员星获取当前时刻接收到的4个基准星的RSS:
根据卫星时钟获取当前时刻t,并在时刻T1,T2,…,Tm中选择与当前时刻t最接近的时刻T*,根据式(17)得到T*时刻的稀疏Radio Map:
将成员星在线获得的RSS与T*时刻的稀疏Radio Map中的RSS进行匹配,根据kNN算法获得成员星的粗定位结果为:
进一步地,所述指纹库在线计算阶段和成员星精定位阶段具体为:
根据成员星粗定位结果确定成员星精定位区域:以成员星粗定位结果为中心,确定边长为l2的正方体区域作为成员星精定位区域;
将精定位区域按照均匀划分为w个边长为l3的网格,将每个边长为l3的网格的中心位置作为精定位区域的参考点;
得到精定位区域参考点位置向量:
根据基准星运行参数得到当前时刻4个基准星的发射功率、在LVLH坐标系下的坐标以及基准星本体坐标系相对于LVLH坐标系旋转角度,从而计算精定位区域在该当前时刻的RSS矩阵:
因此,稠密Radio Map表示为:
将成员星在线获得的RSS向量与当前时刻的稠密Radio Map中的每个参考点的RSS进行匹配,根据最近邻算法获得成员星的精定位结果为:
本发明提出的基于位置指纹的卫星集群分级定位方法首先以较大的参考点间隔计算稀疏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:
其中,(xj,yj,zj),j=1,2,3…n为粗定位区域中第j个参考点在LVLH坐标系下的坐标;
步骤2、将基准星的运动周期T以间隔ΔT划分为m个时间段,每个时间段的中心时刻为:
T1,T2,…,Tm (2)
根据基准星运动方程获取Ti,i=1,2,…m时刻4个基准星位置向量:
步骤3、通过预先设置的参数获取第k个基准星,k=1,2,3,4,在Ti时刻的发射功率第k个基准星在LVLH坐标系下的坐标为(xk(Ti),yk(Ti),zk(Ti)),第k个基准星本体坐标系相对于LVLH坐标系旋转角度为 θk(Ti),ψk(Ti)分别为三个坐标轴相应的旋转角度;
计算在Ti时刻粗定位区域内第j个参考点相对于第k个基准星的位置:
步骤4、计算在Ti时刻LVLH坐标系变换到第k个基准星本体坐标系的旋转矩阵:
其中,
步骤5、计算第j个参考点在Ti时刻第k个基准星本体坐标系下的直角坐标:
为了使用天线方向图计算传输损耗,需要将直角坐标变换为球坐标,如图4所示,从直角坐标到球坐标的变换公式为:
于是,第j个参考点在Ti时刻第k个基准星本体坐标系下的球坐标表示为:
步骤6、根据球坐标计算在Ti时刻天线增益:
计算基准星发射的信号到达成员星的传输损耗:
其中,f是天线发射频率,GR是成员星接收天线增益;
于是,Ti时刻在第j个参考点处接收到第k个基准星的接收信号强度RSS为:
步骤7、重复步骤1至步骤6可得到Ti时刻在所有参考点处接收到4个基准星的RSS,构成如下Ti时刻RSS矩阵:
步骤8、重复步骤1至步骤7可得到T1,T2,…,Tm时刻的RSS矩阵:
R1,R2,…,Rm (16)
于是得到稀疏Radio Map:
M=[L R1 R2 … Rm] (17)
稀疏Radio Map由地面提前计算好并存储至每一个成员星中,成员星进行定位时调用。
所述成员星粗定位阶段具体为:
成员星获取当前时刻接收到的4个基准星的RSS:
根据卫星时钟获取当前时刻t,并在时刻T1,T2,…,Tm中选择与当前时刻t最接近的时刻T*,根据式(17)得到T*时刻的稀疏Radio Map:
将成员星在线获得的RSS与T*时刻的稀疏Radio Map中的RSS进行匹配,根据kNN算法获得成员星的粗定位结果为:
所述指纹库在线计算阶段和成员星精定位阶段具体为:
根据成员星粗定位结果确定成员星精定位区域:以成员星粗定位结果为中心,确定边长为l2的正方体区域作为成员星精定位区域;如图5所示;
将精定位区域按照均匀划分为w个边长为l3的网格,将每个边长为l3的网格的中心位置作为精定位区域的参考点;如图6所示;
得到精定位区域参考点位置向量:
根据基准星运行参数得到当前时刻4个基准星的发射功率、在LVLH坐标系下的坐标以及基准星本体坐标系相对于LVLH坐标系旋转角度,从而计算精定位区域在该当前时刻的RSS矩阵:
因此,稠密Radio Map表示为:
将成员星在线获得的RSS向量与当前时刻的稠密Radio Map中的每个参考点的RSS进行匹配,根据最近邻算法获得成员星的精定位结果为:
仿真实验
在成员星自主定位过程中,成员星所处的时刻和稀疏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本发明所述方法和传统方法的存储空间
固定精定位区域边长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:
其中,(xj,yj,zj),j=1,2,3…n为粗定位区域中第j个参考点在LVLH坐标系下的坐标;
步骤2、将基准星的运动周期T以间隔ΔT划分为m个时间段,每个时间段的中心时刻为:
T1,T2,…,Tm (2)
根据基准星运动方程获取Ti,i=1,2,…m时刻4个基准星位置向量:
步骤3、通过预先设置的参数获取第k个基准星,k=1,2,3,4,在Ti时刻的发射功率第k个基准星在LVLH坐标系下的坐标为(xk(Ti),yk(Ti),zk(Ti)),第k个基准星本体坐标系相对于LVLH坐标系旋转角度为(θk(Ti),ψk(Ti));θk(Ti),ψk(Ti)分别为三个坐标轴相应的旋转角度;
计算在Ti时刻粗定位区域内第j个参考点相对于第k个基准星的位置:
步骤4、计算在Ti时刻LVLH坐标系变换到第k个基准星本体坐标系的旋转矩阵:
其中,
步骤5、计算第j个参考点在Ti时刻第k个基准星本体坐标系下的直角坐标:
为了使用天线方向图计算传输损耗,需要将直角坐标变换为球坐标,从直角坐标到球坐标的变换公式为:
于是,第j个参考点在Ti时刻第k个基准星本体坐标系下的球坐标表示为:
步骤6、根据球坐标计算在Ti时刻天线增益:
计算基准星发射的信号到达成员星的传输损耗:
其中,f是天线发射频率,GR是成员星接收天线增益;
于是,Ti时刻在第j个参考点处接收到第k个基准星的接收信号强度RSS为:
步骤7、重复步骤1至步骤6可得到Ti时刻在所有参考点处接收到4个基准星的RSS,构成如下Ti时刻RSS矩阵:
步骤8、重复步骤1至步骤7可得到T1,T2,…,Tm时刻的RSS矩阵:
R1,R2,…,Rm (16)
于是得到稀疏Radio Map:
M=[L R1 R2 … Rm] (17)
稀疏Radio Map由地面提前计算好并存储至每一个成员星中,成员星进行定位时调用。
5.根据权利要求4所述的方法,其特征在于:所述指纹库在线计算阶段和成员星精定位阶段具体为:
根据成员星粗定位结果确定成员星精定位区域:以成员星粗定位结果为中心,确定边长为l2的正方体区域作为成员星精定位区域;
将精定位区域按照均匀划分为w个边长为l3的网格,将每个边长为l3的网格的中心位置作为精定位区域的参考点;
得到精定位区域参考点位置向量:
根据基准星运行参数得到当前时刻4个基准星的发射功率、在LVLH坐标系下的坐标以及基准星本体坐标系相对于LVLH坐标系旋转角度,从而计算精定位区域在该当前时刻的RSS矩阵:
因此,稠密Radio Map表示为:
将成员星在线获得的RSS向量与当前时刻的稠密Radio Map中的每个参考点的RSS进行匹配,根据最近邻算法获得成员星的精定位结果为:
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)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111836188A (zh) * | 2020-06-17 | 2020-10-27 | 华中科技大学 | 一种基于Wi-Fi RSS的在线协同定位及*** |
Citations (35)
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 | 哈尔滨工程大学 | 一种基于位置指纹筛选的室内定位方法 |
-
2020
- 2020-01-07 CN CN202010014588.1A patent/CN111239777B/zh active Active
Patent Citations (37)
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 (2)
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 |