%升交点经度%根据prn码选择GEO轨道卫星 we=7.292115e-5; if (i >= 1 && i <= 120) || (i >= 1009 && i <= 1056)%选择GEO轨道特殊计算 omegak=OMEGA0+OMEGADOT*tk-we*toe; XGk=xk*cos(omegak)-yk*cos(ik)*sin(omegak); YGk=xk*sin(omegak)+yk*cos(ik)*cos(omegak); ZGk=yk*sin(...