数学建模鱼群竞争

实用文档标准文案2012上海第二工业大学数学建模竞赛承诺书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的,如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。我们参赛选择的题号是（从A/B/C中选择一项填写）：A我们的参赛报名号为（如果赛区设置报名号的话）：所属学校（请填写完整的全名）：上海第二工业大学参赛队员(打印并签名)：1.郑燕燕2.王赞3.指导教师或指导教师组负责人(打印并签名)：日期：2012年5月3日实用文档标准文案竞争鱼群的最优捕获摘要根据目前资源短缺，能源危机的挑战，本文以生态经济着眼，首先用微分方程的思想建立了基本模型。再根据题意，结合Logistic模型，建立了微分方程数学模型，利用Matlab软件编程画出了两种鱼群的数量随时间变化的曲线以及两种鱼群的相轨线，利用图像研究了两鱼群的数量在自然增长、彼此竞争和人类捕捞的情况下，随时间变化的情况。再结合鱼群数量的变化，考虑到经济效益和人与自然的和谐发展，通过改变人类的捕捞方式，来保证鱼群数量的可持续增长，并且在此基础之上，使经济利益达到最大化。关键字鱼群最优捕获Logistic数学模型微分方程Matlab实用文档标准文案目录一、问题重述……………………………………………………………………（3）二、模型假设……………………………………………………………………（3）三、符号说明……………………………………………………………………（4）四、模型建立与分析求解………………………………………………………（4）4.1问题1模型的建立与分析……………………………………………（4）4.2问题2模型的建立与分析……………………………………………（6）4.3问题3模型的建立与分析……………………………………………（7）五、结果分析……………………………………………………………………（10）六、模型评价与改进……………………………………………………………（11）七、参考文献……………………………………………………………………（11）八、附录…………………………………………………………………………（11）8.1问题1的程序代码……………………………………………………（11）8.2问题2的程序代码……………………………………………………（12）8.3问题3的程序代码……………………………………………………（12）实用文档标准文案一问题重述捕鱼与人类的生活品质，社会经济的发展以及人与自然的和谐发展息息相关。某水域中生活着两种鱼类，它们的数量分别记为x(t)和y(t)，内禀增长率分别为r1、r2，最大容纳数量为N1、N2。它们的增长速率除了受自身条件的限制外，还与鱼类之间的竞争和人类的捕捞相关。竞争减少的增长率与两种鱼类数量的乘积成正比，比例系数分别为q1、q2。由于人类的捕捞使得鱼减少的增长率与鱼的数量成正比（比例系数为分别为k1、k2）。作出一些简化假设，用微分方程模型描述这种捕捞方式下鱼群数量的变化过程。根据两种鱼群的数量变化过程，结合当前形式，考虑社会经济的发展效益和协调发展，通过捕捞方式的改变，达到最终目的。1、假设水域中两种鱼类的数据如下：x(0)=1500，y(0)=800，r1=20，r2=30，N1=2500，N2=1000，q1=0.03，q2=0.015。捕鱼的数据如下：k1=0.0001，k2=0.0003。用以上数据编程计算，画出两种鱼类的数量随时间变化的图形和两种鱼类数量的相轨线，用图形分析两种鱼类的数量随时间的变化趋势。2、假设鱼类的数据和问题1相同，捕鱼用如下的数据：k1=0.0002，k2=0.0006。编程计算，并画出两种鱼类的数量随时间变化的图形和两种鱼类数量的相轨线，用图形分析两种鱼类数量随时间的变化趋势。3、试分析如何改变捕鱼方式，既能提高经济效益，又使得人与自然能够和谐发展。二模型假设1、该水域为非开放式的，且不与其他水域发生关系（即鱼群之间不发生大规模迁徙），从而构成一个独立的生态部落；2、鱼群总量的增加虽然是离散的，但对大规模鱼群而言，我们可以认为鱼群的总量变化是随时间连续的；3、持续捕获使各年龄组的鱼群数量呈周期变化，周期为1年，可以只考虑鱼群数量在1年内的情况；4、假设同种鱼群个体之间是同质的，只考虑平均水平，不考虑个体差异；5、假设连续捕获是各年龄的鱼群数量呈周期性变化，可以只考虑鱼群数量在短期内的变化情况；6、捕捞也是一个连续的过程，不是在某一时刻突然发生。实用文档标准文案三符号说明()xt：甲种鱼类随时间变化的增长率；()yt：乙种鱼类随时间变化的增长率；11()()(1)xtrxtN：甲种鱼类满足Logistic模型的自然增长率；22()()(1)ytrxtN：乙种鱼类满足Logistic模型的自然增长率；1()()qxtyt：由题意可得的甲种鱼类的竞争减少的增长率；2()()qxtyt：由题意可得的乙种鱼类的竞争减少的增长率；1()kxt：由题意可得的甲种鱼类因人类捕捞减少的增长率；2()kxt：由题意可得的乙种鱼类因人类捕捞减少的增长率；:ixt在t时刻i龄鱼条数，1,2,3,4i；:k4龄鱼捕捞强度系数；:n每年产卵量；ia每年初i龄鱼的数量，1,2,3,4i；四模型建立与分析求解4.1问题一模型的建立与分析：1、问题分析：根据题意可知，在满足假设的条件下，影响两种鱼群数量的因素有同种鱼群的种类竞争（在内禀增长率之下的数量）、两种鱼群的种间竞争和人类的捕捞，综合以上三种因素，本文以微分方程的方法建立了类似Logistic模型的数学模型。鱼类的增长率由自然增长率、竞争减少的增长率、由于人类捕捞减少的增长率共同组成。实用文档标准文案2、建立模型：11112222()()()(1)()()()()()()(1)()()()xtxtrxtqxtytkxtNytytrxtqxtytkytN3、编程计算：Matlab编程计算可得x和y关于t的图像及x与y的相轨线如图所示：（其中比例系数120.0001,0.0003kk）4、结果分析：①由甲乙两类鱼的数量变化曲线可以看到，甲种鱼群数量的初始值为1500，乙种鱼群数量的初始值为800；②由于甲乙两鱼群的相互竞争以及人类的捕捞，甲鱼群的数量在持续减少，乙鱼群在减少约0.05个时间单位后增长；③在约0.18时间单位时，甲乙两鱼群的数量相等，达到平衡点。④平衡点后甲鱼群数量仍旧减少，直至约0.7个时间单位后减少为零；乙鱼群在约0.05个时间单位后增长，经过平衡点后，仍旧增长，直至增长到最大容纳数量1000后稳定；⑤生态学中有一个竞争排斥原理：若两个种群的单个成员消耗的资源差不实用文档标准文案多相同，而环境能承受的种群甲最大容量比种群乙大，那么种群乙终将灭亡；⑥由相轨线可以看到，随着甲种鱼群数量的增加乙种鱼群的数量一直在减少，直至甲种鱼群的数量从0增加到800时，乙种鱼群的数量由1000减少至约550，而后随着甲种鱼群数量继续增加，乙种鱼群的数量开始增加，当甲种鱼群数量为1500时，乙种鱼群数量达到800。从而可以得到，当甲种鱼群数量在约[0,550]，甲种鱼群数量的增长会抑制乙种鱼群数量，当甲种鱼群数量大于550，甲乙两鱼群会同时增加。4.2问题二模型的建立与分析：1、问题分析：此问题的模型与问题一的模型相同，只是改变了人类捕捞系数的大小。重复问题一的解题步骤，做出新的图象。2、编程计算：Matlab编程计算可得x和y关于t的图像及x与y的相轨线如图所示：（其中比例系数120.0002,0.0006kk）3、结果分析：改变人类捕捞系数后，与问题一相比，甲乙两鱼群的数量变化以及相轨线的变化均不明显。可以得出在此种模型下，由于人类捕捞减少的增长率对鱼群的增实用文档标准文案长率的影响非常小。但也有可能是因为系数过小，并且改变的范围过小，仍旧会没有什么明显的变化。4.3问题三模型的建立与分析：1、问题分析：为了保护人类赖以生存的自然环境，可再生资源（如渔业、林业资源）的开发必须适度。一种合理、简化的策略是，在实现可持续收获的前提下，追求最大产量或最佳效益。我们假设把两种鱼看做一个整体，鱼群之间的竞争为内部竞争，因竞争造成的死亡看作是自然死亡，捕鱼的经济效益与所捕鱼的总产量相关，与鱼的种类无关。考虑鱼群的最优捕捞方式：我们把这两种鱼分4个年龄组，称为1龄鱼，2龄鱼，3龄鱼，4龄鱼。各年龄组每条鱼的平均质量分别为20.50、130.55、517.86、1022.99（g），各年龄组的自然死亡率为0.7（1/年），这种鱼为季节性集中产卵繁殖，平均每条4龄鱼的产卵量为1.309×105（个），3龄鱼的产卵量为这个数的一半，2龄鱼和1龄鱼不产卵，产卵和孵化期为每年的最后4个月，卵孵化并成活为1龄鱼，成活率（1龄鱼条数与产卵量n之比）为1.22×105/（1.22×105+n）。渔业管理部门规定，每年只允许在产卵孵化期前的8个月进行捕捞作业。如果每年投入的捕捞能力（如渔船数、下网次数等）固定不变，这时单位时间捕捞量将与各年龄组鱼群条数成正比，比例系数不妨设为捕捞强度系数，只捕捞3龄鱼和4龄鱼，其两个捕捞系数之比为0.42：1。2、建立模型：①死亡率我们给出鱼的自然死亡率为0.7（1/年），理解为平均死亡率，这是单位时间鱼群死亡数量与现有鱼群数量的比例系数，由假设可知，它是一个与环境等其他因素无关的常数。鱼群的数量是连续变化的，且1龄鱼、2龄鱼在全年及3龄鱼、4龄鱼在后4个月的数量只与死亡率有关，各龄鱼的数量满足()0.7()iidxtxtdt，t∈[0，1]（i=1，2）()0.7()iidxtxtdt，t∈[2/3，1]（i=3，4）②捕捞强度系数单位时间4龄鱼捕捞量与4龄鱼群总数成正比，比例系数即为捕捞强度k，它是一定的，且只在捕捞期内（即每年的前8个月）捕捞3龄鱼，4龄鱼。所以，一方面捕捞强度系数决定了3龄鱼、4龄鱼在捕捞期内的数量，其变实用文档标准文案化规律为33()(0.70.42)()dxtkxtdt，t∈[0，2/3]44()(0.7)()dxtkxtdt，t∈[0，2/3]另一方面也决定了t时刻捕捞3龄鱼、4龄鱼，其数量分别为30.42()kxt和4()kxt。③成活率由于只有3龄鱼、4龄鱼在每年8月底一次产卵，因此可将每年的产卵量n表示为534221.30910[0.5()()]33nxx问题分析中已经说明成活率为551.22101.2210n，所以每年年初的1龄鱼的数量为5151.2210(0)1.2210xnn模型如下：s.t，t∈[0，1]，t∈[0，1]，t∈[0，2/3]，t∈[2/3,1]，t∈[0，2/3]，t∈[2/3,1]22333400max(())517.860.42()1022.99()totalkkxtdtkxtdt11()0.7()dxtxtdt22()0.7()dxtxtdt33()(0.70.42)()dxtkxtdt33()0.7()dxtxtdt44()(0.7)()dxtkxtdt44()0.7()dxtxtdt实用文档标准文案3、编程计算：用Matlab编程可得如下数据和图像：看到k=17.99995188885500,total=1.398478407687244e+007；4、结果分析：由数据可以得到：当4龄鱼捕捞强度系数为k=17.99995188885500时，3龄鱼的捕捞强度系数为0.42k时，最高年收获量为totalT=1.398478407687244