基于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
travelTime	69.1529 69.9666 71.7530 72.3446 73.5044  74.5404 74.9249  76.6978
pseudoranges	1.0e+007 * 2.0732 2.0975 2.1511  2.1688   2.2036  2.2347  2.2462  2.2993