GNSS时间转换程序设计

2024-03-23 11:30

本文主要是介绍GNSS时间转换程序设计,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

GNSS中常见的时间类型:UTC、GPST、BDST、儒略日、年积日 doy(day of year)

​​​​​​本文主要围绕UTC、GPST、BDST、儒略日之间的相互转换,内容包括程序设计思路、预期功能、算例及结果分析的阐述,以及笔者对编程过程中一些常见问题和注意事项的总结。

代码详见:

https://download.csdn.net/download/m0_58307078/88614986?spm=1001.2014.3001.5503icon-default.png?t=N7T8https://download.csdn.net/download/m0_58307078/88614986?spm=1001.2014.3001.5503

1 实验任务及目的

  1. 了解常用的时间系统的定义;
  2. 掌握UTC时间、GPS时、北斗时三者之间相互转换的方法;
  3. 了解并掌握通过儒略日与GPS时、UTC时间之间的转换关系,并将其作为中间桥梁,建立GPS时与UTC时间的相互转换。

2 数据来源及编程测试环境

1.算例数据来源:PPP精密单点定位原理实验平台

2.编程环境:MATLAB 2016a

3.测试环境:MATLAB 2016a及PPP精密单点定位原理实验平台

3 实验内容及成果

3.1 GPS时与UTC的相互转换

3.1.1 GPS时转UTC——gps2utc(t_gps)

1.程序设计

预期功能:将输入参数GPS周和周内秒(week,sec )转换为UTC时(年月日时分秒)并输出。

设计思路:按UTC与GPS时的换算关系编写函数,程序计算流程及步骤如图3-1所示:

图3-1 GPS时转UTC计算流程示意

2.遇到的问题及解决过程

(1)月份计算较为复杂

问题描述:依次计算年于日时分秒(图3-1中step4)中月份计算较为复杂,各种问题层出不穷:1)起初计算思路不明确;2)月份与分钟同用m字母表示导致后期计算错误;3)没有考虑循环遍历时进行的操作会会覆盖原有正确结果。

解决过程

1)由于月份天数不同,导致判断月份较为复杂,但特定的月份内含有的天数固定,因此可通过判断累计天数判断月份,即:若为一月,累计天数不超过31天,若为二月,累计天数在区间[32,59]内,以此类推,设置for循环和if条件判断语句即可确定月份;

2)检查代码逻辑错误的过程中即可发现该错误(细小但致命);

3)排除并修正上述错误后发现依然没有得到正确结果,与其他同学的计算结果比对,定位出错点为月份计算,再次审视月份计算相关代码,发现理论上可以得到正确答案,但由于循环的设置,计算出正确结果后进行下一次循环时的计算结果会覆盖正确结果,因此通过合理使用break语句,在得到正确答案后及时跳出循环,即可保留正确结果。

月份计算代码部分如下图所示:

图3-2 月份计算部分代码截图

(2)有关闰秒调整和GPS周翻转的细节

编写代码时需要注意两个细节:1)闰秒调整要在计算年月日时分秒之前,否则调整后可能会因为进位问题导致结果错误;2)导航电文的GPS周数存储格式导致GPS周数最大为1023,当来到1023周的下一周时,GPS周溢出自动翻转为0周,即:大概每20年发生一次GPS周翻转。目前共发生过两次,第二次也是最近的一次GPS周翻转发生在2019年4月7日0时,因此在原有GPS周的基础上+2048周即为正确的GPS周。

对应代码如图3-3所示:

图3-3 闰秒调整及GPS周翻转代码截图

3.算例结果

GPS时转UTC的算例结果如下表所示:

表3-1 GPS时转UTC算例结果

GPS

Y

M

D

H

Min

S

[55,321011]

2020

4

29

17

9

53

[43,132253]

2020

2

3

12

43

55

[42,479749]

2020

1

31

13

15

31

[8,301098]

2019

6

5

11

38

0

[37,306588]

2019

12

25

13

9

30

4.用儒略日计算

上述GPS时转UTC的算法转换较为复杂,后期发现用儒略日转换更为简单。

算法流程如下图所示:

图3-4 儒略日计算GPS时转UTC流程示意

遇到的问题:日的计算公式中FRAC小数部分计算错误导致没有及时发现该小数部分对应以天为单位的时分秒。

解决方法:通过在程序中合适的位置添加断点,分步排除错误,找出错因。最终发现问题出在FRAC小数部分的计算:编程时打错了一个符号,非常致命的细节错误。

具体实现代码截图见图3-5:

图3-5 儒略日计算GPS时转UTC对应代码截图

3.1.2 UTC转GPS时——utc2gps(t_utc)

1.程序设计

预期功能:将输入参数UTC时(年月日时分秒)转换为GPS周和周内秒(week,sec )并输出。

设计思路:按UTC与GPS时的换算关系编写函数,程序计算流程及步骤如图3-6所示:

图3-6 UTC转GPS时计算流程示意

2.遇到的问题及解决过程

问题描述:对儒略日转GPS时的原理理解不透彻导致编程效率低下。

解决过程:查阅资料+仔细思考,转换时的几点主要理解如下:

  1. 计算GPS周时,由于GPS时的起算时间为1980年1月6日(区分起始时间:1980.1.1),此时的儒略日为2444244.5,因此需要将UTC转换成的儒略日减去2444244.5(如图3-7中第一行代码);
  2. 由于前面的转换都是基于UTC进行的,因此需要对转换后的GPS周内秒进行闰秒调整并进行进位判断(图3-7第三行至最后的代码)。

图3-7 儒略日转GPS时代码截图

3.算例结果

UTC转GPS时的一些算例结果如下表所示:

表3-2 UTC转GPS时算例结果

Y

M

D

H

Min

S

GPS

2019

12

25

13

9

30

9周 177044.00秒

2021

6

30

2

51

26

116周 269504.00秒

2021

10

10

0

29

58

131周 1816.00秒

2019

6

16

20

1

51

10周 72129.00秒

2020

1

12

20

38

9

40周 74307.00秒

3.2 北斗时与UTC的相互转换

1.程序设计

预期功能:基于输入的时间参数,转换为指定时间类型。值得注意的是,输入输出形式均为年月日时分秒。

设计思路:UTC与北斗时之间只相差4秒,即:BDT=UTC+4。如果直接按照两者之间的转换关系编写函数,则需要多次判断进、退位,且日、月之间的进退位情况多种多样,整个算法要逻辑严密、情况复杂,不易编写。基于上述分析,我们可以基于UTC与GPS时两者之间的转换函数进行间接转换:先利用utc2gps(),将UTC与北斗时转换为周和周内秒,通过适当的闰秒调整后,再利用gps2utc()将其转换为年月日时分秒的形式输出。

程序计算流程及步骤如图3-8所示:

图3-8 北斗时与UTC相互转换计算流程示意

2.遇到的问题及解决过程

问题描述:1)即为1.中设计思路所提到的:起初编写程序时直接按照北斗时与UTC之间的转换关系进行,当进位判断进行到月时发现情况过于复杂;2)按照设计思路中的间接转换方式重新编写程序时,对于“适当的闰秒调整”的具体实现容易出错。

解决过程:1)发现直接转换的具体实现较为复杂时立刻止损,思考是否还有其他更简单的转换方法;2)“纯思考”想不通时用笔写写思路会更清晰(见图3-9)。

图3-9 UTC和北斗时转换关系图

3.算例结果

北斗时与UTC之间转换的算例结果如下图所示:

图3-10 北斗时、UTC相互转换算例结果

3.3 GPS时与北斗时的相互转换

1.程序设计

设计思路:与北斗时、UTC相互转换类似,可利用gps2utc()和utc2gps()、基于周和周内秒进行闰秒调整,最终输出想要的时间类型。

需要注意的是:与UTC、GPS时的相互转化区分,虽然都是年月日时分秒与周和周内秒之间的换算,但时间调整不同:±14是GPS时与北斗时之间的闰秒调整,±18是为了补齐gps2utc()、utc2gps()中自带的±18s的闰秒调整。

具体转换关系及程序计算流程如图3-11所示:

图3-11 北斗时、GPS时相互转换计算流程示意

2.算例结果

GPS时与北斗时相互转换的算例结果如图3-12所示:

图3-12 GPS时、北斗时相互转换算例结果

3.4 主程序的设计

1.程序设计

预期功能:1)基于输入类型(type_in)、输出类型(type_out)判断并调用所需的时间转换函数,按特定格式输出转换后的时间;2)基于输入类型给出格式要求语句。

设计思路:通过input()函数输入获取输入、输出类型,设置嵌套if条件语句调用相应的转换函数,通过disp()、sprintf()函数格式化输出转换后的时间。具体判断思路及条件见3-13,输入时交互功能的实现见图3-14。

图3-13 调用所需的时间转换函数判断方法示意

图3-14 输入时交互的代码实现

2.运行效果

主程序运行效果如图3-15所示:

首先输入待转换时间的类型(左上图),输入后程序显示输入格式要求(右上图),按照格式要求输入待转换时间后,提示输入转化的类型(输出类型),输入后显示结果(下图)。

图3-15 输入格式正确时的交互效果

代码详见:

https://download.csdn.net/download/m0_58307078/88614986?spm=1001.2014.3001.5503icon-default.png?t=N7T8https://download.csdn.net/download/m0_58307078/88614986?spm=1001.2014.3001.5503

​​​​​​

这篇关于GNSS时间转换程序设计的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/838184

相关文章

Java中字符串转时间与时间转字符串的操作详解

《Java中字符串转时间与时间转字符串的操作详解》Java的java.time包提供了强大的日期和时间处理功能,通过DateTimeFormatter可以轻松地在日期时间对象和字符串之间进行转换,下面... 目录一、字符串转时间(一)使用预定义格式(二)自定义格式二、时间转字符串(一)使用预定义格式(二)自

python实现svg图片转换为png和gif

《python实现svg图片转换为png和gif》这篇文章主要为大家详细介绍了python如何实现将svg图片格式转换为png和gif,文中的示例代码讲解详细,感兴趣的小伙伴可以跟随小编一起学习一下... 目录python实现svg图片转换为png和gifpython实现图片格式之间的相互转换延展:基于Py

C#实现将Excel表格转换为图片(JPG/ PNG)

《C#实现将Excel表格转换为图片(JPG/PNG)》Excel表格可能会因为不同设备或字体缺失等问题,导致格式错乱或数据显示异常,转换为图片后,能确保数据的排版等保持一致,下面我们看看如何使用C... 目录通过C# 转换Excel工作表到图片通过C# 转换指定单元格区域到图片知识扩展C# 将 Excel

C++使用printf语句实现进制转换的示例代码

《C++使用printf语句实现进制转换的示例代码》在C语言中,printf函数可以直接实现部分进制转换功能,通过格式说明符(formatspecifier)快速输出不同进制的数值,下面给大家分享C+... 目录一、printf 原生支持的进制转换1. 十进制、八进制、十六进制转换2. 显示进制前缀3. 指

使用Python开发一个带EPUB转换功能的Markdown编辑器

《使用Python开发一个带EPUB转换功能的Markdown编辑器》Markdown因其简单易用和强大的格式支持,成为了写作者、开发者及内容创作者的首选格式,本文将通过Python开发一个Markd... 目录应用概览代码结构与核心组件1. 初始化与布局 (__init__)2. 工具栏 (setup_t

Java中Date、LocalDate、LocalDateTime、LocalTime、时间戳之间的相互转换代码

《Java中Date、LocalDate、LocalDateTime、LocalTime、时间戳之间的相互转换代码》:本文主要介绍Java中日期时间转换的多种方法,包括将Date转换为LocalD... 目录一、Date转LocalDateTime二、Date转LocalDate三、LocalDateTim

golang获取当前时间、时间戳和时间字符串及它们之间的相互转换方法

《golang获取当前时间、时间戳和时间字符串及它们之间的相互转换方法》:本文主要介绍golang获取当前时间、时间戳和时间字符串及它们之间的相互转换,本文通过实例代码给大家介绍的非常详细,感兴趣... 目录1、获取当前时间2、获取当前时间戳3、获取当前时间的字符串格式4、它们之间的相互转化上篇文章给大家介

Python实现AVIF图片与其他图片格式间的批量转换

《Python实现AVIF图片与其他图片格式间的批量转换》这篇文章主要为大家详细介绍了如何使用Pillow库实现AVIF与其他格式的相互转换,即将AVIF转换为常见的格式,比如JPG或PNG,需要的小... 目录环境配置1.将单个 AVIF 图片转换为 JPG 和 PNG2.批量转换目录下所有 AVIF 图

详解如何通过Python批量转换图片为PDF

《详解如何通过Python批量转换图片为PDF》:本文主要介绍如何基于Python+Tkinter开发的图片批量转PDF工具,可以支持批量添加图片,拖拽等操作,感兴趣的小伙伴可以参考一下... 目录1. 概述2. 功能亮点2.1 主要功能2.2 界面设计3. 使用指南3.1 运行环境3.2 使用步骤4. 核

Feign Client超时时间设置不生效的解决方法

《FeignClient超时时间设置不生效的解决方法》这篇文章主要为大家详细介绍了FeignClient超时时间设置不生效的原因与解决方法,具有一定的的参考价值,希望对大家有一定的帮助... 在使用Feign Client时,可以通过两种方式来设置超时时间:1.针对整个Feign Client设置超时时间