user_orbit.m

GPS/GALILEO/BEIDOU兼容接收机自主完好性检测程序

function [xut,yut,zut]=user_orbit(t)

%*********************************
% ** 用户轨迹 
% **********************************

X=[];Y=[];Z=[];
for n=1:t
    [long,lat,h]=usercontrail(n);
    we=(7.2921151467e-5)*360/(2*pi);            %因为地球的自转，飞机相对于天球要比相对于地球每秒多走we弧度
    long=long+we*n;
    [xu,yu,zu]=user_ecef(long,lat,h);
    X=[X,xu];
    Y=[Y,yu];
    Z=[Z,zu];
end
xut=X(1,:);
yut=Y(1,:);
zut=Z(1,:);