逆向轨道与原顺行轨道根数关系

---

研究目的与意义

当前空间对抗需求增加，逆向运行的轨道可以带来诸多好处：对于侦查卫星，逆行可以使其在极短时间内飞掠多个侦查目标，对于同轨道仅存在相位差的卫星编队侦查任务尤为适用；对于拦截任务卫星，作为武器以摧毁敌方目标为目的的情况下，巨大的速度差可以使小质量逆行卫星带来对敌毁灭性打击的效果。

---

根数关系

笔者在STK中接触了逆行轨道设计与参数优化相关的研究，当此次进入MATLAB中写代码时，不再有直观的三维信息，面对有偏心率的非圆轨道，也无法轻易利用头脑构建出各根数之间的关系。

因此通过仿真，得出了适用于各类轨道的原轨道与逆向轨道六根数关系：

显然，若出现超出取值范围的情况，自行进行调整即可。

有时在设计案例中，需要给定交会点位置，这涉及到最后一个轨道根数。若是在逆行轨道需要给出原位置，则真近点角需要做如下处理

以上即为逆行轨道与原轨道根数的所有关系。

仿真算例

利用MATLAB内置的随机数函数给出一组目标轨道根数，当然，在随机过程中除了取值范围外唯一需要考虑的约束就是近地点半径约束。利用上方的关系式得出逆行轨道根数，再利用数值积分画出轨道、标注出交会点，通过画图观察仿真结果。

代码如下所示：

miu = 3.986*1e14;
R_Earth = 6378.1*1e3;
pai = 3.1415926535;

% 目标轨道
a = R_Earth+5000*1e3*rand;
e = rand;
while a*(1-e) < R_Earth      % 近地点半径约束
    e = rand;
end
i = pai*rand;
w = pai*rand;
W = pai*rand;
f = 2*pai*rand;

% 逆向轨道
a2 = a;
e2 = e;
i2 = pai-i;
w2 = pai-w;
W2 = pai+W;
f2 = 2*pai-f;

% 计算
data = [a,e,i,w,W,f];
data2 = [a2,e2,i2,w2,W2,f2];

RV = Classic_Coordinate(data);
RV2 = Classic_Coordinate(data2);

T = 2*pai/sqrt(miu)*a^1.5;
tspan = 1:1:round(T);
[t,x] = ode45('dynamic',tspan,RV);
[t2,x2] = ode45('dynamic',tspan,RV2);

% 作图
figure(1)
plot3(RV(1),RV(2),RV(3),'og','LineWidth',2);
hold on;
plot3(RV2(1),RV2(2),RV2(3),'ok','LineWidth',2);
hold on;
plot3(x(:,1),x(:,2),x(:,3),'--r','Linewidth',2);
hold on;
plot3(x2(:,1),x2(:,2),x2(:,3),'-b','Linewidth',2);
hold on;
[Xe,Ye,Ze]=sphere(24);%生成24节点球函数。需乘以地球半径
Xe=R_Earth*Xe;
Ye=R_Earth*Ye;
Ze=R_Earth*Ze;
surf(Xe,Ye,Ze)%画出点构成的表面；
axis equal;
grid on;
legend('交会点1','交会点2','顺行轨道','逆行轨道');
title('轨迹与交会点');

function RV = Classic_Coordinate(data)

miu = 3.986*1e14;

a=data(1);
e=data(2);
i=data(3);
w=data(4);
W=data(5);
f=data(6);
p=abs(a*(1-e^2));
u=w+f;
   
Coordinate=p/(1+e*cos(f))*[cos(W)*cos(u)-sin(W)*sin(u)*...
       cos(i)  sin(W)*cos(u)+cos(W)*sin(u)*cos(i) sin(i)*sin(u)]';
V=(miu/p)^(0.5)*[-cos(W)*(sin(u)+e*sin(w))-sin(W)*(cos(u)...
       +e*cos(w))*cos(i) -sin(W)*(sin(u)+e*sin(w))+cos(W)*(cos(u)+...
       e*cos(w))*cos(i) sin(i)*(cos(u)+e*cos(w))]';

RV = [Coordinate;V];

end

% function dx = dynamic(t,x)
% 
% miu=3.986*1e14;
% 
% rx=x(1);
% ry=x(2);
% rz=x(3);
% vx=x(4);
% vy=x(5);
% vz=x(6);
% 
% r=sqrt(rx^2+ry^2+rz^2);
% 
% drx=vx;
% dry=vy;
% drz=vz;
% dvx=-miu*rx/r^3;
% dvy=-miu*ry/r^3;
% dvz=-miu*rz/r^3;
% 
% dx=[drx;dry;drz;dvx;dvy;dvz];
% 
% end

注：最后被注释掉的为二体动力学微分方程。

通过运行代码，可以得出如下图像：

 当然，也有其他更高可信度的验证方式，此处采取作图方法只是为了更加直观。