gpsvbkf.txt

GPS卫星坐标计算与RENIX数据录入,只是没有进行地面坐标的计算

'''''''''''''GPS卫星位置坐标计算
v = B ^ 2
G = 3.9686005E+15
N1 = Sqr(G / (v ^ 3))
ng = N1 + nx
tk = tx + A0 + a1 * (tx - toc) + a2 * ((tx - toc) ^ 2) - toe
Mk = ng * tk + m1

e1 = Mk
e2 = Mk + e * Sin(e1)
Do Until Abs(e2 - e1) < 0.000001
e1 = e2
e2 = Mk + e * Sin(e1)
Ek = e2
Loop          ''''''''''''''迭代


fk = Atn((Sqr(1 - e ^ 2) * Sin(Ek)) / (Cos(Ek) - e))
uk = fk + w
ux = cuc * Cos(2 * uk) + cus * Sin(2 * uk)
rx = crc * Cos(2 * uk) + crs * Sin(2 * uk)
ix = cic * Cos(2 * uk) + cis * Sin(2 * uk)
ukg = uk + ux
rkg = v * (1 - e * Cos(Ek)) + rx
ikg = i1 + ix + ixx * tk
we = 0.0000729211567
Ωk = Ω + (Ωx - we) * tk - we * toe
X = rkg * Cos(ukg)
Y = rkg * Sin(ukg)
xk = X * Cos(Ωk) - Y * Cos(ikg) * Sin(Ωk)
yk = X * Sin(Ωk) + Y * Cos(ikg) * Cos(Ωk)
zk = Y * Sin(ikg)

Private Sub DYLT_Click()
CommonDialog1.Filter = "text(*.txt)|*.txt|Rich text files(*.rtf)|*.rtf|All files(*.*)|*.*"
CommonDialog1.FilterIndex = 3
CommonDialog1.CancelError = ture
On Error Resume Next
CommonDialog1.Action = 1
If CommonDialog1.filename = "" Then Exit Sub
Dim inputdata As String
i = 0: k = 0
'把文件中的数据付给数组Ａ
Open CommonDialog1.filename For Input As #1
Do While Not EOF(1)
Line Input #1, inputdata
k = k + 1
Loop
Close #1
ZS = Int(k / 8)
Open CommonDialog1.filename For Input As #1
ReDim a(k - 1) As String
Do While Not EOF(1)
Line Input #1, inputdata
a(i) = inputdata
i = i + 1
Loop
Close #1
For i = 1 To 5
For j = 0 To 2
a(i) = LTrim(a(i))
j = j + 4 * (i - 1)
If j < 17 Then
a1 = InStr((a(i)), " ")
Text1(j).Text = Left(a(i), a1)
n = Len(a(i))
m = n - a1
a(i) = Right(a(i), m)
j = j - 4 * (i - 1)
End If
Next j
If j < 17 Then
j = j + 4 * (i - 1)
Text1(j).Text = a(i)
End If
Next i
End Sub