炮弹弹道轨迹仿真：参数设置与代码优化

% 参数设置
T = 300; % 运动时间  
Tt = 0.05; % 仿真时间间隔
P = [0,0]; % 初始发射位置
M = 20; % 炮弹质量
C = 0.45; % 空气阻力系数
rou = 1.29; % 空气密度
D = 0.12; % 炮弹直径
S = pi*D^2/4; % 炮弹迎风面积
k = 0.5*C*rou*S/M; % 空气阻力加速度总系数
g = 9.8; % 重力加速度

% 已知落点距离
dist = 1670;

% 初始化变量
foundMatch = false;
lastTraj = [];
velocities = [];

% 遍历发射角度  
for angle = 0:0.1:20

    V = 100*sqrt(2)*[cos(angle*pi/180), sin(angle*pi/180)]; % 计算飞行速度
    P = [0,0];

    % 存储轨迹点
    traj = []; 
    
    % 存储速度变化
    velocities = [];
    
    for t = 0:Tt:T
        if P(2) < 0
            break;
        end

        af = -k*[V(1)^2, V(2)^2];
        a = af + [0, -g];

        V = V + a*Tt;
        P = P + V*Tt;
        
        % 存储轨迹点
        traj = [traj; P]; 
        
        % 存储速度变化
        velocities = [velocities; norm(V)];
    end

    if abs(P(1) - dist) < 1
        foundMatch = true;
        break
    end
    
    lastTraj = traj;
end  

% 如果未找到匹配的角度，绘制最后一个角度下的轨迹
if ~foundMatch
    lastAngle = angle;
    traj = lastTraj;
end

% 绘制轨迹 
figure;
plot(traj(:,1), traj(:,2), '*');
axis equal; % 设置坐标轴间距相同
title('弹道轨迹');
xlabel('距离/m');
ylabel('高度/m');

% 绘制速度变化曲线
figure;
time = 0:Tt:T;
plot(time(1:length(velocities)), velocities);
title('炮弹速度变化');
xlabel('时间/s');
ylabel('速度/m/s');

% 打印结果
fprintf('落点距离:%d米;
误差:%d米;
对应发射角度:%d度
.', P(1),abs(P(1)-dist),lastAngle);

这里我们在找不到匹配角度时，将轨迹数据保存在 lastTraj 变量中，然后在绘制轨迹和速度变化曲线之前再次使用这些数据来绘制最后一个角度下的轨迹。

请运行修改后的代码，检查是否可以正确生成轨迹和速度变化曲线图。如果您还有其他问题，请随时提问。