数学建模社区-数学中国

标题: 求分析个小程序 [打印本页]

作者: ahy 时间: 2012-9-24 17:46
标题: 求分析个小程序
这是参考别人的程序改动而来，但是总是在小球2运动回初始位置时停止，求指教。
另外我的程序里若有不规范之处，请指出。。
x1=[7.5 7.5 -0.5 -0.5];y1=[4 0 0 4];
x2=0:0.01:7;y2=x2.*(x2-7)/4+4;
x=[x1 x2];y=[y1 y2];%带凹槽的平台边界坐标
fill(x,y,[0,0.5,0]);%对平台补色
axis([-1 10 0 5]);
hold on;
axis off;
title('动量与势能转换');
x1=9;y1=4.1;%小球1坐标
x3=7.1;y3=4.1;%小球2坐标
head1=line(x1,y1,'color','r','marker','.','markersize',40);%定义小球1
head2=line(x3,y3,'color','b','marker','.','markersize',40);%定义小球2
t=0;dt=0.01;%小球1的运动时间
t1=0;dt1=0.01;%小球2的运动时间
while 1
t=t+dt;
x1=9-t;y1=4.1;
if x1<7.3%等小球1运动至此位置，小球2开始运动
      x1=7.3;
      t1=t1+dt1;
      if rem(t1,14)<7
         x3=7-t1;
      else
         x3=t1-7;
      end;%使小球2能够来回运动
      y3=x3*(x3-7)/4+4.1;
end;
set(head2,'xdata',x3,'ydata',y3);
set(head1,'xdata',x1,'ydata',y1);
drawnow;
end;

作者: ahy 时间: 2012-9-24 20:07
把x3的表达式改成x3=abs(7-rem(t1,14));

欢迎光临数学建模社区-数学中国 (http://www.madio.net/)