你需要告诉我们你正在努力实现的微分方程,连同你的J2摄动方程,月亮和太阳。
运动方程的空间是女士在低地球轨道
21日视图(30天)
显示旧的评论
所以我的目标是模拟这项/运动的太空船从低地球轨道初始coditions emuating初始推力。
这里的代码,
这是对的吗?我知道下一步是实现龙格Kutt算法(解决ode),但是我不应该去状态空间表示法(知道ned的目标是实现一个Klaman过滤器),请帮助!
r =√X。Y ^ 2 +。Z ^ ^ 2 + 2);%计算航天器到原点的距离(可能是地球的中心或其他参考点)。
rm =√Xm。^ 2 +,^ 2 + Zm评选^ 2);%计算宇宙飞船到月球的距离。
delta_m =√(X-Xm)。^ 2 + (Y-Ym)。^ 2 + (Z-Zm) ^ 2);%计算宇宙飞船和月球之间的距离。
delta_s =√(X-Xs)。^ 2 +(等号)。^ 2 + (Z-Zs) ^ 2);%计算航天器和太阳之间的距离。
mu_e = 3.986135 e14灯头;%设置为地球重力参数。
mu_m = 4.89820 e12汽油;%设置为月球重力参数。
mu_s = 1.3253 e20;%设置为太阳重力参数。
一个= 6.37826 e6;%地球的赤道半径。
J = 1.6246 e - 3;%集地球的第二个动态形式因素。
(X_2dot =) - (mu_e * X / r ^ 3) * (1 + (J (a / r) ^ 2) * (1 - (5 * Z ^ 2 / r ^ 2)))——(mu_m (X-Xm) / delta_m ^ 3)——(mu_m * Xm / rm ^ 3)——(mu_s * (X-X_s) / delta_s ^ 3)——(mu_s * X / r ^ 3);
Y_2dot = - ((mu_e * Y) / r ^ 3) * (1 + (J (a / r) ^ 2) * (1 - (5 * Z ^ 2 / r ^ 2)))——(mu_m (Y-Ym) / delta_m ^ 3)——(mu_m * Ym / rm ^ 3)——(mu_s * (Y-Y_s) / delta_s ^ 3)——(mu_s * Y / r ^ 3);
(Z_2dot =) - (mu_e * Z / r ^ 3) * (1 + (J (a / r) ^ 2) * (1 - (5 * Z ^ 2 / r ^ 2)))——(mu_m (Z-Zm) / delta_m ^ 3)——(mu_m * Zm评选/ rm ^ 3)——(mu_s * (Z-Z_s) / delta_s ^ 3)——(mu_s * Z / r ^ 3);
答案(1)
詹姆斯Tursa
2023年3月31日
编辑:詹姆斯Tursa
2023年3月31日
比约恩指出,J2的条件都是一样的,但他们应该是不同的,因为赤道隆起影响z轴不同影响x - y轴。例如,你发布方程对Z_2dot相当于3 - (5 * Z ^ 2 / r ^ 2)但代码1 - (5 * Z ^ 2 / r ^ 2)。
另外,我建议您重写代码使用一个向量为你的国家而不是单独的变量。例如,定义状态向量如下:
y = x (1)
(2)= y
y (3) = z
y (4) = xdot
y (5) = ydot
y (6) = zdot
然后实现导数函数具有以下签名:
函数dy = myderiv (t、y,其他东西)
%其他代码
X_2dot =等等。
Y_2dot =等等。
Z_2dot =等等。
dy = [y (4:6); X_2dot Y_2dot; Z_2dot);
返回
结束
其他东西将重力常数,太阳和月亮的位置等。
下游代码将更容易实现,如果你有你的国家在这样一个向量。
似乎rm的评论应该读“地球到月球”而不是“月球飞船”
你失踪一个*用在这些方面:
mu_m * (X-Xm)
mu_m * (Y-Ym)
mu_m * (Z-Zm)
