写一个matlab模拟小球平抛运动过程的代码

在MATLAB中，你可以使用矢量运算和图形绘制功能来模拟小球的平抛运动。下面是一个基本的示例代码，它假设你知道初始位置、初速度、重力加速度以及时间间隔：

% 定义变量
g = 9.8; % 重力加速度，默认值约为9.8 m/s^2
v0 = [0; initial_velocity]; % 初始速度，垂直方向向下，水平方向为0
theta = pi / 4; % 抛射角度
x0 = 0; % 初始水平位置
y0 = 0; % 初始竖直位置
dt = 0.01; % 时间步长

tspan = 0:dt:10; % 计算时间范围直到10秒
numFrames = length(tspan); % 总帧数

% 初始化位置矩阵
[x, y] = zeros(numFrames, 2);

% 模拟平抛运动
for i = 1:numFrames
    % 计算当前时间和位置
    t = tspan(i);
    [x(i), y(i)] = projectile_motion(x0, y0, v0, g, t, theta);
    
    % 更新图形
    if i == 1 || mod(i, round(numFrames/5)) == 0 % 每隔一段时间更新一次图形
        plot(x(1:i), y(1:i), 'ro'); % 用红点表示小球位置
        drawnow;
    end
end

function [px, py] = projectile_motion(x0, y0, v0, g, t, theta)
    % 平抛运动方程
    px = x0 + (v0 * cos(theta)) * t;
    py = y0 - (0.5 * g * t.^2); % 在忽略空气阻力的情况下
end

运行此代码，你会看到一个小球从起点开始沿平抛轨迹移动。每次达到整数倍的时间间隔时，会暂停并显示当前位置。