GPS单点定位实验报告

1、GPS原理与应用实验 题目GPS单点定位 专业：测绘工程 班级：12-01 学号：2012212600 姓名：王威 指导教师：陶庭叶 时间：2014.11 目录 一、实验目的 3 二、实验原理3 三、实验内容3 四、实验效果图 9 五、实验总结 9 实验目的 1. 深入了解单点定位的计算过程； 2. 加强单点定位基本公式和误差方程式，法线方程式的记忆； 3. 通过上机调试程序加强动手能力的培养。 二.实验原理 一个接收机接受三个火三个以上卫星信号，得出卫星坐标和伪距, 利用间接平差计算接收机的坐标。 三实验内容 1. 程序流程图 2、实验数据 某个历元卫星坐标(单位；皿) PRNX G 3,-

2、652183.5613 G 6,-12412058. 2578 G16,-3373128. 6048 G23,10124962.7320 G31,-17262027.6771 Y 2353&896. 6765 22524932. 2093 22004627. 7810 1122T597. 8053 9U899S. 1091 Z -10392867- 6529 -5728436. S4$8 14303822. 4050 21807939. 5847 18295569.1090 伪距观测数扌居 钟差(s) 0.00059981 0. 00051935 -0. 0001026 0. 00035634

3、 -1. 4372e-5 PRN 伪距5) G 323298409. 600 G 622041347.580 G162098932.300 G2323344899.960 G3121424389.-160 3、实验程序代码 Private Sub Comma nd1_Click() Commo nDialogl.Filter = TXT files|*.txt| Common Dialog1.FilterI ndex = 1 Common Dialog1.ShowOpe n Ope n Me.Com mon Dialog1.FileName For In put As #1 Do While

4、Not EOF(1) Line In put #1, Text textbuff = textbuff + Text + vbCrLf Loop Close #1 kk = MSFIexGrid1.Rows - 1 Dim a ReDim a(kk - 1) a = Split(textbuff, vbCrLf) For j = 1 To kk For i = 1 T o 5 MSFlexGrid1.TextMatrix(j, i) = a(j - 1 + 5 * (i - 1) Next i Next j For k = 1 T o kk MSFlexGrid1.TextMatrix(k,

5、0) = 第 & k & 个点 Next k MSFlexGrid1.TextMatrix(0, 1) = X MSFlexGrid1.TextMatrix(0, 2) = Y MSFlexGrid1.TextMatrix(0, 3) = Z MSFlexGrid1.TextMatrix(0, 4) = 伪距 MSFlexGrid1.TextMatrix(0, 5) = 钟差 End Sub Private Sub Command2_Click() kk = MSFlexGrid1.Rows - 1 X0 = 0: Y0 = 0: Z0 = 0 c = 299792458 Dim a() Re

6、Dim a(kk - 1, 3) Dim ll() ReDim ll(kk - 1, 0) For ii = 1 To 100 For i = 1 T o kk l = (MSFlexGrid1.TextMatrix(i, 1) - X0) / Sqr(MSFlexGrid1.TextMatrix(i, 1) - X0) A 2+ (MSFIexGrid1.TextMatrix(i, 2)-Y0) A 2+ (MSFIexGrid1.TextMatrix(i, 3) - ZO) a 2) m = (MSFIexGrid1.T extMatrix(i, 2) - Y0) / Sqr(MSFIex

7、Grid1.TextMatrix(i, 1)- XO) A 2+ (MSFIexGrid1.TextMatrix(i, 2)- YO) A 2+ (MSFIexGrid1.TextMatrix(i, 3) - ZO) A 2) n = (MSFIexGrid1.TextMatrix(i, 3) - ZO) / Sqr(MSFIexGrid1.TextMatrix(i, 1) - XO) A 2+ (MSFIexGrid1.TextMatrix(i, 2)- YO) A 2+ (MSFIexGrid1.TextMatrix(i, 3) - ZO) A 2) a(i - 1, O) = I a(i

8、 - 1, 1) = m a(i - 1, 2) = n a(i - 1, 3) = -1 Ik = MSFIexGrid1.TextMatrix(i, 4) - Sqr(MSFIexGrid1.TextMatrix(i, 1) - XO) A 2 + (MSFIexGrid1.T extMatrix(i, 2) - YO) A 2 + (MSFIexGrid1.TextMatrix(i, 3) - ZO) A 2) + c * MSFIexGrid1.TextMatrix(i, 5) II(i - 1, O) = Ik Next i gzs = xc(qiuni(xc(zz(a), a),

9、xc(zz(a), II) XO = XO - gzs(O, O) Y0 = Y0 - gzs(1, 0) Z0 = Z0 - gzs(2, 0) j = j + 1 Next ii 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 = c X= & c 0 * Sqr(Qx(0

10、, 0) & vbCrLf & vbCrLf & c Y * Sqr(Qx(1, 1) & vbCrLf & vbC rLf & c Z= & c 0 * Sqr(Qx(2, 2) End Sub Private Sub Form_Load() MSFlexGrid1.ColWidth(1) = 1300 MSFlexGrid1.ColWidth(2) = 1300 MSFlexGrid1.ColWidth(3) = 1300 MSFlexGrid1.ColWidth(4) = 1300 Text2.Text = Text3.Text = End Sub 矩阵相减 Public Functio

11、n jian(m, n) Dim i, j As Integer UBound(n, 2) If UBound(m, 1) UBound(n, 1) Or UBound(m, 2) Then MsgBox ( 请确认输入数组是否可以相减！ ) Else Dim c() ReDim c(UBound(m, 1), UBound(n, 2) For i = 0 T o UBound(c, 1) For j = 0 To UBound(c, 2) c(i, j) = m(i, j) - n(i, j) Next j Next i jian = c End If End Function 矩阵的转置

12、Public Function zz(a) Dim i As Integer, j As Integer, t As Integer, b() If UBound(a, 1) = UBound(a, 2) Then For i = 0 T o UBound(a, 1) For j = 0 To UBound(a, 2) If i j Then t = a(i, j) a(i, j) = a(j, i) a(j, i) = t End If Next j Next i zz = a Else ReDim b(UBound(a, 2), UBound(a, 1) For i = 0 T o UBo

13、und(a, 2) For j = 0 To UBound(a, 1) b(i, j) = a(j, i) Next j Next i zz = b End If End Function 两矩阵相乘 Public Function xc(a, b) Dim i As Integer, j As Integer, k As Integer If UBound(a, 2) UBound(b, 1) Then MsgBox （ 这两个矩阵不能够相乘 ） Exit Function End If ReDim sd(UBound(a, 1), UBound(b, 2) For i = 0 T o UB

14、ound(a, 1) For j = 0 To UBound(b, 2) For k = 0 T o UBound(b, 1) sd(i, j) = sd(i, j) + a(i, k) * b(k, j) Next k Next j Next i xc = sd End Function Public Function qiuni(a) Dim c, m%, n%, p#, l%, i%, j%, ab# m = UBound(a, 1) n = UBound(a, 2) If m n Then MsgBox ( 该矩阵不可逆！ ) Exit Function End If ReDim c(

15、m, 2 * n + 1) For i = 0 T o m For j = 0 To n c(i, j) = a(i, j) Next j Next i For i = 0 T o m For j = m + 1 To 2 * m + 1 c(i, j) = 0 Next j Next i i = 0 For j = m + 1 To 2 * m + 1 c(i, j) = 1 i = i + 1 Next j For k = 0 T o n If c(k, k) = 0 Then For i = k + 1 T o n If c(i, k) 0 Then GoTo this End If N

16、ext i If i = n + 1 Then MsgBox ( 该矩阵不可逆 !) Exit Function End If this: For j = 0 To 2 * m + 1 p = c(k, j) c(k, j) = c(i, j) c(i, j) = p Next j End If ab = 1# / c(k, k) For j = 0 To 2 * m + 1 c(k, j) = c(k, j) * ab Next j For i = 0 T o n If i k Then For j = 0 To 2 * m + 1 If j k Then c(i, j) = c(i, j)

17、 - c(i, k) * c(k, j) End If Next j c(i, k) = 0 End If Next i Next k For i = 0 T o m For j = 0 To m a(i, j) = c(i, j + n + 1) a(i, j) = Round(a(i, j), 4) Next j Next i qiuni = a End Function 四实验结果图 r三m 厲 Forml n 回 S3 ! 7 L 禅差 |X=-SKDB4 L7556&3 卜邮厉02361电 2=3235573 ?Ll523t4 Q 珂顾 2E3EP60 。卩二即亚翎空咲巧 0I=l.OilGlG44396 宓溺&765 1EKP陆瞬9 述用9 600 (KE9961 22824332 2J93 5TZ3455.9O 1M1347 590 .cocsiatsl 312S.eO