作业--卫星坐标计算.docVIP

GPS卫星导航定位技术与方法作业报告之 卫星坐标计算  1 作业任务 ------------------------------------------------------------------------------------ 3 2 作业思想 --------------------------------------------------------------------------------------- 33 作业条件及数据-------------------------------------------------------------------- 3 4 作业过程 --------------------------------------------------------------------------- 4 5 源程序 ----------------------------------------------------------------------------- 5 6 计算结果 --------------------------------------------------------------------------- 9 7心得体会与建议----------------------------------------------------------------------------- 9 1 作业任务 根据教材所给广播星历参数，编程计算UTC2004年1月30日8点0分00秒—20分00秒，每隔1分钟的卫星号为7的卫星坐标。 2作业思想 根据教材提供广播星历参数，按照卫星坐标计算步骤一步一步求解表示卫星位置的量，最后求出卫星在地固坐标系中的空间直角坐标。根据历元间隔，计算出共21个历元的卫星坐标。 此次作业的已知条件及计算步骤均为教材提供且思路清晰，简明易懂，只要在编程过程中注意个已知量、中间量和待求量的表示及数据对应正确，涉及角度的计算时小心避免错误即可。 3作业条件及数据 由GPS卫星导航文件（表3.1）提取出7号卫星的轨道参数及其摄动改正项 表3.1 提取及整理后的广播星历参数： 星历参数 参数值 参数含义 Δn 4.85305928777E-09 给定的星历参数计算值的平均速度no之改正数(弧度) √a 5+03 卫星轨道长半轴的平方根 t0e 4.60800000000E+05 星历的参考时刻，从星期日零时开始量度，以秒计 M0 -2.14992476872E+00 参考时刻t0e的平近点角（弧度） e 1.25958691351E-02 轨道偏心率 ω -1.85137234909E+00 近地点角距（弧度） Cuc 2.75485217571E-06 升交点角距的余弦调和改正项的振幅(弧度) Cus 9.03755426407E-06 升交点角距的正弦调和改正项的振幅(弧度) Crc 1.91875000000E+02 轨道半径的余弦调和改正项的振幅（弧度） Crs 5.37812500000E+01 轨道半径的正弦调和改正项的振幅（弧度） Cic 2.06753611565E-07 轨道倾角的余弦调和改正项的振幅（弧度） Cis -7.26431608200E-08 轨道倾角的正弦调和改正项的振幅（弧度） IDOT -1.31434044204E-10 轨道倾角的变率 i0 9.37690094716E-01 参考时刻的轨道倾角（弧度） Ω -7.91890109042E-09 升交点赤经的变率（弧度/s） Ω0 1.93198315097E+00 GPS周开始时刻的升交点经度（弧度） 4作业过程 （10）升交点经度：λt =Ω0 +（Ω-ωe）（t - t0e）-ωe t0e 式中：ωe为地球自转角速度，ωe =7.2921151467×10-5rad/s；Ω为升交点赤经变化率，Ω0为GPS周开始时刻的升交点经度，Ω和Ω0均由星历文件给出；t0为参考时刻，一般取t0 = 0，为一周的开始（星期日的子夜）；t0e为广播星历的参考历元时刻，从t0开始起算； （11）卫星在地固坐标系中的空间直角坐标为： 式中：，为旋转矩阵，将其代入并展开后得 5 源程序 #include iomanip.h #include stdlib.h #include math.h #includefstream.h #includeiostream.h const int k=20;//历元数 doubl