Precision_ephemeris_interpolation.c
精密星历内插示例
用于readrinexobs2.rar的修正,也是用VS2017编写,加上了每颗卫星不是8个数据的检验,但是到第75554秒以后又出现新问题了,他的秒头介绍行从一行变两行了,我懒得重新写代码,所以目前只能到75554秒了,不过还好数据量够了,其实37681行也够了,毕竟一天才86400秒
用于RINEX观测文件的读取,将RINEX转化为2个txt文件,epoch.txt文件是数据文件每秒具有的卫星数,obsdata.txt文件是每秒的8个数据(C1 L1 D1 S1 P2 L2 D2 S2),按照卫星数,每8个一循环。读obsdata.txt的时候用matlab的textread()一个一个的读,然后根据epoch.txt知道每秒多少个卫星,由于每个卫星有8个数据,所以按照8个一循环读取。KMIT2000.20o文件到第37681秒数据之前都是好的,该行之后文件格式有变化,数据文件有缺失,不是每颗卫星都是8个数据文件。所以本程序只能读到37681秒数据,后续会把有不是8个数据的卫星的秒数去除掉。下一个文件就是做这个功能,加上了检验每颗卫星是不是8个数据的功能。 本程序在VS2017上实现
该文件可以将CSV数据转化为.bin二进制数据,我们用射频采集的数据是CSV格式的,不利于数据的读取,可以转化为.bin文件
gps L1 射频信号,1.85S,可用来gps基带捕获跟踪算法的检验,用matlab fopen fread以int8形式打开,l类似C++的char形式存储,一个字节存储一个数据,+1 +3 -1 -3,2bit,采样频率16.368MHz,理论中频4.092MHz