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

集群动物运动的研究和模拟仿真 摘 要 在自然界里 我们经常能够看到某些动物的集群运动行为 比如鱼群的觅 食 躲避危险 鸟群的迁徙等这些高度一致性的行为 这些群体当中的个体的 行为都是相对比较简单的 但是每个个体只需要遵循某种规则后 整个群体就 涌现出高度的群体智能行为 本文主要为了探讨其中的运动机理与规则建立了 相应的模型 对于模拟动物集群运动 我们先抛开具体的物种和运动形式 并把连续运 动进行离散化 构造了某时刻群体的状态矩阵用来表示所有个体的位置和速度 1 11 2 22 1 11 3 S t n nn n nn n xyv xyv xyv xyv 通过建立个体的距离约束方程 速度约束方程 位置约束方程和状态转换条件进 而求出任意时刻的群体状态矩阵 根据状态矩阵就可得到群体的运动规律 鲨鱼捕食鱼群 是一个无领导者的模型 在模拟鲨鱼捕食一问中 经过对 视频的分析 我们将问题归结为小鱼选择最佳的躲避速度 大小和方向 引用 最优化思想建立目标函数 1 1 mincos 1 1 1 i di i P tD t vtv t P tD t 从而确定躲避速度 再根据鲨鱼和小鱼的初始状态以及速度确定了鲨鱼和鱼群 的运动规律 根据该规律进行 Matlab 编程模拟 我们模拟出了 鲨鱼被包围 的情形 在群体中含有信息丰富者情况下 可以将信息丰富者看做群体的领导者 建立了 leader follower 模型 根据条件 1 Q wt t wa v v Q 1 N Q i i ia v v NQ 1 wnwaiavvv cos waiavv 确定了 leader 的运动 进而确定整个群体的运动 最后我们通过编写相关的程序来仿真群体的运动 从而检验模型的正确性 关键词 离散化 状态矩阵 最优化 leader follower 模型 1 1 问题重述问题重述 在动物界 大量集结成群进行移动或者觅食的例子并不少见 这种现象在 食草动物 鸟 鱼和昆虫中都存在 这些动物群在运动过程中具有很明显的特 征 群中的个体聚集性很强 运动方向 速度具有一致性 通过数学模型来模 拟动物群的集群运动行为以及探索动物群中的信息传递机制一直是仿生学领域 的一项重要内容 1 建立数学模型模拟动物的集群运动 2 建立数学模型刻画鱼群躲避黑鳍礁鲨鱼的运动行为 3 假定动物群中有一部分个体是信息丰富者 如掌握食物源位置信息 掌握 迁徙路线信息 请建模分析它们对于群运动行为的影响 解释群运动方向决策 如何达成 2 2 模型假设模型假设 1 假设每个个体的感知范围极限是相同的 2 部信息的强度与传播的距离成反比 3 在短时间内个体的运动是匀速直线的 4 每个个体接受到信息到做出相应的动作的时间是相等的 5 每个个体的最大移动速度是相等的 6 在时间不长内黑鳍礁鲨鱼的速度基本上是不变的 t 3 3 符号说明及意义符号说明及意义 表示个体 i 其中表示横坐标 表示纵坐标 表示 iiii Px y v i x i y iv 速度 表示个体 i 感知范围内个体 k 的速度 kv 表示个体之间能够容忍的最小距离 min r 表示个体能够感知的极限距离 L R 表示 x 方向的位置变动因子 x 表示 y 方向的位置变动因子 y 表示运动速度最大偏差因子 max 表示个体的最大移动速度 max v 表示危险源的速度 比如黑鳍礁鲨鱼 dv 表示速度大小状态因子 1 表示速度方向状态因子 2 个体 t 的自主速度 wtv 表示核心群体下一步速度 wnv 表示与的方向相近程度 wav iav 4 4 问题分析及模型的建立问题分析及模型的建立 4 14 1 模型一的建立模型一的建立 4 1 14 1 1 问题分析问题分析 动物的集群运动会因为不同的物种 种群和环境而存在不同点 但是从本 质上集群运动是有共同点的 不考虑具体的哪种动物和具体的集群运动 从整 体上分析 动物的集群运动可以分为三种 1 觅食运动 2 躲避捕食运动 3 随机运动 没有明确目的的集群运动 群中的每个个体在综合外部信息和内部信息后自己决定如何进行反应动作 外部信息是让他们产生集群运动转变的因素 比如食物位置 捕食者出现和运 动规律等 内部信息是群体中各个个体之间进行传递交流的信息 包含信息发 出者的位置 运动状态 发出的声音 气味和特殊物质等 在集群运动中 每个个体都遵循相同的动作反应规律 因为个体的感知范 围有限 为了保持群体运动的整体性 个体的运动还要遵循下三个约束条件 1 与邻近的个体之间保持必要的距离 距离太近影响个体的运动 距离太远 导致个体容易脱离群体 2 感知范围内的个体运动速度具有一致性 3 感知范围内的个体位置具有中心一致性 4 1 24 1 2 集群运动的离散化处理集群运动的离散化处理 集群运动在时间和空间上是连续的 为了简化描述集群运动 从时间的角 度把该运动进行离散化 每一时刻 该群都有一种状态 把状态按照时间的顺 序排列就可以近似描述集群的运动 在每一时刻 群中的个体必然要满足上面 的三个约束条件 所以我们也把上面两个条件成为状态条件 根据状态条件可 以建立描述集群运动的状态方程 我们首先只靠虑集群在二维平面的运动 在此我们定义了这 iiii Px y v 样一个变量来描述运动的个体 其中表示个体 i 的在世界坐标系中的横坐标 i x 表示个体 i 在世界坐标系的纵坐标 表示个体 i 的运动的方向 进一步简化 i y iv 问题 我们将运动的个体投影到地面就可以把个体看作是点的运动 为此我们 模拟地面上点的运动进而可以推测出个体的运动 个体投影到地面的表示如下图 现在讨论有 n 个个体的集群的运动情况 于是我们建立了一个由 n 个个体的 t 时刻状态组成的状态矩阵 S t 1 11 2 22 1 11 3 S t n nn n nn n xyv xyv xyv xyv 矩阵中表示 i 个体的速度 iv 依据状态条件建立状态方程 距离条件 1 2 22 minikikL rdxxyyR 1 2mk 3 位置中心一致性 1 1 1 1 m k k ix m k k iy x x m y y m 1 2mk 3 速度一致性 1 1 m k k i v v m max 0 其中表示个体之间能够容忍的最小距离 表示个体能够感知的极限 min r L R 距离 m 表示 i 个体感知范围内相邻个体的数量 表示 i 个体运动速度 iv 表示个体 k 的运动速度 表示 x 方向的位置变动因子 表示 y 方向的 kv x y 位置变动因子 表示运动速度最大偏差因子 max t 1 时刻的状态矩阵与 t 时刻具有相同的形式 为此我们只需要知道 i 个 体的速度变化就可以求解 t 1 时刻的状态矩阵 位置变换表达式 速度大小更新式 1 i i tiii tii Px yPx yvt 1 iii v tv tv 1max ii vvv t 1 0 1 加速取正 速度方向更新式 1 iii tt 2 2 i 2 0 1 逆时针改变方向取正 上式中是 t 时刻与 t 1 时刻的间隔时间 表示个体的最大移动速度 t max v 表示个体 i 速度大小的变化值 表示个体 i 速度方向的变化值 表示 i v i 1 速度大小状态因子 表示速度方向状态因子 2 在随机运动过程的中 都是 0 1 范围内的一个随机数 在觅食 1 2 过程中个体的速度方向改变不大 主要是速度大小的改变 于是我们取 1 0 6 在躲避危险的过程中速度大小 方向均有较大的改变 的取值 2 0 1 1 2 满足使其到达最安全的地方 我们发现不管外界的环境如何 某一状态下集群总是一个整体 所以我们 根据这个规律建立任意一状态的方程 并推导出了 t 状态和 t 1 状态之间的关 系 依此就可以像滚雪球一样得到所有状态的矩阵 进而就可得到集群的运动 情况 模拟仿真结果如下 程序代码见附录程序一 初始位置分布情况 10 个步长后的分布情况 50 个步长后的分布情况 100 个步长后的分布情况 500 个步长后的分布情况 4 2 模型二的建立 4 2 1 问题分析 关于鱼群躲避黑鳍礁鲨鱼 从视频中可以看出 在水表面 鲨鱼进入鱼群 时 最先感觉到捕食者迫近的小鱼 靠外层鲨鱼前方 向两边分散游开 速度 很快 距离鲨鱼稍远的鱼则反应一般 只是跟随旁边一侧的鱼游动 速度一般 鲨鱼两侧的鱼在向鲨鱼的后方游去 速度一般 鲨鱼身后的鱼已经重新聚群 鲨鱼身后的空隙由前方游回来的鱼填充 看起来是鲨鱼在鱼群里游动 鱼群采 取分散向后转移的方式躲避鲨鱼的捕食 鲨鱼好像被鱼群包围 如下图 1 图 1 鲨鱼被包围示意 包围圈周围的鱼的运动是我们研究的主要内容 为了研究这个运动 我们 需要的条件有 1 鲨鱼和鱼群中的小鱼耳的初始位置 2 鲨鱼的运动包含速度大小和方向 3 小鱼的速度变化规律 包含大小和方向 这可以帮助确定小鱼某时刻 的位置 4 小鱼的感知和交流方式 这可以帮助我们确定小鱼如何根据捕食者的 运动进行反应 鲨鱼和鱼群运动之间的影响关系如下图 2 图 2 鲨鱼的运动和小鱼的运动之间的影响关系 说明 1 鲨鱼速度大小决定外部信息信号的强弱 而鲨鱼和小鱼的距离决定小鱼所 接收到的信号的强度 小鱼根据接收到的信号强度决定躲避快慢 2 鲨鱼的运动方向和鲨鱼与小鱼之间的位置关系共同决定了小鱼的躲避方向 即如 鲨鱼被包围示意图 中表示的小鱼向左或向右跑和前进的方向 3 鲨鱼的位置由上一时刻的位置和速度决定 小鱼同理 4 2 2 模型建立 根据上述说明我们建立平面直角坐标系来表示 xoy 图 3 平面坐标系中鲨鱼和小鱼个体的运动关系 问题最终归结为为小鱼找到一个合适的躲避速度 到达一个新的位置 使 自己最为安全 满足以下两个条件 目标函数 即可 躲避方向 mincos 1 1 divtv t 下一时刻位置 max 1 1 i P tD t 这是双目标问题 为此我们采用加权系数化为单目标问题求解 mincos 1 1 1 1 1 di i vtv tP tD t 为权重系数 其中认为在短时间内 t 1 DDvtvt 在时间之后t 小鱼的位置 11111 P x yP x yVt 鲨鱼的位置 0 D x yD x yVt 小鱼的速度变为 鲨鱼的速度变为 1 V 0 V 之后的一个中鲨鱼和小鱼以的初状态状态进t 0111 VVD x yP x y 行运动 并且一直如此循环下去 直到鲨鱼离开鱼群 鱼群不再受到被捕食的 威胁 对于其中的一些参数取值的选取 在进行 matlab 模拟时可以适当选取 如 果模拟结果比较符合实际情况则认为参数值的选取是合适的 模拟仿真结果如下 程序代码见附录程序 2 图 4 鲨鱼捕食模型的 Matlab 模拟结果 说明 我们模拟鲨鱼作匀速直线运动进行捕食 红色代表鲨鱼运动轨迹 每个 点代表一个时刻的鲨鱼位置 并且刻画了六条小鱼的躲避轨迹 六条蓝色轨迹 代表 起点为每条轨迹的点稀疏的一端 点越密表示小鱼的速度越小 变化不 大 从模拟结果图中我们可以看出鱼群的运动基本上符合视频中的真实运动特 点 小鱼起初躲避速度很快 图上表示为小鱼轨迹点稀疏 向鲨鱼前进方向 的两侧躲避 之后慢慢减速 在运动到鲨鱼身后时 基本上就不在躲避运动 表现为小鱼运动轨迹末端点的分布密集 所以我们认为模型的模拟效果是比 较好的 4 3 模型三的建立和求解 在群体中包含信息丰富者时 模拟信息丰富个体对群体运动行为的影响 我们分析认为 若动物群体中有一部分个体是信息丰富者 那么这些个体就可 以看作是整个群体的领导者 因为他们掌握了重要的信息 比如食物的位置信 息 他们本身的行为就表现的很有方向性 目标明确 而其他个体只能通过 他们传递的信息知道有相关的信息 但是并不知道信息所反映的实际情况是什 么 所以他们只会跟随这些信息丰富者的运动或者与自己相邻个体的运动而运 动 在寻找食物或迁移的过程中整个群体是跟随信息丰富个体运动的 而且信 息丰富个体都掌握着信息 所以他们本身的运动具有很高的一致性 在此我们 将这些信息丰富个体看做一个整体 定义为核心群体 其他个体定义为跟随群 体 问题三就可以简化成整个群体随核心群体的运动 于是我们建立一个 leader follower 的模型 重点研究核心群体的运动方向 核心群体的运动方向是 由两方面的影响决定的 一是自主运动方向 就是他们根据自己掌握的信息运 动 二是被动运动方向 就是他们的运动还会受到跟随群体运动的影响 假设群体的总量为 N 核心群体中个体的数量为 Q 他们中个体 t 的自主速 度为 跟随群体中个体的数量为 N Q 他们中个体 i 的速度为 由此我们 wtv iv 建立以下关系 1 Q wt t wa v v Q 1 N Q i i ia v v NQ 1 wnwaia vvv 其中 cos waiavv 式中表示核心群体中的平均速度 表示跟随群体的平均速度 表 wav iav wnv 示核心群体下一步速度 表示与的方向相近程度 1 表示与 wav iav wav 同向 iav 另外整个跟随群体是满足问题一中状态矩阵中的随机过程 5 模型的改进与推广 一 改进 鉴于问题的复杂性我们只考虑了群体在二维平面的运动情况 在此我们可 以根据已经建立好的模型将其延伸到三维空间 这样更加符合实际情况 在考 虑小鱼躲避鲨鱼时应该将危险信号强度与躲避速度的关系考虑进来 二 推广 通过对群体集群运动的研究 在过程方面 生物群体的同步 躲避危险的 机制可以应用到分布式机器人集群 无人驾驶飞行器 卫星群的运动控制 在 信息管理方面 可以从生物体如何形成有效决策的得到启示 为管理机制的改 进和管理效率的提高提供新思想 6 模型的评价 一 模型的优势 1 应用离散化的概念 将集群的连续运动化为间隔时间很短的静止状态影像 2 用矩阵简单明了的描述了每一个静止的状态 3 只考虑两个状态之间的联系 不需要分析其中的转换过程 4 引用了最优化的概念有效的解决了躲避运动速度的确定 二 模型的缺点 1 问题三中将整个群体划分为两类 分别作为整体研究简单 但是忽略了个 体之间的差异 2 模型中某些参数的确定带有较强的主观性 7 参考文献 1 沈继红 施久玉等 数学建模 哈尔冰 哈尔冰工程大学出版社 1996 2 郑家茂 数学建模基础 南京 东南大学出版社 1997 3 赵建 鱼群集群行为的建模与仿真 太原 太原科技大学出版社 2005 4 刘秉泉 基于 Vicsek 模型的导弹集群建模与仿真研究 电子信息与电气 工程学院 5 张仲秋 周应琪等 鱼群群体运动的数学模型 上海 上海海洋大学海 洋学院 2011 6 江铭炎 袁东风等 改进的人工鱼群算法 济南 山东大学信息科学与 工程学院 8 附录 程序一 问题一的模拟仿真代码 format long Visual 10 N 100 Gen number 1000 a1 10 b1 10 a2 10 b2 10 maxstep 0 3 dmin 1 dmax 3 k 0 delta max pi 3 180 位置初始化 X1 rand N 1 b1 a1 a1 X2 rand N 1 b2 a2 a2 X1 next X1 X2 next X2 X X1 X2 方向初始化 alpha rand N 1 2 pi while k Gen number 感应伙伴 m zeros 1 N X X1 X2 for i 1 N XX X i 1 X i 2 d clofrd i 100 friend for j 1 N if i j continue end XXX X j 1 X j 2 if norm XX XXX Visual friend numel friend 1 j m i m i 1 end if norm XX XXX d clofrd i clofrd i j d clofrd i norm XX XXX end end mean alpha i mean alpha friend if m i 0 CX1 i 0 CX2 i 0 else CX1 i mean X1 friend CX2 i mean X2 friend d center i sqrt CX1 i X1 i 2 CX1 i X1 i 2 mean alpha i mean alpha friend end end 行为 for i 1 N if m i 0 alpha next i rand 2 pi step maxstep X1 next i step sin alpha next i X1 i X2 next i step cos alpha next i X2 i end if d clofrd i dmin if d clofrd i 0 01 alpha next i rand 2 pi step rand maxstep dmin dmin X1 next i step sin alpha next i X1 i X2 next i step cos alpha next i