动物集群运动行为模型~6

...wd......wd......wd...鱼类集群运动行为研究摘要群体智能是指生物群体中简单个体在相互作用中表现出复杂智能的行为。个体之间的组织构造、关系和群体行为的涌现机制是其研究的关键要素，因此，探究个体在群运动中遵循的规律是极其重要的。我们以鱼群为例对动物集群运动进展研究，针对这个鱼群集群运动模型，我们查找资料完成了以下问题：针对第一问，题目要求我们分析模拟集群运动，我们认为在系统中，集群、觅食、躲避障碍是一个整体，所以我们假设环境为二维平面，将集群、觅食、躲避障碍放在一起考虑，建设模型确定了集群运动的三个准则：对齐准则、靠近准则、防止碰撞准则。在此根基上我们通过编程模拟出鱼类的集群运动。针对第二问，我们在第一问的根基上，为鱼群增加逃避行为，并且鱼群遵守逃避准则。用MATLAB编程，实现了对鱼群逃避黑鳍礁鲨鱼模型的仿真。针对第三问，我们认为信息丰富者是群体内一局部固定的个体，定义为特殊个体，而其他的个体为一般个体，特殊个体的感知范围为一般个体的5倍，一般个体和特殊个体对伙伴中心的影响系数采用二八法。特殊个体主要承受来自环境的信息并影响群体内其他成员的行为，进而对群决策产生影响。在此根基上，我们通过编程进展模拟，可以发现特殊个体周围总是分布着许多一般个体，可见特殊个体对整个群体的影响力较一般个体大。我们在建模过程中并没有局限于题目要求的只研究集群运动规律，而是将觅食集群躲避障碍联系在一起进展研究，提高了仿真的精度。关键字：鱼群仿真集群运动群体智能1.问题重述在动物界，大量集结成群进展移动或者觅食的例子并不少见，这种现象在食草动物、鸟、鱼和昆虫中都存在。这些动物群在自然界中生活有很大优势，如回避捕食者、增加觅食时机等，这些动物群在运动过程中具有很明显的特征：群中的个体聚集性很强，运动方向、速度具有一致性。通过建设数学模型研究集群动物的行为是仿生学一项重要的内容，科学家通过对鱼群运动行为的研究而创造的人工鱼群算法帮助我们解决了许多实际问题，因此分析研究鱼群的集群运动行为具有十分重要的意义。针对这一问题，我们需要解决以下问题：〔1〕建设数学模型模拟动物的集群运动。〔2〕建设数学模型刻画鱼群躲避黑鳍礁鲨鱼的运动行为。〔3〕假定动物群中有一局部个体是信息丰富者(如掌握食物源位置信息，掌握迁徙路线信息)，请建模分析它们对于群运动行为的影响，解释群运动方向决策如何达成。2．问题分析此题要求对动物集群运动躲避威胁等行为建设数学模型，而动物一般不具有人类所具有的高级智能,所以不能采用自下而上的设计思路。我们应采用自下而上的设计思路,从个体行为出发,通过对个体的简单行为和个体与个体之间的行为规则的研究来建设模型。动物行为具有以下几个特点:1)适应性:动物通过感觉器官来感知外界环境,并应激性的做出各种反响,从而影响环境,表现出与环境交互的能力;2)自治性:动物有其特有的某些行为,在不同的时刻和不同的环境中能够自主的选取某种行为,而不是通过外界的控制或指导;3)盲目性:不像传统的基于知识的智能系统,有着明确的目标;单个个体的行为是独立的,与总目标之间没有直接的关系;4)突现性:总目标的完成是在个体行为的运动过程中突现出来的;5)并行性:各个体的行为是实时的、并行进展的。3.模型假设1.假设具体的环境是一个宽为width、长为length的二维平面，人工鱼只在二维平面内进展运动。2.人工鱼的初始速度为，当个体遭遇危险时速度突变为，速度变化时间不计。3.掌握丰富信息的个体是固定的，它们的感知范围是一般个体的5倍，它们受群成员的影响较一般个体小。4.人工鱼的能量小于初始能量时，就认为它处于饥饿状态，无法进展集群运动，即死亡。4.符号说明序号符号解释1能量值2新陈代谢率3感知范围半径4最大能量5初始能量6安全距离7初速度8最大速度9初始位置10t+1时刻个体的的游动方向11远离捕食者的方向12小于安全距离的伙伴到当前个体方向的平均值13伙伴的平均方向14当前个体到的方向15现阶段个体游动方向16伙伴中小于安全距离的伙伴个数17为伙伴的平均值18单个食物的能量19时间内该个体吃掉食物的个数20该个体被捕食的次数(只能取0或1)21初始时刻的新陈代谢率22系统内一点的食物浓度23人工鱼移动步长的最大值24拥挤度因子25与障碍物之间的距离26最大尝试次数27当前各个伙伴的位置28加权后当前各个伙伴的位置29各个伙伴的方向30感知范围内伙伴的个数5.模型的建设与求解5.1模型一的建设与求解5.1.1模型的准备鱼群集群行为的仿真模型由三局部组成:环境、个体、行为规则。环境是个体的生存空间,个体在环境中活动,依靠环境中的资源(食物)存活,个体的总和构成了研究的群体对象;群体的演化过程由其行为规则集来控制,行为规则决定了个体与个体之间、个体与环境之间相互作用的方式。个体依靠摄取环境中的食物资源和消耗自己的能量来生存,任意时刻个体在环境中都有一个确定的位置并且由其所占据的方位来确定。5.1.2环境的描述我们把具体的环境描述为一个宽为width、长为length的二维平面。这样我们就可以把现实中的鱼映射到虚拟环境中来,所有个体在该区域内觅食、集群。这个二维平面区域以及计算机时钟的运行共同构成了鱼群所处的虚拟环境。环境中包含食物、捕食者、天然障碍物等信息。我们用world[x][y]来表示坐标为〔x,y〕的点处的信息，world[x][y]=1，表示此处为捕食者，world[x][y]=2，表示此处为天然障碍物或伙伴。5.1.3个体描述个体〔鱼〕是建模的核心与关键，每条鱼都可以通过有限的感知能力发现外部环境的变化，并进展自主决策。个体依赖环境的资源而存在，由于鱼类视觉很有限，所以我们认为个体之间、个体与环境之间的的信息交流主要由其他的感知器官(主要是嗅觉)来完成,所以我们定义了个体的感知范围是一个半径为R的圆面。个体的行为主要有以下三种：觅食、逃跑、躲避障碍，下面对鱼的生存状态用函数表示:假设，个体处于觅食状态，假设，个体死亡，其中表示初始能量,表示单位时间消耗的能量即新城代谢率，初值为，与速度成正关系,表示单个食物的能量,表示时间内该个体吃掉食物的个数，与食物浓度成正比，与个体数量成反比,表示该个体被捕食的次数(只能取0或1)。个体在系统中具有以下属性:能量值:个体能量的多少由获得的资源量来衡量，当个体的能量值小于或等于0时,该个体死亡；位置:个体位置由一组坐标(x,y)表示；感知范围:个体所能感知的范围是一个半径为的圆；新陈代谢率:单位时间内个体所消耗的能量值。，其中是初始时刻新陈代谢率；最大能量:个体所能拥有能量的最大值。当个体的能量值到达该值时,停顿觅食；初始能量:仿真开场时,个体所拥有的能量值；安全距离:个体之间的最短距离，当个体间的距离小于该值时,个体将朝着远离的方向游动；初速度:仿真开场时,赋予每个个体的速度值,大小相等,方向不同。在没有遇到危险时将保持不变；最大速度:为了使得仿真形象,规定的个体速度最大值。假设在遇到危险时速度由初速度迅速变为5.1.4群体描述个体在大环境中生活，必然会与其他个体产生诸多联系，个体在集群时有以下三种行为准则：〔1〕靠近准则，即每个个体都有向伙伴中心靠拢的特性,伙伴中心为感知范围内各个体所在位置的平均值。〔2〕对齐准则，即个体会和它的伙伴朝同一个方向游动。〔3〕防止碰撞规则，即当个体和它的伙伴靠的太近时(距离小于安全距离),应自动避开。5.1.5模型的建设首先，我们研究一个个体，它的初始位置为，系统内一点的食物浓度，表示伙伴中心位置的食物浓度，STEP表示人工鱼移动步长的最大值;δ表示拥挤度因子，系统开场运行时，它在自身感知范围内，随机选择一个位置：1.觅食行为的完成：如果且，则该个体移动到，否则该个体继续选择另一点进展比拟，假设尝试try-number次后仍未找到适宜点，它将随机选择一点进展游动。2.躲避障碍的完成：如果与障碍物之间的距离不满足，则继续选择其他位置，直至满足上述条件。3.逃跑行为的完成：如果发现捕食者的位置，则个体速度变为，方向为远离捕食者方向为。其次，我们考虑个体在集群时的原则的实现，1.靠近规则的实现：每个个体都有向伙伴中心靠拢的特性,伙伴中心为观察范围内各个体所在位置的平均值。公式表示为其运动方向可表示为其中为伙伴的平均值，为当前个体的位置，为当前各个伙伴的位置，为当前个体到的方向。2.对齐规则的实现：个体会和它的伙伴朝同一个方向游动。公式表示为其中为各个伙伴的方向，为伙伴的个数，为伙伴的平均方向。3.防止碰撞规则的实现:当个体和它的伙伴靠的太近时(距离小于安全距离),自动避开。公式表示为为小于安全距离的伙伴到当前个体方向的平均值,M为伙伴中小于安全距离的伙伴个数。在此根基上鱼下一时刻的游动方向收到四个方向的共同作用，用公式表示为其中为下一时刻个体的的游动方向,为现阶段该个体的游动方向,为现阶段当前个体到伙伴平均位置的方向,为伙伴的平均方向,为小于安全距离的伙伴到当前个体方向的平均值。考虑到各规则对鱼的影响力不同,我们还需要对各个方向加权,取加权平均值,权重的大小可以根据偏好确定。其中5.1.6模型一的求解与检验5.1.6.1步骤一设定二维平面内食物浓度处处相等且不为零，平面内没有障碍。参数设定：λ1=0.5λ2=0.2λ3=0.2λ4=0.1，鱼的数量50，碰撞的最小距离10，最大速5.0度，初始速度3.0，感知半径150.〔蓝色的圆点代表鱼〕我们取鱼群集群过程的截图如图一、图二、图三，发现随着个体的运动，鱼群开场逐渐集群运动。证明我们的模型确实能模拟鱼的集群运动。图一图二图一图二图三图三5.6.1.2步骤二设定〔1,0〕，〔3,0〕，〔3,2〕，〔1,2〕四个点依次连接成的矩形为障碍物，〔-4,0〕，〔-2,0〕，〔-2,2〕，〔-4,2〕四个点依次连接形成的矩形区域内食物浓度为1，其他区域为0，其他参数不变。其过程如下四图，其中图七中只有19条鱼，说明有一条鱼因能量耗尽而死亡。这与我们建设的模型是相吻合的。图五图四图五图四图六图七图六图七5.2模型二的建设与求解5.2.1.鱼群运动模型的准备鱼群躲避黑鳍礁鲨鱼的运动行为可以分为两局部：一是与鲨鱼距离较远的局部鱼群，将其视为不受鲨鱼直接影响，只遵守靠近准则、对齐准则、防止碰撞准图各种域！则；〔红色以内为危险域，红色以外为安全域。〔红色以内为危险域，红色以外为安全域。黑色以内为感知域，蓝色为死亡域〕响，因此产生了逃跑行为，这局部鱼的运动会对第一局部鱼的运动产生间接影响。我们在模型一的根基上建设了模型二，增加了鱼的逃跑行为，即当天敌接近时，鱼会向某个方向逃跑，相应的，我们也增加了逃跑准则，即当天敌进入鱼的戒备域，鱼会逃跑,直到天敌处于安全域，当天敌进入死亡域时，鱼将被捕食。5.2.2鱼群运动模型的建设首先，我们通过对鱼与鲨鱼距离的判定将鱼群进展分组，与鲨鱼距离超过戒备域半径的鱼群作为鱼群一，视其为不受鲨鱼直接影响，与鲨鱼距离小于戒备域半径的鱼群作为鱼群二，受鲨鱼直接影响。1.鱼群一运动模型与模型一完全一样，在适应生存机制的前提下，鱼群的行为有随机行为、觅食行为、追尾行为、躲避障碍和集群行为，遵守准则为靠近准则、对齐准则、防止碰撞准则。具体步骤分析见模型一。2.鱼群二运动模型我们假设鱼群二中的个体逃跑行为和躲避障碍的优先级最高，即该局部鱼群不会发生觅食行为和追尾行为，只遵守靠近准则、对齐准则、防止碰撞准则和逃跑准则，且正在执行逃跑行为。首先，鱼群要适应生存机制，即其中(具体含义见模型一)其次，鱼群在执行逃跑行为时，依然会执行躲避障碍行为，并且遵守靠近准则、对齐准则、防止碰撞准则。〔1〕.躲避障碍的完成：如果与障碍物之间的距离不满足，则继续选择其他位置，直至满足上述条件。〔2〕1〕靠近准则的实现：每个个体都有向伙伴中心靠拢的特性,伙伴中心为观察范围内各个体所在位置的平均值。公式表示为其运动方向可表示为其中为伙伴的平均值，为当前个体的位置，为当前各个伙伴的位置，为当前个体到的方向。2〕对齐准则的实现：个体会和它的伙伴朝同一个方向游动。公式表示为其中为各个伙伴的方向，为伙伴的个数，为伙伴的平均方向。3〕防止碰撞准则的实现:当个体和它的伙伴靠的太近时(距离小于安全距离),自动避开。公式表示为为小于安全距离的伙伴到当前个体方向的平均值,M为伙伴中小于安全距离的伙伴个数。在此根基上鱼下一时刻的游动方向受到四个方向的共同作用，用公式表示为：其中为下一时刻个体的的游动方向,为现阶段该个体的游动方向,为现阶段当前个体到伙伴平均位置的方向,为伙伴的平均方向,为小于安全距离的伙伴到当前个体方向的平均值。考虑到各规则对鱼的影响力不同,我们还需要对各个方向加权,取加权平均值,权重的大小可以根据偏好确定。其中特别的，鱼群二还会执行逃跑行为，并且遵守逃跑准则。〔1〕逃跑行为的完成：，执行逃跑行为时,如果发现捕食者的位置，则个体速度变为，方向为远离捕食者方向为〔2〕逃跑准则的实现：假设，则被天敌捕食，即β=1,；假设，则执行逃跑行为，直到天敌位于安全域。5.2.3模型的求解参数设置保持不变，图中红色符号代表黑鳍礁鲨鱼，鲨鱼在水域中觅食〔鱼群中的个体〕，我们假设鲨鱼的速度为4.0，每次鲨鱼都向距离它最近的鱼游去。运动过程如以下列图：图九图八图九图八5.3模型三的建设与求解5.3.1 模型的描述在第一问的根基上，我们假设在鱼群中存在一局部信息比拟丰富的个体，它们在这类群体中经历比拟丰富，掌握的信息比拟多，在这里我们认为它们是一局部比拟固定的个体，即我们认为这类个体的感知能力比其他成员强，即我们不考虑某一个体偶然一次的获取丰富的信息。5.3.2环境的描述在这里我们仍然把具体的环境描述为一个宽为width、长为length的二维平面，与上述环境一样，不再赘述。5.3.3个体的描述在这里我们将个体分为一般个体和特殊个体，定义了一般个体的感知范围是一个半径为R的圆面，而特殊个体的感知范围为一般个体的5倍。特殊个体对来自于环境的信息较来自群体的信息更为敏感，也就是说特殊个体主要承受来自环境的信息并影响其他成员，而一般个体主要受特殊个体的影响。信息传递的主要内容是寻找食物源和躲避危险。5.3.4群体描述在这里我们主要考虑群体中的信息传递，所以我们认为只有个体生活在群体中才是我们应该考虑的对象，即我们认为当个体的现有能量小于零时，个体死亡。用函数表示为：假设个体死亡。5.3.5模型的建设同样道理在这里，我们还是先研究一个个体，该个体满足生存条件即其次，我们考虑个体在集群时的原则的实现，由于我们在之前已经定义了一般个体和特殊个体，所以在计算伙伴中心时也必须考虑两种个体的差异。我们在这里为每一个体的位置坐标赋权值，一般个体的为为0.2，而特殊个体的为0.8，计算时的坐标为加权后的坐标值。1.靠近规则的实现：每个个体都有向伙伴中心靠拢的特性,伙伴中心为观察范围内各个体所在位置的平均值。公式表示为其运动方向可表示为其中为伙伴的平均值，为当前个体的位置，为当前各个伙伴的位置，为当前个体到的方向。2.对齐规则的实现：个体会和它的伙伴朝同一个方向游动。公式表示为其中为各个伙伴的方向，为伙伴的个数，为伙伴的平均方向。3.防止碰撞规则的实现:当个体和它的伙伴靠的太近时(距离小于安全距离),自动避开。公式表示为为小于安全距离的伙伴到当前个体方向的平均值,M为伙伴中小于安全距离的伙伴个数。在此根基上鱼下一时刻的游动方向收到四个方向的共同作用，用公式表示为其中为下一时刻个体的的游动方向,为现阶段该个体的游动方向,为现阶段当前个体到伙伴平均位置的方向,为伙伴的平均方向,为小于安全距离的伙伴到当前个体方向的平均值。考虑到各规则对鱼的影响力不同,我们还需要对各个方向加权,取加权平均值,权重的大小可以根据偏好确定。其中5.3.6模型的求解参数设定依然保持不变，图中红色符号代表信息丰富的鱼，该鱼食物密集地，即〔-2,0〕，〔-2,2〕，〔-4,2〕四个点依次连接形成的矩形区域。运动过程如以下列图所示：图十一图十图十一图十6.模型的改良与推广该模型只是在二维平面内讨论鱼的集群运动，具有局限性，应该进一步分析三维立体空间内的运动。在事实上，鱼类遇到危险时，速度变化需要时间的，所以我们在改良时，在将二维变空间变为三维空间的同时还应该定义一个加速度来描述鱼类的速度变化。除此之外，还有一局部鱼类的信息传递是依靠光线等，在这种情况下，鱼类的感知范围就不再是一个半径为的圆面，而是一个扇形，所以模型还有很大的改良空间。参考文献[1]赵建,曾建潮.《鱼群集群行为的建模与仿真》北京：太原，2009.[2]班晓娟等，《人工鱼群高级自组织行为研究》北京，2007.[3]于飞等，《根本人工鱼群算法的研究与改良》河北：陕西，2005.附录n=50;fish_setp=50;dim=2;c2=1.2;c1=0.12;w=0.9;fitness=0*ones(n,fish_setp);R1=rand(dim,n);R2=rand(dim,n);current_fitness=0*ones(n,1);current_position=10*(rand(dim,n)-.5);velocity=.3*randn(dim,n);local_best_position=current_position;fori=1:ncurrent_fitness(i)=Live_fn(current_position(:,i));endlocal_best_fitness=current_fitness;[global_best_fitness,g]=min(local_best_fitness);fori=1:nglobl_best_position(:,i)=local_best_position(:,g);endvelocity=w*velocity+c1*(R1.*(local_best_position-current_position))+c2*(R2.*(globl_best_position-current_position));current_positi