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

计算GPS卫星位置实验任务书题目:计算GPS卫星位置专业:测绘工程班级:12-01学号:2012212600姓名:王威指导教师:庭叶时间:....................................................2实验过程流程图........................................2实验过程代码............................................3实验效果图...............................................8实验心得...................................................8实验目的:加深理解对在课堂上对GPS卫星坐标的计算,并复****以前测量学的相关公式和算法,提高我们的程序操作能力和知识转化能力。实验流程图计算卫星平均角速度计算归化时间计算平近点角计算偏近点角计算真近点角计算升交距角计算改正项计算改正后的升交距角、卫星矢径、mand1_Click()μ=398600500000000#ωe='卫星运行平均角速度n0=Sqr(μ)/Val()/Val()/Val()n=n0+Val()'归化时间Δt=Val()+Val()*(Val()-Val())+Val()*(Val()-Val())^2t=Val()-Δttk=t-Val()'偏近点角Mk=Val()+n*tkEk=MkForj=1To100Ek=Mk+Val()*Sin(Ek)NextjEk=-'真近点角Vk=Atn((Sqr(1-Val()^2)*Sin(Ek))/(Cos(Ek)-Val()))'升角距角Φk=Vk+Val()=Φk'摄动改正项的计算δu=Val()*Cos(2*Φk)+Val()*Sin(2*Φk)δr=Val()*Cos(2*Φk)+Val()*Sin(2*Φk)δi=Val()*Cos(2*Φk)+Val()*Sin(2*Φk)'改正后的升交距角、卫星矢径、轨道倾角uk=Φk+δurk=Val()^2*(1-Val()*Cos(Ek))+δrik=Val()+Val()*tk+δi'卫星在轨道平面上的坐标xk=rk*Cos(uk)yk=rk*Sin(uk)'升交点经度ΩkΩk=Val()+(Val()-ωe)*tk-ωe*Val()'w卫星在地固坐标系的坐标