动物集群运动行为模型-3

1动物集群运动行为模型摘要通过观看大量集结成群进行移动或者觅食的动物行为视频和探究动物集群运动的机理，我们建立了鱼群模型模拟动物的集群运动，建立微分方程模型研究鱼群躲避黑鳍礁鲨鱼的运动行为，建立/AR模型分析动物群中有一部分信息丰富者对于群运动行为的影响，并且解释群运动方向决策如何形成。针对问题一，通过个体与个体之间以及个体与环境之间的相互作用来推导模拟整个鱼群的运动。个体鱼具有一定的感知能力以及遵循下列三个行为规则：（1）避免与相邻的鱼发生碰撞冲突；（2）尽量与自己周围的鱼在运动方向上保持协调和一致；（3）向鱼自己周围的邻居的位置中心运动。建立出112341234tttttDDDDD从而通过matlab编程得出模拟动物的集群行为图。（见图1.1）针对问题二，通过对鱼群轨迹和鲨鱼轨迹的分析，在鲨鱼追踪鱼群的任何时刻都要朝向鱼群的运动，我们建立微分方程模型来模拟鲨鱼的追踪和鱼群的躲避的运动过程。鱼群的位置121,+11,12PPjjikvtik鲨鱼的位置2,1,2,12,2,1,PPjjjjjjPPVtPP从而得出鲨鱼的追踪和鱼群躲避图（见图2.1）针对问题三，假设鱼群中有一部分领导者，它们掌握着丰富信息，根据掌握信息的多少，我们将之分为领导者和次领导者。通过建立A/R模型分析发现，次领导者的个数,和预测步长（领导者和次领导者间的距离）是影响集群信息传递的两大因素。领导者将重要信息传递给次领导者，次领导者然后传递给鱼群中的跟随者。领导者和各个次领导者间的距离不宜过大，同时次领导者数量应维持在一定数目，过多的次领导者反而影响信息传递，成为多余。2关键词：鱼群模型集群运动/AR模型微分方程模型目录一、问题重述………………………………………………………………51.1问题背景…………………………………………………………51.2问题提出…………………………………………………………5二、模型假设……………………………………………………………5三、符号说明………………………………………………………………5四、问题分析………………………………………………………………6五、模型的建立与求解……………………………………………………75.1问题一……………………………………………………………75.1.1鱼群模型的建立…………………………………………75.1.2鱼群模型的求解…………………………………………85.2问题二……………………………………………………………115.2.1微分方程模型的建立……………………………………115.2.2微分方程模型的求解……………………………………125.3问题三……………………………………………………………135.3.1/AR模型的建立…………………………………………135.3.2/AR模型的求解…………………………………………14六、模型的评价与推广……………………………………………………196.1模型的优缺点……………………………………………………196.2模型的推广………………………………………………………19参考文献……………………………………………………………………213附录…………………………………………………………………………22一、问题重述1.1问题背景在动物界，大量集结成群进行移动或者觅食的例子并不少见，这种现象在食草动物、鸟、鱼和昆虫中都存在。这些动物群在运动过程中具有很明显的特征：群中的个体聚集性很强，运动方向、速度具有一致性。通过数学模型来模拟动物群的集群运动行为以及探索动物群中的信息传递机制一直是仿生学领域的一项重要内容。1.2问题提出通过观察附件中给出的图片和视频资料和在网上搜索相关资料观察，思考动物集群运动的机理，需要建立数学模型刻画动物集群运动、躲避威胁等行为；从而明确了我们需要建立数学模型解决以下问题:(1)建立数学模型模拟动物的集群运动,即建立集群运动的正常运动的模型。(2)建立数学模型刻画鱼群躲避黑鳍礁鲨鱼的运动行为，即建立鱼群应激反应的模型。(3)假定动物群中有一部分个体是信息丰富者，建立数学模型分析它们对于集群运动行为的影响，并解释集群运动方向决策如何达成的。二、模型假设1.集群中的每个个体都有形同的特性且健康正常。2.每个个体都有向邻居中心靠拢的特性。3.每个个体会和它的邻居朝同一个方向游动。4.当个体和它的邻居靠的太近时(距离小于碰撞距离)都自动避开。5.鲨鱼追踪鱼群的任意时刻,鲨鱼始终朝向鱼群运动。6.集群的领导者是信息丰富的个体。7.模型开始运行时，环境中任意分布一定数量的个体,每个个体具有自己的状态属性。三、符号说明1tD周期t中个体的游动方向2tD周期t个体到邻居平均位置的方向3tD个体的邻居的平均方向,4tD小于碰撞距离的邻居到当前个体方向的平均值Di个体各个邻居的方向iP个体各个邻居的位置4ijd个体i和j的物理距离ptx次领导者的位置状态itx跟随者的位置状态PH次领导者的预测步数plN次领导者的个数pJ有预测机制的有领导者集群的群体性能指标vJ无预测机制的有领导者集群的群体性能指标1,jP鱼群在第j时刻的位置2,jP鲨鱼在第j时刻的位置R鲨鱼的感知范围r鱼群的感知范围1惯性系数2应激系数d鲨鱼与鱼群的距离四、问题分析针对问题一，以鱼群为研究对象，采用鱼群集群运动模型，先从个体鱼出发研究其运动的一般特点，再通过个体与个体之间以及个体与环境之间的相互作用推导来模拟整个鱼群的运动。个体鱼具有一定的感知能力以及遵循下列三个行为规则：（1）避免与相邻的鱼发生碰撞冲突。（2）尽量与自己周围的鱼在运动方向上保持协调和一致。（3）尽量试图向自己所认为的群体中靠近。针对问题二，通过对鱼群轨迹和鲨鱼轨迹的分析，在鲨鱼追踪鱼群的任何时刻都要朝向鱼群的运动，当鲨鱼运动轨迹和鱼群运动轨迹相交时，即鲨鱼捕食成功。要求二者运动轨迹，建立微分方程模型，根据给定的条件，对该模型求解，运用MATLAB软件求出解。针对问题三，动物群中有一部分个体是信息丰富者,因此把经验丰富的个体5视为领导者，把经验较为丰富的个体视为次领导者，其他的个体视为跟随者。为了更清楚的说明集群同步预测机制的优势,我们将对加入预测机制的有领者集群进行统计仿真同时分析和比较无预测机制的有领导者集群和有预测机制的有领导者集群的群体性能指标vJ和pJ。五、模型的建立与求解5.1问题一5.1.1鱼群模型的建立每个个体鱼运动要遵循以下三个行为规则：避免碰撞指向邻居中心邻居方向个体的游动方向图1.1鱼游动方向（1）每个个体鱼有向邻居中心靠近：邻居中心为在观察范围内每个个体鱼所位置的平均值：iPiNNP020arctantyyxxDPxy（,）为邻居平均值，000(,)Pxy为当前位置，iP为当前各个邻居的位置，2tD为当前个体到P的方向。（2）与邻近鱼保持方向上一致：Di为各个邻居的方向,N为邻居的个数,3Dt为邻居的平均方向。（3）避免与相邻的鱼发生碰撞：040arctanitiyyMiMxxD3itDiNND64Dt为小于碰撞距离的邻居到当前个体方向的平均值,M为邻居中小于碰撞距离的邻居个数。另外鱼在运动的过程中，要想改变方向，还需考虑在前一时刻运动方向上的惯性作用。综合考虑四个因素的影响，且每个因素影响的程度不一样，设各个影响因素的权重为1、2、3、4。112341234tttttDDDDD且12341。5.1.2鱼群模型的求解鱼游动方向的确定：112341234tttttDDDDD需要对1、2、3、4的大小进行在matlab程序调试中改进，从而得出了模拟动物的集群行为图。（程序见附录1）图1.1动物的集群行为图（5张）7895.2问题二5.2.1微分方程模型的建立模型模拟是鲨鱼追踪鱼群的过程。鲨鱼以速度V向鱼群追来，鲨鱼始终以朝向鱼群的方向运动，即沿着鲨鱼与鱼群的连线方向，其运动的单位方向向量为：2,1,2,1,PPjjjjPkP然而鱼群当感知到鲨鱼的追踪时，也会调节其运动方向，但鱼群的运动方向要受到多方面徳影响：鱼群本身的运动方向、鲨鱼的追踪方向、水流的方向、以及其他因素造成徳影响。在这里我们主要考虑两方面的影响：（1）鱼群自身的运动方向，其方向向量的初值为：1,1,11,1,1PPPPjjjji（2）鲨鱼的追踪方向，其方向向量即为k。考虑到这两方面因素对鱼群运动徳影响重要性不一样，我们对这两个因素设置权重系数：惯性系数1，应激系数2。则鱼群运动方向单位向量为：1212ikcik10运动时间T离散化:由于鱼群与鲨鱼的运动是随着时间随时发生变化的，因此我们将他们运动的时间离散化。即将时间t等分为为n份，其中一份时间长度定为t。运动规律分析：表示鱼群在第j时刻的位置，2,2,2,(,)jjjPxy表示鲨鱼在第j时刻的位置；则在+1j时刻鱼群的位置向量坐标表示为：121,+11,12PPjjikvtik+1j时刻鲨鱼的位置向量坐标表示为2,1,2,12,2,1,PPjjjjjjPPVtPP5.2.2微分方程模型的求解（程序见附录2）在前面问题的分析中，我们给鱼群的运动方向设置两个影响的参数：惯性参数1和应激参数2。二者满足：121我们在研究其运动时考虑鱼群运动主要是运动方向的两方面参数：初始运动方向以及两个参数1、2。在研究他们作用时，我们分别取不同的值来观察鲨鱼的追踪效果图，如下图所示：图2.1鲨鱼的追踪和鱼群躲避效果图122,0,0.8,0.2r;122,0,0.5,0.5r11122,0,0.2,0.8r;1252,,0.2,0.86r(1)由图知当鱼群受到鲨鱼攻击时，惯性系数1基本不发生作用，起决定性作用的是应激系数2。(2)由后两个图可知，当鱼群受到攻击时，不论其初始时刻运动方向如何，当鱼群发现受攻击时会逐渐沿着追击方向远离。5.3问题三5.3.1/AR模型的建立为了更清楚的说明集群同步预测机制的优势,我们将对加入预测机制的有领者集群进行统计仿真同时分析和比较无预测机制的有领导者集群和有预测机制的有领导者集群的群体性能指标vJ和pJ。首先，我们引入A/R模型对有领导者集群进行仿真，我们用L表示领导者，P表示次领导者，F表示跟随者。其中A/R函数为expijijijdgddabc,,abc为三个自由常量，ijd表示个体i和j的物理距离。则次领导者的位置状态ptx由如下方程确定：,pjPppLpjjLrtgtgtdxddH式中PpLgtdH代表了次领导者与领导者的长边，,pjpjjLrgtdd表示周围邻居对次领导者的作用。跟随者的位置状态itx由如下方程确定：,ijijjtrgtdd12PpLtdH是此领导者与领导者未来PH步时位置之间的距离。5.3.2/AR模型的求解我们考虑一个有M=50个节点的网络，PH表示次领导者的预测步数，PH=10，p-leader表示次领导者，plN表示次领导者个数，2plN，A/R函数参数为a=1，b=2.2，c=0.2。领导者的轨迹函数yx，领导者的速度0.02xkv，0.0210.02xkkky，k表示采样时间。其中◇表示领导者，★表示次领导者，○表示跟随者，蓝线为领导者的轨迹。图3.3有领导者集群的初始状态在经过65步后。当0PH时，未加入预测机制的集群的运动状态如图3.4所示，可以看出群体的个体形成一个环，由于次领导者无法预测领导者未来的状态，群体的个体慢慢地掉队，偏离领导者轨迹。13图3.4无预测机制的集群的稳定状态另一方面我们对有预测机制的集群进行仿真,图3.5为Hp=10时有预测机制的群体经过65步后的运动状态