模拟轨道中卫星间的相对运动
Clohessy–Wiltshire 方程模拟围绕中央天体轨道运动的两颗卫星间的相对运动. 这是一颗卫星相对于另一颗的笛卡尔位置的线性模型. 这里进行两个模拟,每个模拟使用不同的初始条件.
sat = StateSpaceModel[{(| | | | | | |
| ----- | ---- | - | --- | ---- | - |
| 0 | 1 | 0 | 0 | 0 | 0 |
| 3 n^2 | 0 | 0 | 2 n | 0 | 0 |
| 0 | 0 | 0 | 1 | 0 | 0 |
| 0 | -2 n | 0 | 0 | 0 | 0 |
| 0 | 0 | 0 | 0 | 0 | 1 |
| 0 | 0 | 0 | 0 | -n^2 | 0 |), (| | | |
| - | - | - |
| 0 | 0 | 0 |
| 1 | 0 | 0 |
| 0 | 0 | 0 |
| 0 | 1 | 0 |
| 0 | 0 | 0 |
| 0 | 0 | 1 |), (| | | | | | |
| - | - | - | - | - | - |
| 1 | 0 | 0 | 0 | 0 | 0 |
| 0 | 1 | 0 | 0 | 0 | 0 |
| 0 | 0 | 1 | 0 | 0 | 0 |)}];使用 StateResponse 从一组特定的发射条件中获取封闭的相对轨道:
n = 0.001;Subscript[x, 0] = 1;Subscript[xd, 0] = 0.0005;
x = StateResponse[{sat, {Subscript[x, 0], Subscript[xd, 0], (2Subscript[xd, 0]/n), -2n Subscript[x, 0], Sequence@@##}}, {0, 0, 0}, {t, 0, 5 (2Pi/n)}]& /@ {{-2 Subscript[x, 0], -2 Subscript[xd, 0]}, {2 Subscript[x, 0], 2 Subscript[xd, 0]}};ParametricPlot3D[{x[[1, 1 ;; 5 ;; 2]], x[[2, 1 ;; 5 ;; 2]]}, {t, 0, 5 (2Pi/n)}, BoxRatios -> {1, 1, 1}, AxesLabel -> {"X", "Y", "Z"}]