衛(wèi)星導(dǎo)航gps典型例題編程報(bào)告

1、GPSMatlab編程報(bào)告編程報(bào)告一、通過(guò)星歷計(jì)算衛(wèi)星坐標(biāo)（修正）1、算法流程算法流程（1）計(jì)算GPS衛(wèi)星運(yùn)行的平均速度n0nnn?????203snaaa???（2）計(jì)算歸化時(shí)間t?oettt???（3）計(jì)算計(jì)算觀測(cè)歷元的平近點(diǎn)角tM0MMnt???（4）計(jì)算偏近點(diǎn)角E利用不動(dòng)點(diǎn)迭代法求解此方程：得出sinEMeE??E不動(dòng)點(diǎn)迭代法的迭代公式為：??1012kkxxk?????（5）計(jì)算衛(wèi)星的地心矢徑0r??01cosraeE??（6

2、）計(jì)算真近點(diǎn)角f12arctan(tan)12eEfe?????????（7）計(jì)算升交點(diǎn)角距0?0f????（8）計(jì)算攝動(dòng)改正項(xiàng)：uri???000000sin2cos2sin2cos2sin2cos2uusucrrsrciisicCCCCCC???????????????（9）計(jì)算經(jīng)過(guò)攝動(dòng)改正的升交點(diǎn)角距，衛(wèi)星矢徑，和軌道面傾角?riToe=172800%星歷參考時(shí)間IODC=10%鐘差數(shù)據(jù)齡期URAI=0%用戶距離精度標(biāo)志IODE=

3、10%星歷數(shù)據(jù)齡期Toc=172800%本時(shí)段鐘差參數(shù)參考時(shí)間a0=2.72893e005%衛(wèi)星鐘差改正0階多項(xiàng)式系數(shù)a1=6.7531e014%衛(wèi)星鐘差改正1階多項(xiàng)式系數(shù)a2=5.787e018%衛(wèi)星鐘差改正2階多項(xiàng)式系數(shù)%定義常量c=2.99792458e8%光速omegae=7.2921151467e5%地球自轉(zhuǎn)角速度mu=3.986004418e14%地球引力常數(shù)GM%1、計(jì)算GPS衛(wèi)星運(yùn)行的平均速度na=SqrtASqrtAn

4、=sqrt(mu(a^3))Delta_n%2、計(jì)算歸化時(shí)間Delta_tDelta_t=tToe%3、計(jì)算觀測(cè)歷元t的平近點(diǎn)角MM=Mu0nDelta_t%4、計(jì)算偏近點(diǎn)角Eeps=1e20E=Mtol=1while(toleps)%不動(dòng)點(diǎn)迭代法E0=EE=M(Ecc)sin(E0)tol=abs(EE0)end%5、計(jì)算衛(wèi)星的地心矢徑r0r0=a(1Ecccos(E))%6、計(jì)算真近點(diǎn)角f%f=2atan(sqrt((1Ecc)(1