matlab小球碰撞

在MATLAB中，可以通过模拟小球碰撞来研究物体之间的碰撞行为。下面是一个简单的示例代码，用于模拟两个小球的碰撞过程：

% 定义小球的初始位置、速度和半径
pos1 = [0, 0];  % 小球1的初始位置
vel1 = [1, 0];  % 小球1的初始速度
radius1 = 1;    % 小球1的半径

pos2 = [5, 0];  % 小球2的初始位置
vel2 = [-1, 0]; % 小球2的初始速度
radius2 = 1;    % 小球2的半径

% 模拟碰撞过程
numSteps = 100; % 模拟的步数
dt = 0.01;      % 时间步长

for step = 1:numSteps
    % 计算小球之间的距离和相对速度
    distance = norm(pos2 - pos1);
    relativeVel = vel2 - vel1;
    
    % 如果小球之间有碰撞发生
    if distance <= radius1 + radius2
        % 计算碰撞后的速度
        normal = (pos2 - pos1) / distance; % 碰撞法线方向
        impulse = dot(relativeVel, normal) * normal; % 碰撞冲量
        vel1 = vel1 + impulse; % 更新小球1的速度
        vel2 = vel2 - impulse; % 更新小球2的速度
    end
    
    % 更新小球的位置
    pos1 = pos1 + vel1 * dt;
    pos2 = pos2 + vel2 * dt;
    
    % 绘制小球的位置
    plot(pos1(1), pos1(2), 'ro', 'MarkerSize', radius1*10);
    hold on;
    plot(pos2(1), pos2(2), 'bo', 'MarkerSize', radius2*10);
    xlim([-10, 10]);
    ylim([-10, 10]);
    axis equal;
    hold off;
    pause(0.01);
end

这段代码模拟了两个小球的碰撞过程。首先定义了两个小球的初始位置、速度和半径。然后通过循环模拟碰撞过程，计算小球之间的距离和相对速度，如果小球之间有碰撞发生，则根据碰撞法线方向和碰撞冲量更新小球的速度。最后更新小球的位置，并通过绘图函数将小球的位置可视化。