compsatpos.m

matlab版GPS卫星模拟器。给出了详细的设计说明

function R =  CompSatpos(omeg0,M0,t)
%利用星历计算卫星的位置
%omeg0  初始升交点赤经
%Mo 
%t  计算卫星位置的时刻
a = 26609000;
i = 55*pi/180;
n0 = 30.085*pi/180/3600;
wie = 15*pi/180/3600;
t0=0;
omeg0 = omeg0*pi/180;
M0 = M0*pi/180;
Mj = M0 + n0*( t - t0 );
omegG = wie*( t - t0 );
R(1) = a*(cos(Mj)*cos(omeg0 - omegG) - sin(Mj)*sin(omeg0 - omegG)*cos(i));
R(2) = a*(cos(Mj)*sin(omeg0 - omegG) + sin(Mj)*cos(omeg0 - omegG)*cos(i));
R(3) = a*sin(Mj)*sin(i);