免费项目网站,定陶网站建设,失信被执行人查询系统,北京专业做网站2018年认证杯SPSSPRO杯数学建模
基于聚类分析的海豚捕食合作策略
A题 海豚与沙丁鱼
原题再现#xff1a; 沙丁鱼以聚成大群的方式来对抗海豚的捕食。由于水下光线很暗#xff0c;所以在距离较远时#xff0c;海豚只能使用回声定位方法来判断鱼群的整体位置#xff0c;难…2018年认证杯SPSSPRO杯数学建模
基于聚类分析的海豚捕食合作策略
A题 海豚与沙丁鱼
原题再现 沙丁鱼以聚成大群的方式来对抗海豚的捕食。由于水下光线很暗所以在距离较远时海豚只能使用回声定位方法来判断鱼群的整体位置难以分辨每个个体。鱼群的行动是有协调性的在没有外部威胁或障碍物时鱼群常常会聚成接近球形的形态。而当海豚接触甚至冲进鱼群鱼群则会进行协同的躲避所以不易在大鱼群中追踪一个目标。沙丁鱼的这种群体行为降低了其被海豚捕食的概率。 第二阶段问题 智商较高的海豚会配合行动以图提高捕食的效率。请你建立合理的数学模型研究多条海豚配合捕食时应当采取怎样的合作策略。
整体求解过程概述(摘要) 生态系统中无论作为捕食者还是被捕食者好多生物出于生存、觅食、避险和繁殖等原因通常采用集群运动行为。沙丁鱼群和海豚之间的对抗就属于典型的群体行为博弈。第一阶段建立的基于沙丁鱼个体之间疏密度控制的沙丁鱼群逃逸模型被运用到第二阶段中沙丁鱼群的逃逸策略不仅要避免被海豚捕食同时要拉开与沙丁鱼彼此的间距采用了基于最远逃逸距离和最佳逃逸角度相结合的原则进行疏密度控制使得逃逸的沙丁鱼群分布比较分散这增大了海豚的捕食难度。在第二阶段必须制定合适的海豚捕食策略才能实现对分散逃逸、稀疏分布的沙丁鱼群的高效捕食。第二阶段模型主要工作内容在以下两个方面 其一、研究海豚群体行为模式模拟海豚集群协作捕食的场景设定了三种协作阵型圆形包围阵型、三角进攻阵型和混合围追堵截阵型。圆形包围阵型实施对沙丁鱼向中心驱赶沙丁鱼群的最远逃逸距离原则反而使得沙丁鱼群由四周向内聚拢进一步被海豚合围压制因此在圆形包围阵型在捕食初期能制约鱼群逃逸获得很好的捕食效果。三角进攻阵型是从一个方向朝鱼群进攻由于鱼群的较好的自适应逃逸策略通过优化最远逃逸距离和最佳逃逸方向总能避开海豚的追击。混合型阵型结合了两种阵型的特点一部分海豚在一侧追击另一部分在令一侧堵截。通过仿真对三种阵型捕食效果进行对比在同等条件下圆圈包围阵型对沙丁鱼平均捕获率为 16%三角进攻阵型适合驱赶沙丁鱼群平均捕获率为 0混合围追堵截阵型的平均捕获率为 9%。因此圆形包围阵型的效果最好。 其二研究海豚个体行为模式制定了两种海豚个体的捕食策略沙丁鱼近邻策略和沙丁鱼群密度聚类分析策略。近邻策略是指在海豚感知范围内其捕食目标和行动路线始终朝向最近的沙丁鱼聚类分析策略是指利用基于密度的聚类算法DBSCAN对在海豚的有效感知区域内散布的沙丁鱼进行聚类从而判断沙丁鱼聚集的高密度区使其捕食目标和行动路线朝向鱼群高密度区的中心。仿真结果显示聚类法捕获的沙丁鱼数量比近邻法高出 80%在捕食效果上基于密度的聚类法优于近邻策略。通过对海豚多组合作策略对比分析本文得出结论基于聚类法的海豚个体捕食策略结合圆形包围的海豚集群阵型是较好的合作捕食策略该合作策略对进一步建立类似的高效捕食体系仿真模型具有参考价值。同时为三维空间海豚配合捕食运动建模奠定基础。
问题分析 第二阶段在沿用第一阶段背景的基础上讨论针对沙丁鱼群协同避险行为和鱼群的逃逸运动规律研究多条海豚配合捕食时的合作策略。针对多条海豚捕食合作一般分为以下几个阶段 第一阶段当海豚发现一定距离范围有沙丁鱼群提前实施相应的围攻策略形成一定阵型为进一步进攻捕食做足准备。 第二阶段,海豚已经靠近目标鱼群针对海豚进攻时刻周围鱼群的分布情况采取的进攻捕食方向策略。 基于上述两个阶段的海豚捕食过程的分析大致可以将海豚的合作策略分为围攻策略和进攻捕食策略。首先通过查找大量的相关文献本文针对不同围攻策略建立三种排布阵型1圆形阵型2三角阵型3组合阵型即圆形和三角组合。其次考虑到海豚周围鱼群不同的分布情况海豚的进攻捕食方向有两种可能选择:1)向距离最近的沙丁鱼进攻二、向密度较大的逃逸的鱼群展开追逐捕食。为了使问题清晰化本文认为两种进攻过程中海豚应当具有对处在其必定捕杀距离范围内的沙丁鱼一击必杀的能力并且其猎杀过程的时间与整个追逐沙丁鱼群的时间相比可忽略。同时海豚能够利用回声定位对沙丁鱼群保持时刻的整体感知其感知所耗费的时滞几乎不计。当海豚和沙丁鱼距离小于海豚捕食步长时捕杀必定成功。通过分析对比海豚不同的捕食规律所对应的捕食沙丁鱼群的成功率来优化海豚捕食的合作策略。
模型假设 1、鱼群和海豚都活动在足够大的水域中不受地形和海平面的影响 2、海豚的感知范围为以自身为中心半径为 Sense 的圆域简化处理将沙丁鱼和海豚视为质点处理; 3、海豚对沙丁鱼群的每一次进攻捕食都可以简化为直线运动 4、所有的海豚都具有相同的物理特性包括运动能力和感知能力 5、海豚与沙丁鱼距离小于海豚运动步长时海豚必定捕杀成功 6、海豚一次只能吃掉一条沙丁鱼此假设是基于俄勒冈大学和夏威夷大学共同研究成果制定因为海豚不能像蓝鲸一样张大嘴巴一次吞进大量的食物它们需要以单条鱼为目标 7、海豚初始按照阵型规律前进
论文缩略图 全部论文请见下方“ 只会建模 QQ名片” 点击QQ名片即可
部分程序代码(代码和文档not free)
%鲨鱼四周围捕一起进攻
clear all;
clc
M50;N100;ls20; Ns4; h20; Sense150;%ls 鲨鱼步长;h 鲨鱼开始走直线的距离,Ns 鲨鱼
数量
d0[99 99 99 99];
S(1,1)-20;S(2,1)-20;
S(3,1)120;S(4,1)120;
S(5,1)120;S(6,1)-20;
S(7,1)-20;S(8,1)120;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%初始小鱼的聚集球状%%%%%%%%%%%%%%%%%%%%%%%%%%%
x050;y050;O[x0;y0];%圆心坐标
R 10;%圆半径
deta 360.*rand(N,1);%随机角度
r R.*rand(N,1);%随机半径
dx r.*cosd(deta);%计算 x 坐标增量
dy r.*sind(deta);%计算 y 坐标增量
nx x0 dx;%计算随机点 x 坐标
ny x0 dy;%计算随机点 y 坐标
%%%%%%%%%%%%%%%%%%%%%%%%%%%% 鲨 鱼 的 移 动 %%%% 四条集中行动 %%%%% 近 邻 捕
食%%%%%%%%%%%%%%%%%%%%%%%%%
F(1,:)nx;
E(1,:)ny;
FF;EE;
for i1:M;d0(1)((O(1,1)-S(1,i))^2(O(2,1)-S(2,i))^2)^0.5;d0(2)((O(1,1)-S(3,i))^2(O(2,1)-S(4,i))^2)^0.5;d0(3)((O(1,1)-S(5,i))^2(O(2,1)-S(6,i))^2)^0.5;d0(4)((O(1,1)-S(7,i))^2(O(2,1)-S(8,i))^2)^0.5; Dst1sdy2_ss_Dist(S(1,i),S(2,i),F(:,i),E(:,i));Ds(:,1)Dst1;[Dsmin(1),idx(1)]min(Ds(:,1));Dst2sdy2_ss_Dist(S(3,i),S(4,i),F(:,i),E(:,i));Ds(:,2)Dst2;[Dsmin(2),idx(2)]min(Ds(:,2));Dst3sdy2_ss_Dist(S(5,i),S(6,i),F(:,i),E(:,i));Ds(:,3)Dst3;[Dsmin(3),idx(3)]min(Ds(:,3));Dst4sdy2_ss_Dist(S(7,i),S(8,i),F(:,i),E(:,i));Ds(:,4)Dst4;[Dsmin(4),idx(4)]min(Ds(:,4));%NxF(idx);NyE(idx);
%thetasdy2_ss_count(S(1,i-1),S(2,i-1),S(1,i),S(2,i),Nx,Ny);if Dsmin(1)lsS(1,i1)F(idx(1));S(2,i1)E(idx(1));[vol,col]size(F(:,1)); if idx(1)vol;F(idx(1),:)[];E(idx(1),:)[];endelseif Dsmin(1)ls Dsmin(1)Sense[px1,py1]sdy2_point(S(1,i),S(2,i),F(idx(1)),E(idx(1)),Dsmin(1),ls); S(1,i1)px1;S(2,i1)py1;elseS(1,i1)S(1,i)ls*abs((O(1,1)-S(1,i)))/(sqrt((O(1,1)-S(1,i))^2(O(2,1)-
S(2,i))^2));S(2,i1)S(2,i)ls*abs((O(2,1)-S(2,i)))/(sqrt((O(1,1)-S(1,i))^2(O(2,1)-
S(2,i))^2));endif Dsmin(2)lsS(3,i1)F(idx(2));S(4,i1)E(idx(2));[vol,col]size(F(:,1)); if idx(2)~idx(1) idx(2) idx(2)vol;F(idx(2),:)[];E(idx(2),:)[];endelseif Dsmin(2)ls Dsmin(2)Sense[px2,py2]sdy2_point(S(3,i),S(4,i),F(idx(2)),E(idx(2)),Dsmin(2),ls); S(3,i1)px2;S(4,i1)py2;elseS(3,i1)S(3,i)-ls*abs((O(1,1)-S(3,i)))/(sqrt((O(1,1)-S(3,i))^2(O(2,1)-
S(4,i))^2));S(4,i1)S(4,i)-ls*abs((O(2,1)-S(4,i)))/(sqrt((O(1,1)-S(3,i))^2(O(2,1)-
S(4,i))^2));endif Dsmin(3)lsS(5,i1)F(idx(3));S(6,i1)E(idx(3));[vol,col]size(F(:,1)); if idx(3)~idx(2) idx(3)~idx(1) idx(3)vol;F(idx(3),:)[];E(idx(3),:)[];endelseif Dsmin(3)ls Dsmin(3)Sense[px3,py3]sdy2_point(S(5,i),S(6,i),F(idx(3)),E(idx(3)),Dsmin(3),ls);
S(5,i1)px3;S(6,i1)py3;elseS(5,i1)S(5,i)-ls*abs((O(1,1)-S(5,i)))/(sqrt((O(1,1)-S(5,i))^2(O(2,1)-
S(6,i))^2));S(6,i1)S(6,i)ls*abs((O(2,1)-S(6,i)))/(sqrt((O(1,1)-S(5,i))^2(O(2,1)-
S(6,i))^2));endif Dsmin(4)lsS(7,i1)F(idx(4));S(8,i1)E(idx(4)); [vol,col]size(F(:,1)); if idx(4)~idx(3) idx(4)~idx(2) idx(4)~idx(1) idx(4)vol;F(idx(4),:)[];E(idx(4),:)[];endelseif Dsmin(4)ls Dsmin(4)Sense[px4,py4]sdy2_point(S(7,i),S(8,i),F(idx(4)),E(idx(4)),Dsmin(4),ls); S(7,i1)px4;S(8,i1)py4;elseS(7,i1)S(7,i)ls*abs((O(1,1)-S(7,i)))/(sqrt((O(1,1)-S(7,i))^2(O(2,1)-
S(8,i))^2));S(8,i1)S(8,i)-ls*abs((O(2,1)-S(8,i)))/(sqrt((O(1,1)-S(7,i))^2(O(2,1)-
S(8,i))^2));end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%沙丁鱼逃逸%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
[N,nl]size(F(:,1)); numzeros(N,M);
Dszeros(N,Ns);Dst1zeros(N,1);Dst2zeros(N,1);Dst3zeros(N,1);Dst4zeros(N,1);
Df12; Dm7; ohmzeros(N,M);dnzeros(N,M); SSzeros(2,M);%df 沙丁鱼步长for j1:Nmode(1)((F(j,i)-S(1,i))^2(E(j,i)-S(2,i))^2)^0.5;mode(2)((F(j,i)-S(3,i))^2(E(j,i)-S(4,i))^2)^0.5;mode(3)((F(j,i)-S(5,i))^2(E(j,i)-S(6,i))^2)^0.5;mode(4)((F(j,i)-S(7,i))^2(E(j,i)-S(8,i))^2)^0.5;[mm,l]min(mode);dn(j,i)mm;SS(1,i)S((2*l-1),i);SS(2,i)S((2*l-0),i);x(F(j,i)-SS(1,i))/mm*Df;y(E(j,i)-SS(2,i))/mm*Df;x0F(j,i);y0E(j,i);if F(j,i)E(j,i);F(j,i1)x*cos(ohm(j,i))y*sin(ohm(j,i))x0;E(j,i1)-x*sin(ohm(j,i))y*cos(ohm(j,i))y0;elseF(j,i1)x*cos(ohm(j,i))-y*sin(ohm(j,i))x0;E(j,i1)x*sin(ohm(j,i))y*cos(ohm(j,i))y0;
endend
FxF(:,i); FyE(:,i);
numtsdy_Dm(N,M,Fx,Fy,Dm);
num(:,i)numt;
Pnumnum(:,i)/N;
ohmtsdy2_judge(Pnum,dn);
ohm(:,i)ohmt;
end
for i1:5:M
figure(i)
plot(S(1,i),S(2,i),ro);grid on;hold on;plot(S(3,i),S(4,i),ro);hold
on;plot(S(5,i),S(6,i),ro);hold on;plot(S(7,i),S(8,i),ro);hold on;axis([-400,600,-
400,600]);
plot(F(:,i),E(:,i),b.);
title([第 num2str(i) 步]);
end;全部论文及程序请见下方“ 只会建模 QQ名片” 点击QQ名片即可
