GPS_DataFromat

subframe 1

第1子帧、第2子帧、第3子帧的内容通常每两小时更新一次，并且通常发生在整小时交界处。

遥测字(TLM)

每一个子帧的第一个字均为遥测字，因而它在导航电文中每6秒出现一次。下图中显示了遥测字内部码位的分布情况，

固定二进制值10001011：同步码，1：8bit，GPS接收机用它来匹配接收到的数据码，进而搜索、锁定子帧起始边沿，为接下来按照响应格式正确解译二进制数据提供了必要条件，如果接收机找到了一个卫星信号的子帧边沿，那么我们称该接收机对此卫星信号进入了子帧同步状态（或者称为帧同步状态）；

遥测码：9:22bit，提供特许用户所需要的信息；

保留位：23:24 bit

奇偶检验码：25:30 bit

交接字 (HOW)

31:47bit：从Z计数器上得到的截短的周内时计数值。

48bit：交接字的第18bit是警告标志，当标志为1时，提醒非特许用户自己承担使用该卫星信号的风险，因为该卫星第一数据块所提供的URA值有可能比其真实值还要大。

49bit：A-S标志，其值为1时表示对该卫星实施了反电子欺骗措施。

50:52bit：子帧识别标志，共有5个有效二进制值（即：001表示该子帧是第1子帧，010表示第2子帧，011—3，100—4，101—5）。

53:54bit：保留位；

55:60bit：奇偶检验位；

其他字

星期数：WN，61：70bit，

用户测距精度：URA（accuracy），73：76bit，用户测距精度是对所有由GPS地面监控部分和空间星座部分引起的测距误差大小的一个统计值，他是通过导航电文中的一个4bit表示的用户测距精度因子N而提供给非特许用户的。用户测距精度因子N的值在0至15之间，每个值对应于一个用户测距精度URA，而用户可以根据N值用一下公式估计URA值，该值越大，则表示从该卫星信号中得到的GPS距离测量值的精度越低。当N等于15时，URA的估计值缺省，此时用户要自己承担使用该卫星的风险。

卫星健康状况：heltay，77：82 bit，若共计6bit的卫星健康状况的最高位是0，则表示导航电文全部正确。若它的最高位是1，则表示导航电文出错，而低5位又具体指出信号各部分的出错情况。

群波延时校正值：$T_{GD}$，97：204 bit，该值只适合于单频（L1或L2）接收机，而双频接收机则无需此项校正。

时钟数据期号：IODC，83：84bit和197:204 bit，用10bit表示的IODC是时钟数据块的“期刊号”，一个IODC值对应一套时钟校正参数。因为IODC的值在七天内不会出现重复，它可以用来帮助用户接收机快速监测时钟校正参数是否发生了变化。

参考时间$t_{oc}$：t_oc，219:234bit，以秒为单位的时钟数据参考时间

时钟校正参数：${\alpha }_{f0}$（271：292 bit），${\alpha }_{f1}$（249：264 bit），${\alpha }_{f2}$（241：248 bit）。卫星时钟校正模型方程中的三个系数。另外，参数$t_{oc}$称为第一数据块的参考时间，它在时钟校正模型中被用作时间参考点。

eph.weekNumber  = bin2dec(subframe(61:70)) + 1024; % GPS星期数，从1980开始计算
eph.accuracy    = bin2dec(subframe(73:76));        % 
eph.health      = bin2dec(subframe(77:82));
eph.T_GD        = twosComp2dec(subframe(197:204)) * 2^(-31);
eph.IODC        = bin2dec([subframe(83:84) subframe(197:204)]);
eph.t_oc        = bin2dec(subframe(219:234)) * 2^4;
eph.a_f2        = twosComp2dec(subframe(241:248)) * 2^(-55);
eph.a_f1        = twosComp2dec(subframe(249:264)) * 2^(-43);
eph.a_f0        = twosComp2dec(subframe(271:292)) * 2^(-31);

subframe 2

第2子帧和第3子帧属于GPS中的第2数据块。

IODE：IODE_sf2，61：68 bit，星历数据期号，用来标记一套星历参数。

$C_{rs}$：C_rs， 69：84 bit，轨道半径正弦调和校正振幅；

$\Delta n$：deltan，91：106 bit，平均运动角速度校正值；

$M_0$：M_0， 107：114 bit 和121：144 bit，$t_{oe}$时的平近点角；

$C_{uc}$：C_uc， 151：166bit ，升交点角距余弦调和校正振幅；

$e_s$ ：e， 167：174 bit和181:204，轨道偏心率

$C_{us}$：C_us， 211：226 bit，升交点角距正弦调和校正振幅；

$\sqrt{a_s}$：sqrtA，227:234 bit和241:264bit，卫星轨道长半轴as的平方根；

$t_{oe}$ ：t_oe， 276：286 bit，星历参考时间；

eph.IODE_sf2    = bin2dec(subframe(61:68));
eph.C_rs        = twosComp2dec(subframe(69: 84)) * 2^(-5);
eph.deltan      = ...
                 twosComp2dec(subframe(91:106)) * 2^(-43)...
                 * gpsPi;
eph.M_0         = ...
                 twosComp2dec([subframe(107:114) subframe(121:144)]) ...
                 * 2^(-31) * gpsPi;
eph.C_uc        = twosComp2dec(subframe(151:166)) * 2^(-29);
eph.e           = ...
                 bin2dec([subframe(167:174) subframe(181:204)]) ...
                 * 2^(-33);
eph.C_us        = twosComp2dec(subframe(211:226)) * 2^(-29);
eph.sqrtA       = ...
                 bin2dec([subframe(227:234) subframe(241:264)]) ...
                 * 2^(-19);
eph.t_oe        = bin2dec(subframe(271:286)) * 2^4;

subframe 3

$C_{ic}$：C_ic， 61:76bit，轨道倾角余弦调和校正振幅；

${\Omega }_o$：omega_0，77:84bit和91:114bit，周内时等于0时的轨道升交点赤经；

$C_{is}$：C_is， 121:136bit，轨道倾角正弦调和校正振幅；

$i_0$：i_0， 137:144bit和151:174bit，t_oe时的轨道倾角；

$C_{rc}$：C_rc， 181:196bit，轨道半径余弦调和校正振幅；

$\omega$：omega， 197:204bit和211:234bit，轨道近地角距；

$\dot{\Omega }$：omegaDot，241:264bit，轨道升交点赤经对时间的变化率；

IODE：IODE_Sf3，271:278 bit，星历数据期号，用来标记一套星历参数。

$\dot{i}$：iDot， 279:292bit，轨道倾角对时间的变化率；

eph.C_ic        = twosComp2dec(subframe(61:76)) * 2^(-29);
eph.omega_0     = ...
                 twosComp2dec([subframe(77:84) subframe(91:114)]) ...
                 * 2^(-31) * gpsPi;
eph.C_is        = twosComp2dec(subframe(121:136)) * 2^(-29);
eph.i_0         = ...
                 twosComp2dec([subframe(137:144) subframe(151:174)]) ...
                 * 2^(-31) * gpsPi;
eph.C_rc        = twosComp2dec(subframe(181:196)) * 2^(-5);
eph.omega       = ...
                 twosComp2dec([subframe(197:204) subframe(211:234)]) ...
                 * 2^(-31) * gpsPi;
eph.omegaDot    = twosComp2dec(subframe(241:264)) * 2^(-43) * gpsPi;
eph.IODE_sf3    = bin2dec(subframe(271:278));
eph.iDot        = twosComp2dec(subframe(279:292)) * 2^(-43) * gpsPi;