基于Volterra理论的黄鼠狼-老鼠捕食模型

  • 摘要
    • 问题重述
    • 模型假设
    • 符号说明
    • 问题分析
    • 模型建立求解
    • 模型的评价
    • 参考文献

摘要

自然界中存在着一种有趣的基本捕食系统,两个种群之间以既相互依存又相互制约的方式生活,在生态学中我们称其为食饵——捕食者系统。通过研究种群数量的变化对于合理利用资源及保护生态系统有重要意义,因此本文分两步对Volterra食饵——捕食者模型所描述的现象进行分析。首先,利用MATLAB求微分方程的数值解,通过对数值结果和图形的观察,猜测它的解析解的构造;然后,从理论上研究其平衡点及相轨线的形状,验证前面的猜测。

问题重述

假设某区域有黄鼠狼和老鼠两个物种,且老鼠有足够的食物,黄鼠狼捕食老鼠,试建立黄鼠狼和老鼠数量之间关系的微分方程。
注:建议考虑两个物种的出生和死亡率、黄鼠狼对老鼠的依赖关系,以及同一个物种之间的竞争关系等。
这是关于黄鼠狼和老鼠数量的两个因变量与时间t的函数关系构成的两个一阶微分方程的方程组。

模型假设

1.假设是一个包括黄鼠狼和老鼠两个种群生态系统,黄鼠狼靠老鼠而生存,系统与外界没有种群交换关系。
2.假设资源丰富,环境适宜,老鼠独立生存时数量成指数增长。
3.假设黄鼠狼只有老鼠这一种食物来源。

符号说明

符号说明
r 老鼠独立生存时增长率
x(t) 老鼠数量
y(t) 黄鼠狼数量
a 反应黄鼠狼掠取老鼠能力的比例系数
d 黄鼠狼独立生存的死亡率
b 老鼠对黄鼠狼的供养能力

问题分析

老鼠和黄鼠狼在时刻t的数量分别记作x(t),y(t),假设当老鼠独立生存时以指数规律增长,(相对)增长率为r,即而黄鼠狼的存在使老鼠的增长率减少,设减少率与黄鼠狼数量成正比,于是x(t)满足方程

比例系数a反映黄鼠狼掠取老鼠的能力。
假设黄鼠狼只有老鼠这一种食物来源,黄鼠狼离开老鼠将无法生存,设它独自存在时死亡率为d,即,而老鼠的存在为黄鼠狼提供了食物,相当于使黄鼠狼的死亡率降低,且促使其增长,设黄鼠狼增长率与食饵数量成正比,于是y(t)满足

比例系数b反映食饵对捕食者的供养能力
方程(1),(2)是在自然环境中食饵和捕食者之间依存和制约的关系,这里没有考
虑种群自身的阻滞增长作用,是Volterra提出的最简单的模型

模型建立求解

记食饵和捕食者的初始数量分别为:
x(0)=x0,y(0)=y0 (3)
为求微分方程(1),(2)满足初始条件(3)的数值解x(t),y(t)(并作图)及相轨线y(x),设r=1,d=0.5,a=0.1,b=0.02,x0=25,y0=2,用MATLAB软件编程计算,可得x(t) y(t)及相轨线y(x)如图1、图2(数值结果从略)。可以猜测,x(t),y(t)是周期函数,与此相应地,相轨线y(x)是封闭曲线。从数值解近似地定出周期为10.7,x的最大、最小值分别为99.3和2.0,y的最大、最小值分别为28.4和2.0,并且用数值积分容易算出x(t),y(t)在一个周期的平均值为=23, 。


模型的评价

尽管Volterra模型理论上可以解决一些捕食者与被捕食者数量关系的问题,但他作为近似反应现实对象的数学模型仍然存在局限性。
许多生态学家都指出,多数捕食者和被捕食者都观察不vloterra模型显示的周期性震荡,需要在vloterra模型中考虑自身阻滞作用的logistic项目。
另外,一些生态学家认为,自然界里长期存在的呈周期变化的生态系应该是结构稳定的,生态系统受到不可避免的干扰和后,其内部制约作用会使系统自动恢复原来的状态,而volterra模型描述的周期变化状态却不是结构稳定的,而为了得到反映周期变化的结构稳定的模型,要用到极限环的概念

参考文献

[1] 姜启源,谢金星.实用数学数学建模(基础篇).北京:高等教育出版社,2014.
[2] 唐静波,胡智渊,张彦琼,任力锋.确定Lotka-Volterra生态系统模型高精度参数的研究[J].数学的实践与认识,2007(14):77-82.

【数学模型】基于Volterra理论的捕食模型相关推荐

  1. 基于排队理论的云计算模型的MATLAB仿真

    1.简介与仿真结论 用matlab模拟,基本上可以说是一个排队理论.输入一堆包,经过buffer(一个或者几个都行)传给server,这些包会在buffer里面进行排队,采取FIFO的排序方法. 部分 ...

  2. matlab 摄动波浪理论,基于MATLAB的三维海浪模型数值仿真_齐宁.pdf

    ISSN1009-3044 E-mail:eduf@ 第9卷第25期 (2013年09月) ComputerKnowledgeandTechnology电脑知识与技术 ComputerKnowledg ...

  3. 基于博弈理论的入侵检测与响应模型综述

    基于博弈理论的入侵检测与响应模型综述 作者:中国保密协会科学技术分会 2020-05-07 20:52:01 安全 应用安全 当前网络规模急剧增加,各类入侵过程也逐渐向复杂化,多样化和分布式的趋势发展 ...

  4. 《模式识别原理及工程应用》——第2章 基于贝叶斯决策理论的分类器 2.1 分类器的描述方法...

    第2章 基于贝叶斯决策理论的分类器 2.1 分类器的描述方法 分类器是一种计算机程序,它的设计目标是在通过学习后,可自动将数据分到已知类别.分类器常应用在搜索引擎以及各种检索程序中,同时也大量应用于数 ...

  5. 基于matlab的光伏电池通用数学模型,基于MATLAB的光伏电池通用数学模型.doc

    本文由qpadm贡献 pdf文档可能在WAP端浏览体验不佳.建议您优先选择TXT,或下载源文件到本机查看. 第 25 卷第 4 期 2009 年 4 月 电 力 科 学 与 工 程 Vol.25, N ...

  6. 最优传输理论和生成模型的几何观点

    ©PaperWeekly 原创 · 作者|尹娟 学校|北京理工大学博士生 研究方向|随机过程.复杂网络 论文标题:A Geometric View of Optimal Transportation ...

  7. NeurIPS 2021 | CyGen:基于概率论理论的生成式建模新模式!

    关注公众号,发现CV技术之美 本文转自微软研究院AI头条. 编者按:在概率论中,两随机变量的一个联合分布可由一个变量的边缘分布和对应条件分布确定,也可对称地由另一变量的边缘分布和另一方向的条件分布确定 ...

  8. 基于角色的权限控制模型RBAC

    本文来说下基于角色的权限控制模型RBAC 文章目录 概述 RBAC权限模型简介 RBAC的演化进程 用户与权限直接关联 一个用户拥有一个角色 一个用户一个或多个角色 页面访问权限与操作权限 数据权限 ...

  9. 基于MatlabSimulin的微电网模型及光伏电池建模仿真分析

    目录 第一章 微电网 1 1.1概念 1 1.2技术应用 1 1.3 研究意义 1 第二章 微电网组成元件 2 2.1 开关器件 2 2.1.1 断路器 2 2.1.2 静态开关 2 2.2 分布式电 ...

  10. 055B ENMTools教程-基于R语言对MaxEnt模型优化-MaxEnt调参教程--更新日期2021-9

    055B-1 视频附带资料下载和密码:软件-数据-文献下载-持续更新 055B-2 ENMTools软件下载安装 055B-3 R软件和工具包安装 055B-4 生物气候因子的精度说明与选择方法(理论 ...

最新文章

  1. c语言判断数组里的数据先递增后递减,查找两个有序序数组(一个递增、一个递减)中第K大的数...
  2. No module named ‘tensorflow_hub‘
  3. [LeetCode 001] Two Sum
  4. struts导入Excel进行解析
  5. C# winform bin文件夹、obj文件夹、Properties文件夹下分别存放什么文件?
  6. IEDA中彻底删除项目
  7. C#(WinForm)上传图片保存到数据库和从数据库读取图片显示到窗体
  8. Azkaban-solo模式-安装
  9. python 读shell
  10. 装饰器python详解_python装饰器详解
  11. centos5编译内核
  12. Dynamics 365 Online-Unified User Interface
  13. Kubernetes总结一:初识Kubernetes
  14. vue3如何去掉控制台的warn信息
  15. java spring mvc json转对象,SpringMVC中使用@RequestBody,@ResponseBody注解实现Java对象和XML/JSON数据自动转换(上)......
  16. 网桥 网卡网桥有什么区别
  17. 我的梦想是成为一名计算机程序员英语怎么说,英语作文。我的梦想、我的梦想是成为电脑程序员。。80词...
  18. My python learning
  19. Your hostname, xxx resolves to a loopback address: 127.0.1.1; using x.x.x.x instead(on interface xx)
  20. 台式计算机32位和64位的区别,电脑操作系统中32位和64位到底有哪些区别?

热门文章

  1. echarts 地图自定义图标_echarts自定义图标的点击事件怎么添加
  2. coffeescript html5,深入浅出CoffeeScript
  3. 乔纳森-弗莱彻:被遗忘的搜索引擎之父
  4. ANU COMP1100 Lab1简介
  5. 每年10万被动收入目标分解 | 进击
  6. 【死磕 Spring】----- IOC 之深入理解 Spring IoC
  7. -1 转换成二进制数是多少,负数左移右移。
  8. ESXI(VSphere)主机补丁离线升级(多图)
  9. 计算机视觉可解释性——卷积神经网络中间层的可视化
  10. LTP(Linux Test Project)学习(五)——LTP代码学习