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，地月距离。