Universal positioning based on gravity
-
摘要: 万有引力由质量引起,不容易受到干扰. 在地球及其周围空间中,对物体的万有引力影响最大的是地球,其次是地球附近的天体. 因为不同位置所受到的万有引力的大小和方向是随时间规律变化的,而物体受到的万有引力的大小和方向是可以被测量的. 因此,存在通过物体所受万有引力反推其所在位置的可能. 本文探讨了通过测量物体所受万有引力变化进行普适定位的方法,通过正演和反演仿真,证明通过万有引力进行普适定位是可行的. 并且验证了基于万有引力定位的定位误差与观测误差之间的数值关系.Abstract: Gravity is caused by mass and is not susceptible to interference. In the vicinity of the Earth, the greatest influence on the gravitational force of an object is the Earth, followed by the celestial bodies in the vicinity of the Earth. The magnitude and direction of the gravitational force varies with time from and it is regular, while the magnitude and direction of the gravitational force on an object’s location can be measured. Therefore, there is the possibility of positioning inversion by gravity. This paper discusses the universal positioning method by measuring gravity, through forward and inverse simulation, consider it is feasible to positioning through gravity. Then verify the relationship between the positioning error and the observation error.
-
Key words:
- gravity /
- universal positioning /
- inclinometer /
- gravimeter /
- underwater positioning /
- underground positioning
-
表 1 7个观测量,不同观测间隔时间的定位迭代次数和定位误差
Table 1. The number of positioning iterations and positioning errors of seven observations at different observation intervals
间隔时间/
min迭代
次数水平定位误差
(NS)/(°)水平定位误差
(EW)/(°)竖直定位
误差/m3 683 0.027 288 97 0.007 576 43 2 093.994 5 289 0.010 477 68 0.002 833 73 784.955 7 242 0.001 398 91 0.000 368 09 102.327 9 39 0.002 025 03 0.000 516 17 144.329 11 39 0.000 615 99 0.000 152 16 42.826 13 216 0.000 674 43 0.000 160 97 45.728 15 141 0.000 298 18 0.000 068 23 19.655 17 75 0.000 000 02 0.000 000 00 0.001 19 44 0.000 000 07 0.000 000 02 0.005 21 38 0.000 000 04 0.000 000 01 0.003 表 2 不同级别的观测噪声对定位结果的影响对比
Table 2. Comparison of influence of different levels of observation noise on positioning results
噪声
方差迭代
次数水平定位误差
(NS)/(°)水平定位误差
(EW)/(°)竖直定位
误差/m0 118 0.000 000 17 0.000 000 03 0.010 10−16 93 0.000 007 32 0.000 001 41 0.447 10−15 120 0.000 010 25 0.000 002 64 0.698 10−14 95 0.000 154 65 0.000 030 09 9.444 10−13 98 0.003 068 97 0.000 701 79 199.458 10−12 61 0.005 795 25 0.000 992 08 351.853 表 3 2022年12月1日不同观测时刻对定位结果的影响对比
Table 3. Comparison of influence of different observation times on positioning results
观测时刻
(UTC)迭代
次数水平定位误差
(NS)/(°)水平定位误差
(EW)/(°)竖直定位
误差/m00:00 118 0.000 000 17 0.000 000 03 0.010 02:00 54 0.000 000 02 0.000 000 00 0.001 04:00 42 0.000 000 39 0.000 000 06 0.031 06:00 44 0.000 831 12 0.000 122 98 80.545 08:00 30 0.000 000 01 0.000 000 00 0.001 10:00 28 0.000 326 02 0.000 114 48 36.780 12:00 33 0.000 005 29 0.000 002 20 0.939 14:00 61 0.000 003 48 0.000 000 61 0.780 16:00 26 0.000 123 04 0.000 005 56 24.473 18:00 25 0.000 004 89 0.000 000 28 0.890 20:00 29 0.000 003 74 0.000 000 56 0.719 22:00 29 0.000 015 17 0.000 004 78 2.512 -
[1] 李鹏, 魏宗康, 李海兵. 自主式水下无人潜航器导航系统的应用与展望[J]. 导航与控制, 2023, 22(3): 8-20. DOI: 10.3969/j.issn.1674-5558.2023.03.003 [2] 席梦寒, 武凛, 李倩倩, 等. 重力匹配导航基准图特征提取及适配性分析[J/OL]. [2023-08-24]. 武汉大学学报(信息科学版). http://kns.cnki.net/kcms/detail/42.1676.TN.20230407.1152.002.html [3] 王博, 李天姣, 李晓平. 水下惯性/重力梯度匹配导航综述[J/OL]. (2023-03-20)[2023-08-24]. 战术导弹技术. DOI: 10.16358/j.issn.1009-1300.20230020 [4] 胡双贵, 汤井田. 基于重力及其梯度张量的目标在线定位与追踪算法[J/OL]. (2023-02-15)[2023-08-24]. 地球物理学进展. http://kns.cnki.net/kcms/detail/11.2982.P.20230613.1639.042.html [5] 乔中坤, 张家俊, 张宗宇, 等. 重力梯度测量在城市空间探测中的应用[J]. 地质论评, 2023, 69(S1): 419-422. DOI: 10.16509/j.georeview.2023.s1.183 [6] 宋宏伟. 冷原子重力梯度仪发展现状及趋势[J]. 光学与光电技术, 2023, 21(3): 1-14. [7] 李晓平, 陈仕奇, 周贤高, 等. 航行过程中实时校正的重力匹配算法[J]. 中国惯性技术学报, 2023, 31(4): 352-358. [8] 黄祖珂. 潮汐原理与计算[M]. 青岛 : 中国海洋大学出版社, 2005. [9] 雷伟伟, 张捍卫, 孙茜. 精密引潮力的计算及其影响因素分析[J]. 地球物理学进展, 2016, 31(5): 1917-1926. DOI: 10.6038/pg20160505 [10] 吴立恒, 李宏, 陈征. CBT型钻孔倾斜传感器振动时效预老化系统试验研究[J]. 震灾防御技术, 2020, 15(1): 216-223. DOI: 10.11899/zzfy20200122 [11] 贾剑钢. 精密重力测量技术及其应用研究[D]. 武汉: 武汉大学, 2019. [12] NEUMEYER J, DITTFELD H J, PFLUG H, et al. Superconducting gravimeter data from potsdam - Level 1. V. 001[Z/OL]. [2023-08-24]. GFZ Data Services. https://doi.org/10.5880/igets.po.l1.001