0

Space Exploration Beta 的好心人一直在帮助我计算轨道元素,但我一直坚持在 C# 中应用一个公式一段时间,需要一些建议。

该公式适用于偏心矢量,使用:

$$\vec e = \frac{\vec v \times \vec h}{\mu} \ - \frac{\vec r}{||r||} $$

(e = v * h /mu - r / ||r||)

我的公式代码如下;

double V = 7770;
double R = 6771000;
double Mu = 4e+14;

Vector3D velocity = new Vector3D(0, V, 0);
Vector3D Radias = new Vector3D(R, 0, 0);

Vector3D H = Vector3D.CrossProduct(Radias, velocity);//Angular Momentum = r * v
 //e = v * h /mu - r / ||r||
Vector3D rh = Vector3D.CrossProduct(velocity, H);
double t = Radias.Length;
Vector3D e = Vector3D.Divide(rh, Mu) - R / t;

现在这提出了

运算符“-”不能应用于“Vector3D”和“double”类型的操作数

现在我现在说不能将向量除以双倍,但我不确定如何解决这个问题,我想铸造但我不知道该怎么做?

这真的是我的数学技能和 C# 数学编程的边缘。

任何帮助将不胜感激。

编辑这是我正在使用的源材料之一,因此您可以正确查看公式。

编辑 2这是工作代码: Vector3D e = Vector3D.Divide(rh, Mu) - Radias / Radias.Length; double e1 = e.Length; Console.WriteLine(e1);

带来的答案不仅与用户通过 Space Exploration beta 提供的答案相匹配,而且与我尝试过的公式的不同版本相匹配。所以代码和公式都很好。感谢你的协助。

4

2 回答 2

1

r / ||r||是坐标的单位向量,即Radias / Radias.Length= (1, 0, 0)

Vector3D e = rh / Mu - Radias / Radias.Length

您只是将 ,R声明中的数字RadiasRadias它本身混淆了。

[请注意,我们可以简单地使用重载的除法运算/来代替Vector3D.Divide]

于 2018-03-05T09:43:14.160 回答
0

它并不是说你不能分裂。它说你不能减去。如果要为向量的每个分量减去一个双精度值,则需要用它创建一个向量。在您的情况下,问题是您在此处显示的最后一行:

Vector3D e = Vector3D.Divide(rh, Mu) - R / t;

在这一行中,Vector3D.Divide(rh, Mu)是一个向量并且 R / t是一个双精度数。为此,请创建一个向量,其中三个分量设置为R / t

Vector3D rt = new Vector3D(R / t, R / t, R / t);
Vector3D e = Vector3D.Divide(rh, Mu) - rt;

(假设您想从向量的所有分量中减去 R/t)

编辑:再看看你的公式,你确定那R不应该是一个向量吗?似乎它是基本公式中的向量,不在您的代码中

于 2018-03-05T09:37:12.767 回答