matlab代码，已知航天器从轨道半径为7000km的轨道向轨道半径为8000km的轨道进行霍曼转移，要求绘制出航天器的轨道转移的运动轨迹图像

以下是一个简单的 MATLAB 代码，用于绘制航天器从7000km到8000km轨道的霍曼转移轨迹图像：

% 常数定义
mu = 398600; % 地球引力常数
r1 = 7000; % 初始轨道半径
r2 = 8000; % 目标轨道半径
a = (r1+r2)/2; % 中间椭圆轨道长半轴
v1 = sqrt(mu/r1); % 初始轨道速度
v2 = sqrt(mu/r2); % 目标轨道速度
v_inf = v2 - v1; % 转移时无限远处相对速度

% 计算转移半椭圆参数
e = r2/a - 1;
p = a*(1-e^2);
theta = linspace(0,pi,100);
r = p./(1+e*cos(theta));
x = r.*cos(theta);
y = r.*sin(theta);

% 计算转移时刻的速度和位置
delta_t = pi*sqrt(a^3/mu)/2; % 转移时间
delta_v = sqrt(mu/p)*e*sin(pi); % 转移时刻速度增量
v1_trans = [0, v1] + [-delta_v, v_inf];
v2_trans = [0, v2] + [delta_v, v_inf];
[r1_trans, ~] = rv_from_r0v0(r1, v1_trans, mu);
[r2_trans, ~] = rv_from_r0v0(r2, v2_trans, mu);

% 绘制图像
figure; hold on;
title('航天器的轨道转移的运动轨迹');
xlabel('x (km)'); ylabel('y (km)');
axis equal;
plot(0,0,'o','MarkerSize',5,'MarkerFaceColor','k');
plot(r1*cos(theta), r1*sin(theta),'r--');
plot(r2*cos(theta), r2*sin(theta),'b--');
plot(r.*cos(theta), r.*sin(theta),'k-');
plot(r1_trans(1), r1_trans(2),'ro');
plot(r2_trans(1), r2_trans(2),'bo');
legend('地球','初始轨道','目标轨道','转移轨道','转移时刻位置','到达时刻位置');

这个代码将绘制一个包含地球、初始轨道、目标轨道、转移轨道以及转移时刻和到达时刻的位置的图像。你可以通过修改代码中的常数来尝试不同的情况。