资源描述：
《GPS卫星位置计算(C语言)》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、#include#include#include#definebGM843.986005e14#definebOMEGAE847.2921151467e-5voidmain(){longdoubleroota=0.515365263176E+04;//轨道长半轴的平方根（根号a）longdoubletoe=0.720000000000E+04;//观测时刻toelongdoublem0=-0.290282040486E+00;//参考时刻toe的
2、平近点角longdoublee=0.678421219345E-02;//轨道偏心率elongdoubledelta_n=0.451411660250E-08;//卫星的摄动改正数△nlongdoublesmallomega=-0.258419417299E+01;//近地点角距ωlongdoublecus=0.912137329578E-05;//纬度幅角正弦调和项改正的振幅（弧度）longdoublecuc=0.189989805222E-06;//纬度幅角余弦调和项改正的振幅（弧度）longd
3、oublecrs=0.406250000000E+01;//轨道半径的余弦调和项改正的振幅（m）longdoublecrc=0.201875000000E+03;//轨道半径的正弦调和项改正的振幅（m）longdoublecis=0.949949026108E-07;//轨道倾角的余弦调和项改正的振幅（弧度）longdoublecic=0.130385160446E-07;//轨道倾角的正弦调和项改正的振幅（弧度）longdoubleidot=-0.253939149013E-09;//轨道倾角变化
4、率Ilongdoublei0=0.958512160302E+00;//轨道倾角（弧度）longdoublebigomega0=-0.137835982556E+01;//升交点赤经longdoubleearthrate=bOMEGAE84;//地球自转的速率welongdoublebigomegadot=-0.856928551657e-08;longdoublet=0.720000000000E+04;//加入卫星钟差改正的归化时间longdoubleA;longdoublen0=0,n,tk;
5、longdoublemk,ek,tak,ik,omegak,phik,uk,rk;longdoublecorr_u,corr_r,corr_i;longdoublexpk,ypk,xk,yk,zk;inti;printf("输入的数据：");printf("√a=%e",roota);printf("toe=%e",toe);printf("e=%e",e);printf("i0=%e",i0);printf("ω=%e",smallomega);printf("△n=%e
6、n",delta_n);printf("Ω0=%e",bigomega0);printf("I=%e",idot);printf("Cuc=%e",cuc);printf("Cus=%e",cus);printf("Crc=%e",crc);printf("Crs=%e",crs);printf("Cic=%e",cic);printf("Cis=%e",cis);printf("输出的结果为：",e);A=roota*roota;n0=sqrt(bGM8
7、4/(A*A*A));//平均角速度n0printf("n0=%.10lf",n0);tk=t-toe;//相对于参考时刻toe的归化时间tkprintf("tk=%.10lf",tk);n=n0+delta_n;//加摄动改正后的卫星平均角速度printf("n=%.10lf",n);mk=m0+n*tk;//卫星平近点角printf("mk=%.10lf",mk);ek=mk;for(i=0;i<10;i++)ek=mk+e*sin(ek);//利用迭代法求偏近点角ekprint
8、f("ek=%.10lf",ek);tak=atan2(sqrt(1.0-e*e)*sin(ek),cos(ek)-e);//真近点角Vk的计算printf("Vk=%.10lf",tak);phik=tak+smallomega;//升交距角φk的计算printf("φk=%.10lf",phik);corr_u=cus*sin(2.0*phik)+cuc*cos(2.0*phik);//升交距角u的摄动改正δuprintf("δu=%.10lf",
