鱼群运动行为模型 精品

1、鱼群运动行为模型摘要本文研究了鱼群运动时受环境及邻近同族的影响而改变速度方向的机制，并 以此为基础分析了鱼群在躲避捕食者和觅食时的信息传递和转移路线。对于问题一，本文考虑平衡状态时，即没有捕食者威胁也无觅食和迁移的需 求时，个体鱼的游动规律。本文假设个体鱼在二维平面内游动时能够感知到一定 范围（R）内的同族的位置和游动方向，并遵循四个规则：惯性规则、靠近规则、 对齐规则、规避规则，个体鱼的运动方向由这四个规则对鱼的影响大小决定，X = X + V - cos PP 人P +人P +人P +人P，t+11 1t 2 2t3 3t4 4tV +1 - 了 + V. sin P 。由此可对每一条鱼的

2、 y y v sin 1运动状态进行迭代更新。对于问题二，本文考虑在二维平面中引入捕食者，并假设捕食者将游向其感 知范围（R0）内距离其最近的个体鱼，同时受其自身游动惯性的影响，则其游动 方向P -人P +人P。由此可对捕食者的游动状态进行迭代更新。当捕食者靠 t+1 1 1t 2 2t近个体鱼，出现在个体鱼的感知范围内时，小鱼将产生避险意识，避险方向为捕 食者到个体鱼的方向，同时向其感知范围内的个体鱼发送告警信号，接受到告警 信号的个体鱼将产生离散意识，离散方向为其感知到的避险个体鱼游动方向的平均方向。则此时小鱼的游动方向P -人P +人P +人P +人P +人P +人P。由 t+1 1 1

3、t 2 2t 3 3t 4 4t 5 5t 6 6t此可对捕食者和个体鱼的运动状态进行迭代更新。对于问题三，本文仅考虑掌握食物源位置信息的信息丰富者，它们在遵循问 题一中提出的游动规则条件下，将主动靠近食物源，并且把它向食物源游去的信 息告知邻居，召集其它个体鱼共同觅食。对于非信息丰富者来说，它能受到其感 知范围内信息丰富者的召集信息，并趋向这些信息丰富者的实际游动方向的平均 方向，追随它们共同觅食。此时个体鱼 的游动方向：P 人P +人P +人P +人P +人P +人P。对于信息丰富者，受到召集作用的t+1 1 1t 2 2t 3 3t 4 4t 5 5t 6 6t权重X6 0。对于非信息丰

4、富者，游向食物源的权重七-0。由此可得鱼群觅食 的集群运动情况。关键词：个体运动集群运动运动规则一、问题重述1.1问题背景在动物界，大量集结成群进行移动或者觅食的例子并不少见，这种现象在食 草动物、鸟、鱼和昆虫中都存在。这些动物群在运动过程中具有很明显的特征： 群中的个体聚集性很强，运动方向、速度具有一致性。在生态系统中，动物个体 的行为相对简单,集群后却能表现出复杂的群体行为。个体行为是构成群体行为 的基础，个体之间的组织结构、个体行为之间的关系和群体行为的涌现机制是研 究群体行为的关键要素。通过数学模型来模拟动物群的集群运动行为以及探索动 物群中的信息传递机制一直是仿生学领域的一项重要内容

5、。目前主要研究有仿生 的群体优化算法，群体组织内部的通信机制及其应用方面，如微粒群算法、蚁群 算法、群体机器人等。1.2目标任务题目要求查阅相关资料，思考动物集群运动的机理，建立数学模型刻画动物 集群运动、躲避威胁等行为，主要针对以下问题分析建模：建立数学模型模拟动物的集群运动。建立数学模型刻画鱼群躲避黑鳍礁鲨鱼的运动行为。假定动物群中有一部分个体是信息丰富者(如掌握食物源位置信息，掌握 迁徙路线信息)，请建模分析它们对于群运动行为的影响，解释群运动方向决策 如何达成。二、模型假设假设所有个体鱼具有相同的物理特性；假设所有个体鱼具有有限感知能力以及遵循共同的行为规则；假设所有个体鱼之间的相互作

6、用及信息的交互不受外界因素影响;假设捕食者和个体的运动速度相等并且保持不变；假设所有信息丰富者掌握的是同一食物源的信息。三、符号说明符号符号说明A.1R鱼群中的个体鱼个体鱼感知范围的半径Pt 1下一个时间步长周期个体鱼的运动方向Pit本周期t中的游动方向P2t周期t当前个体到邻居平均位置的方向P3t邻居的平均方向P4t小于碰撞距离的邻居到当前个体方向的平均值个体鱼决策方向是的权重值P5t个体鱼为了躲避威胁而以最快方式逃逸的方向P6t该个体鱼受其它个体发出的告警信号的平均方向rR0个体间避免碰撞的最小距离 捕食者的威胁范围四、问题分析在自然界中，动物出于生存、避险、觅食、求偶、繁殖等原因往往选择

7、采取 群体行为的方式。某些动物个体的行为相对简单，集群后却能表现出复杂的群体 行为。这些动物群在运动过程中具有聚集性很强，运动方向、速度保持一致性等 明显的特征。在生态系统中，个体行为是构成群体行为的基础。本文认为由个体 简单运动到群体复杂行为是存在一种映射关系的。个体之间的组织结构、个体行 为之间的关系和群体行为的涌现机制是研究群体行为的关键要素。文章的总体思 路是要从个体的行为、个体与个体的相互作用、个体在群体中的作用等角度出发， 通过数学方法模拟动物的群体运动行为，并利用Matlab软件编程实现模型的仿 真，探索动物群的集群运动行为以及动物群中的信息传递机制。对问题一：题目要求建立合理的

8、数学模型来模拟动物的集群运动。为了能够 建立简洁易懂的模型，在本文中选取鱼群为研究对象，研究鱼群的形成和行为， 希望通过研究个体鱼的运动机制科学阐释集群行为的内在变化原因。针对鱼群的 形成、结构和行为，很多研究者从不同角度提出了一些理论和模型。本文通过研 究个体鱼在群体中位置和速度随时间的不断变化的规律，采取从个体到局部，从 局部到整体，自下向上的建模思想，对个体鱼进行建模，进而通过个体遵循一定 的行为规则来研究个体之间以及个体和环境之间的相互作用，最终探讨出鱼群集 群运动的形成机理，并合理推广至对动物群的研究。对问题二：题目要求建立数学模型刻画鱼群躲避黑鳍礁鲨鱼的运动行为。针 对这个问题，本

9、文在问题一研究的基础上，考虑仍然沿用其个体鱼的运动模型。 当有某些个体受到鲨鱼威胁时，它们一方面选择最快的逃逸方式躲避，一方面发 出告警信号与其感知范围内的其它个体进行信息交互。这样，这些受威胁的个体 鱼和得到告警信号的个体鱼在下一时刻的运动方向就较其它个体及其上一时刻 有了较大不同。基于此，可以采用数学公式计算其不同范围内的个体鱼的运动方 向和位置坐标的变化。对问题三：假设动物群中的部分信息丰富者是提前已知的，并且其不一定是 集聚的。为了简化问题的研究，本文仍然以鱼群为例，探讨鱼群觅食的集群行为。 当有一部分个体掌握食物源位置信息，这部分个体直接向食物源运动，同时通过 个体间的信息交互将信息

10、传递给其它个体，引起其它个体运动状态和位置坐标的 改变，从而形成集群的觅食行为。五、模型的建立5.1问题一的模型5.1.1模型的分析个体行为是构成群体行为的基础，在一个集群中所有个体行为的匹配结合就 映射一种集群的运动行为。为了研究集群的运动行为就必须研究清楚集群中每一 个个体的运动行为和由个体构成的局部的运动行为。对于研究的鱼群行为来说， 本文采取自下而上的建模方法，根据对问题的分析，文章对这种模型分为三层：最下层分析个体鱼（设为A1.A.气）的运动模型，并且假定个体鱼的运动 满足设定的游动规则5，模型具有普遍地适用性。中间层通过分析个体鱼之间 的相互作用，个体鱼的运动信息所能影响到的范围，

11、形成局部（局部1局部n） 的行为。在最上层分析局部行为扩大到集群就构成集群运动的行为模型。从个体 到集群的行为关系略图如下图。图1个体到集群的行为关系图5.1.2模型的建立5.1.2.1个体鱼运动模型个体鱼运动的的简化本文描述的是个数为N的一群可视为质点的个体鱼在L*L的二维周期边界 条件的平面运动的情况，并且每一时刻每个个体鱼运动的速度大小均相同。对个体鱼的定义每一个个体鱼是一个自驱动的个体，形状大小相同，具有一定的感知能力。 它能感知的范围是以其质心为圆点的半径为R的圆形区域，且它能感知这一区 域内其他个体鱼的所有动态信息（包括速度大小和方向，有没有发出告警信号等）。运动规则描述（1）惯性

12、规则，个体鱼在得到需要改变运动的信号时，鱼游动的方向不可 能立刻改变，这时表现出一种惯性的作用。（2）靠近规则，为了不脱离鱼群，需要尽量靠近邻居的中心。如图a.（3）对齐规则，为了保持鱼群运动的连贯性，每个个体鱼尽量与邻居的方 向一致。如图b.（4）规避规则，为了保持鱼群运动的一致性，个体鱼运动尽量避免碰撞。 如图c.图a图b图c图2运动规则描述示意图运动模型的建立假设每一个体鱼在邛寸刻具有相同大小的速度，而运动速度的方向是任意的 或随机的，并且位置坐标在给定的平面内是已知的。建立个体鱼的运动模型就是 要研究个体鱼在下一个时刻t+1（1表示一个时间步长）时刻运动的方向和在坐 标平面内的位置。3

13、中的四条规则对改变鱼下一时刻游动方向都起一定作用，那 么在本文中把这四个方向的平均方向作为鱼下一时刻游动的方向。取个体鱼A1研究其运动。由于方向就是与水平方向的夹角，因此仅仅需要 对上述四个方向与水平方向的夹角进行平均，即为下一时刻该鱼的游动方向。用公式表示为:P = P + P + P + Pt+11t 2t 3t 4t（1）其中Pt+1为下一个时间步长周期个体鱼的运动方向，Pt为本周期t中A1 的游动方向，P2t为周期t当前个体到邻居平均位置的方向，Pt为邻居的平均方 向，乌为小于碰撞距离的邻居到当前个体方向的平均值（见图3）。考虑到各规则对鱼的影响力不同，个体鱼作为自主决策的自驱体在决策

14、时 考虑的先后级不同，所以还需要对各个方向加权，取加权平均值，权重的大小文 中可以根据偏好确定。则公式转化为：P =XP +X P +X P +X P（2）t+1 1 1t2 2t3 3t 4 4t其中七+气+七+人4 =1。下面给出四个规则所代表的四个方向的具体实现方法:（1）惯性规则的实现：P1t为本周期t中A的游动方向，本周期内的运动方向由上一时刻的运动方向所确定，即Pt = P。（2）靠近规则的实现：每个个体都有向邻居中心靠拢的特性，邻居中心为观察范围内各个体所在位置的平均值。假设当前A1所处的位置为D （x , y ），D （x , y）为当前各个邻居的位置， 则邻居平均位置 000

15、i i iz d( x, yP2t为周期t当前个体到邻居平均位置的方向，则d g,y=咨一I (& n N（3）P = arctan Z%2tx - x0（3）对齐规则的实现：个体会和它的邻居朝同一个方向游动。公式表示P3t为：（4）其中P为各个邻居的方向,N为邻居的个数，P3t为邻居的平均方向。（4）避免碰撞规则的实现:当个体和它的邻居靠的太近时（距离小于碰撞距 离），应自动避开。公式表示为：z arctan 0_P广x0 一京（D0D ）其中P4t为小于碰撞距离的邻居到当前个体方向的平均值,M为邻居中小于碰撞 距离的邻居个数。最后根据以上计算可得出个体的位置迭代规则为：x = x + V

16、- cos Py = y + v sin Pt + 1tt5.1.2.2局部运动模型（1）个体间的感知感知问题前以述及，个体鱼在鱼群中相互感知，二维空间内，感知范围为半 径为R的圆面，通过个体鱼之间的通信，传递位置和运动信息以及其他信号。（2）通信机制通信范围：与感知范围一致。通信对象：位于4的通信范围内的所有个体鱼。通信过程：当4的位置坐标发生变化后，4立即更新其通信范围内的通信 对象；能及时将所有通信对象；所有通信对象能将它们的位置、运动等信息传递 给先.。通过先.与通信对象之间的通信，实现个体鱼在群体中的局部交互。（3）鱼群的局部交互鱼在群体中排列十分紧密，但有时拥挤会造成身体遮挡了部分

17、感知范围，但 是离散的鱼却能组合成连贯的鱼群结构，实现全局的统一。但是，如果假设每条 鱼能够与其他所有个体进行交互，获取全局信息，那么所有鱼都将集聚在鱼群的 质心，而不是一种均匀的分布。因而可以推断在鱼群中只存在局部交互。局部的 交互，经过反复协调，最终达成全局的和谐结构。鱼群运动是一个和谐的整体，而局部与局部是在整体中相交互联的。而且个 体之间的交互和局部之间的交互都是在不断地更新中，所以鱼群的整体运动建立 在局部交互的基础之上。5.1.2.3鱼群运动模型根据前面两点分析，求解出在给定初始状态情况下，每个个体鱼在下一时刻 的运动方向能够根据个体鱼运动模型求得，其所有的位置也能通过迭代给出。按

18、 照上述的方法，所有鱼群中的个体趋于同步运动，最后形成鱼群的一种集群运动。图4个体、局部与集群之间的关系图5.2问题二的模型5.2.1模型的建立认为鲨鱼在攻击时有它的威胁范围R0,假定所有在它根据对资料的分析攻击范围内的个体鱼均能感受到威胁的存在，并立即在它的感知范围内向其它个 体发出告警交互信号。此时对于集群里的个体鱼来说可以分成三类：第一类是在 鲨鱼的攻击范围之内，它要躲避威胁，并且在感知范围内发出告警交互信息同时 它能收到威胁范围内其它个体发出的告警交互信号。第二类是不在威胁范围之内 但是在第一类鱼的感知范围内。第三类鱼既不在威胁范围内又不在第一类鱼的感 知范围内，具体用图形表示如下：图

19、5危险信号的传递考虑鲨鱼的运动：鲨鱼的运动受其运动惯性的影响，并游向其感知范围内与其距离最近的小鱼。则其运动方向P =xP +人P，其中，P = P，t+1 1 1t 2 2t1t t-1cv - VP = arctan m&m 0%g t + 1 。其位置vlt+1=% + v - cos P =y + v - sin P。由此可得鲨鱼的位置更替。在问题一模型的基础上，考虑对于第一类鱼，它在下一时刻的运动方向将受 到六个因素的影响，如下图：鲨鱼 指向邻居中心图6第一类鱼的方向确定用公式可以表示为：（7）P =XP +X P +X P +X P +X P +X Pt+1 1 1t 2 2t 3

20、 3t 4 4t 5 5t 6 6t其中P , P , P , P的含义同问题一，P表示个体鱼为了躲避威胁而以最快方式逃 1t 2t 3t 4t5t逸的方向，表示该个体鱼受其它个体发出的告警信号的平均方向;人，人，人，人，人，人表示鱼在决策时对它们的偏好权重，x+x +人+人+人+人=1123456123456对于第二类鱼来说，它在下一时刻的运动方向受五个因素的影响，如下图:指向邻居中心A1运动方向邻居的方向唇警交互图7第二类鱼的方向确定用公式可以表示为：P =XP +X P +X P +X P +X Pt+11 1t22t33t44t66t（8）各符号含义同上，其中人+人+人+人+人=1。1

21、2346对于第三类鱼来说，它在下一时刻的运动方向仍然受问题一中的四个因素影 响。5.2.2参数解算算法个体鱼在下一时刻的运动方向确定：P = arctan 乙一2（ 9）5t* -七其中（七,七）为鲨鱼当前的位置坐标。Z arctan Zo_k. /P6 =Xo - XyK（尾 K）（10）其中（气,七）表示感知范围内第k条个体鱼的位置坐标，K为该个体鱼受其它个 体发出的告警信号的个数。个体鱼在下一时刻的位置坐标确定：x = x + V - cos Py = y + v sin PI11）5.3问题三的模型5.3.1模型的建立假设鱼群中有一部分个体是信息丰富者，本文仅考虑掌握食物源位置信息 的

22、个体鱼Mi，它们在遵循问题一中提出的游动规则条件下，将主动靠近食物源， 并且把它向食物源游去的信息告知邻居，召集其它个体鱼共同觅食。在问题一个体鱼运动的模型中把Mi向食物源游动的趋势加入其实际游动方向的方程模型 中，用公式表示：(12)P =XP +X P +X P +X P +X Pt+111t22t33t44t7 71入+人+人+人+人=1。12347其中P , P , P , P的含义同问题一，P表示个体鱼向食物源游动的方向， 1t 2t 3t 4t7t人,人，人，人，人表示鱼在决策时对它们的偏好权重 12347图8掌握食物源信息的个体的方向确定对于其它没有食物源位置信息的个体鱼N .来

23、说，它能受到其感知范围内信 息丰富者的召集信息，并趋向这些信息丰富者的实际游动方向的平均方向，追随 它们共同觅食。则它在下一时刻的运动方向受五个因素的影响，如下图：图9未掌握食物源信息的个体的方向确定用公式可以表示为：(8)P =XP +X P +X P +X P +X Pt+11 1t 2 2t 3 3t4 4t8 8t人,人，人，人各符号含义同上，力表示个体鱼受召集信息而跟随信息丰富者的偏12348好。其中X+X +人+人+人=1。123485.3.2参数解算算法信息丰富者向食物源游去的方向:(9)P = arctan XX71ys - yi其中（x , y ）为食物源当前的位置坐标。（x

24、 , y ）为信息丰富者的位置坐标。 c ci i没有掌握食物源信息的个体鱼受召集信息影响的游向：arctanP = i=1(mN. 2211K鲨鱼0123 x0a鲨鱼靠近鱼群的初始状态b部分个体鱼受威胁时的反应运行t s后鱼群的集群效果图12受威胁时鱼群运动的效果6.2.2模型的结果仿真利用Matlab软件对模型进行仿真，结果如下图:(2)部分个体鱼受威胁时的反应(1)鲨鱼靠近鱼群的初始状态(3)运行ts后鱼群的集群效果图13受威胁时鱼群运动的仿真效果6.2.3模型的结果分析当鲨鱼入侵鱼群时，一部分个体鱼首先出现在鲨鱼的威胁范围并立刻 感受危险存在。它们会选择最快的方式逃逸并在感知范围内发出

25、告警信号与其它 所有在其感知范围内的个体鱼进行信息交互。随着部分个体的快速逃逸和个体间的不断信息交互，得知危险的鱼群 数量不断扩大，选择逃逸的个体鱼越来越多，但是它们逃逸时仍然受游动规则的 影响。虽然鱼群在不断地逃逸，但要受到游动规则的约束，在保证安全时仍 然表现出集群运动行为，结果与事实吻合。七、模型的评价7.1模型的优点：1、本文将复杂的动物集群运动行为巧妙地转化为在一个二维坐标系内质点 的运动方向和坐标位置不断变化的数学问题，使得问题简单明了，易于解决。2、文章在建立问题二、三的模型时充分借鉴了问题一的模型，充分研究了 三个问题之间的关联性。3、文章对于模型的结果作了合理的预测，结果能够

26、反映模型的真实情况， 与实际比较接近。7.2模型的不足：1、本文在模型建立之初，为了简化问题，假设鱼群运动的周围环境是一个 二维平面，这与现实是不相吻合的。动物群体的运动实际上是在空间和时间上的 不断变化的复杂的量。文中没有考虑模型在三维立体空间的建立过程，具有一定 的缺陷。2、模型的仿真结果比较粗糙，不是非常满意。八、模型的推广1、将二维平面环境变为三维的立体空间可使模型具有更加广泛的适用性。 在向三维推广的过程中，建立三维坐标轴，模型中运动方向的改变和位置坐标的 迭代均应在三维空间上的横、纵、竖坐标上进行。依此，模型在三维的仿真结果 会更加逼真。图13模型一的三维仿真图2、集群运动行为模型

27、的推广：本文研究的集群运动行为，不存在着集中控 制或领导者，只是个体依据简单的行为规则相互交互促使整个群体突现出一种复 杂的协调和适应能力。可见，这种群智能的涌现，正是由于群体内部存在着大量 的正负反馈机制。将这种模型应用到控制领域，建立多样性控制器，对于存在多 种作动器的系统有很好的借鉴。参考文献1姜启源、谢金星、叶俊，数学模型，北京：高等教育出版社，2003。2周义仓、赫孝良，数学建模实验，西安：西安交通大学出版社，2007。3韩中庚，数学建模方法及其应用，北京：高等教育出版社，2005。4王晓红，基于多Agent的人工鱼群自组织行为研究D，北京：北京科技大 学，2006。5赵 健、曾建潮

28、，鱼群集群行为的建模与仿真，太原：太原科技大学，2009。附录：1、 n=100;x(1,:)=20*rand(1,n);y(1,:)=20*rand(1,n);drc(1,:)=(rand(1,n)-0.5)*2火pi;d=;for k=1:50;%EaaAefor i=1:nfor j=1:nd(i,j)=sqrt(x(i)-x(j)八2+(y(i)-y(j)八2);endendfor i=1:nd(i,i)=inf;end%EUiE-IopAiuBAfor i=1:na=0,b=1,c=0,e=0,h=0,l=0;for j=1:nif d(i,j)5a=a+1;c=c+x(j),e=e+y(j);m=atan(y(j)/x(j);l=l+m;endif d(i,j)80 x(k+1,i)=x(k+1,i)-8