GPS报告卫星轨道程序报告.doc

《GPS报告卫星轨道程序报告.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、.计算GPS卫星位置实验任务书题目：计算GPS卫星位置专业：测绘工程班级：12-01学号：2012212600姓名：王威指导教师：庭叶时间：2014.10..目录一、实验目的....................................................2二、实验过程流程图........................................2三、实验过程代码............................................3四、实验效果图........

2、.......................................8五、实验心得...................................................8六、..一．实验目的：加深理解对在课堂上对GPS卫星坐标的计算，并复习以前测量学的相关公式和算法，提高我们的程序操作能力和知识转化能力。二．实验流程图计算卫星平均角速度计算归化时间计算平近点角计算偏近点角计算真近点角计算升交距角计算改正项计算改正后的升交距角、卫星矢径、轨道倾角计算卫星在平面坐标系的坐标升交角点经度的计

3、算得出卫星在地固坐标系下的坐标..三、实验代码PrivateSubCommand1_Click()μ=398600500000000#ωe=0.0000729211567'卫星运行平均角速度n0=Sqr(μ)/Val(Text5.Text)/Val(Text5.Text)/Val(Text5.Text)n=n0+Val(Text11.Text)'归化时间Δt=Val(Text1.Text)+Val(Text2.Text)*(Val(Text20.Text)-Val(Text4.Text))+Val(Text3.

4、Text)*(Val(Text20.Text)-Val(Text4.Text))^2t=Val(Text20.Text)-Δttk=t-Val(Text4.Text)'偏近点角Mk=Val(Text10.Text)+n*tkEk=MkForj=1To100Ek=Mk+Val(Text6.Text)*Sin(Ek)NextjEk=-0.2922365358'真近点角Vk=Atn((Sqr(1-Val(Text6.Text)^2)*Sin(Ek))/(Cos(Ek)-Val(Text6.Text)))'升角距角Φk

5、=Vk+Val(Text8.Text)Text21.Text=Φk'摄动改正项的计算δu=Val(Text15.Text)*Cos(2*Φk)+Val(Text14.Text)*Sin(2*Φk)δr=Val(Text19.Text)*Cos(2*Φk)+Val(Text18.Text)*Sin(2*Φk)δi=Val(Text17.Text)*Cos(2*Φk)+Val(Text16.Text)*Sin(2*Φk)'改正后的升交距角、卫星矢径、轨道倾角uk=Φk+δurk=Val(Text5.Text)^2*

6、(1-Val(Text6.Text)*Cos(Ek))+δrik=Val(Text7.Text)+Val(Text13.Text)*tk+δi'卫星在轨道平面上的坐标xk=rk*Cos(uk)yk=rk*Sin(uk)'升交点经度ΩkΩk=Val(Text9.Text)+(Val(Text12.Text)-ωe)*tk-ωe*Val(Text4.Text)'w卫星在地固坐标系的坐标XXK=xk*Cos(Ωk)-yk*Cos(ik)*Sin(Ωk)YYK=xk*Sin(Ωk)+yk*Cos(ik)*Cos(Ωk)

7、ZZK=yk*Sin(ik)..Text21.Text=n0Text22.Text=nText23.Text=tkText24.Text=MkText25.Text=EkText26.Text=VkText27.Text=ΦkText28.Text=δuText29.Text=δrText30.Text=δiText31.Text=ukText32.Text=rkText33.Text=ikText34.Text=xkText35.Text=ykText36.Text=ΩkText37.Text=XXKText

8、38.Text=YYKText39.Text=ZZKEndSubPrivateSubCommand2_Click()Text1.Text=""Text2.Text=""Text3.Text=""Text4.Text=""Text5.Text=""Text6.Text=""Text7.Text=""Text8.Text=""Text9.Text=""Text10.Text=""Text1