GNSS学习笔记：PPPH源码阅读笔记

PPPH文件存储结构

1. 原始观测数据

%max=(start_time-end_time)/interval, sno是卫星数
p1s  = NaN(max,sno); %频率1的伪距观测值
p2s  = NaN(max,sno); %频率2的伪距观测值
l1s  = NaN(max,sno); %频率1的载波观测值
l2s  = NaN(max,sno); %频率2的载波观测值
eps  = NaN(max, 1); %历元数
st   = zeros(max,sno); %可用卫星

all = p1s + p2s + l1s + l2s; %把p1s/p2s/l1s/l2s相加，得到all矩阵
st(~isnan(all)) = 1; %将all中所有不是NAN值的位置在st中设为1，则得到可观测卫星信息。

对于5个NAN矩阵，矩阵行表示历元，列表示卫星（1-32 GPS，33-58 GLO, 59-88 GAL, 89-105 BDS）。如果GPS有DCB，则在p1s上加上DCB，l1s和l2s需要用原始观测值乘以对应波长。

%最后，将6个矩阵组成结构体，命名obs。
obs.p1 = p1s;
obs.p2 = p2s;
obs.l1 = l1s;
obs.l2 = l2s;
obs.ep = eps;
obs.st = st;

2. 卫星精密星历sp3文件

sp3=NAN[历元, 4, 卫星数]; %矩阵第二维包括卫星X/Y/Z和钟差。

3. 卫星精密钟差clk文件

clk=NAN[历元,卫星数]

4. 卫星相位改正atx文件

sat_neu = NaN(sn,3,2); %卫星天线相位中心偏差值改正
%sn卫星编号 3个方向 两种频率
rcv_neu = NaN( 1,3,4); %接收机天线相位中心偏差值改正
%3个方向 四种频率