1、-. z.GPS原理與應用實驗題 目： GPS單點定位 專 業： 測繪工程 班 級： 12-01 學 號： 2012212600 姓 名： 王威 指導教師： 陶庭葉 時間：2014.11目錄實驗目的.3實驗原理.3實驗內容.3實驗效果圖.9實驗總結.9一 實驗目的1.深入了解單點定位的計算過程；2.加強單點定位根本公式和誤差方程式，法線方程式的記憶；3.通過上機調試程序加強動手能力的培養。二實驗原理 一個接收機承受三個火三個以上衛星信號，得出衛星坐標和偽距，利用間接平差計算接收機的坐標。 三實驗內容1.程序流程圖輸入原始數據打印輸出精度評定列出偽距觀測方程循環屢次間接平差得出新的接收機坐標*,

2、Y,Z組成誤差方程組Ai*+Li=Vi得出各個坐標的l，m，n偽距方程線性化2、實驗數據3、實驗程序代碼Private Sub mand1_Click()monDialog1.Filter = T*T files|*.t*t|monDialog1.FilterInde* = 1monDialog1.ShowOpenOpen Me.monDialog1.FileName For Input As #1Do While Not EOF(1)Line Input #1, Te*tte*tbuff = te*tbuff + Te*t + vbCrLfLoopClose #1kk = MSFle*Gri

3、d1.Rows - 1Dim aReDim a(kk - 1)a = Split(te*tbuff, vbCrLf)For j = 1 To kk For i = 1 To 5 MSFle*Grid1.Te*tMatri*(j, i) = a(j - 1 + 5 * (i - 1)Ne*t iNe*t jFor k = 1 To kkMSFle*Grid1.Te*tMatri*(k, 0) = 第 & k & 個點Ne*t kMSFle*Grid1.Te*tMatri*(0, 1) = *MSFle*Grid1.Te*tMatri*(0, 2) = YMSFle*Grid1.Te*tMatri

4、*(0, 3) = ZMSFle*Grid1.Te*tMatri*(0, 4) = 偽距MSFle*Grid1.Te*tMatri*(0, 5) = 鐘差End SubPrivate Sub mand2_Click()kk = MSFle*Grid1.Rows - 1*0 = 0: Y0 = 0: Z0 = 0c = 299792458Dim a()ReDim a(kk - 1, 3)Dim ll()ReDim ll(kk - 1, 0)For ii = 1 To 100For i = 1 To kk l = (MSFle*Grid1.Te*tMatri*(i, 1) - *0) / Sqr(

5、MSFle*Grid1.Te*tMatri*(i, 1) - *0) 2 + (MSFle*Grid1.Te*tMatri*(i, 2) - Y0) 2 + (MSFle*Grid1.Te*tMatri*(i, 3) - Z0) 2) m = (MSFle*Grid1.Te*tMatri*(i, 2) - Y0) / Sqr(MSFle*Grid1.Te*tMatri*(i, 1) - *0) 2 + (MSFle*Grid1.Te*tMatri*(i, 2) - Y0) 2 + (MSFle*Grid1.Te*tMatri*(i, 3) - Z0) 2) n = (MSFle*Grid1.T

6、e*tMatri*(i, 3) - Z0) / Sqr(MSFle*Grid1.Te*tMatri*(i, 1) - *0) 2 + (MSFle*Grid1.Te*tMatri*(i, 2) - Y0) 2 + (MSFle*Grid1.Te*tMatri*(i, 3) - Z0) 2) a(i - 1, 0) = l a(i - 1, 1) = m a(i - 1, 2) = n a(i - 1, 3) = -1lk = MSFle*Grid1.Te*tMatri*(i, 4) - Sqr(MSFle*Grid1.Te*tMatri*(i, 1) - *0) 2 + (MSFle*Grid

7、1.Te*tMatri*(i, 2) - Y0) 2 + (MSFle*Grid1.Te*tMatri*(i, 3) - Z0) 2) + c * MSFle*Grid1.Te*tMatri*(i, 5)ll(i - 1, 0) = lkNe*t igzs = *c(qiuni(*c(zz(a), a), *c(zz(a), ll)*0 = *0 - gzs(0, 0)Y0 = Y0 - gzs(1, 0)Z0 = Z0 - gzs(2, 0)j = j + 1Ne*t iiTe*t2.Te*t = *= & *0 & vbCrLf & vbCrLf & Y= & Y0 & vbCrLf &

8、vbCrLf & Z= & Z0V = jian(ll, *c(a, gzs)zjl = *c(zz(V), V)0 = Sqr(zjl(0, 0) / (kk - 3)Q* = qiuni(*c(zz(a), a)Te*t3.Te*t = *= & 0 * Sqr(Q*(0, 0) & vbCrLf & vbCrLf & Y= & 0 * Sqr(Q*(1, 1) & vbCrLf & vbCrLf & Z= & 0 * Sqr(Q*(2, 2)End SubPrivate Sub Form_Load()MSFle*Grid1.ColWidth(1) = 1300MSFle*Grid1.Co

9、lWidth(2) = 1300MSFle*Grid1.ColWidth(3) = 1300MSFle*Grid1.ColWidth(4) = 1300Te*t2.Te*t = Te*t3.Te*t = End Sub矩陣相減Public Function jian(m, n)Dim i, j As IntegerIf UBound(m, 1) UBound(n, 1) Or UBound(m, 2) UBound(n, 2) ThenMsgBo* (請確認輸入數組是否可以相減！)ElseDim c()ReDim c(UBound(m, 1), UBound(n, 2)For i = 0 To

10、 UBound(c, 1)For j = 0 To UBound(c, 2)c(i, j) = m(i, j) - n(i, j)Ne*t jNe*t ijian = cEnd IfEnd Function矩陣的轉置Public Function zz(a)Dim i As Integer, j As Integer, t As Integer, b()If UBound(a, 1) = UBound(a, 2) ThenFor i = 0 To UBound(a, 1)For j = 0 To UBound(a, 2)If i j Thent = a(i, j)a(i, j) = a(j,

11、i)a(j, i) = tEnd IfNe*t jNe*t izz = aElseReDim b(UBound(a, 2), UBound(a, 1)For i = 0 To UBound(a, 2)For j = 0 To UBound(a, 1)b(i, j) = a(j, i)Ne*t jNe*t izz = bEnd IfEnd Function兩矩陣相乘Public Function *c(a, b)Dim i As Integer, j As Integer, k As IntegerIf UBound(a, 2) UBound(b, 1) ThenMsgBo* (這兩個矩陣不能夠

12、相乘)E*it FunctionEnd IfReDim sd(UBound(a, 1), UBound(b, 2)For i = 0 To UBound(a, 1)For j = 0 To UBound(b, 2)For k = 0 To UBound(b, 1)sd(i, j) = sd(i, j) + a(i, k) * b(k, j)Ne*t kNe*t jNe*t i*c = sdEnd FunctionPublic Function qiuni(a)Dim c, m%, n%, p#, l%, i%, j%, ab#m = UBound(a, 1)n = UBound(a, 2)If

13、 m n ThenMsgBo* (該矩陣不可逆！)E*it FunctionEnd IfReDim c(m, 2 * n + 1)For i = 0 To m For j = 0 To n c(i, j) = a(i, j) Ne*t jNe*t iFor i = 0 To m For j = m + 1 To 2 * m + 1 c(i, j) = 0 Ne*t jNe*t ii = 0For j = m + 1 To 2 * m + 1 c(i, j) = 1 i = i + 1 Ne*t jFor k = 0 To n If c(k, k) = 0 Then For i = k + 1

14、To n If c(i, k) 0 Then GoTo this End If Ne*t i If i = n + 1 Then MsgBo* (該矩陣不可逆!) E*it Function End Ifthis: For j = 0 To 2 * m + 1 p = c(k, j) c(k, j) = c(i, j) c(i, j) = p Ne*t j End Ifab = 1# / c(k, k)For j = 0 To 2 * m + 1 c(k, j) = c(k, j) * abNe*t jFor i = 0 To n If i k Then For j = 0 To 2 * m + 1 If j k Then c(i, j) = c(i, j) - c(i, k) * c(k, j) End If Ne*t j c(i, k) = 0 End IfNe*t iN