用matlab实现碰撞模型程序代码

用m a t l a b实现碰撞模型程序代码

标准化工作室编码[XX968T-XX89628-XJ668-XT689N]

c l c;

clear;

fill([6,7,7,6],[5,5,0,0],[0,0.5,0]);%右边竖条的填充

holdon;%保持当前图形及轴系的所有特性

fill([2,6,6,2],[3,3,0,0],[0,0.5,0]);%左边竖条的填充

holdon;%保持当前图形及轴系的所有特性

t1=0:pi/60:pi;

plot(4-2*sin(t1-pi/2),5-2*cos(t1-pi/2));%绘制中间的凹弧图形gridon;%添加网格线

axis([0,9,0,9]);%定义坐标轴的比例%

axis('off');%关闭所有轴标注，标记，背景

fill([1,2,2,1],[5,5,0,0],[0,0.5,0]);%中间长方形的填充

holdon;%保持当前图形及轴系的所有特性

title('碰撞');%定义图题

x0=6;

y0=5;

head1=line(x0,y0,'color','r','linestyle','.','erasemode','xor','marke rsize',30);

head2=line(x0,y0,'color','r','linestyle','.','erasemode','xor','marke rsize',50);%设置小球颜色,大小,线条的擦拭方式

t=0;%设置小球的初始值

dt=0.001;%设置运动周期

t1=0;%设置大球的初始值

dt1=0.001;

while1%条件表达式

t=t+dt;

x1=9-1*t;

y1=5;

x3=6;

y3=5;

ift>0

x2=6;

y2=5;%设置小球的运动轨迹

end

ift>2.8

t=t+dt;

a=sin(t-3);

x1=6.1;

y1=5.1;

x3=4-2*sin(1.5*a);

y3=5-2*cos(1.5*a);%设置大球的运动轨迹

end

set(head1,'xdata',x1,'ydata',y1);%设置球的运动set(head2,'xdata',x3,'ydata',y3);

drawnow;

end