高中/碰撞

Posted on

MATLAB模拟碰撞

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
m=[3,1,2];%小球质量
ball=[5 5;20 10;11 15];%球的坐标
r=[3,4,5];%半径
color=['b','y','r'];%球的颜色
v= rand(3,2)*0.8;%随机生成速度
disp(v)
dt=0.55;
figure
xlim([-5 45])
ylim([-5 35])
while true%逐帧绘制
    yes=[1,1,1];%未处理,防止重复计算
    for i=1:3
        for j=1:3
            if yes(i)
                if i==j
                    if ball(j,1)<r(j) || ball(j,1)>=40-r(j)%x边界
                        v(j,1)=-v(j,1);
                    end
                    if ball(j,2)<r(j) || ball(j,2)>=30-r(j)%y边界
                        v(j,2)=-v(j,2);
                    end
                else
                    distance=sqrt((ball(i,1)-ball(j,1)).^2+(ball(i,2)-ball(j,2)).^2);
                    if distance<=(r(i)+r(j))%速度计算
                        vix=((m(i)-m(j))*v(i,1)+2*m(j)*v(j,1))/(m(i)+m(j));
                        viy=((m(i)-m(j))*v(i,2)+2*m(j)*v(j,2))/(m(i)+m(j));
                        vjx=((m(j)-m(i))*v(j,1)+2*m(i)*v(i,1))/(m(i)+m(j));
                        vjy=((m(j)-m(i))*v(j,2)+2*m(i)*v(i,2))/(m(i)+m(j));
                        v(i,1)=vix;
                        v(i,2)=viy;
                        v(j,1)=vjx;
                        v(j,2)=vjy;
                        yes(j)=0;
                    end    
                end
            end
        end 
    end
    rectangle('Position',[0,0,40,30])
    for i=1:3
        ball(i,:)=ball(i,:)+v(i,:)*dt;%近似匀速
        rectangle('Position',[ball(i,1)-r(i),ball(i,2)-r(i),2*r(i),2*r(i)],'Curvature',[1,1],  'FaceColor',color(i))
    end
    drawnow
    cla%清空
end
Shiwei Pan wechat