getepochobsinfo.m

差分载波定位

%读取一个历元观测数据
function [obsInfo, isEnd] = getEpochObsInfo(obsTypeNum,fid)
%obsTypeNum:观测值的类型个数
%obsInfo：存放观测信息的结构体
%Satprns：当前观测到的卫星PRN号数组
%Values： 观测值信息
%GPST:    观测历元的GPS时
isEnd = 0;
%读取第一行信息
line = fgetl(fid);  
if line == -1
    isEnd = 1;
end
char = 'G';
[prns,ind] = strSplit(line,char); %字符串解析,prns:‘G’后面的字符；ind：第一个‘G’在哪一列出现。
                                  % strSplit的返回值是数值类型
svnum = length(prns);
char = ' ';
info = line(1:ind-2);
time = strSplit(info,char);
%观测历元的GPS时: 星期六\星期日的零时到当前的时间长度
obsInfo.GPST = UTC2GPST(time(1:6));
%按行读取,读取了svnum个观测类型的观测数据
format long;
obsInfo.Values = zeros(svnum,obsTypeNum);
k = 1;
for i = 1 : svnum
    line = fgetl(fid);
    val = strSplit(line,char);
    if (length(val) == obsTypeNum)
        obsInfo.Values(k,:) = val;
        obsInfo.SatPrns(k) = prns(i);
        k = k + 1;
    end
end