GPS单点定位实验报告.docx
- 文档编号:436255
- 上传时间:2022-10-10
- 格式:DOCX
- 页数:11
- 大小:48.28KB
GPS单点定位实验报告.docx
《GPS单点定位实验报告.docx》由会员分享,可在线阅读,更多相关《GPS单点定位实验报告.docx(11页珍藏版)》请在冰豆网上搜索。
GPS单点定位实验报告
GPS原理与应用实验
题目GPS单点定位
专业:
测绘工程
班级:
12-01
学号:
2012212600
姓名:
王威
指导教师:
陶庭叶
时间:
2014.11
一、实验目的3
二、实验原理3
三、实验内容……3
四、实验效果图9
五、实验总结9
实验目的
1.深入了解单点定位的计算过程;
2.加强单点定位基本公式和误差方程式,法线方程式的记忆;
3.通过上机调试程序加强动手能力的培养。
二.实验原理
一个接收机接受三个火三个以上卫星信号,得出卫星坐标和伪距,利用间接平差计算接收机的坐标。
三•实验内容
1.程序流程图
2、实验数据
某个历元卫星坐标(单位;皿)
PRNX
G3,-652^183.5613
G6,-12412058.2578
G16,-3373128.6048
G23,10124962.7320
G31,-17262027.6771
Y
2353&896.6765
22524932.2093
22004627.78101122T597.80539U899S.1091
Z
-10392867-6529-5728436.S4$814303822.405021807939.5847
18295569.1090
伪距观测数扌居
钟差(s)
0.00059981
0.00051935
-0.0001026
0.00035634
-1.4372e-5
PRN伪距5)
G323298409.600
G622041347.580
G1620^98932.300
G2323344899.960
G3121424389.-160
3、实验程序代码
PrivateSubCommand1_Click()
CommonDialogl.Filter="TXTfiles|*.txt|"
CommonDialog1.FilterIndex=1
CommonDialog1.ShowOpen
OpenMe.CommonDialog1.FileNameForInputAs#1
DoWhileNotEOF
(1)
LineInput#1,Texttextbuff=textbuff+Text+vbCrLf
Loop
Close#1kk=MSFIexGrid1.Rows-1
Dima
ReDima(kk-1)a=Split(textbuff,vbCrLf)
Forj=1Tokk
Fori=1To5
MSFlexGrid1.TextMatrix(j,i)=a(j-1+5*(i-1))
Nexti
Nextj
Fork=1Tokk
MSFlexGrid1.TextMatrix(k,0)="第"&k&"个点"
Nextk
MSFlexGrid1.TextMatrix(0,1)="X"
MSFlexGrid1.TextMatrix(0,2)="Y"
MSFlexGrid1.TextMatrix(0,3)="Z"
MSFlexGrid1.TextMatrix(0,4)="伪距"
MSFlexGrid1.TextMatrix(0,5)="钟差"
EndSub
PrivateSubCommand2_Click()
kk=MSFlexGrid1.Rows-1
X0=0:
Y0=0:
Z0=0c=299792458
Dima()
ReDima(kk-1,3)
Dimll()
ReDimll(kk-1,0)
Forii=1To100
Fori=1Tokk
l=(MSFlexGrid1.TextMatrix(i,1)-X0)/Sqr((MSFlexGrid1.TextMatrix(i,1)
-X0)A2+(MSFIexGrid1.TextMatrix(i,2)-Y0)A2+
(MSFIexGrid1.TextMatrix(i,3)-ZO)a2)
m=(MSFIexGrid1.TextMatrix(i,2)-Y0)/Sqr((MSFIexGrid1.TextMatrix(i,
1)-XO)A2+(MSFIexGrid1.TextMatrix(i,2)-YO)A2+
(MSFIexGrid1.TextMatrix(i,3)-ZO)A2)
n=(MSFIexGrid1.TextMatrix(i,3)-ZO)/Sqr((MSFIexGrid1.TextMatrix(i,1)
-XO)A2+(MSFIexGrid1.TextMatrix(i,2)-YO)A2+
(MSFIexGrid1.TextMatrix(i,3)-ZO)A2)
a(i-1,O)=I
a(i-1,1)=m
a(i-1,2)=n
a(i-1,3)=-1
Ik=MSFIexGrid1.TextMatrix(i,4)-Sqr((MSFIexGrid1.TextMatrix(i,1)-XO)
A2+(MSFIexGrid1.TextMatrix(i,2)-YO)A2+(MSFIexGrid1.TextMatrix(i,
3)-ZO)A2)+c*MSFIexGrid1.TextMatrix(i,5)
II(i-1,O)=Ik
Nexti
gzs=xc(qiuni(xc(zz(a),a)),xc(zz(a),II))
XO=XO-gzs(O,O)
Y0=Y0-gzs(1,0)
Z0=Z0-gzs(2,0)j=j+1
Nextii
Text2.Text="X="&X0&vbCrLf&vbCrLf&"Y="&Y0&vbCrLf&
vbCrLf&"Z="&Z0
V=jian(ll,xc(a,gzs))
zjl=xc(zz(V),V)
o0=Sqr(zjl(0,0))/(kk-3)
Qx=qiuni(xc(zz(a),a))
Text3.Text="cX="&c0*Sqr(Qx(0,0))&vbCrLf&vbCrLf&"cY
*Sqr(Qx(1,1))&vbCrLf&vbCrLf&"cZ="&c0*Sqr(Qx(2,2))
EndSub
PrivateSubForm_Load()
MSFlexGrid1.ColWidth
(1)=1300
MSFlexGrid1.ColWidth
(2)=1300
MSFlexGrid1.ColWidth(3)=1300
MSFlexGrid1.ColWidth(4)=1300
Text2.Text=""
Text3.Text=""
EndSub'矩阵相减
PublicFunctionjian(m,n)
Dimi,jAsInteger
UBound(n,2)
IfUBound(m,1)<>UBound(n,1)OrUBound(m,2)<>
Then
MsgBox("请确认输入数组是否可以相减!
")
Else
Dimc()
ReDimc(UBound(m,1),UBound(n,2))
Fori=0ToUBound(c,1)
Forj=0ToUBound(c,2)
c(i,j)=m(i,j)-n(i,j)
Nextj
Nexti
jian=c
EndIf
EndFunction
'矩阵的转置
PublicFunctionzz(a)
DimiAsInteger,jAsInteger,tAsInteger,b()
IfUBound(a,1)=UBound(a,2)Then
Fori=0ToUBound(a,1)
Forj=0ToUBound(a,2)
Ifi a(i,j)=a(j,i) a(j,i)=t EndIf Nextj Nextizz=a Else ReDimb(UBound(a,2),UBound(a,1)) Fori=0ToUBound(a,2) Forj=0ToUBound(a,1) b(i,j)=a(j,i) Nextj Nexti zz=b EndIf EndFunction '两矩阵相乘 PublicFunctionxc(a,b) DimiAsInteger,jAsInteger,kAsIntegerIfUBound(a,2)<>UBound(b,1)Then MsgBox("这两个矩阵不能够相乘") ExitFunction EndIf ReDimsd(UBound(a,1),UBound(b,2)) Fori=0ToUBound(a,1) Forj=0ToUBound(b,2) Fork=0ToUBound(b,1) sd(i,j)=sd(i,j)+a(i,k)*b(k,j) Nextk Nextj Nexti xc=sd EndFunction PublicFunctionqiuni(a) Dimc,m%,n%,p#,l%,i%,j%,ab# m=UBound(a,1) n=UBound(a,2) Ifm<>nThen MsgBox("该矩阵不可逆! ! ! ") ExitFunction EndIf ReDimc(m,2*n+1) Fori=0Tom Forj=0Ton c(i,j)=a(i,j) Nextj Nexti Fori=0Tom Forj=m+1To2*m+1 c(i,j)=0 Nextj Nexti i=0 Forj=m+1To2*m+1 c(i,j)=1 i=i+1 Nextj Fork=0Ton Ifc(k,k)=0Then Fori=k+1Ton Ifc(i,k)<>0Then GoTothis EndIf Nexti Ifi=n+1Then MsgBox("该矩阵不可逆! ! ! ") ExitFunction EndIf this: Forj=0To2*m+1 p=c(k,j) c(k,j)=c(i,j)c(i,j)=p
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- GPS 单点 定位 实验 报告
![提示](https://static.bdocx.com/images/bang_tan.gif)