clear all
syms aa bb
R=5;i=7; ? ? ? ? ? ? ? ? %設(shè)置小球半徑以及小球數(shù)量
X=randi([10 90],1,i); ? ? ? %在[10，90]間取整數(shù)作為每個(gè)小球圓心的X坐標(biāo)
Y=randi([10 90],1,i); ? ? ? %在[10，90]間取整數(shù)作為每個(gè)小球圓心的Y坐標(biāo)
t=0:pi/100:2*pi;
x=2.*rand(1,i)+1; ? ? ? ? ?%在[1，3]間取隨機(jī)數(shù)作為每個(gè)小球x軸方向速度
y=2.*rand(1,i)+1; ? ? ? ? ?%在[1，3]間取隨機(jī)數(shù)作為每個(gè)小球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畫出每個(gè)小球
? ? 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) ? ? ? ? ? ? ? ? ? ? ? ? %將存好的畫面以更快的幀數(shù)播放一遍

總結(jié)

以上是生活随笔為你收集整理的MATLAB——小球碰撞的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。