GNSS时间转换

算法部分

GitHub完整源程序地址
https://github.com/brucema49/GNSStime_conversion.git
代码部分：

#include <iostream>
#include<string>
#include <iomanip>
#include<math.h>
#include<stdio.h> 
using namespace std;
// 格里高利日期转换为儒略日 
double def1(double a[6])
{double s1,s2,s3,jd;s1=(int ) (1461*(a[0]+4800+(int) ((a[1]-14)/12 ))/4 );s2= (int) (367*(a[1]-2-(int) ((a[1]-14)/12 )*12)/12 );s3=(int) (3*(int)((a[0]+4900+(int) ((a[1]-14)/12) )/100 )/4 );jd= a[2]-32075+s1+s2-s3-0.5+a[3]/24+a[4]/1440+a[5]/86400;return jd; }
//儒略日转换为格里高利日期 
double *def2(double jd)
{double  j,n,l1,l2,l3,y1,m1,t,t1;
//定义静态变量，防止返回地址在函数结束时被销毁 
static    double a[6]; j=(int) (jd+0.5 ),n=(int) (4*(j+68569)/146097);l1=j+68569-(int) ((n*146097+3)/4 ), y1=(int) ((4000*(l1+1)/1461001) );l2=l1-(int) ((1461*y1/4))+31,m1=(int) (80*l2/2447 );a[2]=l2-(int) (2447*m1/80 ),l3= (int) (m1/11 );a[1]=m1+2-12*l3,a[0]=(int) (100*(n-49)+y1+l3 );t=fmod(((jd+0.5-(int) jd )*24),24 ); a[3]=(int) t;t1=fmod((t-a[3])*60,60),a[4]=(int) t1,a[5]=fmod((t1-a[1])*60,60 );return a;
}
//年+年积日+天内秒转换为儒略日 
void  def3(double a[3])
{double jd1,jd2,jd,b[6];
b[0]=a[0],b[1]=1,b[2]=1,b[3]=b[4]=b[5]=0;jd1=def1(b);jd2=a[1]+jd1-1,jd=jd2+a[2]/86400;printf("儒略日为 %lf",jd);}
//儒略日转换为年+年积日+天内秒 
void  def4(double jd)
{double *a,b[6],c1;a=def2(jd);b[0]=a[0],b[1]=b[2]=1,b[3]=b[4]=b[5]=0;c1=def1(b),a[1]=jd-c1+1,a[2]=a[3]*3600+a[4]*60+a[5];cout<<a[0]<<"年 "<<a[1]<<"年积日 "<<a[2]<<" 天内秒"; 
}//GPS周+周内秒转换为儒略日 
void def5()
{cout<<"请输入GPS周+周内秒";double a[2],jd;cin>>a[0]>>a[1];jd=2444244.5+a[0]*7+a[1]/86400;printf("儒略日为%lf",jd);
}
//儒略日转换为GPS周+周内秒 
void  def6()
{double m,jd;cout<<"请输入儒略日"<<endl;cin>>jd;int z;z=(int) ((jd-2444244.5)/7);m=( (jd-2444244.5)/7-z)*604800;cout<<z<<"周"<<m<<"周内秒"; 
}
/*根据BD周+周内秒得到的儒略日具有如下约束条件：jd对应的儒略日大于(1980,1,6,0,0,0.0)对应的儒略日，即jd≥2444244.5。
//BD周+周内秒转为儒略日 */
void  def7()
{cout<<"请输入BD周+周内秒";double a[2],jd;cin>>a[0]>>a[1];jd=2453736.5+a[0]*7+a[1]/86400; printf("儒略日为%lf",jd);
}//儒略日转为BD周+周内秒 
void  def8()
{cout<<"请输入儒略日"<<endl;double jd,m;int z;cin>>jd;
z=int((jd-2453736.5)/7);
m=((jd-2453736.5)/7-z)*604800;    
cout<<z<<"BD周"<<m<<"周内秒";
}//儒略日转换为Galileo周+周内秒 
void  def9()
{cout<<"请输入Galileo周+周内秒";double a[2],jd;cin>>a[0]>>a[1];jd=2451412.5+a[0]*7+a[1]/86400;printf("儒略日为%lf",jd);
}
//儒略日转换为Galileo周+周内秒 
void  def10()
{cout<<"请输入儒略日";
double jd,m;
int z;
cin>>jd;
z= ((jd-2451412.5)/7);
m=( (jd-2451412.50 )/7-z)*604800;
printf("%dGalileo周+%lf周内秒",z,m); 
}
/* 
Galileo周+周内秒时间标示法具有如下约束条件
周:z≥0.
周内秒(Time of Week)：0≤m<604800
根据GPS周+周内秒得到的儒略日具有如下约束条件：jd对应的儒略日大于(1999,8,22,0,0,0.0)对应的儒略日，即jd≥2451412.5
*///glonasst计时转换为儒略日
void  def11()
{cout<<"请输入glonasst计时N4+Nt+h+m+s\n";double Nt,N4,a[6],Jdn,jd,h,m,s;cin>>N4>>Nt>>h>>m>>s;a[0]=1996+4*(N4-1),a[1]=a[2]=1,a[3]=a[4]=a[5]=0;Jdn=def1(a);jd=Nt-1+Jdn+h/24+m/1440+s/86400;printf("儒略日为%lf",jd);
}
//儒略日转换为glonasst 
void  def12()
{cout<<"请输入儒略日";double jd,*a,b[6],JDn,Nt,c;int h,m,s; int N4;cin>>jd;a=def2(jd);N4=(int) ((a[0]-1996)/4)+1;b[0]=1996+4*(N4-1),b[1]=b[2]=1,b[3]=b[4]=b[5]=0;JDn=def1(b);Nt=jd-JDn+1;c=Nt-(int) Nt;Nt=(int) Nt;h=(int)(c*24);m=(int)((c*24-h)*60);s=((c*24-h)*60-m)*60;cout<<"GLONASST计时为"<<"N4:"<<N4<<" Nt:"<<Nt<<" 时:"<<h<<" 分:"<<m<<" 秒:"<<s;
}int main()
{
cout<<"格里高利历日期转换为儒略日请输入1        "<<"儒略日转换成格里高利日期请输入2"<<endl;
cout<<"年+年积日+天内秒转换为儒略日请输入3      "<<"儒略日转换为年+年积日+天内秒请输入4"<<endl; 
cout<<"GPS周+周内秒转为儒略日请输入5            "<<"儒略日转换到GPS周+周内秒请输入6"<<endl;
cout<<"BD周+周内秒转为儒略日请输入7             "<<"儒略日转换到BD周+周内秒请输入8"<<endl;
cout<<"Galileo周+周内秒转为儒略日请输入9        "<<"儒略日转换到Galileo周+周内秒请输入 10"<<endl;
cout<<"GLONASST计时转换到儒略日请输入11         "<<"儒略日转换到GLONASST计时请输入12"<<endl<<"请输入选项（数字）";int i;
cin>>i;switch(i)
{ case 1:int i;cout<<"请输入格里高利日期年+月+日+时+分+秒";double a[6],jd1;for(i=0;i<6;i++){cin>>a[i];}
jd1=def1(a);
printf("儒略日日期为%lf",jd1);break;case 2:double jd;double *p;cout<<"请输入儒略日 ";cin>>jd;p= def2(jd);
cout<<"格里高利日期为"<< (int) p[0]<<"year"<<(int) p[1]<<"month"<<(int) p[2]<<"day"<<(int) p[3]<<"hour"<<(int) p[4]<<"minute"<<(int) p[5]<<"second";break;case 3:int j;double b[3];cout<<"请输入年+年积日+天内秒  "<<endl;for(j=0;j<3;j++)  cin>>b[j];def3(b);break;case 4:double jd2;cout<<"请输入儒略日  "<<endl;cin>>jd2;def4(jd2);break;case 5:def5();break;case 6:def6();break;case 7:def7();         break;    case 8:def8();break;case 9:def9();break;case 10:def10();break;case 11:def11();break;case 12 :def12();break;            
}return 0;
}
禁止一切商业转载，网络非商业转载请注明出处