资源描述：
《GPS单点定位实验报告.doc》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。
1、GPS原理与应用实验题目：GPS单点定位专业：测绘工程班级：12-01学号：2012212600姓名：王威指导教师：陶庭叶时间：2014.11目录一、实验目的......................................................3二、实验原理………..........................................3三、实验内容……..............................................3四、实验效果图.......................
2、...........................9五、实验总结......................................................9一．实验目的1.深入了解单点定位的计算过程；2.加强单点定位基本公式和误差方程式，法线方程式的记忆；3.通过上机调试程序加强动手能力的培养。二．实验原理一个接收机接受三个火三个以上卫星信号，得出卫星坐标和伪距，利用间接平差计算接收机的坐标。三．实验内容1.程序流程图输入原始数据打印输出精度评定列出伪距观测方程循环多次间接平差得出新的接收机坐标X,Y,Z组成误差方
3、程组AiX+Li=Vi得出各个坐标的l，m，n伪距方程线性化2、实验数据3、实验程序代码PrivateSubCommand1_Click()CommonDialog1.Filter="TXTfiles.txt"CommonDialog1.FilterIndex=1CommonDialog1.ShowOpenOpenMe.CommonDialog1.FileNameForInputAs#1DoWhileNotEOF(1)LineInput#1,Texttextbuff=textbuff+Text+vbCrLfLoopClose#1kk=M
4、SFlexGrid1.Rows-1DimaReDima(kk-1)a=Split(textbuff,vbCrLf)Forj=1TokkFori=1To5MSFlexGrid1.TextMatrix(j,i)=a(j-1+5(i-1))NextiNextjFork=1TokkMSFlexGrid1.TextMatrix(k,0)="第"&k&"个点"NextkMSFlexGrid1.TextMatrix(0,1)="X"MSFlexGrid1.TextMatrix(0,2)="Y"MSFlexGrid1.TextMatrix(0,3)="
5、Z"MSFlexGrid1.TextMatrix(0,4)="伪距"MSFlexGrid1.TextMatrix(0,5)="钟差"EndSubPrivateSubCommand2_Click()kk=MSFlexGrid1.Rows-1X0=0:Y0=0:Z0=0c=299792458Dima()ReDima(kk-1,3)Dimll()ReDimll(kk-1,0)Forii=1To100Fori=1Tokkl=(MSFlexGrid1.TextMatrix(i,1)-X0)/Sqr((MSFlexGrid1.TextMatrix(
6、i,1)-X0)^2+(MSFlexGrid1.TextMatrix(i,2)-Y0)^2+(MSFlexGrid1.TextMatrix(i,3)-Z0)^2)m=(MSFlexGrid1.TextMatrix(i,2)-Y0)/Sqr((MSFlexGrid1.TextMatrix(i,1)-X0)^2+(MSFlexGrid1.TextMatrix(i,2)-Y0)^2+(MSFlexGrid1.TextMatrix(i,3)-Z0)^2)n=(MSFlexGrid1.TextMatrix(i,3)-Z0)/Sqr((MSFlex
7、Grid1.TextMatrix(i,1)-X0)^2+(MSFlexGrid1.TextMatrix(i,2)-Y0)^2+(MSFlexGrid1.TextMatrix(i,3)-Z0)^2)a(i-1,0)=la(i-1,1)=ma(i-1,2)=na(i-1,3)=-1lk=MSFlexGrid1.TextMatrix(i,4)-Sqr((MSFlexGrid1.TextMatrix(i,1)-X0)^2+(MSFlexGrid1.TextMatrix(i,2)-Y0)^2+(MSFlexGrid1.TextMatrix(i,3
8、)-Z0)^2)+cMSFlexGrid1.TextMatrix(i,5)ll(i-1,0)=lkNextigzs=xc(qiuni(xc(zz(a),a)),xc(zz(a),ll))X0=X0-gzs(
