動物集群運動行為模型-(編輯修改稿)
2025-08-18 00:50 本頁面
【文章內(nèi)容簡介】 s(1,3)。 endfor t=1:1000for 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。 %內(nèi)聚性和排列性 direc8{i}=(direc8{i}+direc1{r})./r。 end if (norm(pos{i}pos{j})=5)amp。amp。(norm(pos{i}pos{j})=) %向心性和同向性 m=m+1。k=k+norm(direc1{j})。l=l+direc1{j}。 direc2{i}= (direc2{i}+(pos{j}pos{i}))./m。 direc3{i}=(direc3{i}+direc1{j})./m。 cons(i)=norm(l)/k。 elseif norm(pos{i}pos{j}) n=n+1。 direc4{i}=(direc4{i}+pos{i}pos{j})./n。 %排斥性 end end end if (pos{i}(1)=10)||(pos{i}(2)=10)||(pos{i}(1)=0)||(pos{i}(2)=0)||(pos{i}(3)=10)||(pos{i}(3)=0) direc1{i}=direc1{i}*(1)。 %限定個體活動范圍 else direc1{i}=*direc1{i}+*direc2{i}+*direc3{i}+*direc4{i}+*direc7{i}+*direc8{i}。 end pos{i}=pos{i}+direc1{i}/norm(direc1{i})*v*exp(5*(cons(i)1))。endfor i=1:30 x(i)=pos{i}(1)。 y(i)=pos{i}(2)。 z(i)=pos{i}(3)。endplot3(x,y,z,39。.39。)axis([0 10 0 10 0 10],39。square39。,39。manual39。)grid on。pause()end、問題二的程序代碼：pos=cell(1,100)。for i=1:100 pos{i}=[rand(1,1)*5,rand(1,1)*5,rand(1,1)*5]。 %縮小活動范圍至5*5*5endpre=[rand(1,1)*5,rand(1,1)*5,rand(1,1)*5]。direcpre=pos{1}pre。for i=2:100 %增加個體數(shù)目至100 if norm(pos{i}pre)norm(direcpre) direcpre=pos{i}pre。 endend x=zeros(1,100)。y=zeros(1,100)。z=zeros(1,100)。v=。direc1=cell(1,100)。direc2=cell(1,100)。direc3=cell(1,100)。direc4=cell(1,100)。direc5=cell(1,100)。direc6=cell(1,100)。direc7=cell(1,100)。direc8=cell(1,100)。cons=ones(1,100)。for (i=1:100) direc1{i}=[rand(1,1)*5,rand(1,1)*5,rand(1,1)*5]pos{i}。 direc2{i}=zeros(1,3)。 direc3{i}=zeros(1,3)。 direc4{i}=zeros(1,3)。 direc7{i}=zeros(1,3)。 direc8{i}=zeros(1,3)。endfor t=1:2000 pre=pre+direcpre/norm(direcpre)**v/10。for i=1:100 m=0。 n=0。 k=0。 l=[0,0,0]。 direc5{i}=zeros(1,3)。 direc6{i}=zeros(1,3)。 for j=1
范文總結(jié)相關(guān)推薦
鄂ICP備17016276號-1