2021年美赛A题思路详解

2021年美赛A题思路详解2021年数模美赛A题思路详解题目分析思路详解由于和队友思路不一致,导致最后我的思路只算了前两问,而后几问用了我认为离题的PCA(主成分分析)的方法,我的建模思路没有得到完全实现,总体情况很不满意,特此写下这篇文章。题目分析从题目前面所提供的背景知识可以看出,C指出分解速率与菌丝伸长速率成正相关关系,我队友认为是线性关系而我认为是对数近似的关系。第二长图给了一个正比的关系,但是坐标却很容易理解错。这个moisturetrde-off不是湿度耐受性(moisturenichewidth),更

大家好,又见面了,我是你们的朋友全栈君。

2021年数模美赛A题思路详解

由于和队友思路不一致,导致最后我的思路只算了前两问,而后几问用了我认为离题的PCA(主成分分析)的方法,我的建模思路没有得到完全实现,总体情况很不满意,特此写下这篇文章。

题目分析

在这里插入图片描述

在这里插入图片描述
从题目前面所提供的背景知识可以看出,C指出分解速率与菌丝伸长速率成正相关关系,我队友认为是线性关系而我认为是对数近似的关系。第二长图给了一个正比的关系,但是坐标却很容易理解错。这个moisture trde-off 不是湿度耐受性(moisture niche width),更不是湿度的忍耐性。图表中描述这个是物种的竞争能力分数(competitive ranking)与湿度耐受范围(moisture niche width)的差值。如果还不能理解的话,向大家推荐一篇研究内容近似的nature的文章(也是我建模数据的主要来源) Consistent trade-offs in fungal trait expression across broad spatial scales 在这里我们能看到相比于原文引用的文章对湿度耐受性和各种trade-off更准确的定义:
在这里插入图片描述

观察上面这个图标,我们还发现温度和湿度的niche width是有上下界的,这说明这个系数并不是单一的对高温或者对湿度的忍耐度,而是对温度、湿度容忍力的合适的范围。
而对于竞争分数competitive rank,虽无明确定义,但nature的文章有这样一段表述;在这里插入图片描述
这说明,这个ranking跟各个条件下菌丝的延展速率成正相关的关系。

因此,总体上来看,题目给的这张图和nature的文章都在说明逻辑相同的,也是符合我们常识认知的一件事请:菌类会为了适应更加极端的气候条件而“进化”,放弃部分生长速率。

接下来看题目具体要求:在这里插入图片描述
放眼一瞄,看到了唯一一个加粗的biodiversity,还要说明他的重要性,这才是这个题目的最终目的。再看看前面递进的几问,判断文章重点应该在于论证生物多样性对环境变化的优势,而非环境对生物有怎样的影响

思路详解

第一问,说明多物种对分解率的影响。首先不可避免我就想到了物种竞争,接着想到了论文里说的分解率和菌丝延申率的关系。在空间有界、自然资源分布均匀、竞争均匀发生的空间内,我们可以引入Lokta-Volterra模型(不懂的戳进百度百科康康),由二种物种竞争的式子我们可以容易推出多种物种竞争的式子(以物种1为例,其余物种式子省略):在这里插入图片描述
我们发现我们需要物种内禀增长率,环境容纳量,以及对其他物种的竞争系数。’

(1)内禀增长率关于温度和湿度的函数可由实验测得的查找菌丝生长速度随温度变化曲线拟合获得,实际选用了顶点和边边界点的进行双线性插值,拟合出一个简单二元函数:在这里插入图片描述
(2)单物种环境容纳量也会随着气候变化而变化,这里需要假设环境容纳量与菌丝生长速率成正比(菌丝环境容纳量与气候的关系无相关资料),即菌丝生长速率需受气候影响。否则根据模型,无论气候怎么变化后一项多项式在平衡条件下将恒为0,不会产生随环境变化的波动。
(3)物种间的互竞争系数互为倒数,如a对b是2,表示每个a个体占据2个b个体的位置,则b对a为0.5。在此可以物种间内禀增长率的比值近似作为竞争系数。

将三个系数都化作内禀增长率r1的,构建出实际增长率与物种数量的常微分方程,作为菌类系统,这样基本模型就完成建立了。接下来假设温度与湿度的函数做输入信号,以菌丝生长率/分解率作为输出信号迭代计算即可:
(1)设定温度、湿度函数:三角函数、锯齿波函数、方波函数(占空比<10%或>90%表急剧天气变化)等均可。周期可以一年为准(热带为半年),振幅贴近实际情况选取。如要表示全球气候变暖,则可以选取多几年的周期,加上表示气温升高的一次函数。
(2)设定各物种数量初值:设定为输入温度湿度函数的数学期望时的稳定值域,使波动在开始时不会太剧烈。
(3)根据各物种数量值与气候条件,通过Lokta-Volterra模型计算实际增长率。
(4)各物种数量叠加上实际增长率,得到下一刻(第二天)物种数量。
(5)循环到(3),迭代计算,最终得到增长率曲线。
(6)由增长率曲线通过换算得到分解率曲线。

第三问短期和长期的天气变化,通过改变温度、湿度的函数即可获得多物种/单物种情况下分解率的变化趋势。
第四问要求对每一类物种和能共存的物种的组合在不同的气候下进行分析需要考虑该物种的最佳湿度、最佳温度是不是在干旱~热带雨林某个区域的条件中,如果在的话就说明他适合在这个地方生存。

关于主成分分析:虽然在包括题目引用的多篇真菌研究的论文中,对物种的评估都用到了环境因子主成分分析,但个人认为在实际建模种不需要再做。主成分分析的目的在于评估那种环境对物种的影响更大,但是却完全无法满足题目中对论述生物多样性的要求,只能描绘环境对生物的影响,而不能描绘生物对环境的适应力。

至于对于生物多样性重要的论证,首先需要挑选适合的物种组合,设置温度湿度值在这些物种的最佳生长率对应的温度和湿度之间游走波动。如此则易计算出组合的总体分解率是一定会比绝大部分单个物种的情况下来的高的。因此论证生物多样性能提升总体的生物效率。

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。

发布者:全栈程序员-用户IM,转载请注明出处:https://javaforall.cn/133086.html原文链接:https://javaforall.cn

【正版授权,激活自己账号】: Jetbrains全家桶Ide使用,1年售后保障,每天仅需1毛

【官方授权 正版激活】: 官方授权 正版激活 支持Jetbrains家族下所有IDE 使用个人JB账号...

(1)


相关推荐

  • 【解决方案】“‘javac’不是内部或外部命令,也不是可运行的程序或批处理文件”「建议收藏」

    【解决方案】“‘javac’不是内部或外部命令,也不是可运行的程序或批处理文件”「建议收藏」简述在学着使用Java的命令行来编译java文件的时候,遇到了这个问题Windows操作系统“‘javac’不是内部或外部命令,也不是可运行的程序或批处理文件”下面就讲讲解决方法解决过程根据探究的情况来看,这里其实是没有安装好(或者是安装了之后但是在环境变量上却没有设置好)第一步:在http://www.oracle.com/technetwork/java/ja…

  • Python 判断文件是否存在的三种方法

    Python 判断文件是否存在的三种方法这里将介绍三种判断文件或文件夹是否存在的方法,分别使用os模块、Try语句、pathlib模块。1.使用os模块os模块中的os.path.exists()方法用于检验文件是否存在。判断文件是否存在import osos.path.exists(test_file.txt)#True

  • Pycharm中三种注释[通俗易懂]

    Pycharm中三种注释[通俗易懂]”””注释””””’注释”’\#注释前两种不是真正的注释,知识字符串的另一种表现形式,例如会被用于写__doc__等长字符串的时候起到保留格式的目的,虽然能够完成注释功能,但实际上还是字符串。按ctrl+/即可快速注释代码块ps.转载出处作者已注销账号…

  • 链路层链路发现协议(LLDP)讲解「建议收藏」

    链路层链路发现协议(LLDP)讲解「建议收藏」一、LLDP协议概述 随着网络技术的发展,接入网络的设备的种类越来越多,配置越来越复杂,来自不同设备厂商的设备也往往会增加自己特有的功能,这就导致在一个网络中往往会有很多具有不同特性的、来自不同厂商的设备,为了方便对这样的网络进行管理,就需要使得不同厂商的设备能够在网络中相互发现并交互各自的系统及配置信息。 LLDP(LinkLayerDiscoveryProtocol,链路层发现协

  • 数据挖掘十大算法–Apriori算法

    数据挖掘十大算法–Apriori算法一、Apriori算法概述Apriori算法是一种最有影响力的挖掘布尔关联规则的频繁项集的算法,它是由RakeshAgrawal和RamakrishnanSkrikant提出的。它使用一种称作逐层搜索的迭代方法,k-项集用于探索(k+1)-项集。首先,找出频繁1-项集的集合。该集合记作L1。L1用于找频繁2-项集的集合L2,而L2用于找L2,如此下去,直到不能找到

  • mysql主从和mycat读写分离的安装及验证「建议收藏」

    mysql主从和mycat读写分离的安装及验证「建议收藏」mysql主从和mycat读写分离的安装及验证

发表回复

您的电子邮箱地址不会被公开。

关注全栈程序员社区公众号