写一段matlab代码，巡航导弹打击舰船目标，要包括导弹的各个飞行阶段，并画图轨迹图

以下是一个简单的巡航导弹打击舰船目标的 MATLAB 代码，其中包括导弹的各个飞行阶段和绘制轨迹图：

% 巡航导弹打击舰船目标 MATLAB 代码

% 初始化导弹和目标的位置和速度
missile_pos = [0; 0; 0]; % 导弹的初始位置
missile_vel = [200; 0; 0]; % 导弹的初始速度
target_pos = [5000; 500; 0]; % 目标的位置
target_vel = [0; 0; 0]; % 目标的速度

% 设定模拟时间和时间步长
T = 60; % 模拟时间60秒
dt = 0.1; % 时间步长0.1秒
t = 0:dt:T; % 时间向量

% 初始化导弹的加速度和速度向量
missile_acc = zeros(3, length(t));
missile_vel_vec = zeros(3, length(t));
missile_vel_vec(:, 1) = missile_vel;

% 计算导弹的加速度向量
for i = 1:length(t)-1
    % 计算导弹与目标的距离和方向向量
    dist = target_pos - missile_pos;
    dist_norm = norm(dist);
    dir_vec = dist ./ dist_norm;
    
    % 计算导弹的加速度
    missile_acc(:, i) = 1000 * dir_vec - 0.1 * missile_vel_vec(:, i);
    
    % 计算导弹的速度和位置
    missile_vel_vec(:, i+1) = missile_vel_vec(:, i) + dt * missile_acc(:, i);
    missile_pos = missile_pos + dt * missile_vel_vec(:, i);
    
    % 判断导弹是否已经击中目标
    if dist_norm < 50
        fprintf('导弹已经击中目标！\n');
        break;
    end
end

% 绘制导弹和目标的轨迹图
figure;
plot3(missile_pos(1,:), missile_pos(2,:), missile_pos(3,:), 'r-', 'LineWidth', 2);
hold on;
plot3(target_pos(1,:), target_pos(2,:), target_pos(3,:), 'g-', 'LineWidth', 2);
grid on;
xlabel('X');
ylabel('Y');
zlabel('Z');
legend('导弹', '目标');
title('导弹和目标的轨迹图');

这段代码模拟了一个巡航导弹打击舰船目标的过程，其中导弹在飞行过程中会进行加速度计算，然后根据加速度来更新速度和位置。当导弹与目标的距离小于50时，程序会输出“导弹已经击中目标！”并退出循环。

最后，程序会绘制导弹和目标的轨迹图，可以清晰地看到导弹在飞行过程中如何追踪目标并最终击中目标。可以使用 `plot3` 函数绘制三维图形，并在图例中显示导弹和目标的轨迹。