麦蚜生态系统燕尾突变的模型划分.docxVIP

麦蚜生态系统燕尾突变的模型划分 小麦是世界上最重要的粮食。为了确保小麦的稳定产量，防治病虫害已成为小麦产量提高的重要措施之一。麦蚜作为麦区的主要害虫, 常见种类有麦二叉蚜、麦长管蚜和禾谷缢管蚜, 其防治问题引起了一些学者的关注。许多学者通过建立麦蚜生态系统模型, 解决了生态领域的诸多问题。Matis等用尖点突变模型探讨了蚜虫的内禀增长率和气候生态因子的关系;赵慧燕等研究了农田生态系统中的麦蚜种群动态模型, 通过实验数据建立了折叠突变模型, 给出了蚜虫防治的关键期及防治强度;赵慧燕等又建立了尖角突变模型, 确立突变区域, 并给出了病虫害发生的趋势及防治对策;魏雪莲等建立燕尾突变模型, 通过对分歧点集的分析, 说明了害虫种群数量发生突变的条件和机理, 为实际情况提供理论依据;Piyaratne等利用实际数据, 并通过计算机编程实现了对燕尾突变模型的模拟和分析。因为种群数量的变化往往不受单一因素的影响, 所以用多个控制变量来研究种群数量更全面, 然而由于缺乏燕尾突变模型及更复杂初等突变模型的研究依据, 因此突变模型在生态学中的研究者相对较少, 研究比较困难。本文根据燕尾突变模型分歧点集的截线图形来分析小麦蚜虫种群动态的变化, 蚜虫种群动态的变化取决于当控制点穿过分歧点集不同区域时, 它们所在区域平衡点的性质是否发生变化。通过不用势函数而直接研究平衡曲线方程的办法, 得到平衡点个数的变化情况。燕尾突变模型的平衡曲面是一元四次方程, 本文主要研究可否用一元四次方程根的判别式来划分燕尾突变模型的突变区域。 1 麦自身分歧点集的建立及其实现 考虑受天敌和气象环境因素及植被状况影响的害虫种群动态模型为 式中:N为害虫密度;r为害虫种群的内禀增长率;t为时间;e为气象环境对害虫种群密度的影响因子;K为环境容纳量;P为天敌种群密度;k为天敌的捕食率;Nm为常数, 是捕食行为可实现时害虫的最小量;d为常数, 是天敌的饱和捕食量;r为模型参数, 其值与具体的害虫种群有关;P, e, K为控制变量;N-Nm表示天敌的实际捕食量。 当时, 种群达到平衡状态, 模型 (1) 的平衡曲面为 根据式 (3) 得模型 (1) 的奇点集为 将式 (3) 和式 (4) 联立消去x, 即可得到模型 (1) 的分歧点集方程, 但是这个过程较繁琐。文献通过固定参数u, 找到v-w截面曲线, 若对所有u值能画出该曲线, 用“平衡截面法”就可得到分歧点集的图形, 文献给出了通过联立消去x得到的分歧点集方程。本文利用文献中一元四次方程实根的判别依据, 对模型 (3) , 通过根的变化情况推出v-w面的截线方程及分歧点集的表达式, 并做出分歧点集的三维几何图形来确定麦蚜爆发区域, 为实际应用提供理论依据。 定理1对四次函数 注:这个定理是不完全的判据, 因为Δ0包括了两种不同的情况, 故定理2给出了四次函数实根的完全判据。 定理2对四次函数 有下面结论成立: Ⅰ.当a0a2a12时: (1) 若Δ0, 或 (q, p) = (1/4, 0) 时, f (x) 无实数根; (2) 若Δ=0, 且 (q, p) ≠ (1/4, 0) 时, f (x) 有一个二重实根; (3) 若Δ0, f (x) 有两个不相等的实根。 Ⅱ.当a0a2=a12时: (1) 若Δ0, f (x) 无实数根; (2) 若Δ=0, 分2种情况: (1) 当c2≠0时, f (x) 有一个二重实根, (2) 当c2=0时, f (x) 有一个四重实根; (3) 若Δ0, f (x) 有两个不相等的实根。 Ⅲ.当a0a2a12时: (1) 若Δ0, 分2种情况: (1) q1/4时, f (x) 无实数根, (2) q1/4时, f (x) 有4个不相等的实根; (2) 若Δ=0, 分4种情况: (1) q1/4, f (x) 有一个二重实根, (2) q=1/4, 即 (q, p) = (1/4, 0) , f (x) 有两个二重实根, (3) q1/4, f (x) 有一个二重实根和两个不同实根, (4) q=-1/12, f (x) 有一个三重实根和一个单实根; (3) 若Δ0, f (x) 有两个不相等的实根。 2 =0的模型结构 利用上述定理划分麦蚜生态系统的突变区域, 设 对于g (x) , 其根有下面三种情况。 情况1当u0时, 根据定理2中的结论Ⅰ有 (3) 若Δ0, 方程 (3) 有两个不相等的实根, 其中正值为稳定平衡点, 负值为不稳定平衡点。 由此可知, Δ从负到零再到正这个变化过程, 模型 (1) 的平衡点个数从两个变到一个再到零, 固定u值, Δ=0对应模型分歧点集的截线 (见图1a) , 此时麦蚜生态系统可能发生蚜虫爆发现象。 情况2当u=0时, 根据定理2中的结论Ⅱ有 相应判别式