不同氮水平下大豆与玉米单、间作氮素吸收的微分方程建模.pdf

Advances in Applied Mathematics 应用数学进展应用数学进展,2024,13(1),313-330 Published Online January 2024 in Hans.https:/www.hanspub.org/journal/aam https:/doi.org/10.12677/aam.2024.131034 文章引用文章引用:普雯瑞,化存才.不同氮水平下大豆与玉米单、间作氮素吸收的微分方程建模J.应用数学进展,2024,13(1):313-330.DOI:10.12677/aam.2024.131034 不同氮水平下大豆与玉米单、间作氮素吸收的不同氮水平下大豆与玉米单、间作氮素吸收的微分方程建模微分方程建模 普雯瑞普雯瑞，化存才化存才*云南师范大学数学学院，云南 昆明 收稿日期：2023年12月19日；录用日期：2024年1月13日；发布日期：2024年1月23日 摘摘 要要 大豆大豆/玉米间作具有一定的养分利用优势，但是不同氮水平下的间作体系对氮素吸收的响应不同。在本文玉米间作具有一定的养分利用优势，但是不同氮水平下的间作体系对氮素吸收的响应不同。在本文中，首先分别建立了自然条件下大豆和玉米单作氮素吸收的非自治中，首先分别建立了自然条件下大豆和玉米单作氮素吸收的非自治Logistic微分方程模型，探讨了不同微分方程模型，探讨了不同氮水平下大豆和玉米单作时氮素吸收的非自治氮水平下大豆和玉米单作时氮素吸收的非自治Logistic微分方程模型，得到单作大豆和玉米的最佳施氮微分方程模型，得到单作大豆和玉米的最佳施氮量分别为量分别为0.4808 gplant1、2.6199 gplant1。然后，进一步建立了不同氮水平下大豆和玉米间作氮素。然后，进一步建立了不同氮水平下大豆和玉米间作氮素吸收的二维非自治微分方程模型，对模型拟合的相关系数均在吸收的二维非自治微分方程模型，对模型拟合的相关系数均在98%以上，数值计算得到间作大豆和玉米以上，数值计算得到间作大豆和玉米的最佳施氮量分别为的最佳施氮量分别为0.4821 gplant1、2.6201 gplant1，与单作结果差别不大，然而，两作物的氮素，与单作结果差别不大，然而，两作物的氮素吸收量在不同氮水平交叉种植模式下却有着显著差异。在最佳施氮量下，间作玉米的氮素吸收量最大为吸收量在不同氮水平交叉种植模式下却有着显著差异。在最佳施氮量下，间作玉米的氮素吸收量最大为4.48 gplant1，较自然间作提升了，较自然间作提升了59.43%，而大豆的氮素吸收量为，而大豆的氮素吸收量为1.31 gplant1，较自然间作也提升，较自然间作也提升了了32.32%；就氮肥利用率而言，最佳施氮量下间作大豆相比最佳单作大豆下降了；就氮肥利用率而言，最佳施氮量下间作大豆相比最佳单作大豆下降了24.01%，而最佳间作，而最佳间作玉米相比最佳单作玉米则提高了玉米相比最佳单作玉米则提高了30.45%。关键词关键词 大豆大豆/玉米间作，氮素吸收，非自治微分方程，最佳施氮量，氮利用率玉米间作，氮素吸收，非自治微分方程，最佳施氮量，氮利用率 Differential Equation Modeling of Nitrogen Uptake in Soybean and Maize Monoculture and Intercropping under Different Nitrogen Levels Wenrui Pu,Cuncai Hua*School of Mathematics,Yunnan Normal University,Kunming Yunnan *通讯作者。普雯瑞，化存才 DOI:10.12677/aam.2024.131034 314 应用数学进展 Received:Dec.19th,2023;accepted:Jan.13th,2024;published:Jan.23rd,2024 Abstract Soybean/corn intercropping has certain advantages in nutrient utilization,but the response of in-tercropping systems to nitrogen uptake is different under different nitrogen levels.In this paper,non-autonomous Logistic differential equation models of nitrogen absorption of monoculture soy-bean and corn under natural conditions were established respectively,and non-autonomous Logis-tic differential equation models of nitrogen absorption of soybean and corn under different nitrogen levels were discussed at first.The optimal nitrogen application rates for the soybean and maize were 0.4808 gplant1 and 2.6199 gplant1,respectively.Then,a two-dimensional non-autonomous dif-ferential equation model was established further for the nitrogen uptake between soybean and corn at different nitrogen levels.The correlation coefficients of the model were all above 98%.The optimal nitrogen application rates for intercropping soybean and corn were 0.4821 gplant1 and 2.6201 gplant1 respectively,which was no significant difference between the intercrops and mo-noculture results.However,the nitrogen uptake of the two crops was significantly different under different nitrogen levels of cross-cropping.Under the optimal nitrogen application rate,the maxi-mum nitrogen uptake of intercropping maize was 4.48 gplant1,which was 59.43%higher than that of natural intercropping,while the nitrogen uptake of soybean was 1.31 gplant1,which was 32.32%higher than that of natural intercropping.In terms of the nitrogen utilization ratios,the best inter-cropping soybean decreased by 24.01%by comparing with the best monoculture soybean,while the best intercropping corn increased by 30.45%by comparing with the best monoculture corn.Keywords Soybean/Maize Intercropping,Nitrogen Uptake,Non-Autonomous Differential Equation,Optimal Nitrogen Application Rate,Nitrogen Utilization Rate Copyright 2024 by author(s)and Hans Publishers Inc.This work is licensed under the Creative Commons Attribution International License(CC BY 4.0).http:/creativecommons.org/licenses/by/4.0/1.引言引言 大豆、玉米是我国极其重要的粮油饲作物。氮素是植物生长过程中所必需的第一元素。肥料的利用在很大程度上解决了植物生长过程中的需氮问题，然而，氮肥施用量大且利用率偏低已给我国带来了严重的环境问题1。研究表明，豆科与非豆科作物的间作已成为一种农业可持续发展模式，它不仅减少了化学肥料的投入，减轻农业生产对环境的负面影响，而且还提高了各作物的氮素吸收量及其产量2，比如大豆/玉米间作3、玉米/花生间作4、玉米/毛苕子间作5等都使得间作玉米的氮素吸收积累量较单作时有了不同程度的提升，是一种“共赢”的种植模式。但是，在不同氮水平下间作体系的养分利用效率表现不同。当水稻/花生间作6时，在不同氮肥施用量下，花生体内有 4%10%的氮素能够准确转移到水稻植株体内。大豆/玉米间作研究表明7，在施氮量180 kghm2水平下，单位面积间作玉米的氮素最大吸收量相比于单作提高了 18.4%，间作大豆的则比单作处理降低了 15.9%，而整个间作系统的氮素则提高了 0.4%。总之，禾本科与豆科间作能够显著增加作物氮吸收优势，减轻了豆科作物自身的“氮遏制”，提高了豆科作物的固氮效率和间作体系的氮素利用率，进而增加了其干物质的累积量和产量。Open AccessOpen Access普雯瑞，化存才 DOI:10.12677/aam.2024.131034 315 应用数学进展 在数学交叉生物应用的研究中，将植物或农作物生长发育上升为微分方程建模研究具有极其重要的理论和应用意义。早在 1926 年，意大利数学家 Volterra 就两种群提出了 Volterra 食饵捕食者微分方程模型8。2012 年，张林、化存才和张继业等考虑到 UV-B(波长为 280320 nm 的紫外光)辐射对烤烟株高的影响，建立了非自治 Logistic 模型9。2021 年，梁仕军考虑到花生间作可缓解连作障碍问题，建立了花生与玉米间作生长的二维微分方程模型10，数值模拟显示：花生和玉米总产量增加了 11.15%，验证了间作的增产优势。2023 年，段国梅考虑到作物的氮素吸收是农作物提升产量的关键机制，选取大豆和玉米这两种作物，建立了在自然条件间作模式下氮素吸收量的二维微分方程模型11，数值模拟显示：大豆和玉米两者的氮素吸收总量比起两者单作时的总量增加了 5.33%。综上所述，关于大豆、玉米间作时氮素吸收的微分方程建模才刚开始，没有从数学本质上阐明间作减少供氮水平或者有效提高氮利用率的生态机理问题。为此，本文在考虑到大豆本身的固氮作用和不同供氮水平对植物的氮素吸收产生显著影响的基础上，主要地建立在不同氮水平下大豆、玉米间作时氮素吸收的微分方程组模型，数值计算阐明间作减少供氮水平或有效提高氮利用率等问题。2.不同氮水平下大豆单作氮素吸收的非自治不同氮水平下大豆单作氮素吸收的非自治 Logistic 微分方程模型微分方程模型 大豆是我国食用油和蛋白粕饲料的主要来源，其需氮量是其他谷类作物的 4 倍12，而在实际中施加于土壤中的氮养分情况对于大豆的氮素吸收具有直接的重要影响。在本节中，将基于农学实验的一些数据，结合施氮量对于大豆氮素吸收积累量的影响，把大豆在每个生长阶段的氮素吸收积累量与其生长的时间节点联系起来，建立大豆在不同氮水平下的非自治 Logistic 微分方程模型，分析单作大豆在不同氮水平下的氮素吸收机制，以揭示大豆的需肥规律，为实现大豆的合理施肥、提高其氮素吸收积累量及氮肥利用率提供理论依据。2.1.模型的建立模型的建立 通过分析大豆在整个生育期内的氮素吸收积累量13，我们总结出以下基本事实：(F1)大豆对氮素的吸收积累量是大豆产量提高的基础，也是影响氮肥利用率的主要因素；(F2)大豆的氮素吸收积累过程符合“慢快慢”的“S”型曲线增长生长规律，并且在种植 82102天内氮素积累能力达到最快；(F3)随着时间的推移，大豆的氮素吸收积累量最终达到一个饱和值(最大值)。基于以上事实，我们作出如下合理假设：(H1)单株大豆氮素吸收积累量的饱和值为一个常数；(H2)大豆的氮素吸收积累受自身条件以及土壤肥力水平、土壤 PH 值、灌溉条件等多种因素的综合影响，由此可设其内禀增长率是一个随时间呈“先增后减”变化趋势的二次函数；(H3)在整个生育期内大豆的氮素吸收积累速度与环境中所剩的资源成正比。根据上述基本事实和合理假设，现设t时刻单株大豆的氮素吸收积累量为()S t，单位为 g/plant，最大值为mS，大豆氮素吸收的内禀增长率为()()21192Srtt=+，其中1，1为待定参数且10，于是，在自然条件下，我们建立单株大豆氮素吸收积累的非自治 Logistic 模型的如下初值问题：()()()21100d1921dSmmSSSrt StStSSS tS=+=(1)由于不同施氮量处理下大豆氮素吸收积累量仍然满足Logistic微分方程，故借鉴文献14的建模思想，我们基于模型(1)考虑不同氮水平，建立大豆氮素吸收积累的微分方程模型：普雯瑞，化存才 DOI:10.12677/aam.2024.131034 316 应用数学进展 ()d1dS SmSSh rt StS=(2)其中，Sh为有关施氮量的影响因子，最大值mSmSh S=，其他参数同模型(1)。2.2.参数拟合值参数拟合值 由文献13易得，大豆在苗期、开花期、结荚期、鼓粒期、成熟期各个生育阶段中的氮素吸收积累量的数据，见表 1：Table 1.Nitrogen uptake accumulation data of soybean monoculture under natural conditions 表表 1.自然条件下单作大豆氮素吸收积累量数据 t/day 48 54 76 93 109 128 132 S/gplant1 0.08 0.12 0.38 0.74 1.07 1.31 1.34 由表 1 可知，模型(1)的初值条件为()400.05S=。利用 OpenLu 编程软件拟合得其中的参数值为：1=0.000001019，1=0.0627，1.45mS=。根据文献13整理得到不同施氮量下大豆最大氮素吸收积累量的数据，见表 2。Table 2.Nitrogen application rate and maximum nitrogen uptake accumula-tion in soybean 表表 2.大豆施氮量与其最大氮素吸收积累量 施氮量：x/gplant1 0 0.20 0.40 0.61 最大氮素吸收积累量：mS/gplant1 1.45 1.63 1.80 1.75 由表 2 得：在一定施氮量范围内，大豆的氮素吸收积累量随着施氮量的增加而增大；而当施氮量超过一定的范围后，大豆的氮素吸收积累量随着施氮量的增加而减少。利用 MATLAB 编程对表 2 中的数据进行拟合，得到如下的三次拟合函数：323.87462.14290.62641.45mSxxx=+(3)拟合函数的曲线见图 1。对式(3)求导，当取施氮量为 x=0.4808 gplant1时，大豆的氮素吸收积累量能够达到最大值，此时的施氮量即为影响单作大豆氮素吸收积累量的最佳施氮量。2.3.数值模拟结论数值模拟结论 拟合效果：根据文献13的相关数据，我们整理得到了不同供氮水平下，单作大豆的氮素吸收积累量的相关数据见表 3。Table 3.Nitrogen uptake and accumulation of monocropping soybean at different nitrogen levels(unit:gplant1)表表 3.不同氮水平下单作大豆的氮素吸收积累量(单位：gplant1)时间 day 施氮量 40 48 54 76 93 109 128 132 0 gplant1 0.05 0.08 0.12 0.38 0.74 1.07 1.31 1.34 0.2 gplant1 0.08 0.12 0.16 0.47 0.85 1.20 1.46 1.49 普雯瑞，化存才 DOI:10.12677/aam.2024.131034 317 应用数学进展 续表 0.4 gplant1 0.08 0.13 0.18 0.53 0.96 1.34 1.62 1.66 0.61 gplant1 0.08 0.13 0.17 0.48 0.86 1.23 1.53 1.57 借助上述表 3 中的数据，通过 MATLAB 编程计算分别得到不同氮水平下单作大豆的实验数据与模拟数据之间的拟合程度见表 4。Table 4.Table of fitting effects under different nitrogen application rates 表表 4.不同施氮量下的拟合效果表 施氮量(x/gplant1)0 0.2 0.4 0.61 相关系数(R2)0.9998 0.9994 0.9974 0.9960 数值计算模拟：利用 MATLAB 编程对模型(2)进行数值计算模拟，得到在不同氮水平下单作大豆氮素吸收动态变化的曲线图，见图 2。由图 2 我们可得如下的结论：1)不同施氮量下，大豆整株的氮素吸收积累量均符合 S 型曲线生长特征，且其氮素吸收积累依然主要集中在结荚鼓粒期。2)当施氮量为 0.4808 gplant1时，大豆的氮素吸收积累速度比施氮量为 0 gplant1，0.2 gplant1，0.4 gplant1和 0.61 gplant1时快，说明此时的施氮量即为影响单作大豆氮素吸收积累量的最佳施氮量，与前面的理论推导是一致的。3)单作大豆在生长到 200 天时的氮素吸收积累量已接近各不同施氮量下的最大值，见表 5。Figure 1.Relationship between nitrogen fertilization rate and maximum nitrogen uptake accumulation in soybean 图图 1.大豆施氮量与其氮素吸收积累量最大值的关系 普雯瑞，化存才 DOI:10.12677/aam.2024.131034 318 应用数学进展 Figure 2.Changes in nitrogen uptake and accumulation of monoculture soybean under different nitrogen levels 图图 2.不同氮水平下单作大豆氮素吸收积累量变化 Table 5.Nitrogen uptake and accumulation in 200 days of soybean with different nitrogen application rates and monoculture 表表 5.不同施氮量与单作大豆 200 天时氮素吸收积累量 施氮量(x/gplant1)0 0.2 0.4 0.4808 0.61 200 天时氮素吸收积累量(Sm/gplant1)1.45 1.63 1.80 1.8155 1.75 3.不同氮水平下玉米单作氮素吸收的非自治不同氮水平下玉米单作氮素吸收的非自治 Logistic 微分方程模型微分方程模型 玉米作为我国的第一大粮食作物，保证其产量是我国粮食安全问题的重中之重问题15。氮元素作为植物生长过程中的第一元素，对玉米的产量有着举足轻重的意义。在本节中，将类似于前面第 2 节，基于一些农学实验数据，结合施氮量对玉米氮素吸收积累量的影响，把玉米每个阶段的氮素吸收积累量与其生长的时间节点联系起来，建立玉米在不同氮水平下的非自治 Logistic 微分方程模型，分析单作玉米在不同氮水平下的氮素吸收机制，揭示玉米的需肥规律，为实现玉米的合理施肥、提高其氮素吸收积累量及氮肥利用率提供理论依据。3.1.模型的建立模型的建立 通过分析玉米在整个生育期内的氮素吸收积累量13，我们总结出有关玉米的一些基本事实与大豆的大致相同，而其在种植 6484 天内积累能力达到最快。基于上述事实，我们作出如下合理假设：(H1)单株玉米氮素吸收积累量的饱和值为一个常数；(H2)玉米的氮素吸收积累受自身条件以及土壤肥力水平、土壤 PH 值、灌溉条件等多种因素的综合影响，由此可设其内禀增长率是一个随时间呈“先增后减”变化趋势的二次函数；(H3)在整个生育期内玉米的氮素吸收积累速度与环境中所剩的资源成正比。普雯瑞，化存才 DOI:10.12677/aam.2024.131034 319 应用数学进展 根据上述基本事实和合理假设，现设t时刻单株玉米的氮素吸收积累量为()M t，单位为 g/plant，最大值为mM，大豆氮素吸收的内禀增长率为()()22274Mrtt=+，其中2，2为待定参数且20，于是，我们建立单株玉米在自然条件下氮素吸收积累的非自治 Logistic 模型的如下初值问题：()()()22200d1741dMmmMMMrt MtMtMMM tM=+=(4)在不同施氮量处理下，玉米氮素吸收积累量的 Logistic 微分方程模型为：()d1dM MmMMh rt MtM=(5)其中Mh为有关施氮量的影响因子，最大值为mMmMh M=，其他参数同模型(4)。3.2.参数拟合值参数拟合值 王雪蓉在文献13中给出了试验研究得到的玉米氮素吸收积累量的实验数据。现将玉米在各生育期时氮素吸收的特点和实验数据进行比对，选取苗期、拔节期、抽雄期、灌浆期、成熟期等各个阶段相应的数据，见表 6：Table 6.Nitrogen uptake accumulation of maize monoculture under natural conditions 表表 6.自然条件下单作玉米氮素吸收积累量数据 t/day 34 48 70 85 99 109 126 140 M/gplant1 0.06 0.20 0.89 1.60 2.01 2.15 2.23 2.25 从表 6 中得到模型(4)的初值条件为()340.06M=，再利用 OpenLu 软件编程拟合得到了各个参数的最佳值为：2=0.000001615，2=0.0873，2.26mM=。根据文献13整理得到不同施氮量下玉米最大氮素吸收积累量的数据见表 7。Table 7.Data of nitrogen application rate and maximum nitrogen uptake ac-cumulation of maize 表表 7.玉米施氮量与其最大氮素吸收积累量数据 施氮量(x/gplant 1)0 1.84 2.45 3.06 最大氮素吸收积累量(mM/gplant1)2.26 2.97 3.5 3.26 由表 7 得：在一定施氮量范围内，玉米的氮素吸收积累量随着施氮量的增加而增大；而当施氮量超过一定的范围后，玉米的氮素吸收积累量随着施氮量的增加而减少。利用 MATLAB 编程对表 7 中的数据进行拟合，得到三次拟合函数为：320.40261.92411.79162.26mMxxx=+(6)拟合函数的曲线见图 3。对(6)求导得，当取施氮量为 x=2.6199 gplant1时，玉米的氮素吸收积累量达到最大值，此时的施氮量即为影响单作玉米氮素吸收积累量的最佳施氮量。3.3.数值模拟结论数值模拟结论 拟合效果：根据文献内容，整理得到了不同供氮水平下，单作玉米的氮素吸收积累量的相关数据，见表 8。普雯瑞，化存才 DOI:10.12677/aam.2024.131034 320 应用数学进展 Table 8.Nitrogen uptake and accumulation of monocropping soybean under different nitrogen levels(unit:gplant1)表表 8.不同氮水平下单作大豆的氮素吸收积累量(单位：gplant1)时间 day 施氮量 34 48 70 85 99 109 126 140 0 gplant1 0.06 0.20 0.89 1.60 2.01 2.15 2.23 2.25 1.84 gplant1 0.11 0.30 1.17 2.02 2.57 2.78 2.92 2.95 2.45 gplant1 0.14 0.38 1.40 2.37 3.01 3.26 3.43 3.48 3.06 gplant1 0.11 0.31 1.22 2.17 2.79 3.02 3.17 3.21 利用上述表 8 数据，通过 MATLAB 编程计算，分别得到不同氮水平下单作玉米的实验数据与模拟数据之间的拟合程度见表 9。Table 9.Table of fitting effects under different nitrogen application rates 表表 9.不同施氮量下的拟合效果表 施氮量(x/gplant1)0 0.2 0.4 0.61 相关系数(R2)0.9999 0.9821 0.9571 0.9665 数值计算模拟：利用 MATLAB 编程对模型(5)进行数值计算模拟，我们得到在不同氮水平下，单作玉米氮素吸收积累的曲线图，见图 4。由图 4 我们可得如下的结论：1)不同施氮量下，玉米整株的氮素吸收积累量均符合 S 型曲线生长特征；抽雄期灌浆期依然是玉米氮素吸收积累的主要集中区。2)当施氮量为 2.6199 gplant1时，玉米的氮素吸收积累速度比施氮量为 0 gplant1，1.84 gplant1，2.45 gplant1和 3.06 gplant1时快，说明此时的施氮量为影响单作玉米氮素吸收积累量的最佳施氮量，与前面的理论推导是一致的。单作玉米140天时的氮素吸收积累量已接近各不同施氮量下的最大值，见表10。Figure 3.Relationship between nitrogen fertilization rate and maximum nitrogen uptake accumulation in maize 图图 3.玉米施氮量与其氮素吸收积累量最大值的关系 普雯瑞，化存才 DOI:10.12677/aam.2024.131034 321 应用数学进展 Figure 4.Changes in nitrogen uptake and accumulation of monoculture maize under different nitrogen levels 图图 4.不同氮水平下单作玉米氮素吸收积累量变化 Table 10.Nitrogen uptake and accumulation of 140 days of maize with dif-ferent nitrogen application rates and monoculture 表表 10.不同施氮量与单作玉米 140 天时氮素吸收积累量 施氮量 x/gplant1 0 1.84 2.45 2.6199 3.06 140 天时氮素吸收积累量(Mm/gplant1)2.25 2.97 3.50 3.5330 3.26 4.不同氮水平下大豆与玉米间作氮素吸收的非自治微分方程模型不同氮水平下大豆与玉米间作氮素吸收的非自治微分方程模型 研究表明，豆科与非豆科作物间作已成为一种农业可持续发展模式，而与豆科作物根系共生的根瘤菌可固定大气中的氮，在满足自身需求的前提下，还可以通过各种途径为伴生的非豆科植物提供氮源16 17。但是，无论在何种施氮量下，大豆和玉米间作时的氮素各自主体吸收积累量仍然具有 Logistic 微分方程中的项，再加上间作时共享土壤中营养的相互作用项。因此，我们可以建立在不同氮水平下的大豆与玉米间作时氮素吸收的非自治微分方程模型。讨论在不同氮水平下，间作对于大豆和玉米各自的氮素吸收积累量的影响，以及最佳供氮量相比单作时是否减少，氮素利用率是否提升等相关问题。4.1.模型的建立模型的建立 从参考文献13中，我们分析提取如下重要事实信息：(F1)在不同氮水平下，间作玉米和大豆各自的氮素吸收积累过程依然符合“S”型增长规律。相比于单作，大豆与玉米间作能够有效的发挥大豆的固氮作用，为间作系统提供除了土壤氮和肥料氮以外的第三氮源，由此不仅能使玉米氮素吸收量增加，而且还能有效地抑制大豆自身的“氮遏制”作用；(F2)由于不同株型的作物占据不同的生态位，玉米的根系分布较浅，大豆根系较深且为直根系，因此，间作模式充分利用了空间生态位的差异，进而提高了整个系统的养分利用率；(F3)玉米的种间竞争能力要强于大豆，说明了玉米更具有间作优势；普雯瑞，化存才 DOI:10.12677/aam.2024.131034 322 应用数学进展 为此，我们可以作出如下合理假设：(H1)在间作系统下，可将大豆和玉米从根部吸收到氮素的土壤空间分为两个：自主空间和共享空间；(H2)在自主土壤空间里，大豆和玉米各自主体对氮素的吸收都符合“S”型的 Logistic 增长模式；(H3)在共享土壤空间，与大豆共生的固氮根瘤菌释放了一定量的氮养分，以供给系统中的其他伴生植物，故而使得大豆的氮素吸收速率减缓，处于弱劣势；同时，玉米则获得更多的氮养分供给，处于竞争优势；因此，可设大豆的氮素吸收积累速率相对于玉米的成反比，而玉米的氮素吸收积累速率相对于大豆的则成正比；(H4)只考虑交错行单株大豆和单株玉米的氮素吸收积累量情况；(H5)忽略间作对于大豆和玉米这两种作物的内禀增长率的微弱影响，内禀增长率依旧取单作时的关系式。根据以上事实和假设，并参考11的建模思想，我们基于单作自主吸氮模型(2)和(5)，建立间作共享吸氮空间模型。设单株大豆与玉米间作时在t时刻的氮素吸收积累量分别为()S t、()M t初值分别为()340.03S=，()340.06M=，氮素吸收的增长率分别为()S Sh rt、()M Mh rt，最大氮素吸收积累量分别为0Sh S、0Mh M，抑制作用影响因子为，促进作用影响因子为，且01。于是，建立在不同氮水平下大豆与玉米间作氮素吸收积累的二维非自治微分方程模型如下：()()d1dd1dS SmmM MmmSSMh rt StSMMMSh rt MtMS=+(7)其中Sh、Mh分别为间作大豆和玉米有关施氮量的影响因子，最大值分别为mSmSh S=，mMmMh M=，其他参数同模型(1)、(4)。4.2.数值模拟及结论数值模拟及结论 根据文献11中的结论，可选定微分方程组模型(7)的参数为0.1994=，0.3804=，满足条件01 2。拟合效果：根据文献13内容，我们整理得到不同供氮水平下，间作大豆及玉米的氮素吸收积累量的相关数据见表 11 及表 12。Table 11.Nitrogen uptake accumulation of intercropping soybean under different nitrogen levels 表表 11.不同氮水平下间作大豆的氮素吸收积累量 时间 day 施氮量 40 48 54 76 93 109 128 132 0 gplant1 0.06 0.09 0.12 0.30 0.53 0.76 0.95 0.98 0.2 gplant1 0.07 0.11 0.14 0.34 0.59 0.85 1.11 1.14 0.4 gplant1 0.09 0.13 0.18 0.46 0.77 1.04 1.24 1.26 0.61 gplant1 0.08 0.12 0.16 0.40 0.68 0.93 1.12 1.15 借助于上述数据，通过 matlab 编程计算分别得到不同氮水平下间作大豆及玉米的实验数据与模拟数据之间的拟合程度见表 13 及表 14。普雯瑞，化存才 DOI:10.12677/aam.2024.131034 323 应用数学进展 Table 12.Nitrogen uptake accumulation of intercropping maize at different nitrogen levels 表表 12.不同氮水平下间作玉米的氮素吸收积累量 时间 day 施氮量 34 48 70 85 99 109 126 140 0 gplant1 0.09 0.21 1.19 2.10 2.50 2.60 2.65 2.65 1.84 gplant1 0.07 0.28 1.53 2.73 3.28 3.42 3.50 3.51 2.45 gplant1 0.18 0.50 1.81 3.01 3.76 4.03 4.23 4.28 3.06 gplant1 0.12 0.37 1.60 2.77 3.43 3.64 3.77 3.80 Table 13.Table of fitting effect of intercropping soybean under different ni-trogen application rates 表表 13.不同施氮量下间作大豆的拟合效果表 施氮量(x/gplant1)0 0.2 0.4 0.61 相关系数(R2)0.9984 0.9932 0.9985 0.9978 Table 14.Table of fitting effect of intercropping maize under different nitro-gen application rates 表表 14.不同施氮量下间作玉米的拟合效果表 施氮量(x/gplant1)0 1.84 2.45 3.06 相关系数(R2)0.9878 0.9986 0.9846 0.9926 结论一：表 13 及表 14 中的相关系数均在 98%以上说明，我们的模型拟合效果很好，误差相对较小，可以较精确地预测不同氮水平下间作大豆及玉米在各生育期内的氮素吸收积累量。数值模拟：对模型(7)进行数值模拟，得到在不同氮水平下间作大豆和玉米氮素吸收积累量变化的相关曲线图，见图 5 及图 6。Figure 5.Changes in nitrogen uptake in intercropping soybean under different nitrogen levels 图图 5.不同氮水平下间作大豆氮素吸收变化 普雯瑞，化存才 DOI:10.12677/aam.2024.131034 324 应用数学进展 Figure 6.Changes in nitrogen uptake in intercropping maize under different nitrogen levels 图图 6.不同氮水平下间作玉米氮素吸收变化 结论二：从图 5 和图 6 中可得知，不同氮水平下间作大豆和玉米的氮素吸收积累量都符合 S 型曲线生长特征，且间作大豆的氮素吸收积累主要集中在结荚鼓粒期，间作玉米的则是抽雄期灌浆期，与单作时保持一致。最大氮素吸收积累量模拟：根据模型(5)，我们分别模拟得到了对应不同施氮量时，间作大豆及玉米的最大氮素吸收积累量，见表 15。对表 15 中的数据进行拟合，得到不同施氮量与大豆和玉米的最大氮素吸收积累量之间的关系式分别为：32322.56041.286