用Matlab模拟四人追逐问题

qibbxxt

看到有人用mathematica模拟，我也写了一个Matlab的模拟程序：

1. 
   
2. clear;clc;close all
   
3. rectangle('Position',[1,1,1,1],...
   
4.     'Curvature',[0,0],...
   
5.     'LineWidth',2,'LineStyle','-')
   
6. axis([0.9 2.1 0.9 2.1])
   
7. daspect([1,1,1])
   
8. x=[1 1 2 2];
   
9. y=[1 2 2 1];
   
10. n=10000;
    
11. axis off
    
12. hold on
    
13. h=zeros(length(x),1);
    
14. for i=1:4
    
15.     h(i)=plot(x(i),y(i),'r.','markersize',10,'erasemode','none');
    
16. end
    
17. title('\copyright:qibbxxt')
    
18. t=0.001;
    
19. dirx=[0,t,0,-t];
    
20. diry=[t,0,-t,0];
    
21. for i=1:n
    
22.     x=x+t*dirx*i/10;
    
23.     y=y+t*diry*i/10;
    
24.     dirx=diff([x,x(1)]);
    
25.     diry=diff([y,y(1)]);
    
26.     if dist([x(1),y(1)],[x(2);y(2)])<15*t
    
27.         break
    
28.     end
    
29.     arrayfun(@set,h,repmat({'xdata'},4,1),num2cell(x'),repmat({'ydata'},4,1),num2cell(y'));
    
30.     pause(0.01);
    
31. end
    

复制代码

mathematica程序见：http://forum.simwe.com/thread-949737-1-1.html