基于MATLAB的GPS信号解调验证（3）

计算转换

postNavigation

通过load testfs38.mat, 加载前面的跟踪结果trackResults。

八 查找子帧   findPreambles 寻找遥测字和交接字

     需要navPartyChk 计算奇偶校验位

     需要trackResults的结果status/i_p/   

     通过前导字和子帧帧号找子帧

九 解码星历表 ephemerris  解码星历表和第一子帧的周内时

       需要生成5帧数据

       length(bits)=1500

      需要checkPhase  如果数据流中的相位反相，那么将相位翻转  需要invert

      需要twosComp2dec

     将从子帧1、2、3得到的星历转换成十进制格式

十 求伪距观测量 calculatePseudoranges

    跟踪变量、CA码起点和初值

十一 求卫星位置和卫星钟差 satpos

      求卫星参量、轨道变化、ECEF坐标

      需要check_t

十二 最小二乘法求位置  leastSqoarePos计算接收机位置

    卫星位置、钟差和伪距求位置和时钟误差。

     需要e_r_corr、

     需要topocent、需要togeod

     需要tropo

计算的结果是北纬40.008度，西经105.2627度，对应科罗拉多州的丹佛市。

另一个数据计算出来的值为

对应为意大利的都灵

计算的数值结果如下：

<

firstSubframe	4524   4529   4525   4525   4530   4536  4536     4534
track.absoluteSample	51595  82671  150897  173491  217786  257353  272039  339747
travelTime	1.3509  2.1646  3.9510  4.5426  5.7024  6.7384  7.1229  8.8958