c 计算delta_生活丨由广播星历(导航电文)简单计算卫星位置

//写这个花了一天的时间 主要是那一块费了很久

//于是专门把时间背景知识加进来做一个总结

//经常感觉很口渴 思路也很闷 就像上面的停了水的喷灌器

 一、时间背景知识 

---

 二、不同时间系统之间的转换 

 利用导航电文计算卫星位置时 需要将目标计算时刻换算到GPS时间 

//由于导航电文中的是GPS周内的秒数

//进行时间转换时：需要先将(通常是格里高利历)转换到，再将<儒略历下时刻>转换到

//再进行T_GPS与TOE做差时，还要考虑加减604800s，以保证(T_GPS-TOE)在区间[-302400, 302400]，即不跨越周

---

 三、导航电文的数据结构 

格式说明采用的是Fortran程序设计语言中的格式说明方式, 一个格式说明项通常具有如下形式: [r]fw.[m] 其中:
r : 重复因子, 表示后面的内容将重复的次数, 该部分是可选的;
f: 数据类型符, 在RINEX 格式的说明中, 用到了如下数据类型: X, 空格
    A, 字符型
    I , 整型
    F, 单精度浮点型
    D, 双精度浮点型
w: 字段宽度;
m: 在字段中最少的数字或字符数, 当数据类型为单精度浮点型或双精度浮点型时, 表示小数位数, 该部分是可选的。 例如:
格式说明符“ F9.2, 11X, A1, 19X”表示
  这一行的内容从第1 列开始依次是宽度为9 位、小数点后有2 位的单精度浮点数;
  11 个空格;
  宽度为1 的字符串;
  19 个空格;
格式说明符“3F14.4”表示这一行内容从第1列开始依次是3个宽度为14、小数点后有4位的单精度浮点数；而格式说明符“7( 3X, A1, I2) ”则表示这一行内容从第1 列开始, 将3 个空格、宽度为1的字符串和宽度为2 的整型这些内容重复7 次。

---

 四、导航电文的计算流程 

//首先计算卫星在 轨道平面坐标系 下的坐标 //然后将卫星在轨道坐标系下的坐标分别 绕X轴旋转-i角、绕Z轴旋转-Wk(升交点经度)角 ，求出卫星在地心地固系下的坐标