zhuizhu.m

来自「本程序是用MATLAB编写」· M 代码 · 共 43 行

   function zhuizhu

    %取v=1,t=12,A,B,C,D点的坐标分另为(0,10),（10,10）,（10,0）,(0, 0)

    v=1;

    dt=0.05;

    d=20;

    x=[0  0  0   10   10   10   10  0];

    x(9)=x(1);

    x(10)=x(2);

    hold

    %axis('equal')

    axis([0 10 0 10]);
    
     while(d>0.1)

          for  i=1:2:7

               d=sqrt((x(i)-x(i+1))^2+(x(i+1)-x(i+3))^2);

               x(i)=x(i)+v*dt*(x(i+2)-x(i))/d;

               x(i+1)=x(i+1)+v*dt*(x(i+3)-x(i+1))/d;

               plot(x(i),x(i+1),'ob')

          end

         x(9)= x(1);x(10)= x(2);

      end
  
   hold