2024年美国大学生数学建模竞赛A题资源可用性和性别比例解题全过程文档及程序

2024年美国大学生数学建模竞赛

A题 资源可用性和性别比例

原题再现：

  虽然有些动物物种不存在通常的雄性或雌性，但大多数物种基本上都是雄性或雌性。虽然许多物种在出生时的性别比例为1:1,但其他物种的性别比例却偏离了平均值。这就是所谓的适应性性别比例变异。例如，美洲鳄孵卵巢的温度会影响其出生时的性别比例。
  七鳃鳗的作用十分复杂。在一些湖泊栖息地，七鳃鳗被视为对生态系统有重大影响的寄生虫,而在世界上的一些地区，如斯堪的纳维亚半岛、波罗的海地区，以及北美西北太平洋地区的一些土著居民，七鳃鳗也是他们的食物来源。
  海洋七鳃鳗的性别比例会因外部环境而变化。海洋七鳃鳗变成雄性或雌性取决于它们在幼虫阶段的生长速度。幼体的生长速度受食物供应的影响。在食物供应较少的环境中，生长速度会较低，雄性海洋七鳃鳗的比例可达到约78%。在食物比较容易获得的环境中，雄性的比例据观
察约占种群的56%。
  我们的研究重点是性别比例及其对当地条件的依赖性，特别是海灯鱼的性别比例。海洋七鳃鳗生活在湖泊或海洋栖息地，并溯流而上产卵。我们的任务是研究一个物种能够根据资源可用性改变性别比例的利弊。您的团队应开发并研究一个模型，以便深入了解生态系统中由此产生的相互作用。
  需要研究的问题包括以下方面
  ●当七鳃鳗种群的性别比例发生变化时，会对更大的生态系统产生什么影响?
  ●七鳃鳗的数量有哪些优势和劣势?
  ●七鳃鳗性别比例的变化对生态系统的稳定性有何影响?
  ●七鳃鳗种群性别比例可变的生态系统能否为生态系统中的其他生物(如寄生虫) 带来优势?

整体求解过程概述(摘要)

  本文对七鳃鳗生态系统进行了建模，分析了七鳃鳗性比变化对种群和生态系统的影响。首先，我们建立了一个具有年龄组的离散种群动力学模型。根据七鳃鳗的生物学特性，将七鳃鳗种群划分为7个年龄组，并考虑七鳃鳗所生活的生态系统，包括幼虫饵料、幼虫竞争物、成虫饵料、成虫竞争物和人工渔业。该模型考虑了性别比例随外部环境变化的事实。
  任务1主要研究七鳃鳗性别比变化对大生态系统的影响。离散种群动力学模型的结果表明，七鳃鳗性别比的变化有利于种群数量的稳定，使其在生态系统中处于优势地位。通过比较七鳃鳗入侵生态系统前后种群数量的变化，得出七鳃鳗作为入侵物种的性别比例变化对生态系统的本地种群是不利的，但同时，七鳃鳗的幼虫可以净化水质，促进底栖环境中的养分循环。
  在任务2中，我们应用我们的模型，重点研究了七鳃鳗变性比、变态、一次性繁殖和迁移等生物学特性的优缺点。通过离散种群动力学模型计算环境资源变化和外部环境突变的情况，得出变性和变性别比的七鳃鳗种群更稳定，变性七鳃鳗种群更能实现环境资源的充分利用。通过模拟人类活动对七鳃鳗洄游产卵过程的影响和一年内产卵数的急剧减少，得出七鳃鳗的洄游特征和一次性繁殖的缺点是洄游过程可能是危险的，并且七鳃鳗的幼虫依赖于生境，并且给定世代的低存活率使整个种群处于危险之中。
  在任务3中，我们研究了七鳃鳗性别比例的变化对生态系统稳定性的影响。以生态系统的抗性稳定性、恢复力稳定性和可持续性为评价指标。离散种群动力学模型的数值模拟结果表明：当不同性别比的七鳃鳗和浮游生物种群数量发生突变时，生态系统波动度分别为23.28%和33.19%，恢复时间分别为18年和28年。与雌雄比不可变的数值模拟结果相比，雌雄比可变的七鳃鳗种群使生态系统具有较强的抗性稳定性和恢复稳定性。引入人工干预后，不同性别比例的七鳃鳗种群的生态系统是可持续的。
  在任务4中，我们主要利用我们的模型分析了七鳃鳗不同性别比例对其他物种的影响。将七鳃鳗的寄生种群添加到离散种群动力学模型中，模型结果表明七鳃鳗的变性比有利于其他寄生生物和种群。综上所述，七鳃鳗性别比的变化有利于在恶劣环境下保证种群数量的稳定。这反过来又有助于更大生态系统的稳定和可持续性。

模型假设：

  •假设1：我们假设所有七鳃鳗都是在第七年繁殖的，并且父母在繁殖结束时死亡。
  理由：七鳃鳗是一次性繁殖生物，寿命为7年，为了便于计算七鳃鳗繁殖产生的新生幼体数量，我们假设它们全部在第7年繁殖。

  •假设2：我们假设七鳃鳗在进入第五个年头时就完成了性别的改变。
  理由：在我们的模型中，我们按年龄对七鳃鳗进行了分组，并按年度间隔计算了每一年龄的七鳃鳗的种群规模，因此为了便于计算数量，我们考虑了改变七鳃鳗性别的行为是瞬时进行的。

  •假设3：我们假设七鳃鳗没有天敌，即我们不将七鳃鳗视为诱饵，但主要受到人类活动的威胁。
  理由：幼虫七鳃鳗生活在相对隐秘的环境中，在很大程度上受到定期捕食的保护，而成年七鳃鳗由于其体型和习性相对较少受到天敌的威胁。因此，我们不考虑七鳃鳗的捕食者种群，而是增加人工干预项来模拟七鳃鳗种群对人类活动的影响。

  •假设4：我们假设七鳃鳗的寄生虫种群只寄生于七鳃鳗的成熟体。
  理由：现有的七鳃鳗寄生虫研究表明，80%以上的寄生虫存在于七鳃鳗的成熟体种群中，只有11%的寄生虫存在于七鳃鳗的幼体种群中，为了简化模型，我们主要考虑寄生虫对七鳃鳗成熟体种群的影响。

问题重述：

  对于七鳃鳗的性别比例及其对当地条件的依赖性问题，我们开发了合理的数学模型，以深入了解生态系统中出现的相互作用，并回答以下问题：

  •分析七鳃鳗性别比例与其种群规模相关时对生态系统的影响。

  •使用相应开发的模型分析并证明七鳃鳗种群的优势和劣势。

  •分析七鳃鳗性别比例对生态系统稳定性的影响，考虑到其随外部环境而变化。

  •根据七鳃鳗种群性别比的特异性，解释七鳃鳗是否能为生态系统中的其他物种（例如寄生虫）提供优势。

模型的建立与求解整体论文缩略图

全部论文请见下方“ 只会建模 QQ名片” 点击QQ名片即可

代码：(文档not free)

%厄农映射
%x(j+1)=1+yj-a*xj*xj
%y(j+1)=b*xj
clc
clear all
close all
a=0.4;b=0.3;%系统参数设置
xjL=[];yjL=[];
for k=1:140%迭代次数
    a=0.5+(k-1)*0.01;
xj=[];yj=[];
    x=0.1;y=0.1;
for j=1:200
        x=1+y-a*x*x;y=b*x;%系统函数
xj=[xj;j x];  
yj=[yj;j y];
if j>=120
xjL=[xjL;a b x];    
yjL=[yjL;a b y]; 
end
end
end
plot(xjL(:,1),xjL(:,3),'k.')
%坐标轴设置
xlabel('\itt','FontSize',20,'FontName','Times New Roman');
ylabel('\itx','FontSize',20,'FontName','Times New Roman');
grid on;


clc, clear
a=0:0.01:1.6; n=140; S=[];
for j=1:length(a)
    x=0.2; y=0.3; %x,y的初值
    for i=1:n
        x2=1-a(j)*x^2+y; y=0.3*x; x=x2; %首先进行序列迭代
    end
    if x(end)>-100 & x(end)<100 %若不发散，再计算指数
        x=0.2; y=0.3; %设立二阶单位矩阵
        JJ=eye(2);
        for i=1:n
            x2=1-a(j)*x^2+y; y=0.3*x; x=x2;
            J=[-2*a(j)*x, 1; 0.3, 0]; %分别对两迭代式对于x、y求偏导数
            JJ=JJ*J;
        end
        L=eigs(JJ,1); %求模最大的特征值
        S=[S,[a(j);log(abs(L))/n]]; %把a值及指数值保存
    end
end
plot(S(1,:),S(2,:),'.-') %画出a对应的最大Lyapunov指数
hold on, plot([a(1),a(end)],[0,0],'k')
xlabel('\it a'), ylabel('最大Lyapunov指数')

全部论文请见下方“ 只会建模 QQ名片” 点击QQ名片即可