[发明专利]一种基于位置指纹的卫星集群分级定位方法有效
申请号: | 202010014588.1 | 申请日: | 2020-01-07 |
公开(公告)号: | CN111239777B | 公开(公告)日: | 2023-07-25 |
发明(设计)人: | 马琳;黄鹏飞;谭学治;王孝 | 申请(专利权)人: | 哈尔滨工业大学 |
主分类号: | G01S19/25 | 分类号: | G01S19/25;G06F17/16 |
代理公司: | 哈尔滨奥博专利代理事务所(普通合伙) 23220 | 代理人: | 叶以方 |
地址: | 150001 黑龙*** | 国省代码: | 黑龙江;23 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 一种 基于 位置 指纹 卫星 集群 分级 定位 方法 | ||
1.一种基于位置指纹的卫星集群分级定位方法,其特征在于:所述方法分为三个阶段:构建坐标系阶段、离线阶段和在线阶段;在构建坐标系阶段建立LVLH坐标系和基准星本体坐标系;在离线阶段,由地面计算出稀疏Radio Map,并储存在成员星中;在线阶段分为三个子阶段:成员星粗定位阶段、指纹库在线计算阶段和成员星精定位阶段;在成员星粗定位阶段,成员星根据当前时刻和稀疏Radio Map选取最近时刻的RSS矩阵,并结合实时获取到的RSS向量进行粗定位;在指纹库在线计算阶段,成员星根据粗定位结果计算稠密Radio Map;在成员星精定位阶段,成员星根据实时获取到的RSS向量和计算出的稠密Radio Map进行精定位,得到成员星的最终估计位置;
所述离线阶段具体步骤为:
步骤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个基准星本体坐标系下的球坐标表示为:
其中,rk(Ti)是第j个参考点在Ti时刻与第k个基准星之间的距离;θG,k(Ti)为第j个参考点在Ti时刻相对于第k个基准星的天顶角;为第j个参考点在Ti时刻相对于第k个基准星的方位角;
步骤6、根据球坐标计算在Ti时刻天线增益:
其中,Gmax为发射天线最大增益,为天线方向性函数,fmax是天线方向性函数最大值;
计算基准星发射的信号到达成员星的传输损耗:
其中,f是天线发射频率,GR是成员星接收天线增益;
于是,Ti时刻在第j个参考点处接收到第k个基准星的接收信号强度RSS为:
步骤7、重复步骤1至步骤6可得到Ti时刻在所有参考点处接收到4个基准星的RSS,构成如下Ti时刻RSS矩阵:
其中,为Ti时刻在第j个参考点处接收到4个基准星的RSS构成的向量;
步骤8、重复步骤1至步骤7可得到T1,T2,…,Tm时刻的RSS矩阵:
R1,R2,…,Rm(16)
于是得到稀疏RadioMap:
M=[L R1 R2 … Rm] (17)
稀疏RadioMap由地面提前计算好并存储至每一个成员星中,成员星进行定位时调用;
所述成员星粗定位阶段具体为:
成员星获取当前时刻接收到的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进行匹配,根据最近邻算法获得成员星的精定位结果为:
2.根据权利要求1所述的方法,其特征在于:所述构建坐标系阶段具体为:
在卫星集群中基准星的数量为4,其余卫星为成员星,基准星成对进行绕飞,绕飞的中心称为“中心星”,中心星为一虚拟的点,实际上并不是卫星;成员星则在基准星附近进行伴飞;
LVLH坐标系原点OL位于卫星集群中的“中心星”,坐标轴之间满足右手定则;OLXL轴方向为中心星运动的切线方向,OLYL轴由中心星指向地球质心,OLZL轴垂直于XLYL平面;
基准星本体坐标系原点OB位于基准星的质心,坐标轴方向与基准星姿态有关;OBXB轴为天线方向图0度方向,OBYB轴为天线方向图90度方向,OBZB轴垂直于XBYB平面,OBXB轴和OBYB轴满足右手定则;并且需要对4个基准星都建立对应的基准星本体坐标系。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于哈尔滨工业大学,未经哈尔滨工业大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/202010014588.1/1.html,转载请声明来源钻瓜专利网。
- 同类专利
- 专利分类
G01S 无线电定向;无线电导航;采用无线电波测距或测速;采用无线电波的反射或再辐射的定位或存在检测;采用其他波的类似装置
G01S19-00 卫星无线电信标定位系统;利用这种系统传输的信号确定位置、速度或姿态
G01S19-01 .传输时间戳信息的卫星无线电信标定位系统,例如,GPS [全球定位系统]、GLONASS[全球导航卫星系统]或GALILEO
G01S19-38 .利用卫星无线电信标定位系统传输的信号来确定导航方案
G01S19-39 ..传输带有时间戳信息的卫星无线电信标定位系统,例如GPS [全球定位系统], GLONASS [全球导航卫星系统]或GALILEO
G01S19-40 ...校正位置、速度或姿态
G01S19-42 ...确定位置