关于GLONASS卫星位置的计算，由广播星历进行计算，精度达到0.5m，我的原创

利用广播星历进行计算GLONASS卫星位置，精度达到0.5m。利用matlab编写的，有详细说明。

function y=f（xiyizivxvyvzxlsylszls）
format long e;
GM=398600.44E9;%-m3/s2
ae = 6378136; %地球半径；单位m
J0_2= 1082625.7E-9;%地球重力系数；
we = 7.292115e-5;%地球自转速度；
r = sqrt（xi*xi+yi*yi+zi*zi）;
f1 = -GM*xi/r^3-1.5*J0_2*GM*ae*ae*xi*（1-5*zi^2/r^2）/r^5+we*we*xi+2*we*vy+xls;
f2 = -GM*yi/r^3-1.5*J0_2*GM*ae*ae*yi*（1-5*zi^2/r^2）/r^5+we*we*yi-2*we*vx+yls;
f3 = -GM*zi/r^3-1.5*J0_2*GM*ae*ae*zi*（3-5*zi^2/r^2）/r^5+zls;
y = [f1f2f3];
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        455  2010-03-23 14:55  glonasspvt\f.asv

     文件        457  2010-03-24 11:23  glonasspvt\f.m

     文件       1980  2010-03-23 17:10  glonasspvt\pvtcal.asv

     文件       2083  2010-03-24 11:55  glonasspvt\pvtcal.m

     目录          0  2010-03-23 14:55  glonasspvt

----------- ---------  ---------- -----  ----

                 4975                    5