（二）基于MATLAB的RINEX3.02格式星历匹配

以历元中所接收到的PRN数目作为循环次数去匹配该历元卫星PRN

function  [date_0]=match(Obs,Nav)
C=2.997925*10^8;
for i=1:size(Obs.GPS,2) 
    date_0(i).epoch=Obs.GPS(i).t;
    m=0;    
    for a=1:length(Obs.GPS(i).PRN)   
        for j=1:length(Nav)    
            if (Obs.GPS(i).PRN(a)==Nav(j).prn)
               tk=Obs.GPS(i).t-Nav(j).toe;   
                if abs(tk)<7200
                   m=m+1;
                   PseudoRange1=Obs.GPS(i).PR1(a);
                   PseudoRange2=Obs.GPS(i).PR2(a);
                   date_0(i).time(:,m)=[j;Nav(j).prn;Nav(j).toe;tk;PR1;PR2];
                   break;
                end     
            end  
        end
    end
end