MATLAB——小球碰撞

clear all
syms aa bb
R=5;i=7;                 %设置小球半径以及小球数量
X=randi([10 90],1,i);       %在[10，90]间取整数作为每个小球圆心的X坐标
Y=randi([10 90],1,i);       %在[10，90]间取整数作为每个小球圆心的Y坐标
t=0:pi/100:2*pi;
x=2.*rand(1,i)+1;          %在[1，3]间取随机数作为每个小球x轴方向速度
y=2.*rand(1,i)+1;          %在[1，3]间取随机数作为每个小球y轴方向速度
for n=1:10000
    for v=1:i
        if 100-X(v)<R          %小球与右器壁相碰的速度判定
            x(v)=-abs(x(v));
        end
        if X(v)<R              %小球与左器壁相碰的速度判定
            x(v)=abs(x(v));
        end
        if Y(v)<R              %小球与下器壁相碰的速度判定
            y(v)=abs(y(v));
        end
        if 100-Y(v)<R          %小球与上器壁相碰的速度判定
            y(v)=-abs(y(v));
        end
        for b=v:i
            if sqrt((X(v)-X(b))^2+(Y(v)-Y(b))^2)<=2*R   %小球间相碰速度交换
                aa=x(v);bb=y(v);
                x(v)=x(b);y(v)=y(b);x(b)=aa;y(b)=bb;
            end
        end
        X(v)=X(v)+x(v);Y(v)=Y(v)+y(v);            %小球在X与Y上的位置变化
    end
    x1=X(1)+R*cos(t);y1=Y(1)+R*sin(t);            %在圆心处以半径R画出每个小球
    x2=X(2)+R*cos(t);y2=Y(2)+R*sin(t);
    x3=X(3)+R*cos(t);y3=Y(3)+R*sin(t);
    x4=X(4)+R*cos(t);y4=Y(4)+R*sin(t);
    x5=X(5)+R*cos(t);y5=Y(5)+R*sin(t);
    x6=X(6)+R*cos(t);y6=Y(6)+R*sin(t);
    x7=X(7)+R*cos(t);y7=Y(7)+R*sin(t);
    fill(x1,y1,'r',x2,y2,'b',x3,y3,'g',x4,y4,'y',x5,y5,'m',x6,y6,'w',x7,y7,'k')
    axis equal                         %定义容器为100×100的正方形
    axis([0 100 0 100])
    M(n)=getframe;                    %将每帧画面存入M
end
movie(M,1,48)                         %将存好的画面以更快的帧数播放一遍