基于Biod模型的动物集群运动分析

基于Biod模型的动物集群运动分析附录问题一的程序：clear;n=50;eta=1;%设定参数range=0.5;%间距speed=0.03;TIME=600;density=1;RUNS=1;fieldsize=sqrt(n/density);%边长pos=fieldsize*rand(n,2);%随机生成n个坐标（位置）heading=2*pi*rand(n,1);%方向relposX=zeros(n);%水平距离relposY=zeros(n);%垂直距离neighbours=zeros(n);%标记0-1矩阵RelHead=zeros(n,1);%方向差newHeading=zeros(n,1);newPos=zeros(n,2);meanHeading=zeros(TIME,1);%任意时刻的平均方向deviationMean=zeros(n,TIME);forrun=1:RUNS%第一层循环Run次数fortime=1:TIME%第二层循环时间fori=1:nforj=1:nrelposX(i,j)=abs(pos(i,1)-pos(j,1));%水平距离relposY(i,j)=abs(pos(i,2)-pos(j,2));%垂直距离if(sqrt(relposX(i,j)^2+relposY(i,j)^2)<=range)neighbours(i,j)=1;%判断是否位于排斥区endendendfori=1:nforj=1:nif(neighbours(i,j)==1)%如果位于排斥区relheading(i,j)=heading(j,1)-heading(i,1);%方向差endendendRelHead=sum(relheading,2);%每行求和，结果为列向量fori=1:nwhileRelHead(i,1)<-piRelHead(i,1)=RelHead(i,1)+pi;endwhileRelHead(i,1)>piRelHead(i,1)=RelHead(i,1)-pi;endendnoise=(rand(n,1)*eta)-eta/2;%噪声newHeading(:,1)=heading(:,1)+RelHead(:,1)./sum(neighbours,2)+noise(:,1);newPos(:,1)=pos(:,1)+cos(newHeading(:,1))*speed;newPos(:,2)=pos(:,2)+sin(newHeading(:,1))*speed;fork=1:nwhilenewPos(k,1)<0newPos(k,1)=newPos(k,1)+fieldsize;endwhilenewPos(k,1)>=fieldsizenewPos(k,1)=newPos(k,1)-fieldsize;endwhilenewPos(k,2)<0newPos(k,2)=newPos(k,2)+fieldsize;endwhilenewPos(k,2)>=fieldsizenewPos(k,2)=newPos(k,2)-fieldsize;endendheading=newHeading;pos=newPos;meanHeading(time,1)=mean(heading,1);deviationMean(time,1)=mean(abs(meanHeading(time,1)-heading(:,1)));scatter(pos(:,1),pos(:,2),'xr');%画图axis([0fieldsize0fieldsize]);xlabel(time);%画图控制轴线范围M(time)=getframe;%为输出视频建立矩阵endrunDe