卫星导航中的时间标示转换

最新推荐文章于 2025-01-29 23:24:20 发布

活着就要测量

最新推荐文章于 2025-01-29 23:24:20 发布

阅读量1.1k

点赞数 3

CC 4.0 BY-SA版权

文章标签： c++ 学习

本文链接：https://blog.youkuaiyun.com/2201_75704169/article/details/130046157

此程序用C++代码实现卫星导航中常用的计时方法之间的相互转化，主要包括的计时方法有：格里高利日期，儒略日，年+年积日+周内秒，GPS周+周内秒，BD周+周内秒，Galileo周+周内秒，GLONASS计时。

具体实现算法如下：

C++代码实现如下：

#include<iostream>
#include<math.h>
#include<iomanip>
using namespace std;

double fun1_1(float Y, float M, float D, float h, float m, float s)//此函数为格里高利日转儒略日
{
	if (M <= 2)
	{
		Y = Y - 1;
		M = M + 12;
	}
	double JD = floor(365.25 * Y) + floor(30.6001 * (M + 1)) + D + 1720981.5 + h / 24 + m / 1440 + s / 86400;
	return JD;
}

float* fun1_2(double JD) //此函数为儒略日转格里高利日
{
	static float arr[6];
	float L1, L2, L3, Y1, M1, T1, J, N, T;
	J = floor(JD + 0.5);
	N = floor(4 * (J + 68569) / 146097);
	L1 = J + 68569 - floor((N * 146097 + 3) / 4);
	Y1 = floor(4000 * (L1 + 1) / 1461001);
	L2 = L1 - floor(1461 * Y1 / 4) + 31;
	M1 = floor(80 * L2 / 2447.0);
	arr[2] = L2 - floor(2447 * M1 / 80);
	L3 = floor(M1 / 11);
	arr[1] = M1 + 2 - 12 * L3;
	arr[0] = floor(100 * (N - 49) + Y1 + L3);
	T = fmod((JD + 0.5 - floor(JD)) * 24, 24);
	arr[3] = floor(T);
	T1 = (T - arr[3]) * 60;
	arr[4] = floor(T1);
	arr[5] = (T1 - arr[4]) * 60;

	return arr;
}

double fun2_1(float year, float day, float sec) //此函数为年+年积日+天内秒转儒略日
{
	double JD1 = fun1_1(year, 1, 1, 0, 0, 0);
	double JD2 = day + JD1 - 1;
	double JD = JD2 + sec / 86400.0;
	return JD;
}

float* fun2_2(double JD)  //此函数为儒略日转年+年积日+天内秒
{
	static float arr1[3];
	float* arr=fun1_2(JD);
	arr1[0] = arr[0];
	double JD1