动物集群运动行为模型-12(5)

end

if x(k+1,i)>40

x(k+1,i)=x(k+1,i)-40; elseif x(k+1,i)<0 x(k+1,i)=x(k+1,i)+40; end

if y(k+1,i)>40

y(k+1,i)=y(k+1,i)-40; elseif y(k+1,i)<0 y(k+1,i)=y(k+1,i)+40; end end end

for i=1:k pause(0.2)

plot(x(i,:),y(i,:),'o',xs(i,:),ys(i,:),'r*') axis([0 40 0 40]) getframe; end

20