三维航迹动态绘制plot3

clc
figure;

% 飞机1的初始位置和速度
x1 = 0;
y1 = 0;
z1 = 0;
vx1 = 1; % x 方向速度
vy1 = 0.5; % y 方向速度
vz1 = 0.2; % z 方向速度

% 飞机2的初始位置和速度
x2 = 0;
y2 = 0;
z2 = 0;
vx2 = 0.7; % x 方向速度
vy2 = -0.3; % y 方向速度
vz2 = -0.1; % z 方向速度

% 创建一个空的图形对象
% axis ([0 1000 -500 500 -100 200]);
h = plot3(x1, y1, z1, 'b-', x2, y2, z2, 'r-');
grid on;
axis ([0 1000 -500 500 -100 200]);
xlabel('X');
ylabel('Y');
zlabel('Z');
title('飞机运动轨迹');
legend('飞机1', '飞机2');
view(-30, 20);

% 模拟飞机的动态运动
for t = 1:1000 % 循环次数可以根据需要调整
    % 更新飞机1的位置
    x1 = x1 + vx1;
    y1 = y1 + vy1;
    z1 = z1 + vz1;
    
    % 更新飞机2的位置
    x2 = x2 + vx2;
    y2 = y2 + vy2;
    z2 = z2 + vz2;
    
    % 更新图形对象的数据
    set(h(1), 'XData', [get(h(1), 'XData'), x1], 'YData', [get(h(1), 'YData'), y1], 'ZData', [get(h(1), 'ZData'), z1]);
    set(h(2), 'XData', [get(h(2), 'XData'), x2], 'YData', [get(h(2), 'YData'), y2], 'ZData', [get(h(2), 'ZData'), z2]);
    
    % 设置绘图的延迟，可根据需要调整速度
    pause(0.01);
end