基于突变现象的麦蚜种群广义系统模型构建与解析：理论、实践与展望

基于突变现象的麦蚜种群广义系统模型构建与解析：理论、实践与展望一、引言1.1研究背景与意义小麦作为全球重要的粮食作物之一，在保障粮食安全方面发挥着不可替代的关键作用。据联合国粮食及农业组织（FAO）的统计数据显示，全球小麦的种植面积广泛，年产量巨大，是人类主要的食物来源之一，为数十亿人口提供了基本的能量和营养支持。在中国，小麦同样占据着举足轻重的地位，是北方地区的主要粮食作物，种植历史悠久，种植区域覆盖了从东北平原到华北平原，再到西北干旱地区等广袤的土地。其产量和质量直接关系到国家的粮食供应稳定以及人民的生活福祉，对于维护社会稳定和促进经济发展意义深远。然而，在小麦的生长过程中，麦蚜的危害成为了影响小麦产量和质量的重要因素。麦蚜，作为小麦生产上的重大害虫，种类繁多，常见的有麦长管蚜、麦二叉蚜、禾缢管蚜等。这些麦蚜以其刺吸式口器刺入小麦植株，贪婪地吸食汁液，给小麦的生长发育带来了严重的负面影响。在叶片上，麦蚜的取食会导致叶片出现黄斑，随着危害的加重，黄斑逐渐扩大并连接成片，最终使叶片全部枯黄，严重影响了叶片的光合作用，导致小麦无法正常制造和积累养分。在茎秆部位，麦蚜吸食汁液会阻碍水分和养分的运输，使小麦生长停滞，分蘖减少，植株变得矮小瘦弱，无法形成健壮的个体，从而影响后续的抽穗和结实。当麦蚜大量聚集在麦穗上时，会造成籽粒饥瘦或不能结实，严重降低了小麦的产量和品质。更为严重的是，麦蚜还是多种病毒病的传播媒介，如小麦黄矮病等。当麦蚜在感染病毒的植株上取食后，再迁移到健康植株上时，就会将病毒传播给健康植株，引发病毒病的大面积流行，进一步加剧了小麦的减产损失。据相关研究表明，受麦蚜危害的小麦田块，一般减产幅度在10%-20%之间；在麦蚜爆发严重的年份，减产幅度可达50%以上，甚至导致绝收，给农业生产带来了巨大的经济损失。为了有效地控制麦蚜的危害，保障小麦的安全生产，众多学者从不同角度开展了广泛而深入的研究。在化学防治方面，开发了一系列高效、低毒的化学农药，如吡虫啉、啶虫脒、吡蚜酮等，这些农药在一定程度上能够快速有效地降低麦蚜的种群数量。然而，长期大量使用化学农药也带来了一系列严重的问题。一方面，化学农药的残留会对土壤、水源和空气等环境要素造成污染，破坏生态平衡，影响生物多样性；另一方面，麦蚜对化学农药的抗药性逐渐增强，使得农药的防治效果逐年下降，需要不断加大用药量和用药频率，进一步加剧了环境污染和农业生产成本的增加。在生物防治领域，研究人员致力于挖掘和利用麦蚜的天敌资源，如瓢虫、草蛉、食蚜蝇、蚜茧蜂等，通过释放这些天敌昆虫来控制麦蚜的种群增长。同时，微生物农药如白僵菌、绿僵菌等也被应用于麦蚜的防治，这些微生物能够感染麦蚜并导致其死亡，具有环境友好、不易产生抗药性等优点。此外，农业防治措施如合理密植、科学施肥、及时清除田间杂草等，通过改善麦田的生态环境，减少麦蚜的滋生和繁殖场所，也在一定程度上起到了控制麦蚜危害的作用。然而，单一的防治措施往往难以达到理想的防治效果，综合运用多种防治手段成为了当前麦蚜防治的发展趋势。在这样的背景下，建立麦蚜种群广义系统模型具有重要的理论和现实意义。从理论层面来看，麦蚜种群的动态变化受到多种复杂因素的相互作用，包括生物因素（如作物、天敌等）、环境因素（如温度、湿度、光照等）以及人为因素（如农业生产活动、农药使用等）。传统的研究方法往往难以全面、深入地揭示这些因素之间的内在联系和作用机制。而广义系统理论作为一种强大的分析工具，能够将麦蚜种群与其他相关因素视为一个有机的整体，从系统的角度出发，综合考虑各种因素的影响，为研究麦蚜种群的动态变化提供了一个全新的视角。通过建立麦蚜种群广义系统模型，可以更加准确地描述麦蚜种群的增长、波动和突变等现象，深入剖析麦蚜种群爆发的生态学机制，丰富和完善害虫种群生态学的理论体系，为进一步研究其他害虫种群的动态变化提供有益的借鉴和参考。从现实应用角度而言，准确预测麦蚜种群的动态变化对于制定科学合理的防治策略至关重要。通过建立的广义系统模型，可以根据不同的环境条件和防治措施，对麦蚜种群的未来发展趋势进行模拟和预测，为农业生产者提供及时、准确的预警信息。例如，模型可以预测在不同的气候条件下，麦蚜种群的爆发时间、规模和危害程度，从而指导生产者提前做好防治准备，合理安排防治时间和措施。同时，模型还可以评估不同防治策略的效果，如不同农药的使用剂量和时机、天敌昆虫的释放数量和时间等，帮助生产者选择最优的防治方案，实现精准防治。这样不仅可以提高防治效果，有效降低麦蚜的危害，减少小麦的减产损失，还可以避免盲目使用农药和过度防治带来的环境污染和资源浪费问题，促进农业的可持续发展。此外，麦蚜种群广义系统模型的建立和应用，还可以为政府部门制定农业政策和规划提供科学依据，有助于加强对农业害虫的宏观管理和控制，保障国家的粮食安全和生态安全。1.2国内外研究现状在麦蚜种群动态研究领域，国内外学者开展了大量的工作。张德昌、张广林等对麦蚜种群动态进行分析，揭示出其种群动态存在6个状态集，明确了防治重点应放在初盛期，为制定防治策略提供了一定依据。通过对麦长管蚜和禾缢管蚜的研究发现，二者在发生时间和数量变化上存在差异，麦长管蚜发生相对较早，数量占主导，禾缢管蚜发生较晚且数量较少。同时，利用逻辑斯蒂模型对麦蚜混合种群动态增长进行拟合，得出环境负荷量和种群增长拐点，为预测麦蚜种群数量变化提供了方法。在突变现象研究方面，赵惠燕、汪世泽等运用突变论对麦蚜种群数量变动中的突变机理进行证实，指出防治害虫强度的不当增加可能导致害虫密度加大，只有使防治强度达到一定水平（如74.38％）才能控制害虫在经济阈值以下。韩立红从广义系统理论角度出发，在麦蚜种群折叠突变模型基础上，考虑作物种群对麦蚜种群的消长作用，建立两种群广义麦蚜生态系统模型，利用非线性广义系统稳定性理论和奇异值理论，分析平衡点稳定性和麦蚜种群突变机理。还进一步考虑作物种群、天敌种群对麦蚜种群的影响，建立三种群广义麦蚜生态系统模型，给出麦蚜种群发生突变的充分条件，为研究麦蚜种群突变现象提供了新的视角和方法。在广义系统模型应用于生物种群研究领域，虽然取得了一些进展，但仍存在不足。在模型构建方面，目前的研究大多侧重于考虑单一或少数几种因素对麦蚜种群的影响，而实际生态系统中，麦蚜种群受到生物因素（如多种天敌、不同作物品种等）、环境因素（如温度、湿度、光照等）以及人为因素（如农业生产活动、农药使用等）的综合作用，如何全面、准确地将这些复杂因素纳入广义系统模型是亟待解决的问题。在模型参数确定上，许多参数的获取依赖于实验数据和经验估计，存在一定的不确定性，这可能影响模型的准确性和可靠性。在模型验证和应用方面，目前的研究在实际田间验证和大规模推广应用方面还存在欠缺，如何将模型与实际农业生产相结合，为麦蚜防治提供切实可行的指导，还需要进一步的研究和实践。1.3研究目标与内容本研究旨在基于突变现象，构建麦蚜种群广义系统模型，深入分析麦蚜种群突变机理，为麦蚜的科学防治提供理论依据和技术支持。具体研究内容如下：麦蚜种群广义系统模型构建：在现有研究基础上，全面考虑生物因素（如作物种群、天敌种群等）、环境因素（如温度、湿度、光照等）以及人为因素（如农药使用、农事操作等）对麦蚜种群的综合影响，运用广义系统理论，构建更加完善的麦蚜种群广义系统模型。通过对模型的建立和参数设定，准确描述麦蚜种群在不同条件下的动态变化过程，为后续的分析和预测提供基础。麦蚜种群突变现象分析：结合突变理论和广义系统理论，对构建的麦蚜种群广义系统模型进行深入分析。研究麦蚜种群在各种因素作用下发生突变的条件和机制，确定影响麦蚜种群突变的关键因素。通过数值模拟和实例分析，验证模型的准确性和可靠性，揭示麦蚜种群突变的规律和特点，为理解麦蚜种群的爆发机制提供理论依据。模型验证与参数优化：收集实际田间数据，对构建的麦蚜种群广义系统模型进行验证和校准。通过对比模型预测结果与实际观测数据，评估模型的准确性和适用性。针对模型存在的偏差和不足，运用参数估计和优化算法，对模型参数进行调整和优化，提高模型的预测精度和可靠性。同时，考虑不同地区、不同年份的环境差异和数据特点，对模型进行适应性调整，使其能够更好地应用于实际生产中。防治策略制定与效果评估：根据麦蚜种群广义系统模型的分析结果，结合实际农业生产情况，制定科学合理的麦蚜防治策略。包括生物防治措施（如释放天敌昆虫、利用微生物农药等）、化学防治措施（如合理选择农药种类和使用剂量、优化施药时间和方法等）以及农业防治措施（如合理密植、科学施肥、及时清除田间杂草等）的综合应用。利用模型对不同防治策略的效果进行模拟和评估，比较不同策略下麦蚜种群的动态变化和防治成本，筛选出最优的防治方案，为农业生产者提供决策支持。不确定性分析与风险评估：考虑到生态系统的复杂性和不确定性，对麦蚜种群广义系统模型中的参数和输入变量进行不确定性分析。评估参数不确定性对模型预测结果的影响程度，确定模型的不确定性范围。在此基础上，开展麦蚜种群爆发的风险评估，分析不同环境条件和防治措施下麦蚜种群爆发的可能性和风险程度，为制定风险管理策略提供依据。通过不确定性分析和风险评估，提高麦蚜防治决策的科学性和稳健性，降低农业生产风险。1.4研究方法与技术路线本研究综合运用多种研究方法，从理论分析、模型构建、数值模拟到实际验证，全面深入地开展对麦蚜种群广义系统模型的研究。生态学原理：深入研究麦蚜与周围生物和环境的相互关系，包括麦蚜与小麦作物之间的取食与防御关系，麦蚜与天敌之间的捕食与寄生关系，以及温度、湿度、光照等环境因素对麦蚜生长发育、繁殖和种群动态的影响。通过实地观察、田间试验和文献调研，获取麦蚜生态系统中各生物和环境因子的详细信息，为后续的模型构建提供坚实的生态学基础。例如，通过在不同生态区域设置长期监测样地，定期调查麦蚜及其天敌的种群数量、分布情况，以及小麦的生长状况和环境因子的变化，分析它们之间的相互作用规律。数学建模：运用广义系统理论，将麦蚜种群与作物种群、天敌种群以及环境因素等视为一个相互关联的广义系统。建立状态方程来描述系统中各变量（如麦蚜数量、作物生物量、天敌数量等）随时间的变化，考虑各因素之间的非线性相互作用，引入合适的参数来刻画这些关系。例如，对于麦蚜与作物的关系，通过建立麦蚜取食对作物生长影响的函数，以及作物对麦蚜繁殖和生存的反馈函数，来准确描述两者之间的动态关系；对于麦蚜与天敌的关系，考虑天敌的捕食效率、搜索半径等因素，建立相应的捕食模型。在建模过程中，充分参考前人的研究成果，并结合实际观测数据，对模型进行不断的完善和优化。稳定性分析：利用非线性广义系统稳定性理论，对构建的麦蚜种群广义系统模型进行稳定性分析。确定系统的平衡点，通过计算雅可比矩阵的特征值，判断平衡点的稳定性。分析不同参数变化对系统稳定性的影响，找出影响麦蚜种群动态稳定性的关键因素。例如，研究天敌数量的变化、农药使用强度的改变等因素对麦蚜种群平衡点稳定性的影响，确定在何种条件下麦蚜种群能够保持稳定，何种条件下会发生突变。通过稳定性分析，深入理解麦蚜种群动态变化的内在机制，为制定有效的防治策略提供理论依据。数值模拟：运用Matlab、Python等数学软件，对麦蚜种群广义系统模型进行数值模拟。设置不同的初始条件和参数值，模拟在不同环境条件和防治措施下麦蚜种群的动态变化过程。通过数值模拟，直观地展示麦蚜种群数量的增长、波动和突变等现象，预测麦蚜种群的未来发展趋势。例如，模拟在不同温度、湿度条件下麦蚜种群的增长情况，以及在不同农药使用剂量和释放天敌数量的情况下麦蚜种群的控制效果。将数值模拟结果与实际观测数据进行对比，验证模型的准确性和可靠性，对模型进行进一步的优化和改进。实验验证：在实际麦田中开展实验，收集麦蚜种群数量、作物生长状况、天敌数量以及环境因素等数据。将实验数据与模型预测结果进行对比分析，评估模型的准确性和适用性。根据实验结果，对模型参数进行调整和优化，提高模型的精度和可靠性。例如，在不同的麦田区域设置对照实验，分别采用不同的防治措施，定期监测麦蚜种群数量的变化，并将这些数据用于模型的验证和参数优化。通过实验验证，确保模型能够真实反映实际麦蚜生态系统的动态变化，为实际应用提供有力支持。本研究的技术路线如图1-1所示：首先，通过广泛的文献调研和实地观察，收集麦蚜生态系统的相关资料，包括麦蚜的生物学特性、生态习性，以及与作物、天敌和环境的相互关系等。在此基础上，运用生态学原理和广义系统理论，构建麦蚜种群广义系统模型，并对模型进行稳定性分析和参数估计。然后，利用数值模拟方法，对模型进行仿真研究，分析不同因素对麦蚜种群动态的影响。同时，开展实地实验，收集实验数据，对模型进行验证和优化。最后，根据模型分析和实验结果，制定科学合理的麦蚜防治策略，并对防治效果进行评估和反馈，为麦蚜的有效防治提供科学依据和技术支持。\begin{figure}[h]\centering\includegraphics[width=0.8\textwidth]{技术路线图.jpg}\caption{技术路线图}\end{figure}\begin{figure}[h]\centering\includegraphics[width=0.8\textwidth]{技术路线图.jpg}\caption{技术路线图}\end{figure}\centering\includegraphics[width=0.8\textwidth]{技术路线图.jpg}\caption{技术路线图}\end{figure}\includegraphics[width=0.8\textwidth]{技术路线图.jpg}\caption{技术路线图}\end{figure}\caption{技术路线图}\end{figure}\end{figure}二、麦蚜种群动态及突变现象分析2.1麦蚜种群概述麦蚜隶属同翅目蚜科，是一类对小麦生产造成严重危害的害虫，其种类繁多，在全球范围内已被记录的种类超过百种。在我国，对小麦危害较为严重且普遍发生的麦蚜种类主要有麦长管蚜（Macrosiphumavenae）、麦二叉蚜（Schizaphisgraminum）和禾缢管蚜（Rhopalosiphumpadi）。这三种麦蚜在形态特征、生活习性以及对小麦的危害方式上既有相似之处，又存在一定的差异。麦长管蚜无翅孤雌蚜体长卵形，颜色从草绿色至橙红色不等，头部略显灰色，腹侧具灰绿色斑；有翅孤雌蚜体呈圆形，绿色，触角黑色，腹管长圆筒形且为黑色，尾片长圆锥状。麦二叉蚜无翅孤雌蚜体卵圆形，淡绿色，背中线深绿色，腹管浅绿色，顶端黑色，中胸腹部具短柄，触角6节，尾片长圆锥形；有翅孤雌蚜体长卵形，体绿色，背中线深绿色，头、胸黑色，腹部色浅，触角黑色共6节，前翅中脉二叉状。禾缢管蚜无翅孤雌蚜体宽卵形，体表颜色从绿色至橙红色，常被薄粉，头部光滑，胸腹部背面有清楚网纹，腹管黑色，长圆角形，端部略凹缢，有瓦纹，触角6节；有翅孤雌蚜体长卵形，头胸黑色，腹部绿色至深绿色，触角黑色6节，短于体长。在生活习性方面，麦蚜的越冬虫态及场所因各地气候条件而异。在南方温暖地区，麦蚜无明显的越冬期，可全年进行繁殖和为害；而在北方麦区以及黄河流域麦区，麦蚜多以无翅胎生雌蚜的形态在麦株基部叶丛或土缝内越冬；在北部较为寒冷的麦区，麦蚜则多以卵的形式在麦苗枯叶上、杂草上、茬管中或土缝内越冬，且越往北，以卵越冬的比例越高。从发生时间来看，麦二叉蚜通常发生时间较早，在小麦苗期即可大量出现；麦长管蚜一般在小麦拔节后，随着气温的升高和小麦植株的生长，其种群数量逐渐增加并加重为害；禾缢管蚜的发生时间相对较晚，多在小麦生长后期，特别是灌浆期前后，其种群数量会有所上升。麦蚜的繁殖能力极强，在适宜的环境条件下，如温度、湿度、食物资源等满足其生长需求时，麦蚜能够以孤雌生殖的方式迅速繁殖后代，短时间内即可形成庞大的种群数量。据研究表明，麦蚜在适宜条件下，每头雌蚜每天可产若蚜数头，一生可产若蚜数十头甚至上百头，且若蚜发育至成蚜的时间较短，一般只需数天至十几天不等，这使得麦蚜种群能够在短时间内迅速增长并对小麦造成严重危害。麦蚜的分布范围极为广泛，几乎遍及世界各个产麦国家和地区。在我国，麦蚜的分布也十分普遍，从北方的东北平原、华北平原到南方的长江流域、华南地区，以及西部的干旱半干旱地区，只要有小麦种植的地方，就可能有麦蚜的存在。不同种类的麦蚜在分布上也存在一定的偏好性。麦长管蚜适应能力较强，在我国南北方麦区的密度均相对较大，但在偏北方地区，由于其对温度和湿度的适应性特点，发生更为严重；麦二叉蚜主要分布于长江以北各省，在比较少雨的西北冬春麦区，因其耐旱怕湿的特性，发生频率最高；禾缢管蚜在各地均有分布，但其对湿度和光照的特殊需求，使其在某些特定环境条件下更容易大量繁殖和为害。在麦田生态系统中，麦蚜占据着重要的生态位，是麦田生态系统中的关键生物因子之一。麦蚜以小麦植株为食，通过其刺吸式口器刺入小麦的叶片、茎秆、嫩穗等组织，吸食汁液，直接影响小麦的生长发育、光合作用和营养物质的合成与运输。在小麦苗期，麦蚜的为害会导致叶片枯黄，生长停滞，分蘖减少，使小麦植株无法形成健壮的个体，影响后续的生长和产量；在小麦生长后期，尤其是灌浆期，麦蚜集中在穗部为害，会造成麦粒不饱满，形成秕粒，严重降低小麦的千粒重和品质，导致减产甚至绝收。此外，麦蚜还会分泌蜜露，这些蜜露覆盖在小麦叶片和穗部表面，不仅影响小麦的光合作用和呼吸作用，还容易诱发煤污病等病害，进一步加重对小麦的危害。更为严重的是，麦蚜是多种病毒病的传播媒介，如小麦黄矮病等。当麦蚜在感染病毒的小麦植株上取食后，病毒会在其体内增殖，随后麦蚜再迁移到健康的小麦植株上取食时，就会将病毒传播给健康植株，引发病毒病的大面积流行，导致小麦减产甚至绝收。据统计，在麦蚜与病毒病共同为害的情况下，小麦的减产幅度可达30%-50%，严重威胁着小麦的安全生产和粮食安全。同时，麦蚜的种群动态变化也会对麦田生态系统中的其他生物产生影响。作为众多天敌昆虫的食物来源，麦蚜种群数量的波动会直接影响到瓢虫、草蛉、食蚜蝇、蚜茧蜂等天敌昆虫的种群数量和分布。当麦蚜种群数量增加时，天敌昆虫的食物资源丰富，其种群数量也会相应增加；反之，当麦蚜种群数量减少时，天敌昆虫的食物短缺，其种群数量也会随之下降。这种相互作用关系维持着麦田生态系统的相对平衡，对生态系统的稳定性和生物多样性具有重要意义。2.2麦蚜种群动态特征麦蚜种群数量在时间尺度上呈现出复杂的变化规律，受到多种因素的综合影响。从季节消长来看，在冬麦区，麦蚜的发生通常可分为几个明显的阶段。秋季冬麦出苗后，麦蚜从夏寄主上迁入麦田，此时种群数量开始逐渐增加，形成一个小高峰，但由于气温逐渐下降，麦蚜的繁殖速度受到一定限制，整体为害程度相对较轻。进入冬季后，麦蚜以不同的虫态在麦田中越冬，种群数量处于相对较低的水平。春季随着气温回升，小麦返青，当气温高于6℃时，麦蚜开始繁殖，种群数量缓慢上升。随着气温进一步升高，麦苗抽穗，麦蚜的繁殖速度明显加快，大量麦蚜迁移至穗部为害，虫口数量迅速上升。到了灌浆和乳熟期，麦蚜种群数量达到高峰，对小麦的产量和品质造成严重威胁。此后，随着气温继续升高，当气温高于22℃时，麦蚜会产生大量有翅蚜，迁飞到阴凉地带越夏，种群数量急剧下降。在春麦区，麦蚜的发生规律也类似，春小麦出苗后，麦蚜开始迁入并繁殖，随着小麦生长发育，在抽穗至灌浆期达到种群数量高峰，之后随着小麦成熟和环境变化，麦蚜种群数量减少。麦蚜种群数量的年际波动也较为明显。据相关研究表明，在不同年份，麦蚜的发生程度和种群数量存在显著差异。例如，通过对某地区连续多年的麦蚜种群数量监测发现，在某些年份，麦蚜的发生量较小，对小麦的危害较轻；而在另一些年份，麦蚜可能会大规模爆发，导致小麦严重减产。这种年际波动与多种因素密切相关，其中气候条件是一个重要的影响因素。气温、降水、光照等气候因子的年际变化会直接影响麦蚜的生长发育、繁殖和生存。在冬暖春早的年份，有利于麦蚜的越冬和早春繁殖，使得麦蚜种群数量迅速增加，容易引发麦蚜的猖獗发生；而在气温偏低、降水较多的年份，麦蚜的繁殖和活动可能会受到抑制，种群数量相对较低。农业生产活动的变化也会对麦蚜种群数量的年际波动产生影响。种植制度的调整、农药使用的频率和剂量、农事操作的时间和方式等，都会改变麦田的生态环境，进而影响麦蚜的种群动态。若某一年份大面积种植了对麦蚜具有抗性的小麦品种，或者加强了田间管理，及时清除杂草，合理施肥灌溉，都可能降低麦蚜的发生数量；相反，若不合理使用农药，导致麦蚜的天敌数量减少，或者种植制度不合理，为麦蚜提供了适宜的生存和繁殖条件，就可能使得麦蚜种群数量在该年份增加。在小麦的不同生长阶段，麦蚜种群数量也呈现出不同的变化特征。在小麦苗期，麦二叉蚜通常发生较早，其种群数量会在一定程度上逐渐增加，主要集中在叶片背面为害，吸食叶片汁液，导致叶片出现黄斑、枯黄等症状。随着小麦生长进入拔节期，麦长管蚜的种群数量开始逐渐上升，此时麦长管蚜和麦二叉蚜共同为害小麦，它们在叶片、茎秆等部位取食，影响小麦的正常生长发育，导致小麦分蘖减少，植株生长缓慢。当小麦进入抽穗期后，麦长管蚜的种群数量迅速增长，成为优势种群，大量聚集在穗部为害，严重影响小麦的灌浆和结实。禾缢管蚜在小麦生长后期，特别是灌浆期前后，种群数量会有所上升，与麦长管蚜一起对小麦穗部造成危害，导致麦粒不饱满，千粒重降低。2.3麦蚜种群突变现象表现麦蚜种群在特定条件下会呈现出明显的突变现象，其中种群数量的急剧增长是最为显著的表现形式之一。以2017年河南省部分麦区的麦蚜爆发为例，在当年4月中旬，随着气温的迅速回升，降水相对减少，气候条件变得极为适宜麦蚜的繁殖。据当地农业部门的调查数据显示，在4月15日时，麦蚜的平均百株蚜量仅为50头左右，处于相对较低的水平。然而，仅仅过了一周，到4月22日，平均百株蚜量就急剧攀升至200头以上，增长幅度超过了300%。此后，麦蚜种群数量继续快速增长，到5月上旬，部分麦田的百株蚜量甚至突破了1000头，达到了严重爆发的程度。这种短期内种群数量的爆发式增长，使得小麦植株在短时间内遭受了严重的侵害。大量麦蚜聚集在小麦叶片和穗部，疯狂吸食汁液，导致叶片迅速枯黄卷曲，光合作用无法正常进行，小麦的生长发育受到极大阻碍。麦穗也因麦蚜的为害而无法正常灌浆，籽粒干瘪，严重影响了小麦的产量和品质。据统计，受麦蚜严重危害的麦田，小麦减产幅度达到了20%-30%，给当地的小麦生产带来了巨大的经济损失。麦蚜种群在空间分布上也会出现突变现象。在麦蚜爆发初期，其通常在麦田中呈零星分散分布，密度相对较低。但随着种群数量的迅速增加，麦蚜会在短时间内迅速扩散，呈现出聚集分布的特征。在2020年山东省某地区的麦田中，在麦蚜发生初期，通过5点取样法调查发现，麦蚜主要集中在少数几株小麦上，且分布较为分散，平均每株小麦上的蚜量不足5头。然而，随着时间的推移，在短短一周内，麦蚜迅速向周围植株扩散，原本分散的分布格局被打破，大量麦蚜聚集在部分区域的小麦植株上，形成了明显的聚集斑块。这些聚集斑块内的麦蚜密度极高，每株小麦上的蚜量可达数十头甚至上百头，而周围区域的麦蚜密度相对较低。这种空间分布的突变现象，使得麦蚜对小麦的危害范围迅速扩大，增加了防治的难度。同时，聚集分布的麦蚜也更容易对小麦造成局部严重危害，导致小麦生长的不均匀性加剧，进一步影响了小麦的整体产量和质量。麦蚜的为害程度也会发生突变。在正常情况下，麦蚜对小麦的为害相对较轻，仅会导致小麦叶片出现少量黄斑，生长速度略有减缓。但在某些特殊条件下，如气候异常、天敌数量骤减等，麦蚜的为害程度会突然加重。在2019年河北省某麦区，前期麦蚜的为害程度较轻，小麦生长基本正常。然而，在5月中旬，由于连续的高温干旱天气，使得麦蚜的繁殖速度大幅加快，同时天敌昆虫的活动受到抑制，数量明显减少。在这种情况下，麦蚜的为害程度急剧加重，原本只是叶片出现少量黄斑的小麦，在短短几天内，叶片大面积枯黄，茎秆变得脆弱易折，麦穗发育不良，出现大量秕粒。原本预计的小麦产量因麦蚜为害程度的突变而大幅下降，严重影响了农民的经济收入。2.4影响麦蚜种群突变的因素麦蚜种群的突变受到多种因素的综合影响，这些因素相互作用，共同决定了麦蚜种群的动态变化。气候因素在麦蚜种群突变中起着关键作用，其中温度的影响尤为显著。麦蚜的生长发育和繁殖对温度具有较高的敏感性，不同的温度条件会直接影响麦蚜的生物学特性。麦长管蚜适宜在中温环境下生存，其最适宜的温度范围一般在15℃-20℃之间。在这个温度区间内，麦长管蚜的新陈代谢活动较为活跃，酶的活性较高，能够高效地进行物质和能量的转化，从而促进其生长发育和繁殖。据研究表明，在适宜温度条件下，麦长管蚜的若蚜发育历期较短，一般只需5-7天即可发育为成蚜，且成蚜的繁殖能力较强，平均每头成蚜每天可产若蚜3-5头。然而，当温度过高或过低时，都会对麦长管蚜的生存和繁殖产生不利影响。当温度超过25℃时，麦长管蚜的生长发育会受到抑制，其若蚜发育历期延长，繁殖能力下降，平均每头成蚜每天的产蚜量可降至1-2头。这是因为高温会导致麦长管蚜体内的生理生化过程紊乱，酶的活性降低，影响了其对营养物质的摄取和利用，进而影响了其生长和繁殖。麦二叉蚜则具有较强的耐高温能力，能够在相对较高的温度下生存和繁殖。其适宜的温度范围相对较宽，一般在10℃-30℃之间。在25℃-30℃的温度条件下，麦二叉蚜的繁殖速度较快，种群数量增长迅速。但当温度低于10℃时，麦二叉蚜的活动能力和繁殖能力会明显减弱。低温会使麦二叉蚜的新陈代谢减缓，酶的活性降低，导致其生长发育缓慢，繁殖周期延长。降水对麦蚜种群的影响也不容忽视。降水通过改变麦田的湿度环境，间接影响麦蚜的生存和繁殖。适度的降水能够维持麦田的适宜湿度，为麦蚜提供较为有利的生存条件。在湿度适宜的环境下，麦蚜的体表水分蒸发较慢，能够保持体内的水分平衡，有利于其正常的生理活动。同时，适宜的湿度也有利于小麦植株的生长，为麦蚜提供了充足的食物资源。然而，过多的降水会导致麦田湿度过高，容易引发各种病害，如白粉病、锈病等，这些病害不仅会影响小麦的生长发育，降低其抗蚜能力，还会直接对麦蚜造成危害。高湿度环境还会影响麦蚜的取食和繁殖行为，使其繁殖速度下降。相反，降水过少会导致麦田干旱，小麦生长受到抑制，麦蚜的食物质量和数量都会下降，从而不利于麦蚜的生存和繁殖。干旱会使小麦叶片失水卷曲，营养物质的合成和运输受阻，麦蚜吸食的汁液中营养成分减少，影响其生长和繁殖。生物因素同样对麦蚜种群突变有着重要影响。天敌昆虫作为麦蚜种群的重要调控因子，对麦蚜的数量起着关键的制约作用。常见的麦蚜天敌包括瓢虫、草蛉、食蚜蝇、蚜茧蜂等。这些天敌昆虫通过捕食或寄生麦蚜，有效地控制了麦蚜的种群数量。七星瓢虫是麦蚜的重要捕食性天敌之一，其成虫和幼虫都以麦蚜为食。据观察，一头七星瓢虫成虫每天可捕食麦蚜50-100头，幼虫每天可捕食麦蚜20-50头。七星瓢虫具有敏锐的嗅觉和视觉，能够迅速发现麦蚜的踪迹，并通过其强大的口器将麦蚜捕食。蚜茧蜂则是麦蚜的寄生性天敌，它将卵产在麦蚜体内，卵孵化后，幼虫在麦蚜体内生长发育，最终导致麦蚜死亡。研究表明，当蚜茧蜂的寄生率达到一定程度时，麦蚜的种群数量会显著下降。若蚜茧蜂的寄生率达到30%以上，麦蚜种群的增长速度会明显减缓，甚至出现负增长。如果麦田中这些天敌昆虫的数量充足，麦蚜种群的增长就会受到有效抑制，从而降低麦蚜种群突变的可能性。相反，当天敌昆虫数量减少时，麦蚜失去了有效的自然控制，种群数量可能会迅速增长，进而引发麦蚜种群的突变。小麦品种的抗性差异也会影响麦蚜种群的动态变化。不同小麦品种对麦蚜的抗性存在显著差异，这种抗性主要体现在形态结构和生理生化特性等方面。一些小麦品种具有特殊的形态结构，如叶片表面的蜡质层较厚、叶片茸毛较多等，这些结构特征能够阻碍麦蚜的取食和活动。蜡质层较厚的叶片能够增加麦蚜口器刺入的难度，使其取食效率降低；叶片茸毛较多则会干扰麦蚜的爬行和定位，使其难以找到合适的取食部位。部分小麦品种在生理生化方面具有抗性机制，如能够分泌一些对麦蚜有毒或抑制其生长发育的物质。这些物质可能会影响麦蚜的消化酶活性，使其无法正常消化和吸收食物中的营养成分，从而抑制麦蚜的生长和繁殖。高抗麦蚜的小麦品种能够有效地抵御麦蚜的侵害，减少麦蚜在其上的繁殖和生存，降低麦蚜种群突变的风险。而感蚜品种则容易受到麦蚜的攻击，为麦蚜的繁殖提供了适宜的环境，增加了麦蚜种群突变的可能性。人为因素在麦蚜种群突变中也扮演着重要角色。农药的不合理使用是导致麦蚜种群突变的一个重要人为因素。长期大量使用单一类型的农药，会使麦蚜逐渐产生抗药性。以吡虫啉为例，在一些地区连续多年使用吡虫啉防治麦蚜后，麦蚜对吡虫啉的抗性倍数不断增加。据监测，部分地区麦蚜对吡虫啉的抗性倍数已达到几十倍甚至上百倍。抗药性的产生使得农药的防治效果大幅下降，为了达到相同的防治效果，不得不增加农药的使用剂量和频率。然而，这种做法不仅会进一步污染环境，还会大量杀死麦蚜的天敌昆虫。当麦蚜的天敌数量减少后，麦蚜失去了自然控制，种群数量可能会迅速反弹，甚至出现爆发性增长，从而导致麦蚜种群的突变。种植制度的变化也会对麦蚜种群产生影响。随着农业生产的发展，种植制度不断调整，如小麦与其他作物的轮作、间作模式的改变等。不合理的种植制度可能会为麦蚜提供更多的生存和繁殖机会。在小麦与玉米轮作的农田中，如果玉米收获后未能及时清理田间杂草和残茬，这些杂草和残茬就可能成为麦蚜的越冬场所和食物来源，增加了麦蚜的越冬基数。当第二年小麦种植后，麦蚜就能够迅速繁殖并迁移到小麦上为害，导致麦蚜种群数量增加，增加了麦蚜种群突变的风险。而合理的种植制度，如采用小麦与豆类等作物间作，能够改善农田的生态环境，增加生物多样性，从而抑制麦蚜的生长和繁殖，降低麦蚜种群突变的可能性。豆类作物能够固定空气中的氮素，改善土壤肥力，促进小麦的生长，增强小麦的抗蚜能力。豆类作物还能够吸引一些有益昆虫，如寄生蜂等，这些昆虫能够捕食或寄生麦蚜，从而控制麦蚜的种群数量。三、麦蚜种群广义系统模型的建立3.1建模理论基础突变理论作为一门研究不连续现象的新兴数学分支，运用拓扑学、奇点理论和结构稳定性等数学工具，深入剖析自然界各种形态、结构以及社会经济活动中的非连续性突然变化现象。其核心在于用形象且精确的数学模型来阐释质量互变的过程，为理解复杂系统中的突变现象提供了独特的视角。突变理论起源于19世纪末，早期关于方程定性理论的研究孕育了其基本思想。庞加莱（Poincare）提出常微分方程解的三个要素，为突变理论的发展奠定了一定基础，尽管当时他的思想未被广泛接受。20世纪30年代，莫尔斯（HaroldMarstonMorse）提出临界理论，进一步推动了相关研究。1955年，惠特尼（Whitney）发表论文《曲面到平面的映射》，指出空间曲面到平面的投影仅可能存在折叠和尖点这两种奇异性，这便是突变理论中最低阶的两种突变形式，为突变理论的数学基础构建了关键的一环。1968年，雷内・托姆（ReneThom）在总结和继承前人研究成果的基础上，发表了论文《生物学中的拓朴模型》，并于1972年出版了《稳定性结构与形态形成学》专著，系统地阐述了突变理论，标志着这一理论的正式诞生。此后，经过众多学者的不断改进与完善，突变理论逐渐引起国际上各领域专家的广泛关注以及媒体报道的重视。突变理论具有7种基本突变模型，分别为折叠型、尖点型、燕尾型、蝴蝶型、双曲脐点型、椭圆脐点型和抛物脐点型。以折叠型突变模型为例，它描述了系统在一个控制参数作用下，从一种稳定状态突然跃迁到另一种稳定状态的过程。在实际应用中，当研究麦蚜种群数量的突变时，若将温度作为控制参数，当温度逐渐变化并达到某个临界值时，麦蚜种群数量可能会突然发生急剧增长或减少，这种现象就可以用折叠型突变模型来进行分析和解释。尖点型突变模型则涉及两个控制参数，能够更复杂地描述系统状态的突变情况。比如在考虑麦蚜种群动态时，同时考虑温度和湿度这两个控制参数，当它们在一定范围内变化时，麦蚜种群可能保持相对稳定；但当温度和湿度的变化达到特定的组合时，麦蚜种群就可能发生突变，尖点型突变模型能够很好地刻画这种复杂的关系。广义系统理论是一类形式更为一般化且具有广泛应用背景的动力系统理论，自20世纪70年代以来，其研究取得了长足的进展，许多正常系统的结论被相继推广到广义系统中。广义系统包含了大量带有时滞的广义系统，它比正常时滞系统更为广泛，产生于电机、经济中的投入产出、计量经济、环境污染、宇宙飞船等多种模型中。在麦蚜种群研究中，广义系统理论能够将麦蚜种群与其他相关因素，如作物种群、天敌种群以及环境因素等，视为一个相互关联的整体系统进行分析。从状态空间描述的角度来看，连续广义系统的一般表达形式为E(t)\dot{x}(t)=f(t,x(t),u(t))，y(t)=g(t,x(t),u(t))。其中，t为时间变量，t\geq0；E(t)一般为不可逆矩阵，且E(t)\inR^{n\timesn}；x(t)\inR^{n}是状态向量，代表了麦蚜种群数量、作物生物量、天敌数量等系统状态变量；u(t)\inR^{m}是输入向量，可表示外界对系统的干扰，如农药使用、气候变化等因素；y(t)\inR^{r}是输出向量，例如可以是麦蚜对小麦的危害程度、小麦的产量等可观测的结果变量；f(t,x(t),u(t))和g(t,x(t),u(t))分别为关于t、x(t)、u(t)的向量函数，它们描述了系统状态变量随时间的变化规律以及输入与输出之间的关系。离散广义系统的状态方程与连续广义系统相对应，其一般表达形式为E(k)x(k+1)=f(k,x(k),u(k))，y(k)=g(k,x(k),u(k))。其中，k为时间变量，k\inN；其他量的含义与连续广义系统类似，只是时间变量变为离散的整数形式，更适用于描述时间以离散间隔进行测量的系统，如在定期监测麦蚜种群数量的情况下，就可以采用离散广义系统模型来分析其动态变化。3.2两种群广义麦蚜生态系统模型构建3.2.1模型假设与变量定义为了构建两种群广义麦蚜生态系统模型，做出以下合理假设：麦蚜种群：麦蚜种群数量的增长不仅受到自身密度的制约，还受到作物种群的影响。随着麦蚜种群数量的增加，种内竞争加剧，对食物和生存空间的争夺会限制其增长速度。同时，作物的生长状况和生物量也会直接影响麦蚜的食物资源，进而影响麦蚜种群的增长。假设麦蚜种群的出生率和死亡率是其种群数量和作物生物量的函数，且麦蚜种群在麦田中的分布是均匀的。作物种群：作物种群的生长受到麦蚜取食的影响，麦蚜吸食作物汁液会导致作物生长受到抑制，光合作用减弱，生物量积累减少。同时，作物种群的生长还受到环境因素（如光照、温度、水分、养分等）的影响，假设在模型中环境因素对作物生长的影响是通过一个综合的环境因子来体现的。作物种群具有一定的自我修复和补偿能力，在受到麦蚜一定程度的取食后，能够通过自身的生理调节机制来恢复部分生长能力。基于上述假设，定义以下变量：麦蚜种群数量：用x(t)表示t时刻麦蚜的种群数量，单位为头/平方米。x(t)是模型中的一个关键状态变量，其动态变化受到多种因素的影响，如自身的出生率、死亡率、迁移率以及作物种群的影响等。作物生物量：用y(t)表示t时刻作物的生物量，单位为克/平方米。作物生物量反映了作物的生长状况和资源状况，它不仅是麦蚜的食物来源，其自身的变化也会受到麦蚜取食以及环境因素的影响。麦蚜出生率：b(x,y)表示麦蚜在t时刻的出生率，单位为头/(头・天)。出生率是麦蚜种群增长的重要因素之一，它与麦蚜种群数量x和作物生物量y相关。当作物生物量充足时，麦蚜有足够的食物资源，出生率可能较高；随着麦蚜种群数量的增加，种内竞争加剧，出生率可能会下降。麦蚜死亡率：d(x,y)表示麦蚜在t时刻的死亡率，单位为头/(头・天)。死亡率同样与麦蚜种群数量x和作物生物量y有关。当麦蚜种群数量过多，食物资源匮乏时，死亡率可能会上升；作物生物量的变化也会影响麦蚜的生存环境，从而影响死亡率。作物生长率：g(y)表示作物在t时刻的生长率，单位为克/(克・天)。作物生长率反映了作物在自然状态下的生长速度，它主要取决于环境因素以及作物自身的生物学特性。在适宜的环境条件下，作物生长率较高；而当环境条件不利时，如光照不足、温度不适宜、水分或养分缺乏等，作物生长率会降低。麦蚜对作物的取食率：h(x,y)表示麦蚜在t时刻对作物的取食率，单位为克/(头・天)。取食率体现了麦蚜对作物的危害程度，它与麦蚜种群数量x和作物生物量y密切相关。麦蚜种群数量越多，对作物的取食就越严重；作物生物量的变化也会影响麦蚜的取食行为，当作物生物量较低时，麦蚜可能会更加集中地取食，取食率可能会相应提高。环境因子：用e(t)表示t时刻的环境因子，它是一个综合指标，反映了光照、温度、水分、养分等环境因素对作物生长和麦蚜生存的综合影响。环境因子e(t)的取值范围在0到1之间，1表示最适宜的环境条件，0表示最恶劣的环境条件。例如，当光照充足、温度适宜、水分和养分充足时，e(t)的值接近1；而当出现干旱、高温、病虫害等不利因素时，e(t)的值会降低。3.2.2模型构建过程基于生态学原理和上述假设，构建包含麦蚜种群和作物种群的两种群广义系统模型。对于麦蚜种群，其种群数量的变化率等于出生率减去死亡率，再加上可能的迁移率（在本模型中暂不考虑迁移率），即：\frac{dx(t)}{dt}=b(x(t),y(t))x(t)-d(x(t),y(t))x(t)出生率b(x,y)可以表示为一个与作物生物量y正相关，与麦蚜种群数量x负相关的函数，例如：b(x,y)=\frac{\alphay}{\beta+x}其中，\alpha表示作物生物量对麦蚜出生率的影响系数，\alpha越大，说明作物生物量对麦蚜出生率的促进作用越强；\beta表示麦蚜种群数量对出生率的抑制系数，\beta越大，说明麦蚜种群数量对自身出生率的抑制作用越明显。死亡率d(x,y)可以表示为一个与麦蚜种群数量x正相关，与作物生物量y负相关的函数，例如：d(x,y)=\gammax-\frac{\deltay}{\epsilon+x}其中，\gamma表示麦蚜种群数量对死亡率的影响系数，\gamma越大，说明麦蚜种群数量越多，死亡率越高；\delta表示作物生物量对麦蚜死亡率的抑制系数，\delta越大，说明作物生物量对麦蚜死亡率的降低作用越显著；\epsilon是一个调节参数，用于调整函数的形式。将出生率和死亡率代入麦蚜种群数量变化率的方程中，得到：\frac{dx(t)}{dt}=\frac{\alphay(t)}{\beta+x(t)}x(t)-(\gammax(t)-\frac{\deltay(t)}{\epsilon+x(t)})x(t)对于作物种群，其生物量的变化率等于生长率减去麦蚜的取食率，即：\frac{dy(t)}{dt}=g(y(t))y(t)-h(x(t),y(t))x(t)作物生长率g(y)可以表示为一个与环境因子e(t)正相关的函数，例如：g(y)=\mue(t)(1-\frac{y(t)}{K})其中，\mu表示环境因子对作物生长率的影响系数，\mu越大，说明环境条件对作物生长的促进作用越强；K表示作物在理想环境下的最大生物量，1-\frac{y(t)}{K}表示作物生长的自我限制项，当作物生物量y(t)接近最大生物量K时，生长率会逐渐降低。麦蚜对作物的取食率h(x,y)可以表示为一个与麦蚜种群数量x和作物生物量y都相关的函数，例如：h(x,y)=\frac{\thetax}{\varphi+y}其中，\theta表示麦蚜种群数量对取食率的影响系数，\theta越大，说明麦蚜种群数量越多，对作物的取食越严重；\varphi表示作物生物量对取食率的调节系数，\varphi越大，说明作物生物量越大，麦蚜对其取食的相对强度可能会降低。将作物生长率和麦蚜取食率代入作物生物量变化率的方程中，得到：\frac{dy(t)}{dt}=\mue(t)(1-\frac{y(t)}{K})y(t)-\frac{\thetax(t)}{\varphi+y(t)}x(t)综上，构建的两种群广义麦蚜生态系统模型为：\begin{cases}\frac{dx(t)}{dt}=\frac{\alphay(t)}{\beta+x(t)}x(t)-(\gammax(t)-\frac{\deltay(t)}{\epsilon+x(t)})x(t)\\\frac{dy(t)}{dt}=\mue(t)(1-\frac{y(t)}{K})y(t)-\frac{\thetax(t)}{\varphi+y(t)}x(t)\end{cases}这个模型综合考虑了麦蚜种群和作物种群之间的相互作用，以及环境因素对作物生长的影响，能够较为全面地描述两种群广义麦蚜生态系统的动态变化。通过对该模型的分析，可以深入研究麦蚜种群和作物种群在不同条件下的变化规律，以及它们之间的相互关系，为麦蚜的防治和小麦生产的管理提供理论依据。3.3三种群广义麦蚜生态系统模型构建3.3.1引入天敌种群的考虑在麦田生态系统中，天敌种群对麦蚜种群的调控作用至关重要，是维持生态平衡的关键因素之一。天敌昆虫通过捕食、寄生等方式，直接作用于麦蚜个体，从而有效控制麦蚜种群的数量增长。七星瓢虫作为麦蚜的常见捕食性天敌，其成虫和幼虫均以麦蚜为食，具有较强的捕食能力。在适宜的环境条件下，一头七星瓢虫成虫每天可捕食麦蚜50-100头，幼虫每天也能捕食20-50头。草蛉同样是重要的捕食性天敌，其幼虫和成虫对麦蚜的捕食量也相当可观。食蚜蝇幼虫以麦蚜为主要食物来源，在麦田中大量捕食麦蚜，对麦蚜种群数量的增长起到了明显的抑制作用。蚜茧蜂则是寄生性天敌的代表，它将卵产在麦蚜体内，卵孵化后，幼虫在麦蚜体内生长发育，逐渐消耗麦蚜的营养物质，最终导致麦蚜死亡。据研究表明，当蚜茧蜂的寄生率达到一定程度时，麦蚜的种群增长会受到显著抑制。若蚜茧蜂的寄生率达到30%以上，麦蚜种群的增长率会大幅下降，甚至出现负增长。天敌种群的存在还能够影响麦蚜的行为习性，从而间接控制麦蚜种群数量。麦蚜在感知到天敌存在时，会改变其取食、繁殖和迁移等行为。麦蚜可能会减少在叶片表面的活动时间，降低取食频率，以避免被天敌发现和捕食。这种行为的改变会影响麦蚜的生长发育和繁殖效率，进而对麦蚜种群的增长产生负面影响。天敌的存在还会促使麦蚜产生更多的有翅蚜，这些有翅蚜会迁移到其他区域，从而分散了麦蚜种群的分布，降低了局部地区麦蚜的种群密度。从生态系统的稳定性角度来看，引入天敌种群能够增强麦田生态系统的自我调节能力。当麦蚜种群数量增加时，天敌昆虫的食物资源丰富，其种群数量也会相应增加，从而加大对麦蚜的捕食和寄生压力，抑制麦蚜种群的进一步增长。相反，当麦蚜种群数量减少时，天敌昆虫的食物短缺，其种群数量也会随之下降，减少对麦蚜的控制压力，使得麦蚜种群有一定的恢复空间。这种动态的平衡关系使得麦田生态系统能够保持相对稳定，避免麦蚜种群的过度增长或灭绝。在实际农业生产中，由于不合理使用农药等人为因素，常常导致天敌种群数量减少，从而破坏了这种自然的调控机制。长期大量使用广谱性农药，在杀死麦蚜的同时，也会大量杀死天敌昆虫，使得麦蚜失去了有效的自然控制。当麦蚜种群失去天敌的制约后，容易出现爆发性增长，对小麦造成严重危害。为了恢复和维持麦田生态系统的平衡，有必要引入天敌种群，充分发挥其对麦蚜种群的调控作用。通过人工释放七星瓢虫、草蛉、蚜茧蜂等天敌昆虫，或者采取保护和营造适宜天敌生存的环境等措施，增加天敌种群的数量和多样性，能够有效地控制麦蚜种群数量，减少化学农药的使用，实现农业的可持续发展。3.3.2模型结构与方程建立在构建三种群广义麦蚜生态系统模型时，除了考虑麦蚜种群和作物种群外，还将天敌种群纳入其中，以更全面地描述麦田生态系统的动态变化。假设天敌种群对麦蚜种群的捕食作用符合HollingⅡ型功能反应，即天敌的捕食率随着麦蚜种群数量的增加而增加，但增加的速率逐渐减缓。同时，考虑到天敌种群的繁殖也受到自身密度的制约，以及麦蚜种群数量的影响。定义以下变量：麦蚜种群数量：仍用x(t)表示t时刻麦蚜的种群数量，单位为头/平方米。作物生物量：用y(t)表示t时刻作物的生物量，单位为克/平方米。天敌种群数量：用z(t)表示t时刻天敌的种群数量，单位为头/平方米。麦蚜出生率：b(x,y)表示麦蚜在t时刻的出生率，单位为头/(头・天)，与之前两种群模型中的定义相同，即b(x,y)=\frac{\alphay}{\beta+x}。麦蚜死亡率：d(x,y)表示麦蚜在t时刻的死亡率，单位为头/(头・天)，与之前定义相同，即d(x,y)=\gammax-\frac{\deltay}{\epsilon+x}。作物生长率：g(y)表示作物在t时刻的生长率，单位为克/(克・天)，与之前定义相同，即g(y)=\mue(t)(1-\frac{y(t)}{K})。麦蚜对作物的取食率：h(x,y)表示麦蚜在t时刻对作物的取食率，单位为克/(头・天)，与之前定义相同，即h(x,y)=\frac{\thetax}{\varphi+y}。天敌的捕食率：p(x,z)表示天敌在t时刻对麦蚜的捕食率，单位为头/(头・天)，采用HollingⅡ型功能反应，可表示为p(x,z)=\frac{\lambdaxz}{\omega+x}，其中\lambda表示天敌的捕食能力系数，\lambda越大，说明天敌的捕食能力越强；\omega表示半饱和常数，反映了天敌捕食行为对麦蚜种群数量变化的敏感程度。天敌的出生率：r(x,z)表示天敌在t时刻的出生率，单位为头/(头・天)，可表示为r(x,z)=\frac{\sigmax}{\tau+x}z，其中\sigma表示麦蚜种群数量对天敌出生率的影响系数，\sigma越大，说明麦蚜种群数量越多，越有利于天敌的繁殖；\tau是一个调节参数。天敌的死亡率：s(z)表示天敌在t时刻的死亡率，单位为头/(头・天)，可表示为s(z)=\rhoz，其中\rho表示天敌的自然死亡率系数。基于以上定义和假设，构建三种群广义麦蚜生态系统模型如下：\begin{cases}\frac{dx(t)}{dt}=b(x(t),y(t))x(t)-d(x(t),y(t))x(t)-p(x(t),z(t))x(t)\\\frac{dy(t)}{dt}=g(y(t))y(t)-h(x(t),y(t))x(t)\\\frac{dz(t)}{dt}=r(x(t),z(t))z(t)-s(z(t))z(t)\end{cases}将各函数表达式代入上述方程，得到：\begin{cases}\frac{dx(t)}{dt}=\frac{\alphay(t)}{\beta+x(t)}x(t)-(\gammax(t)-\frac{\deltay(t)}{\epsilon+x(t)})x(t)-\frac{\lambdax(t)z(t)}{\omega+x(t)}x(t)\\\frac{dy(t)}{dt}=\mue(t)(1-\frac{y(t)}{K})y(t)-\frac{\thetax(t)}{\varphi+y(t)}x(t)\\\frac{dz(t)}{dt}=\frac{\sigmax(t)}{\tau+x(t)}z(t)-\rhoz(t)z(t)\end{cases}这个模型综合考虑了麦蚜种群、作物种群和天敌种群之间的相互作用关系，以及环境因素对作物生长的影响。通过对该模型的分析，可以深入研究三种群广义麦蚜生态系统的动态变化规律，以及各因素对麦蚜种群突变的影响，为制定有效的麦蚜防治策略提供更全面、准确的理论依据。四、麦蚜种群广义系统模型的分析4.1模型平衡点分析对于构建的两种群广义麦蚜生态系统模型：\begin{cases}\frac{dx(t)}{dt}=\frac{\alphay(t)}{\beta+x(t)}x(t)-(\gammax(t)-\frac{\deltay(t)}{\epsilon+x(t)})x(t)\\\frac{dy(t)}{dt}=\mue(t)(1-\frac{y(t)}{K})y(t)-\frac{\thetax(t)}{\varphi+y(t)}x(t)\end{cases}令\frac{dx(t)}{dt}=0，\frac{dy(t)}{dt}=0，求解平衡点。由由\frac{dx(t)}{dt}=0可得：\frac{\alphay}{\beta+x}x-(\gammax-\frac{\deltay}{\epsilon+x})x=0即x(\frac{\alphay}{\beta+x}-\gamma+\frac{\deltay}{\epsilon+x})=0则x=0或\frac{\alphay}{\beta+x}-\gamma+\frac{\deltay}{\epsilon+x}=0（1）当x=0时，代入\frac{dy(t)}{dt}=0，即\mue(1-\frac{y}{K})y-\frac{\thetax}{\varphi+y}x=0，可得\mue(1-\frac{y}{K})y=0，解得y=0或y=K。所以得到平衡点所以得到平衡点(0,0)和(0,K)。当\frac{\alphay}{\beta+x}-\gamma+\frac{\deltay}{\epsilon+x}=0时，化简可得：\alphay(\epsilon+x)-\gamma(\beta+x)(\epsilon+x)+\deltay(\beta+x)=0\alphay\epsilon+\alphayx-\gamma(\beta\epsilon+\betax+\epsilonx+x^2)+\deltay\beta+\deltayx=0(\alpha\epsilon+\delta\beta)y+(\alpha+\delta)yx-\gamma(\beta\epsilon+(\beta+\epsilon)x+x^2)=0（2）再由\frac{dy(t)}{dt}=0可得：\mue(1-\frac{y}{K})y-\frac{\thetax}{\varphi+y}x=0\muey-\frac{\muey^2}{K}-\frac{\thetax^2}{\varphi+y}=0\muey(\varphi+y)-\frac{\muey^2(\varphi+y)}{K}-\thetax^2=0\mue\varphiy+\muey^2-\frac{\mue\varphiy^2}{K}-\frac{\muey^3}{K}-\thetax^2=0（3）联立方程（2）和（3）求解，过程较为复杂，可通过数值方法或特定参数取值来求解其他可能的平衡点。对于平衡点的稳定性分析，计算系统在平衡点处的雅可比矩阵J。J=\begin{pmatrix}\frac{\partial\frac{dx}{dt}}{\partialx}&\frac{\partial\frac{dx}{dt}}{\partialy}\\\frac{\partial\frac{dy}{dt}}{\partialx}&\frac{\partial\frac{dy}{dt}}{\partialy}\end{pmatrix}\frac{\partial\frac{dx}{dt}}{\partialx}=\frac{\alphay(\beta+x)-\alphaxy}{(\beta+x)^2}-\gamma+\frac{\deltay(\epsilon+x)-\deltaxy}{(\epsilon+x)^2}\frac{\partial\frac{dx}{dt}}{\partialy}=\frac{\alphax}{\beta+x}+\frac{\deltax}{\epsilon+x}\frac{\partial\frac{dy}{dt}}{\partialx}=-\frac{2\thetax}{\varphi+y}\frac{\partial\frac{dy}{dt}}{\partialy}=\mue(1-\frac{2y}{K})+\frac{\thetax^2}{(\varphi+y)^2}将平衡点(0,0)代入雅可比矩阵J：J_{(0,0)}=\begin{pmatrix}-\gamma&0\\0&\mue\end{pmatrix}其特征值为\lambda_1=-\gamma，\lambda_2=\mue。当当\gamma>0且\mue>0时，(0,0)是鞍点，不稳定。将平衡点(0,K)代入雅可比矩阵J：J_{(0,K)}=\begin{pmatrix}-\gamma+\frac{\deltaK}{\epsilon}&\frac{\alpha+\delta}{\beta+\epsilon}K\\0&-\mue\end{pmatrix}其特征值为\lambda_1=-\gamma+\frac{\deltaK}{\epsilon}，\lambda_2=-\mue。当当-\gamma+\frac{\deltaK}{\epsilon}<0且-\mue<0时，(0,K)是稳定的节点。对于三种群广义麦蚜生态系统模型：\begin{cases}\frac{dx(t)}{dt}=\frac{\alphay(t)}{\beta+x(t)}x(t)-(\gammax(t)-\frac{\deltay(t)}{\epsilon+x(t)})x(t)-\frac{\lambdax(t)z(t)}{\omega+x(t)}x(t)\\\frac{dy(t)}{dt}=\mue(t)(1-\frac{y(t)}{K})y(t)-\frac{\thetax(t)}{\varphi+y(t)}x(t)\\\frac{dz(t)}{dt}=\frac{\sigmax(t)}{\tau+x(t)}z(t)-\rhoz(t)z(t)\end{cases}令\frac{dx(t)}{dt}=0，\frac{dy(t)}{dt}=0，\frac{dz(t)}{dt}=0，求解平衡点。由由\frac{dx(t)}{dt}=0可得：\frac{\alphay}{\beta+x}x-(\gammax-\frac{\deltay}{\epsilon+x})x-\frac{\lambdaxz}{\omega+x}x=0x(\frac{\alphay}{\beta+x}-\gamma+\frac{\deltay}{\epsilon+x}-\frac{\lambdaz}{\omega+x})=0则x=0或\frac{\alphay}{\beta+x}-\gamma+\frac{\deltay}{\epsilon+x}-\frac{\lambdaz}{\omega+x}=0（4）当x=0时，代入\frac{dy(t)}{dt}=0可得\mue(1-\frac{y}{K})y=0，解得y=0或y=K；代入\frac{dz(t)}{dt}=0可得-\rhoz^2=0，解得z=0。所以得到平衡点所以得到平衡点(0,0,0)和(0,K,0)。当\frac{\alphay}{\beta+x}-\gamma+\frac{\deltay}{\epsilon+x}-\frac{\lambdaz}{\omega+x}=0时，结合\frac{dy(t)}{dt}=0和\frac{dz(t)}{dt}=0，通过复杂的代数运算和可能的数值方法求解其他平衡点。对于平衡点稳定性分析，计算系统在平衡点处的雅可比矩阵J：J=\begin{pmatrix}\frac{\partial\frac{dx}{dt}}{\partialx}&\frac{\partial\frac{dx}{dt}}{\partialy}&\frac{\partial\frac{dx}{dt}}{\partialz}\\\frac{\partial\frac{dy}{dt}}{\partialx}&\frac{\partial\frac{dy}{dt}}{\partialy}&\frac{\partial\frac{dy}{dt}}{\partialz}\\\frac{\partial\frac{dz}{dt}}{\partialx}&\frac{\partial\frac{dz}{dt}}{\partialy}&\frac{\partial\frac{dz}{dt}}{\partialz}\end{pmatrix}分别计算各偏导数：\frac{\partial\frac{dx}{dt}}{\partialx}=\frac{\alphay(\beta+x)-\alphaxy}{(\beta+x)^2}-\gamma+\frac{\deltay(\epsilon+x)-\deltaxy}{(\epsilon+x)^2}-\frac{\lambdaz(\omega+x)-\lambdaxz}{(\omega+x)^2}\frac{\partial\frac{dx}{dt}}{\partialy}=\frac{\alphax}{\beta+x}+\frac{\deltax}{\epsilon+x}\frac{\partial\frac{dx}{dt}}{\partialz}=-\frac{\lambdax^2}{\omega+x}\frac{\partial\frac{dy}{dt}}{\partialx}=-\frac{2\thetax}{\varphi+y}\frac{\partial\frac{dy}{dt}}{\partialy}=\mue(1-\frac{2y}{K})+\frac{\thetax^2}{(\varphi+y)^2}\frac{\partial\frac{dy}{dt}}{\partialz}=0\frac{\partial\frac{dz}{dt}}{\partialx}=\frac{\sigmaz(\tau+x)-\sigmaxz}{(\tau+x)^2}\frac{\partial\frac{dz}{dt}}{\partialy}=0\frac{\partial\frac{dz}{dt}}{\partialz}=\frac{\sigmax}{\tau+x}-2\rhoz将平衡点(0,0,0)代入雅可比矩阵J：J_{(0,0,0)}=\begin{pmatrix}-\gamma&0&0\\0&\mue&0\\0&0&0\end{pmatrix}其特征值为\lambda_1=-\gamma，\lambda_2=\mue，\lambda_3=0。由于存在零特征值，(0,0,0)的稳定性需进一步分析。将平衡点(0,K,0)代入雅可比矩阵J，并通过分析特征值的实部来判断其稳定性。通过对两种群和三种群广义麦蚜生态系统模型平衡点的分析，可以了解系统在不同状态下的稳定性，为进一步研究麦蚜种群的动态变化提供重要的理论依据。4.2模型稳定性分析运用非线性广义系统稳定性理论，对构建的麦蚜种群广义系统模型进行稳定性分析，以探究系统在不同条件下的稳定状态，这对于深入理解麦蚜种群动态变化机制至关重要。对于两种群广义麦蚜生态系统模型，首先考虑其平衡点的稳定性。前文已求得平衡点(0,0)和(0,K)，通过计算系统在平衡点处的雅可比矩阵J，并分析其特征值来判断稳定性。在平衡点(0,0)处，雅可比矩阵J_{(0,0)}=\begin{pmatrix}-\gamma&0\\0&\mue\end{pmatrix}，特征值为\lambda_1=-\gamma，\lambda_2=\mue。当\gamma>0且\mue>0时，(0,0)是鞍点，不稳定。这意味着在该平衡点附近，系统的状态容易发生变化，麦蚜种群数量和作物生物量不会稳定在(0,0)状态。在平衡点(0,K)处，雅可比矩阵J_{(0,K)}=\begin{pmatrix}-\gamma+\frac{\deltaK}{\epsilon}&\frac{\alpha+\delta}{\beta+\epsilon}K\\0&-\mue\end{pmatrix}，特征值为\lambda_1=-\gamma+\frac{\deltaK}{\epsilon}，\lambda_2=-\mue。当-\gamma+\frac{\deltaK}{\epsilon}<0且-\mue<0时，(0,K)是稳定的节点，即系统在该平衡点处能够保持相对稳定，麦蚜种群数量趋近于0，作物生物量稳定在K，表明此时麦蚜对作物的影响较小，作物能够在相对稳定的状态下生长。从生物学意义角度深入分析，\gamma表示麦蚜种群数量对死亡率的影响系数，\gamma>0说明随着麦蚜种群数量增加，其死亡率上升，这符合一般的生态规律，过多的麦蚜个体竞争有限的资源，导致生存压力增大，死亡率升高。\mue体现了环境因子对作物生长率的影响，\mue>0表明适宜的环境条件有利于作物生长。在平衡点(0,0)不稳定的情况下，系统会向其他状态发展，这可能是由于环境中存在一些因素促使麦蚜种群或作物种群发生变化，比如适宜的气候条件可能导致麦蚜繁殖，或者有外来的麦蚜迁入，使得麦蚜种群数量不再保持为0。而在平衡点(0,K)稳定时，说明在当前环境和生物因素的综合作用下，麦蚜种群得到了有效控制，作物能够充分利用资源生长到最大生物量K。对于三种群广义麦蚜生态系统模型，同样考虑平衡点的稳定性。前文已得到平衡点(0,0,0)和(0,K,0)，计算系统在平衡点处的雅可比矩阵J。在平衡点(0,0,0)处，雅可比矩阵J_{(0,0,0)}=\begin{pmatrix}-\gamma&0&0\\0&\mue&0\\0&0&0\end{pmatrix}，特征值为\lambda_1=-\gamma，\lambda_2=\mue，\lambda_3=0。由于存在零特征值，其稳定性需进一步分析，零特征值的出现表明系统在该平衡点处可能存在一些特殊的动态行为，可能处于一种临界状态，微小的扰动就可能导致系统状态发生改变。在平衡点(0,K,0)处，将其