compsatvel.m

来自「matlab版GPS卫星模拟器。给出了详细的设计说明」· M 代码 · 共 20 行

function V =  CompSatvel(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 );
v = n0*a;
Vx= v*[cos(Mj+pi/2)*cos(omeg0 - omegG)-sin(Mj+pi/2)*sin(omeg0 - omegG)*cos(i)];
Vy= v*[cos(Mj+pi/2)*sin(omeg0 - omegG)+sin(Mj+pi/2)*cos(omeg0 - omegG)*cos(i)];
Vz= v*sin(Mj+pi/2)*sin(i);
V = [Vx,Vy,Vz];