行为研究,北京大学Р[5]翟超,张海涛,生命群协调行为模型的改进及同步控制研究,华中科技大学Р七、附录Р8.1、问题一的程序代码:Рpos=cell(1,30);Рfor i=1:30Р pos{i}=[rand(1,1)*10,rand(1,1)*10,rand(1,1)*10] %定义30个个体РendРx=zeros(1,30);Рy=zeros(1,30);Рz=zeros(1,30);Рv=0.1; %最大速度Рdirec1=cell(1,30); %定义各个方向向量Рdirec2=cell(1,30); Рdirec3=cell(1,30);Рdirec4=cell(1,30);Рdirec7=cell(1,30);Рdirec8=cell(1,30);Рcons=zeros(1,30);Рfor i=1:30Р direc1{i}=[rand(1,1)*10,rand(1,1)*10,rand(1,1)*10]-pos{i};Р direc2{i}=zeros(1,3);Р direc3{i}=zeros(1,3);Р direc4{i}=zeros(1,3);Р РendРfor t=1:1000Рfor i=1:30Р m=0;n=0;k=0;l=[0,0,0];direc7{i}=zeros(1,3);Р direc8{i}=zeros(1,3);Р for j=1:30Р if j==iР continueР elseР for r=1:30Р direc7{i}=(direc7{i}+pos{r}-pos{i})./r; %内聚性和排列性Р direc8{i}=(direc8{i}+direc1{r})./r;Р endР if (norm(pos{i}-pos{j})<=5)&&(norm(pos{i}-pos{j})>=0.1) %向心性和同向性