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

计算GPS卫星位置 实验任务书题 目： 计算GPS卫星位置 专 业： 测绘工程 班 级： 12-01 学 号： 2012212600 姓 名： 王威 指导教师： 陶庭叶 时间：2014.10 目录1、 实验目的.22、 实验过程流程图.23、 实验过程代码.34、 实验效果图.85、 实验心得.89 1 实验目的：加深理解对在课堂上对GPS卫星坐标的计算，并复习以前大地测量学的相关公式和算法，提高我们的程序操作能力和知识转化能力。2 实验流程图计算卫星平均角速度 计算归化时间 计算平近点角 计算偏近点角 计算真近点角 计算升交距角 计算改正项计算改正后的升交距角、卫星矢径、轨道倾角计算卫星在平面坐标系的坐标升交角点经度的计算得出卫星在地固坐标系下的坐标3、 实验代码Private Sub Command1_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.Text) * (Val(Text20.Text) - Val(Text4.Text) 2t = Val(Text20.Text) - ttk = t - Val(Text4.Text)偏近点角Mk = Val(Text10.Text) + n * tkEk = MkFor j = 1 To 100Ek = Mk + Val(Text6.Text) * Sin(Ek)Next jEk = -0.2922365358真近点角Vk = Atn(Sqr(1 - Val(Text6.Text) 2) * Sin(Ek) / (Cos(Ek) - Val(Text6.Text)升角距角k = 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 * (1 - Val(Text6.Text) * Cos(Ek) + rik = Val(Text7.Text) + Val(Text13.Text) * tk + i卫星在轨道平面上的坐标xk = rk * Cos(uk)yk = rk * Sin(uk)升交点经度kk = 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)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 = XXKText38.Text = YYKText39.Text = ZZKEnd SubPrivate Sub Command2_Click()Text1.Text = Text2.Text = Text3.Text = Text4.Text = Text5.Text = Text6.Text = Text7.Text = Text8.Text = Text9.Text = Text10.Text = Text11.Text = Text12.Text = Text13.Text = Text14.Text = Text15.Text = Text16.Text = Text17.Text = Text18.Text = Text19.Text = Text20.Text = Text21.Text = Text22.Text = Text23.Text = Text24.Text = Text25.Text = Text26.Text = Text27.Text = Text28.Text = Text29.Text = Text30.Text = Text31.Text = Text32.Text = Text33.Text = Text34.Text = Text35.Text = Text36.Text = Text37.Text = Text38.Text = Text39.Text = End Sub读取文件Private Sub Command3_Click()CommonDialog1.ShowOpenOpen Me.CommonDialog1.FileName For Input As #2Do While Not EOF(2)Line Input #2, Texttextbuff = textbuff + Text + vbCrLfLoopClose #2Dim aReDim a(20)a = Split(textbuff, vbCrLf)Text1.Text = a(0)Text2.Text = a(1)Text3.Text = a(2)Text4.Text = a(3)Text5.Text = a(4)Text6.Text = a(5)Text7.Text = a(6)Text8.Text = a(7)Text9.Text = a(8)Text10.Text = a(9)Text11.Text = a(10)Text12.Text = a(11)Text13.Text = a(12)Text14.Text = a(13)Text15.Text = a(14)Text16.Text = a(15)Text17.Text = a(16)Text18.Text = a(17)Text19.Text = a(18)Text20.Text = a(19)End SubPrivate Sub Form_Load()Text1.Text = Text2.Text = Text3.Text = Text4.Text = Text5.Text = Text6.Text = Text7.Text = Text8.Text = Text9.Text = Text10.Text = Text11.Text = Text12.Text = Text13.Text = Text14.Text = Text15.Text = Text16.Text = Text17.Text = Text18.Text = Text19.Text = Text20.Text = Text21.Text = Text22.Text = Text23.Text = Text24.Text = Text25.Text = Text26.Text = Text27.Text = Text28.Text = Text29.Text = Text30.Text = Text31.Text = Text32.Text = Text33.Text = Text34.Text = Text35.Text = Text36.Text = Text37.Text = Text38.Text = Text39.Text = End Sub三、实验效果图四、