GPS单点定位实验报告

GPS原理与应用实验题目：GPS单点定位

专业：测绘工程

班级：12-01

学号：2012212600

姓名：王威

指导教师：陶庭叶

时间：2014.11

目录

一、实验目的 (3)

二、实验原理 (3)

三、实验内容 (3)

四、实验效果图 (9)

五、实验总结 (9)

一．实验目的

1.深入了解单点定位的计算过程；

2.加强单点定位基本公式和误差方程式，法线方程式的记忆；

3.通过上机调试程序加强动手能力的培养。

二．实验原理

一个接收机接受三个火三个以上卫星信号，得出卫星坐标和伪距，利用间接平差计算接收机的坐标。

三．实验内容

1.程序流程图

2、实验数据

3、实验程序代码

Private Sub Command1_Click()

CommonDialog1.Filter = "TXT files|*.txt|" CommonDialog1.FilterIndex = 1 CommonDialog1.ShowOpen

Open https://www.360docs.net/doc/c512678578.html,monDialog1.FileName For Input As #1 Do While Not EOF(1)

Line Input #1, Text

textbuff = textbuff + Text + vbCrLf

Loop

Close #1

kk = MSFlexGrid1.Rows - 1

Dim a

ReDim a(kk - 1)

a = Split(textbuff, vbCrLf)

For j = 1 To kk

MSFlexGrid1.TextMatrix(j, i) = a(j - 1 + 5 * (i - 1)) Next i

Next j

For k = 1 T o kk

MSFlexGrid1.TextMatrix(k, 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()

ReDim a(kk - 1, 3)

Dim ll()

ReDim ll(kk - 1, 0)

For ii = 1 To 100

l = (MSFlexGrid1.TextMatrix(i, 1) - X0) / Sqr((MSFlexGrid1.TextMatrix(i, 1) - X0) ^ 2 + (MSFlexGrid1.TextMatrix(i, 2) - Y0) ^ 2 + (MSFlexGrid1.TextMatrix(i, 3) - Z0) ^ 2)

m = (MSFlexGrid1.T extMatrix(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((MSFlexGrid1.TextMatrix(i, 1) - X0) ^ 2 + (MSFlexGrid1.TextMatrix(i, 2) - Y0) ^ 2 + (MSFlexGrid1.TextMatrix(i, 3) - Z0) ^ 2)

a(i - 1, 0) = l

a(i - 1, 1) = m

a(i - 1, 2) = n

a(i - 1, 3) = -1

lk = MSFlexGrid1.TextMatrix(i, 4) - Sqr((MSFlexGrid1.TextMatrix(i, 1) - X0) ^ 2 + (MSFlexGrid1.T extMatrix(i, 2) - Y0) ^ 2 + (MSFlexGrid1.TextMatrix(i, 3) - Z0) ^ 2) + c * MSFlexGrid1.TextMatrix(i, 5)

ll(i - 1, 0) = lk

Next i

gzs = xc(qiuni(xc(zz(a), a)), xc(zz(a), ll))

X0 = X0 - gzs(0, 0)

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)

σ0 = Sqr(zjl(0, 0)) / (kk - 3)

Qx = qiuni(xc(zz(a), a))

Text3.Text = "σX=" & σ0 * Sqr(Qx(0, 0)) & vbCrLf & vbCrLf & "σY=" & σ0 * Sqr(Qx(1, 1)) & vbCrLf & vbC rLf & "σZ=" & σ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 Function jian(m, n)

Dim i, j As Integer

If UBound(m, 1) <> UBound(n, 1) Or UBound(m, 2) <> UBound(n, 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

'矩阵的转置

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 UBound(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 UBound(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(m, 2 * n + 1)