用Matlab绘制地球与火箭轨道的三维图形

为了绘制地球和火箭轨道的三维图形，我们需要使用Matlab中的三维绘图工具。以下是一个简单的示例代码，可以帮助您开始：

% 定义常数
G = 6.67e-11; % 万有引力常数
M = 5.97e24; % 地球质量
m = 1000; % 火箭质量
r0 = 6.38e6; % 地球半径
h0 = 200e3; % 初始高度
v0 = 7.9e3; % 初始速度

% 计算初始位置和速度
r = r0 + h0;
v = v0;
theta = 0;

% 定义时间步长和总时间
dt = 10;
t = 0:dt:3600;

% 初始化位置和速度向量
r_vec = zeros(length(t), 3);
v_vec = zeros(length(t), 3);

% 计算位置和速度向量
r_vec(1, :) = [r*cosd(theta), r*sind(theta), 0];
v_vec(1, :) = [v*cosd(theta+90), v*sind(theta+90), 0];

for i = 2:length(t)
    % 计算火箭的加速度
    Fg = -G*M*m/(r^2);
    a = Fg/m;
    
    % 计算位置和速度
    r = r + v*dt;
    v = v + a*dt;
    theta = theta + v*dt/r;
    
    % 更新位置和速度向量
    r_vec(i, :) = [r*cosd(theta), r*sind(theta), 0];
    v_vec(i, :) = [v*cosd(theta+90), v*sind(theta+90), 0];
end

% 绘制地球
[x, y, z] = sphere(50);
r_earth = r0;
surf(x*r_earth, y*r_earth, z*r_earth)

hold on

% 绘制火箭轨迹
plot3(r_vec(:,1), r_vec(:,2), r_vec(:,3), 'r')

% 设置坐标轴和标题
axis equal
title('地球与火箭轨迹')
xlabel('X (m)')
ylabel('Y (m)')
zlabel('Z (m)')

hold off

在这个示例中，我们首先定义了常数和初始条件。然后我们使用一个循环来计算火箭的位置和速度，并将它们存储在位置和速度向量中。最后，我们使用Matlab的三维绘图工具来绘制地球和火箭轨迹。

运行代码后，您应该可以看到一个三维图形，其中地球是一个球体，火箭轨迹是一条红色的线。您可以通过旋转图形来查看它们的不同角度和方向。