多个小球碰撞有背景颜色支持更换小球个数项目 - matlab其它源码

下载频道> 资源分类> matlab源码> 其它源码> 多个小球碰撞有背景颜色支持更换小球个数项目

标题：多个小球碰撞有背景颜色支持更换小球个数项目
分享到：

所属分类: 其它源码资源类型：程序源码文件大小: 1.5 KB 上传时间: 2019-07-16 22:24:13 下载次数: 2 资源积分：1分提供者: jiqiren 20190716102448802

内容：
function[]=Billiards(NumberOfBalls)

% Call: Billiards(Number of Balls), e.g.:

% Billiards(10)

if nargin == 0

NumberOfBalls = 10;%规定球的个数重要

end

close all;

hold on;

drawflag=1;

factor=82; %Adjust this factor when scaling 重要

DT=5e-3;

Bound=[-4 4 -2 2];

BallColour=[[1 0 0];[1 0 0.5];[1 0.5 0];[0 1 0];[0 0 1];[1 1 0];[1 1 1];...

[0 0.3 0];[0 0 0];[0.65 0.65 0.65];[0 0.75 0.75];[0.3 0 0.6];[0.95 0.65 0.75];...

[0.5 0.25 0];[0 0.2 0.4];[0.9 0.4 0.7];[0.4 0.2 0.3];[0.65 0.55 0.15];[0.25 0.35 0.25];[0.5 0 0]];

TableColour=[.4 .5 .8];

%Plothandle

axis(Bound);

set(gca,'Color',TableColour,'xcolor',TableColour,'ycolor',...

TableColour,'PlotBoxAspectRatio',[1 abs((Bound(3)-Bound(4))/(Bound(2)-Bound(1))) 1],'xtick',[],'ytick',[])

%===============================================================

r=0.15+0.25*rand(NumberOfBalls,1);

%Radii==========================================================

for count1=1:NumberOfBalls;

for count2=1:count1,

rmatrix(count2,count1)=r(count1)+r(count2);

end;

rmatrix=rmatrix.*triu(abs(-1+eye(size(rmatrix))));

%Mass ==========================================================

massa=pi*r.^2;

X=[(Bound(2)-Bound(1)-2*max(r))*rand(NumberOfBalls,1)+Bound(1)+max(r),...

(Bound(4)-Bound(3)-2*max(r))*rand(NumberOfBalls,1)+Bound(3)+max(r)];

for j=1:NumberOfBalls;

for i=1:j;

distmatrix(i,j)=sqrt((X(j,1)-X(i,1))^2+(X(j,2)-X(i,2))^2);

end;

%Initial Edgedetectionmatrix====================================

Botsmatrix=(distmatrix-rmatrix)+tril(abs(-1+eye(size(distmatrix))))+eye(size(distmatrix));

while find(Botsmatrix<=0);

X=[(Bound(2)-Bound(1)-2*max(r))*rand(NumberOfBalls,1)+Bound(1)+max(r),...

(Bound(4)-Bound(3)-2*max(r))*rand(NumberOfBalls,1)+Bound(3)+max(r)];

for j=1:NumberOfBalls;

for i=1:j;

distmatrix(i,j)=sqrt((X(j,1)-X(i,1))^2+(X(j,2)-X(i,2))^2);

end;

Botsmatrix=(distmatrix-rmatrix)+tril(abs(-1+eye(size(distmatrix))))+eye(size(distmatrix));

end

%Initial velocities

V=5*(-1+2*rand(NumberOfBalls,2));

%Plot startingpositions

for k=1:NumberOfBalls;

h(k) = plot(X(k,1),X(k,2),'o', 'MarkerEdgeColor',BallColour(mod(k-1,length(BallColour))+1,:),...

'MarkerFaceColor',BallColour(mod(k-1,length(BallColour))+1,:),'MarkerSize',factor*r(k));

end

pause(1e-1);

%Loop

while drawflag==1;

for k=1:NumberOfBalls;

delete(h(k));

end

%Edgedetecton positive======================================

for j=1:NumberOfBalls;

for i=1:2;

d=X(j,i)+r(j)-Bound(2*i);

if d>=0

dt=d/V(j,i);

X(j,i)=X(j,i)-V(j,i)*dt;

V(j,i)=-V(j,i);

end

%Edgedetecton negative======================================

for j=1:NumberOfBalls;

for i=1:2;

d=X(j,i)-r(j)-Bound(2*i-1);

if d<=0

dt=d/V(j,i);

X(j,i)=X(j,i)-V(j,i)*dt;

V(j,i)=-V(j,i);

end

%Distancematrix============================================

for j=1:NumberOfBalls;

for i=1:j;

distmatrix(i,j)=sqrt((X(j,1)-X(i,1))^2+(X(j,2)-X(i,2))^2);

end;

%Collisiondetectionmatrix=========================================================

Botsmatrix=(distmatrix-rmatrix)+tril(abs(-1+eye(size(distmatrix))))+eye(size(distmatrix));

%===========================================================================

if find(Botsmatrix<0);

[I,J]=find(Botsmatrix<0);

for i=1:length(I)

normdist=normr([X(I(i),1)-X(J(i),1) X(I(i),2)-X(J(i),2)]);

vaA=(V(I(i),1)*normdist(1)+V(I(i),2)*normdist(2));

vaB=(V(J(i),1)*normdist(1)+V(J(i),2)*normdist(2));

dt=abs(r(I(i))+r(J(i))-distmatrix(I(i),J(i)))/(abs(vaA)+abs(vaB));

X(I(i):J(i),:)=X(I(i):J(i),:)-V(I(i):J(i),:)*dt;

end

for i=1:length(I)

normdist=normr([X(I(i),1)-X(J(i),1) X(I(i),2)-X(J(i),2)]);

vaA=(V(I(i),1)*normdist(1)+V(I(i),2)*normdist(2));

vbA=(-V(I(i),1)*normdist(2)+V(I(i),2)*normdist(1));

vaB=(V(J(i),1)*normdist(1)+V(J(i),2)*normdist(2));

vbB=(-V(J(i),1)*normdist(2)+V(J(i),2)*normdist(1));

va_A=vaA+2*(vaB-vaA)/(1+massa(I(i))/massa(J(i)));

va_B=vaB+2*(vaA-vaB)/(1+massa(J(i))/massa(I(i)));

V(I(i),1)=va_A*normdist(1)-vbA*normdist(2);

V(I(i),2)=va_A*normdist(2)+vbA*normdist(1);

V(J(i),1)=va_B*normdist(1)-vbB*normdist(2);

V(J(i),2)=va_B*normdist(2)+vbB*normdist(1);

X(I(i):J(i),:)=X(I(i):J(i),:)+V(I(i):J(i),:)*dt;

end

%Propagation=================================================

X=X+V*DT;

%Plotting====================================================

for k=1:NumberOfBalls;

h(k) = plot(X(k,1),X(k,2),'o', 'MarkerEdgeColor',BallColour(mod(k-1,length(BallColour))+1,:),...

'MarkerFaceColor',BallColour(mod(k-1,length(BallColour))+1,:),'MarkerSize',factor*r(k));

end

drawnow;

end

文件列表（点击上边下载按钮，如果是垃圾文件请在下面评价差评或者投诉）:

Billiards.m

关键词：多个小球 碰撞 个数 项目

matlab源码下载排行

回到顶部

友情链接

matlab源码网站源码编程语言网络编程手机通讯编程 KUNAG 官网 KUNAG 电商更多

: 联系方式| 版权声明| 招聘信息| 广告服务| 银行汇款| 法律顾问| 兼职技术| 付款方式| 关于我们|; 网站客服程序员兼职招聘 411241577@qq.com; 沪ICP备19040327号-3; 公安备案号：沪公网安备 31011802003874号; 库纳格流体控制系统（上海）有限公司版权所有; Copyright © 1999-2014, GUSUCODE.COM, All Rights Reserved

关键词： 多个小球 碰撞 个数 项目

相关推荐

matlab源码下载排行

友情链接

关键词：多个小球碰撞个数项目