Private Sub CommandButton2_Click()
x = 1000
y = 1000
z = 1000
r1 = 0
l = 0
Jd = Sheet1.Cells(8, 5)
100 If l <= 1000000 Then
r = Sqr(x ^ 2 + y ^ 2 + z ^ 2)
x1 = Sheet1.Cells(1, 6) * r + Sheet1.Cells(8, 1)
y1 = Sheet1.Cells(2, 6) * r + Sheet1.Cells(8, 2)
z1 = Sheet1.Cells(3, 6) * r + Sheet1.Cells(8, 3)
r1 = Sqr(x1 ^ 2 + y1 ^ 2 + z1 ^ 2)
If (Abs((r - r1) / r) > Jd Or Abs((x - x1) / x) > Jd Or Abs((y - y1) / y) > Jd Or Abs((z - z1) / z) > Jd) Then
x = x1
y = y1
z = z1
r = r1
l = l + 1
GoTo 100
End If
End If
Sheet1.Cells(10, 1) = r
Sheet1.Cells(10, 2) = x
Sheet1.Cells(10, 3) = y
Sheet1.Cells(10, 4) = z
Sheet1.Cells(10, 5) = l
End Sub
输入观测点ECEF坐标数据 ：
a a a a x y z
-1724.234， 2869.607451， 5420.518 a
输入目标点的ECEF坐标增量/r数据：
d1,d2,d3
a a 0.437128, 0.228882, 0.869789 a
计算结果：r=38.4wKm，地月距离。
转载本文请联系原作者获取授权，同时请注明本文来自刘志明科学网博客。
链接地址：https://wap.sciencenet.cn/blog-3556836-1499153.html?mobile=1
收藏
