基于biod模型的動物集群運動分析-資料下載頁
2025-06-20 12:28本頁面
【正文】 2*pi*rand(n,1)。%方向relposX = zeros(n)。%水平距離relposY = zeros(n)。%垂直距離neighbours = zeros(n)。%標記01矩陣RelHead = zeros(n,1)。%方向差newHeading = zeros(n,1)。newPos = zeros(n,2)。meanHeading = zeros(TIME,1)。%任意時刻的平均方向deviationMean = zeros(n,TIME)。for run=1:RUNS %第一層循環(huán)Run次數(shù)for time=1:TIME %第二層循環(huán) 時間for i=1:n for j=1:n relposX(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。%判斷是否位于排斥區(qū)endendendfor i=1:nfor j=1:nif(neighbours(i,j)==1)%如果位于排斥區(qū)relheading(i,j) = heading(j,1)heading(i,1)。%方向差endendendRelHead = sum(relheading,2)。%每行求和，結(jié)果為列向量for i=1:nwhile RelHead(i,1) piRelHead(i,1) = RelHead(i,1) + pi。endwhile RelHead(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。for k=1:nwhile newPos(k,1) 0newPos(k,1) = newPos(k,1) + fieldsize。endwhile newPos(k,1) = fieldsizenewPos(k,1) = newPos(k,1) fieldsize。endwhile newPos(k,2) 0newPos(k,2) = newPos(k,2) + fieldsize。endwhile newPos(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), 39。xr39。)。%畫圖axis([0 fieldsize 0 fieldsize])。 xlabel(time )。%畫圖 控制軸線范圍M(time) = getframe。%為輸出視頻建立矩陣endrunDev(:,run) = deviationMean(:,1)。end%plot(runDev)%movie2avi(M,39。39。, 39。quality39。,100)。
物理相關(guān)推薦
鄂ICP備17016276號-1