资源描述:
《GPS卫星位置计算.doc》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、GPS卫星位置计算实验学校:工业大学学号:20104168专业班级:测绘工程10-1班学院:土木与水利工程学院指导教师:陶庭叶一、实验目的1、通过对GPS卫星位置的计算,增强我们对《GPS测量原理及应用》课程的理解,使我们牢固掌握GPS测量的基本原理和过程,熟悉GPS测量数据处理的基本技能和计算方法。2、提高运用计算机语言编程开发能力,运用测量平差知识、数学知识和计算机知识,通过给定的程序算法,编制程序实现GPS卫星位置的计算过程。二、实验容及要求1、通过课堂对GPS测量原理的学习,在课后自主完成GPS卫星位置的计算过程并按照课本上的步
2、骤设计程序完成GPS卫星位置的计算过程;2、要求每位同学独立完成GPS卫星位置的计算过程,编写程序,调试程序,并编写程序设计文档。要求过程和算确、程序运行正确、设计文档完备;三、课程设计工具运用自己熟悉的编程开发语言(C、C#、C++、VC、VS、VB、FORTRAN等)。四、实验总结1、实验简单介绍运行后的主界面比以前做的程序要漂亮些,给界面增加了一副图片,让界面看起来还是比较的漂亮,但由于起始的参数较多,因此,界面整体看起来比较的拥挤,但这个不影响程序的计算过程。整个程序的功能不是很多,但满足了基本的要求,能够进行卫星的位置计算,并
3、且能够达到很好的精度,与课本上给出的结果相差很小,符合实验的要求。为了避免繁琐的数据的输入,在本次程序中我增加了一个的功能按钮,点击后,程序自动给相关的起始数据赋值并显示在界面上,如图所示:这些数据是教材《GPS测量原理与应用》41页的卫星星历数据;用户也可以根据需要更改相应的数据,如上图所示,当点击按钮时,就可以在对话框相应的位置计算出卫星的坐标,如下图所示:其中,计算出的卫星的坐标为Xk=4589210.3510074820,Yk=25371005.6999580190,Zk=-5618292.2998269377,这是卫星在地固坐
4、标系中的坐标,如果还知道极移参数就可以计算出卫星在协议地球坐标系中的坐标了。用户也可以选择清除计算的结果,点击即可完成,如下图所示:用户也可以选择重新输入数据,点击就可以将输入的参数清零,如下图所示。每次的默认数据均相同,用户可以选择输入默认数据后再更改需要更改的地方,来完成计算过程。当用户选择了重新输入后,即所有的起始数据均清零后,此时不能参与计算,若用户点击了按钮,则系统会提示用户输入的数据有误!当运行程序显示主界面时,如果用户同样不小心点击了按钮,则会出现同样的提示信息,以此来提示用户没有输入起始数据,避免产生不必要的错误。当数据
5、清零后点击会提示用户输入的数据有误!实验心得GPS卫星位置的计算实验不是很难,因为实验过程在书上很详细,只要按照步骤进行,计算过程就不会出现问题。我起始是先把程序调试正确,在调试的过程中,采取了简单的直接给参数赋值的方法,因此,很快就把程序调试正确了。在后期的修改工作中,我修改掉了直接赋值的地面观测时刻对应的标准GPS时间t,而是通过计算过程得出的。这样要不了多久就把整个程序调通了,调试过程没有花很多的时间。在查了一些资料后,我给对话框的界面增加了一幅位图,并在图片上打印出字体,调整界面的编辑框的分布后看上去整个程序就变得挺不错的,尽管
6、界面看起来有点儿拥挤,但并不影响程序的计算过程的正确性。整体来说效果不错,但显然需要改进的地方很多,就像刚才说的,界面整体看起来比较拥挤,这个更改起来比较容易,如下图所示:这样看起来就比较的舒服,不会导致眼睛疲劳,但中间的计算过程均没有,调试时我们可以选择前者,后期修改时我们就可以去掉中间的计算过程。这是做的比较好的方面。程序可能还曾在很多我没有发觉的问题,还请老师多多指教。实验源代码(部分)voidCGPSDlg::OnJisuan(){UpdateData();if(!(m_a0
7、
8、m_a1
9、
10、m_a2
11、
12、m_Cic
13、
14、m_Cis
15、
16、
17、m_Crc
18、
19、m_Crs
20、
21、m_Cuc
22、
23、m_Cus
24、
25、m_dn
26、
27、m_e&&m_roota
28、
29、m_I
30、
31、m_i0
32、
33、m_M0
34、
35、m_Q
36、
37、m_Q0
38、
39、m_toc
40、
41、m_toe
42、
43、m_tph
44、
45、m_tpm
46、
47、m_w)){MessageBox("输入数据有误,请重新输入!");SetFocus();}else{//n0,n,tk,Mk,Ek,Vk,cosVk,sinVk,du,dr,di,uk,rk,ik,xk,yk,,Xk,Yk,Zk均为课本上参数//Ek1为趋近计算迭代参数,Fk为φk,Qk为Ωk,ωdoublen0,n,tk,Mk
48、,Ek,Ek1,cosVk,sinVk,Vk,Fk,du,dr,di,uk,rk,ik;doublexk,yk,Qk,Xk,Yk,Zk,t,t1,dn;charx[40],y[40],z[40];//cha