CN105913326B - 约束性知识与精英个体策略遗传算法融合的作物生育期模型品种参数优化方法 - Google Patents
约束性知识与精英个体策略遗传算法融合的作物生育期模型品种参数优化方法 Download PDFInfo
- Publication number
- CN105913326B CN105913326B CN201610220442.6A CN201610220442A CN105913326B CN 105913326 B CN105913326 B CN 105913326B CN 201610220442 A CN201610220442 A CN 201610220442A CN 105913326 B CN105913326 B CN 105913326B
- Authority
- CN
- China
- Prior art keywords
- indicate
- phenological period
- individual
- parameter
- formula
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000004422 calculation algorithm Methods 0.000 title claims abstract description 71
- 238000005457 optimization Methods 0.000 title claims abstract description 48
- 238000000034 method Methods 0.000 title claims abstract description 44
- 230000002068 genetic effect Effects 0.000 title claims abstract description 39
- 230000004927 fusion Effects 0.000 title claims description 18
- 108090000623 proteins and genes Proteins 0.000 claims abstract description 65
- 230000035772 mutation Effects 0.000 claims abstract description 16
- 238000002360 preparation method Methods 0.000 claims abstract description 3
- 230000000694 effects Effects 0.000 claims description 16
- 230000008859 change Effects 0.000 claims description 13
- 238000004088 simulation Methods 0.000 claims description 11
- 238000010276 construction Methods 0.000 claims description 8
- 238000005259 measurement Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 4
- 230000002596 correlated effect Effects 0.000 claims description 4
- 238000011161 development Methods 0.000 claims description 4
- 230000006978 adaptation Effects 0.000 claims description 3
- 229910002056 binary alloy Inorganic materials 0.000 claims description 3
- 239000000463 material Substances 0.000 claims description 3
- 230000012010 growth Effects 0.000 description 18
- 241000209140 Triticum Species 0.000 description 10
- 235000021307 Triticum Nutrition 0.000 description 10
- 238000009826 distribution Methods 0.000 description 10
- 238000012360 testing method Methods 0.000 description 10
- 241000208340 Araliaceae Species 0.000 description 7
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 7
- 235000003140 Panax quinquefolius Nutrition 0.000 description 7
- 235000008434 ginseng Nutrition 0.000 description 7
- 235000013339 cereals Nutrition 0.000 description 6
- 230000008901 benefit Effects 0.000 description 5
- 238000011156 evaluation Methods 0.000 description 5
- 230000006872 improvement Effects 0.000 description 5
- 230000007246 mechanism Effects 0.000 description 5
- 230000018109 developmental process Effects 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 4
- 230000015572 biosynthetic process Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 230000005855 radiation Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000010792 warming Methods 0.000 description 3
- 241000196324 Embryophyta Species 0.000 description 2
- 238000012271 agricultural production Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 2
- 230000000875 corresponding effect Effects 0.000 description 2
- 235000013399 edible fruits Nutrition 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 230000035800 maturation Effects 0.000 description 2
- 230000035479 physiological effects, processes and functions Effects 0.000 description 2
- PXFBZOLANLWPMH-UHFFFAOYSA-N 16-Epiaffinine Natural products C1C(C2=CC=CC=C2N2)=C2C(=O)CC2C(=CC)CN(C)C1C2CO PXFBZOLANLWPMH-UHFFFAOYSA-N 0.000 description 1
- 241001269238 Data Species 0.000 description 1
- 101150017040 I gene Proteins 0.000 description 1
- 240000007594 Oryza sativa Species 0.000 description 1
- 235000007164 Oryza sativa Nutrition 0.000 description 1
- 230000018199 S phase Effects 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 238000009305 arable farming Methods 0.000 description 1
- 238000009395 breeding Methods 0.000 description 1
- 230000001488 breeding effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000008602 contraction Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000010429 evolutionary process Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 229940074869 marquis Drugs 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000010899 nucleation Methods 0.000 description 1
- 235000015816 nutrient absorption Nutrition 0.000 description 1
- 235000015097 nutrients Nutrition 0.000 description 1
- 230000029553 photosynthesis Effects 0.000 description 1
- 238000010672 photosynthesis Methods 0.000 description 1
- 230000000243 photosynthetic effect Effects 0.000 description 1
- 230000001766 physiological effect Effects 0.000 description 1
- 230000008635 plant growth Effects 0.000 description 1
- VBUNOIXRZNJNAD-UHFFFAOYSA-N ponazuril Chemical compound CC1=CC(N2C(N(C)C(=O)NC2=O)=O)=CC=C1OC1=CC=C(S(=O)(=O)C(F)(F)F)C=C1 VBUNOIXRZNJNAD-UHFFFAOYSA-N 0.000 description 1
- 238000004064 recycling Methods 0.000 description 1
- 230000006903 response to temperature Effects 0.000 description 1
- 235000009566 rice Nutrition 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 238000012549 training Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/12—Computing arrangements based on biological models using genetic models
- G06N3/126—Evolutionary algorithms, e.g. genetic algorithms or genetic programming
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B25/00—ICT specially adapted for hybridisation; ICT specially adapted for gene or protein expression
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Biophysics (AREA)
- Theoretical Computer Science (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- General Physics & Mathematics (AREA)
- Genetics & Genomics (AREA)
- Business, Economics & Management (AREA)
- Artificial Intelligence (AREA)
- Computational Linguistics (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Mining & Mineral Resources (AREA)
- Marine Sciences & Fisheries (AREA)
- Animal Husbandry (AREA)
- Physiology (AREA)
- Primary Health Care (AREA)
- Marketing (AREA)
- Biomedical Technology (AREA)
- Human Resources & Organizations (AREA)
- Data Mining & Analysis (AREA)
- Evolutionary Computation (AREA)
- Agronomy & Crop Science (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Biotechnology (AREA)
- Medical Informatics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Breeding Of Plants And Reproduction By Means Of Culturing (AREA)
- Feedback Control In General (AREA)
Abstract
本发明提出一种约束性知识引导精英个体策略遗传算法的作物生育期模型参数优化方法，属于作物模型品种参数优化的最优化方法领域。其步骤：作物生育期模型品种参数优化的数据准备；作物生育期模型品种参数初始范围的等级划分；确定对模型品种参数具有关键约束性的物候期；关键约束性物候期实测数据范围的估算；设置改进型遗传算法控制参数的初始值；生成初始种群；构造适应度函数；构造个体的基因位方向因子；采用轮盘赌方法构造选择算子；构造方向引导交叉算子；构造方向与精英个体引导的变异算子；基于精英个体的局部搜索。本发明克服了效率低、优化结果不符合品种生物学含义等问题，增加了参数优化的准确性和效率，具有科学意义和实用价值。
Description
一、技术领域
本发明属于农业信息技术领域，是领域知识引导进化算法与作物生长模型的交叉领域。涉及一种可用于非线性、不连续、多参数的作物生育期机理模型品种参数估算的自动优化方法，可用于小麦、水稻等禾谷类作物。具体地讲，该方法是一种约束性知识与精英个体策略遗传算法融合的作物生育期模型品种参数优化方法。
二、背景技术
气候变化问题已成为国际社会关注的热门话题，被列为全球十大环境问题之一。20世纪80年代初，科学家们开始发现全球的气候正发生着显著的变化。在2007年气候变化专门委员会(IPCC)发布的第四次评估报告指出：全球气候变暖已是不争的事实。农业生产与气候息息相关，气候变暖导致的极端气候事件趋多、趋强。这势必会对农业及生态环境等领域产生深远的影响。气候变化会导致粮食生产的不稳定。因此，气候变暖对农业生产的影响评价研究近年来受到重视，一方面积极开展作物对气候变化响应的观测试验研究，另一方面通过作物生长模型模拟评估气候变化对作物系统的可能影响，为应对气候变化条件下，农业生产策略与措施变化提供基础数据与决策支持。与观测试验研究相比，利用作物生长模型进行请模拟研究，由于它具有机理强、成本低、时效高、变量易于控制等优点，目前已成为气候变化影响评估的主要研究方法，被世界各地的研究者普遍采用。
作物生长模型是基于作物生理学和生态学原理，运用系统动力学的方法，综合考虑作物品种、栽培措施与气候环境、土壤条件之间的互作关系，量化作物生长发育、光合生产、器官建成及产量形成等生理过程的数学模型。由于作物生长模型充分考虑了作物对土壤、天气、水分、养分等多种影响因子的机理响应，能够模拟不同自然环境和管理措施下作物的生长状况，因此是有效探索气候变化对农业生产的影响，评价作物系统生产力的有力工具。作物生育期模型是作物生长模型的重要组成部分，它通过定量模拟作物对温度和光照的反应，预测和划分作物的发育进程，是作物光合作用、干物质积累与分配、养分吸收与转移、产量品质形成等方面模拟的基础，直接影响作物生长模型的模拟结果。
作物生长模型通过品种参数量化不同品种的遗传与生理特性。由于自然环境复杂性，作物生长模型应用中具有较强的不确定性，需要根据应用品种和应用区域的试验数据，对模型品种参数进行校正。目前作物生长模型参数校正方法主要分为知识驱动法以及数据驱动法。知识驱动法是指作物模型专家以作物模型建模机理、作物栽培经验、作物品种特性等专业知识为引导，采用手工试错或计算机循环迭代的方式，在分阶段、分目标等策略指导下不断试错寻找最优的参数集合。该方法以领域知识为支撑，得到的结果较为可靠，但当应用尺度增大时，为了使模型对全部试验资料均取得较好的模拟效果，调参的数据量及工作量增大，要求调参人员具有较为完备的作物模型知识和丰富的模型调参经验，因而受调参人员的主观因素影响较大；且循环试错的过程需要不断重复迭代，费事费力，参数优化效率较低。数据驱动法是将该问题转化为最优化问题，采用演化计算等全局优化算法，通过大量的随机搜索试验寻找模型参数的最优解。遗传算法是一种常用的随机化、启发式优化算法，具有鲁棒性强、搜索效率高等显著特点，适合求解作物长模型参数优化这样一类结构复杂、非线性程度高、搜索空间巨大且非单峰的问题，已经在ORYZA2000、CropSyst、RiceGrow、SUCROS等模型中取得了应用。现有研究更为关注优化算法本身的改进，如在种群划分时加入动态聚类的思想、采用两级遗传算法的方式进行参数优化等，提升了算法的优化性能。但在这类算法改进中，往往进化寻优的角度增加算法的优化策略，往往控制参数众多，算法时间和空间复杂度变高。此外，作物生育期模型的品种参数反应了作物品种的遗传特性，而基于遗传算法的作物模型参数的校正与优化仅从数值寻优的角度寻找品种参数最优解，往往忽略了作物模型品种参数所表达的生理遗传特性，应用过程中会造成优化后的品种参数出现不符合作物品种生物学意义的现象。
本发明采用两种方法有效的融合的方式，将作物品种及模型知识融合到改进型遗传算法中，设计了一种适用于作物生育期模型品种参数自动校正的方法，对快速而准确的优化作物模型参数具有重要的实际应用价值。
三、发明内容
本发明提出一种约束性知识与精英个体策略遗传算法融合的作物生育期模型参数自动优化方法(Constraints Knowledge and Elite Individual Strategy GeneticAlgorithm Fusion，简称CK-EISGA)。该方法将品种特征知识与改进型精英策略遗传算法有效融合，解决作物生育期模型品种参数自动校正时，由于模型算法多峰非线性以及进化算法优化结果的生物学意义不明确等问题，可快速准确地估算作物生育期模型品种参数。
本发明所述的约束性知识引导精英策略遗传算法的作物生育期模型品种参数优化方法，包括以下步骤：
步骤1：驱动数据准备
作物生育期模型品种参数优化需准备的基础数据包括气象数据(最低温、最高温、日照时数或辐射量)、不同品种在出苗、拔节、抽穗、成熟等物候时期的历史实测数据以及相应的播期、播种深度等实测数据。
步骤2：品种参数初始范围的等级划分
作物生育期模型的品种参数量化了作物生长发育的遗传特性，建模时确定的初始取值范围通常涵盖了某种类型作物的遗传特性，而没有针对一类具体品种遗传特性的强弱进一步细化，会造成优化算法搜索空间的扩大，参数优化后会出现模型模拟值与实测值误差虽小，但品种参数已失去生物学含义的问题。因此，品种参数的初始范围应根据具体品种的遗传特性缩小限制在某一区间内，按公式(1)(2)对品种参数范围进行等级划分。
式中：表示品种参数i原始范围最大值，表示品种参数i原始范围最小值，表示参数i重新划分等级后，第j级特性的最小值，表示参数i重新划分等级后，第j级特性的最大值，N表示参数i品种特性的等级个数。
步骤3：确定对模型品种参数具有关键约束性的物候期
作物生育期模型模拟了多个物候阶段，每个物候阶段通常受到不同发育因子的影响且影响程度不同，造成了各物候阶段对品种参数的约束程度不同。关键物候期是指对待优化品种参数具有明显约束能力的物候期。只有当选用了关键物候期作为模拟目标，并分析与品种参数的数据相关性，才能使优化效果更为准确。
所述步骤3中包括：
步骤3.1定义模型品种参数集合P＝{X1，X2，…，Xn}，模型模拟物候期集合W＝{Y1，Y2，…，Yn}。通过公式(3)衡量物候期Yi与品种参数Xi的约束关系以及正反比关系。若Zi≠0，则Yi约束Xi；若Zi＝0，则Yi不约束Xi。
式中：Xi表示品种参数i的初始值，X′i表示品种参数i改变后的值，Yi，j表示在品种参数i条件下物候期j的日期，Yi，j+1表示品种参数条件i下物候期j+1的日期，Y′i，j表示改变品种参数i后物候期j的日期，Y′i，j+1表示改变品种参数i后物候期j+1的日期，Zi表示品种参数i的敏感度。
步骤3.2定义关键物候期集合K＝{k1，k2，…，kn}，采用以下规则确定作物模型参数优化的关键物候期，：
(1)若某一Xi仅被唯一的Ki约束时，则Ki为关键物候期；
(2)若Xi被Kj至Kk发育阶段中的物候期约束时，则选取Kj和Kk的两个物候期作为关键物候期。
步骤3.3判定品种参数与模型模拟物候期的相关性
根据步骤3.1得到的Z的情况能够确定品种参数Xi与物候期持续天数之间存在一定的相关性，当Z＞0时，Xi与成正相关，当Z＜0时，Xi与成负相关，当Z＝0时，Xi与不相关。
步骤4关键约束性物候期实测数据范围的估算
由于作物生长需要延续若干个月，品种参数优化时所选用的物候期实测数据往往不够完备，会出现缺少步骤2确定的部分关键约束性物候期ki的实测数据，需要利用当地栽培专家的经验知识加以推理估算数据范围。
所述步骤4中包括：
步骤4.1估算关键物候期k1的基础日期
首先判断实测物候期数据是否完全包含关键物候期集合K中全部元素的实测数据，若包含，则跳过该步骤；若缺失某一关键物候期k1的实测数据，则根据公式(4)推理估算该关键约束性物候期k1的物候日期。
式中：Dj表示待估算的物候日期，Di表示已知物候期的日期，dij表示物候期i与物候期j的持续天数，i，j表示物候期索引，i＞j表示物候期j在物候期i之前，其中，dij的值可由文献资料或栽培经验获取。
步骤4.2构造关键物候期实测数据的允许误差范围
由估算得到关键物候期ki的日期具有一定误差。因此，构造关键物候期估算数据的运行误差范围[Dj-E，Dj+E]E为相差范围。
步骤5设置改进型遗传算法控制参数的初始值
主要包括：种群规模BN、最大迭代次数TM、交叉概率Pc、变异概率Pm、待优化参数个数PN、及算法终止条件等；算法终止条件可采用到达最大迭代次数或达到收敛条件，由用户自定义。
步骤6初始种群的生成
种群个体采用实数编码方式，每一个个体对应一套模型品种参数X＝{xi，x2，…，xn}。在步骤2划分的品种参数范围内随机生成BN套个体。
步骤7适应度函数的构造
适应度函数是评价个体优劣的指标，主要由目标函数和惩罚函数两部分构成。
所述步骤7中包括：
步骤7.1利用模型中关键物侯期的模拟值和实际观测值之间的根均方差(RMSE)构造目标函数，公式见(5)。当关键物候期实测数据是由步骤3推算估计获得的范围时，对该范围进行一下处理，公式见(6)。
式中：Pi表示关键物候期i的模拟值，Oi表示关键物候期i的实测值，n表示其物候期实测数据样本数，i表示1，2，…，n，E表示估算物候期允许误差天数，Omax表示误差范围的上限，Omin表示误差范围的下限。
步骤7.2构造惩罚函数，对关键物候期集合K中任一元素ki的模拟效果极差的参数集进行惩罚，降低其传入下一代的概率，惩罚函数见公式(7)。
式中：ftarget表示目标函数值，Pi表示关键物候期ki的模拟值，Oi表示关键物候期ki的物候期实测值，ε表示误差临界值，超过即认为模拟效果极差。
步骤7.3首先采用公式(8)将最小化的目标函数转化为适应度函数，再对各目标的适应度函数值进行加权求和，从而计算个体的适应度函数值。
式中：表示关键物候期j的目标函数值，ffit表示适应度函数值，N表示关键物候期个数，wj表示关键物候期j的权重。
步骤8构造个体的基因位方向因子
传统遗传算法的交叉和变异算子通常缺少方向性。而在交叉和变异过程中，某一基因位的交叉或变异方向与下一代期望调整方向有可能不一致，导致了作物模型参数优化背景下个体更新的无效搜索。由步骤2.3可以获得，作物生育期模型的品种参数的变化规律与各物候期模拟结果存在着正相关或负相关性，评价个体时，能够得到个体物候期模拟结果与实测结果之间的正负关系，利用该正负关系，结合模型参数与物候期间的相关性，能够对个体中每一基因位(即品种参数)下一代更新期望进化的方向加以有效引导，提高全局搜索的效率。
作物模型品种参数个体的基因位期望进化方向可由公式(9)计算。
式中：di表示个体中第i个基因位(品种参数)的期望调整方向，P表示任意一个约束个体中第i基因位的关键物候期模拟值，O表示任意一个约束个体中第i基因位的关键物候期实测值，Zi表示个体中第i个品种参数与物候期的正反比关系，计算公式见公式(3)。
步骤9采用轮盘赌选择算子，首先根据个体适应度值大小，计算个体遗传进入下一代的概率。根据轮盘赌选择原则，从父代中选取部分个体。
步骤10：构造方向引导交叉算子。在模拟二进制交叉算子的基础上，结合个体基因位方向引导因子，对各基因的交叉方向进行引导，如下式所示：
式中，表示两个父代个体的第i位基因，表示两个子代个体的第i位基因，U表示区间[0，1]内的均匀分布随机数，ηc表示交叉参数，表示两个父代个体第i位基因方向因子。
步骤11：构造方向与精英个体引导的变异算子。在非均匀变异算子的基础上，利用个体基因位方向引导因子与种群精英个体综合引导变异算子，从而加快算法收敛到全局最优解的速度，公式见(13)。
式中：X′k表示变异后第k位基因，Xk表示变异前第k位基因，dk表示第k位基因的方向因子，γ表示[0，1]内的均匀分布随机数，t表示当前进化代数，T表示最大进化代数，b表示非均匀度决定系数，c表示加速系数，表示种群中第k位基因的最大值，表示最优个体第k位基因。步骤12：基于精英个体的局部搜索。由于种群精英个体与全局最优解之间的亲和度要大于种群中的其他个体，利用精英个体信息产生多个与精英个体有较大亲和度的个体，可以提高种群优势个体的多样性，增加算法搜索到全局最优解的概率。因此在精英个体附近进行随机扰动，产生一定数量的新个体，提高算法后期局部搜索能力。新个体产生公式如下：
式中：X′k表示新个体第k位基因，表示最优个体第k位基因，表示种群中第k位基因的最大值，)表示种群中第k位基因的最小值，r表示[0，1]内的均匀分布随机数。
步骤13重复步骤7至步骤12，直至算法运行代数达到算法终止条件。
本发明与现有的作物模型参数估算方法相比，其改进主要体现如下：
(1)根据待优化品种的生物学特性对作物生育期模型的品种参数初始范围进行了进一步划分，，缩小遗传算法随机搜索空间，解决了基于遗传算法进行作物模型参数优化可能得到的不符合品种生物学特性的参数集问题。
(2)给出了关键约束性物候期的确定方法和日期范围的估算方法，为作物生育期模型参数优化的目标变量选择提供了有效方法。
(3)构造模型输出物候期与品种参数间的基因位方向引导因子，对遗传算法的交叉变异操作进行引导，提高了进化算法在该类应用中有效搜索能力。
(4)个体变异操作时加入精英个体的信息进行引导，提高了算法向最优解附近靠拢的速度。并对精英个体附近进行局部搜素，提高了种群中优势个体的多样性，增加算法搜索到全局最优解的概率。
有益效果：
本发明提出一种约束性知识与精英个体策略遗传算法融合的作物生育期模型品种参数优化方法，该方法一方面将作物品种特性的相关约束性知识引入到改进型遗传算法优化框架中，得出约束参数的初始范围和关键物候期的实测值范围，使得到的品种参数更符合作物品种的生物学特性。另一方面通过改进具有方向性的精英个体遗传算法的种群更新策略，提高了算法寻找最优解的有效性和寻优效率。
实际应用表明，本发明克服了传统参数优化方法中由于大量的迭代试错而带来的优化效率低下、调参周期长等问题，并解决了由于传统遗传算法的随机搜索机制而带来的优化结果不符合品种生物学含义的问题，增强了模型品种参数优化的准确性，提高了作物模型品种参数的优化效率，有助于作物生长模型的发展与应用，具有科学意义和实用价值。
四、说明书附图
图1本发明的流程示意图
图2基于精英个体策略遗传算法与蜜蜂进化型遗传算法的种群进化分布图
图3各项约束性知识的改进效果(图3-1为作物模型品种范围约束对参数PVT分布结果的影响，图3-2为关键约束性物候期实测数据约束对品种参数IE分布结果的影响。图3-3为方向因子对收敛速率的影响)
图4为本方法在徐州、济宁、潍坊、赣榆四个地点模拟效果的1∶1关系图；
图5为徐州、济宁、潍坊、赣榆四地的历史实测数据(图5-1为开花期和成熟期的实测数据，图5-2为四个地点的平均温度变化图，图5-3四个地点的辐射量变化图)。
五、具体实施方式
以下通过实施案例对本发明进一步描述
实施例
以南京农业大学的WheatGrow小麦生育期模型模拟济南13号小麦品种在山东省济宁市、胶州市、潍坊市和江苏省徐州市四个地点种植的物候期变化情况。结合附图1，具体说明本实施方式：步骤1：作物生育期模型品种参数优化的数据准备
选取济南13号小麦品种在山东省济宁市、胶州市、潍坊市和江苏省徐州市四个地点种植的历史数据，试验地点及年份数据见表1，各套数据均记录开花期和成熟期的实测日期，见附图5-1。气象数据均来自当地气象站，包括日最高温度(℃)、日最低温度(℃)、辐射量(MJ/m2))等，见附图5-2和5-3。播种深度为2.5cm。
表1试验选用的地点及年份数据
步骤2：WheatGrow小麦生育期模型品种参数范围划分。
小麦为冷季作物，其生育期的特性包括感温性、感光性、基本早熟型和春化效应，WheatGrow的品种参数TS、PS、IE和PVT分别表现了以上四种特性，根据公式(1)对品种参数初始范围进行划分，结果见表2。
表2 WheatGrow小麦生育期模型参数初始范围划分及生物学含义
步骤3：确定WheatGrow小麦生育期模型具有关键约束性的物候期
步骤3.1分析模型模拟输出物候期Yi到Yi+1的持续天数受到品种参数Xi变化所带来的影响。根据公式(3)计算获得Zi，具体数据见表3。
表3 WheatGrow小麦生育期模型各物候期持续天数随参数变化的情况
步骤3.2分析物候期Yi和品种参数Xi的约束关系发现，参数PVT被拔节期唯一约束，参数FD被成熟期唯一约束，IE和PS被拔节至抽穗阶段约束，而TS则影响整个出苗至成熟阶段。因此根据调参关键物候期的确定规则，确定WheatGrow小麦生育期模型的拔节期、抽穗期和成熟期为调参关键物候期。
步骤3.3确定品种参数与模型模拟物候期的相关性
根据步骤2.1得到的Zi的情况能够确定品种参数Xi与物候期Yi持续天数之间存在一定的相关性。根据表4中Zi值，确定品种参数TS、PS和PVT与物候期持续天数之间呈正比关系，品种参数IE、FD与物候期持续天数之间呈反比关系。
步骤4关键约束性物候期实测数据范围的估算
徐州等四个地点的实测数据仅包含开花期和成熟期，缺少抽穗期和拔节期的实测数据，需要利用当地栽培专家的经验知识加以推理估算数据范围。
步骤4.1估算关键物候期中抽穗期和拔节期的基础日期
根据文献资料搜集和专家调研，确定济南13号品种在江苏徐州等地种植时，拔节期距抽穗期约35天，抽穗期距开花期约5天。徐州地点实测数据中1982年的开花期为5月2日，则根据公式(4)，估算抽穗期为4月27日，拔节期为3月23日。其他地点和年份的估算方法与上例相同。
步骤4.2构造关键物候期实测数据的允许误差范围
由估算得到关键约束性物候期中拔节期和抽穗期的实测数据具有一定误差。因此，需构造拔节期和抽穗期的实测数据的允许误差范围。根据当地栽培经验，抽穗期和拔节期的误差分别设置为1天和3天。
步骤5设置算法的初始控制参数
设置种群规模BN＝100、最大迭代次数TM＝100、交叉概率Pc＝0.8、变异概率Pm＝0.05、待优化参数个数PN＝5，算法终止条件为迭代次数达到最大迭代次数。
步骤6初始种群的生成
种群个体的编码方式采用实数编码，每一组个体对应一套作物模型品种参数，每个基因位对应一个品种参数。由于济南13号品种的春化特性为半冬性，因此根据表3，将PVT参数的初始范围设置为[30，40]，其他参数为初始参数。在该范围内随机生成BN个个体。
步骤7适应度函数的构造
适应度函数是评价个体优劣的指标，主要由目标函数和惩罚函数两部分构成。
所述步骤7中包括：
步骤7.1目标函数设计
利用模型模拟值和实际观测值之间的根均方差(RMSE)构造目标函数，公式见(5)。当物候期实测数据是由步骤3推算估计获得的范围时，对该范围进行一下处理，公式见(6)。
式中：Pi表示物候期模拟值，Oi表示物候期实测值，n表示物候期实测数据样本数，i表示1，2，…，n，E表示估算物候期允许误差天数，Omax表示估算范围的上限，Omin表示估算范围的下限。
步骤7.2构造惩罚函数，对关键物候期集合K中任一元素ki的模拟效果极差的参数集进行惩罚，降低其传入下一代的概率，惩罚函数见公式(7)。
式中：ftarget表示目标函数值，Pi表示物候期模拟值，Oi表示物候期实测值，ε表示误差临界值，本例取5。
步骤7.3首先采用公式(8)将最小化的目标函数转化为适应度函数，再对各目标的适应度函数值进行加权求和，从而计算个体的适应度函数值。
式中：表示物候期j的目标函数值，ffit表示适应度函数值，N表示关键物候期个数，wj表示物候期j的权重。
步骤8构造方向因子
传统遗传算法的交叉和变异算子通常缺少方向性。而在交叉变异等算子中，某一基因位的交叉或变异方向与其期望调整方向有可能不一致，导致了无效搜索。由于作物生育期模型的品种参数的变化规律与各物候期模拟结果存在一定相关性，主要表现为正相关和负相关，由步骤3.3可以获得。在评价个体时，能够得到个体物候期模拟结果与实测结果之间的正负关系，而利用该正负关系，结合模型参数与物候期间的相关性，能够对个体中每一基因位下一步期望进化的方向进行引导。
作物模型品种参数个体中基因位期望进化的方向可由公式(9)
式中：di表示个体中第i个品种参数的期望调整方向，P表示任意一个约束个体中第i基因位的关键物候期模拟值，O表示任意一个约束个体中第i基因位的关键物候期实测值，Zi表示个体中第i个品种参数与物候期的正反比关系，计算公式见步骤3.3。
步骤9采用轮盘赌选择算子，首先根据个体适应度值大小，计算个体遗传进入下一代的概率。根据轮盘赌选择原则，从父代中选取部分个体。
步骤10：构造方向引导交叉算子。在模拟二进制交叉算子的基础上，结合个体基因位的方向引导因子，对各基因的交叉方向进行引导，如下式所示：
式中，表示两个父代个体的第i位基因，表示两个子代个体的第i位基因，U表示区间[0，1]内的均匀分布随机数，ηc表示交叉参数，一般取值为2，表示两个父代个体第i位基因方向因子。
步骤11：构造方向和精英个体引导的变异算子。在非均匀变异算子的基础上，利用个体基因位的方向引导因子和种群精英个体引导变异算子，从而加快算法收敛到全局最优解的速度，公式见(13)。
式中：X′k表示变异后第k位基因，Xk表示变异前第k位基因，dk表示第k位基因的方向因子，r表示[0，1]内的均匀分布随机数，t表示当前进化代数，T表示最大进化代数，b表示非均匀度决定系数，c表示加速系数，表示最优个体第k位基因，表示种群中第k位基因的最大值。
步骤12：基于精英个体的局部搜索。种群精英个体与全局最优解之间的亲和度要大于种群中的其他个体，利用精英个体信息产生多个与精英个体有较大亲和度的个体，可以提高种群优势个体的多样性，增加算法搜索到全局最优解的概率。因此在精英个体附近进行随机扰动，产生一定数量的新个体，提高算法后期局部搜索能力。新个体产生公式如下：
式中：X′k表示新个体第k位基因，表示最优个体第k位基因，表示种群中第k位基因的最大值，表示种群中第k位基因的最小值，r表示[0，1]内的均匀分布随机数。
步骤13重复步骤7至步骤12，直至算法运行代数达到最大迭代次数。
基于精英个体策略的改进型遗传算法(EISGA)的改进效果见表4和附图2至附图4。表4为基于精英个体策略的改进型遗传算法与其他遗传算法在f1-f5五个测试函数上的平均进化代数及收敛率的比较结果；附图2为EISGA与蜜蜂进化型遗传算法在有关过程中种群的分布情况。试验结果能够发现，基于精英个体策略的改进型遗传算法的收敛速度快且精度更高，表明变异算子及精英个体局部搜索等策略增强了算法的寻优能力；进化过程的种群分布情况表明了EISGA具有更强的后期在最优解処的局部搜索能力，使算法搜索到全局最优解的概率增大。
f1＝100(x1 2-x2)2+(1-x1)2，x1，x2∈[-5，5]
f2＝(-1)(x2+2y2-0.3cos(3x)-0.4cos(3y))+4，x，y∈[-1，1]
f4＝|(1-x)x2sin(200x)|，x∈[0，1]
表4各测试函数的平均进化代数及其收敛率
约束性知识与精英个体策略遗传算法融合的优化方法在作物生育期模型参数优化中的各项改进如表5、附图3和附图4；表5为分别采用关键约束目性物候期的估算日期和估算误差范围构造适应度函数的优化结果对比；图3-1为作物模型品种范围约束对参数PVT分布结果的影响，图3-2为关键约束性物候期实测数据约束对品种参数IE分布结果的影响。图3-3为徐州和济宁两地试验中，方向因子对收敛速率的影响。图4为本方法在徐州、济宁、潍坊、赣榆四个地点模拟效果的1∶1关系图；试验结果能够发现，在作物模型品种参数优化中，在加入PVT品种参数范围约束后，能够有效的约束PVT的优化结果均符合半冬性的品种特性；加入关键约束性物候期的实测数据约束后，品种参数IE的分布也有明显约束，并且采用关键物候期估算范围构造适应度函数的优化效果较好，表明估算关键物候期的实测数据范围对品种参数的分布具有明显的约束作用，并且能够减小经验估算的误差；加入方向因子后，算法的寻优速度有所提高，表明了针对该应用问题方向算子是有效的。算法优化结果在四个地点均取得了良好的优化效果，证明了本方法的有效性。
表5分别采用估算的关键物候期的物候日期和物候范围构造适应度函数的对比结果
Claims (9)
1.约束性知识与精英个体策略遗传算法融合的作物生育期模型品种参数优化方法，其主要包括以下步骤：
1)作物生育期模型品种参数优化的数据准备；
2)作物生育期模型品种参数初始范围的等级划分；
3)确定对模型品种参数具有关键约束性的物候期；
4)关键约束性物候期实测数据范围的估算；
5)设置改进型遗传算法控制参数的初始值；
6)初始种群的生成；
7)适应度函数的构造；
8)构造个体的基因位方向因子；
9)采用轮盘赌方法构造选择算子；
10)构造方向引导交叉算子；
11)构造方向与精英个体引导的变异算子；
12)基于精英个体的局部搜索；
13)重复步骤7至步骤12，直至算法运行代数达到算法终止条件。
2.根据权利要求1所述的约束性知识与精英个体策略遗传算法融合的作物生育期模型品种参数优化方法，其特征在于步骤2)中作物生育期模型品种参数初始范围的等级划分，按公式(1)(2)对品种参数范围进行等级划分；
式中：表示品种参数i原始范围最大值，表示品种参数i原始范围最小值，表示参数i重新划分等级后，第j级特性的最大值，表示参数i重新划分等级后，第j级特性的最小值，n表示参数i品种特性的等级个数。
3.根据权利要求1所述的约束性知识与精英个体策略遗传算法融合的作物生育期模型品种参数优化方法，其特征在于步骤3)中确定对模型品种参数具有关键约束性的物候期，包括：
步骤3.1 定义模型品种参数集合P＝{X1,X2,…,Xn}，模型模拟物候期集合W＝{Y1,Y2,…,Yn}；采用以下方法确定作物模型模拟物候期对品种参数的约束能力：通过公式(3)获取品种参数i的敏感度Zi，Zi的值以衡量物候期Yi与品种参数Xi的约束关系以及正反比关系：若Zi≠0，则Yi约束Xi；若Zi＝0，则Yi不约束Xi；
式中：Xi表示品种参数初始值，X′i表示品种参数改变后的值，Yi表示物候期i日期，Yi+1表示物候期i+1的日期，Y′i表示改变后Xi后物候期i的日期，Y′i+1表示改变后Xi后物候期i+1的日期，Zi表示品种参数i的敏感度；
步骤3.2 定义关键物候期集合K＝{k1，k2，…，kn}，采用以下规则确定作物模型参数优化的关键物候期，若物候期Yi满足以下两个条件之一时认为该物候期为关键物候期，将该物候期Yi加入关键物候期集合K中：
(1)若某一Xi仅被唯一的Ki约束时，则Ki为关键物候期；
(2)若Xi被Kj至Kk发育阶段中的物候期约束时，则选取Kj和Kk的两个物候期作为关键物候期；
步骤3.3 判定品种参数与模型模拟物候期的相关性
根据步骤3.1得到的Z的情况能够确定品种参数Xi与物候期Yi持续天数之间存在一定的相关性，当Z>0时，Xi与Yi成正相关，当Z<0时，Xi与Yi成负相关，当Z＝0时，Xi与Yi不相关。
4.根据权利要求1所述的约束性知识与精英个体策略遗传算法融合的作物生育期模型品种参数优化方法，其特征在于步骤4)中关键约束性物候期实测数据范围的估算，包括：
步骤4.1 估算关键物候期k1基础日期
首先判断实测物候期数据是否完全包含关键物候期集合K中全部元素的实测数据，若包含，则跳过该步骤；若缺失某一关键物候期k1的实测数据，则根据公式(4)推理估算该关键约束性物候期k1的物候日期；
式中：Dj表示待估算的物候日期，Di表示已知物候期的日期，dij表示物候期i与物候期j的持续天数，i，j表示物候期索引，i>j表示物候期j在物候期i之前，其中，dij的值由文献资料或栽培经验获取；
步骤4.2 构造关键物候期实测数据的允许误差范围
构造关键目标物候期估算数据的运行误差范围[Dj-E,Dj+E]，E为相差范围。
5.根据权利要求1所述的约束性知识与精英个体策略遗传算法融合的作物生育期模型品种参数优化方法，其特征在于步骤7)中适应度函数的构造，主要由目标函数和惩罚函数两部分构成，包括：
步骤7.1 利用模型模拟值和实际观测值之间的根均方差构造目标函数ftarget，公式见(5)；当物候期实测数据是由步骤3推算估计获得的范围时，对该范围进行一下处理，公式见(6)；
式中：Pi表示物候期模拟值，Oi表示物候期实测值，n表示物候期实测数据样本数，i表示1,2，…，n，E表示估算物候期允许误差天数，Omax表示估算范围的上限，Omin表示估算范围的下限；
步骤7.2 构造惩罚函数，对关键物候期集合K中任一元素ki的模拟效果极差的参数集进行惩罚，降低其传入下一代的概率，惩罚函数见公式(7)；
式中：ftarget表示目标函数值，Pi表示物候期模拟值，Oi表示物候期实测值，ε表示误差临界值，超过ε即认为模拟效果极差；
步骤7.3首先采用公式(8)将最小化的目标函数转化为适应度函数，再对各目标的适应度函数值进行加权求和，从而计算个体的适应度函数值；
式中：表示物候期j的目标函数值，ffit表示适应度函数值，N表示关键物候期个数，wj表示物候期j的权重。
6.根据权利要求3所述的约束性知识与精英个体策略遗传算法融合的作物生育期模型品种参数优化方法，其特征在于步骤8)中构造个体的基因位方向因子；作物模型品种参数个体的基因位期望进化方向由公式(9)计算：
式中：di表示个体中第i个基因位的期望调整方向，P表示任意一个约束个体中第i基因位的关键物候期模拟值，O表示任意一个约束个体中第i基因位的关键物候期实测值，Zi的值表示个体中第i个品种参数与物候期的正反比关系，计算公式见公式(3)。
7.根据权利要求1所述的约束性知识与精英个体策略遗传算法融合的作物生育期模型品种参数优化方法，其特征在于包括步骤10)中构造方向引导交叉算子；在模拟二进制交叉算子的基础上，结合个体基因位方向引导因子，对各基因的交叉方向进行引导，如下式所示；
式中，表示两个父代个体的第i位基因，表示两个子代个体的第i位基因，
u表示区间[0,1]内的均匀分布随机数，ηc表示交叉参数，取值为2，表示两个父代个体第i位基因方向因子。
8.根据权利要求1所述的约束性知识与精英个体策略遗传算法融合的作物生育期模型品种参数优化方法，其特征在于步骤11)中构造方向与精英个体引导的变异算子；在非均匀变异算子的基础上，利用个体基因位方向引导因子与种群精英个体综合引导变异算子，从而加快算法收敛到全局最优解的速度，公式见(13)：
式中：X′k表示变异后第k位基因，Xk表示变异前第k位基因，dk表示第k位基因的方向因子，r表示[0,1]内的均匀分布随机数，t表示当前进化代数，T表示最大进化代数，b表示非均匀度决定系数，c表示加速系数，dk表示第k位基因的方向引导因子，表示种群中第k位基因的最大值，表示最优个体第k位基因。
9.根据权利要求1所述的约束性知识与精英个体策略遗传算法融合的作物生育期模型品种参数优化方法，其特征在于步骤12)中基于精英个体的局部搜索；在精英个体附近进行随机扰动，产生一定数量的新个体，提高算法后期局部搜索能力；产生新个体的个数由种群大小决定，取种群大小的1/10；新个体产生公式如下：
式中：X′k表示新个体第k位基因，表示最优个体第k位基因，表示种群中第k位基因的最大值，表示种群中第k位基因的最小值，r表示[0,1]内的均匀分布随机数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610220442.6A CN105913326B (zh) | 2016-04-06 | 2016-04-06 | 约束性知识与精英个体策略遗传算法融合的作物生育期模型品种参数优化方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610220442.6A CN105913326B (zh) | 2016-04-06 | 2016-04-06 | 约束性知识与精英个体策略遗传算法融合的作物生育期模型品种参数优化方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105913326A CN105913326A (zh) | 2016-08-31 |
| CN105913326B true CN105913326B (zh) | 2019-08-06 |
Family
ID=56744843
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610220442.6A Expired - Fee Related CN105913326B (zh) | 2016-04-06 | 2016-04-06 | 约束性知识与精英个体策略遗传算法融合的作物生育期模型品种参数优化方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105913326B (zh) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110503287A (zh) * | 2018-05-18 | 2019-11-26 | 南京农业大学 | 基于形态与混合梯度的动态时间弯曲度量水稻物候期多元气象数据相似性方法 |
| CN109472320B (zh) * | 2018-11-30 | 2021-09-07 | 南京农业大学 | 不确定条件下作物生育期模型品种参数自动校正方法 |
| CN114444381A (zh) * | 2021-12-30 | 2022-05-06 | 南京诺源医疗器械有限公司 | 一种适用于光源的控制方法及系统 |
| CN115828438B (zh) * | 2023-02-20 | 2023-05-05 | 中汽研汽车检验中心（天津）有限公司 | 汽车极限性能预测方法、介质、设备 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102393634A (zh) * | 2011-11-28 | 2012-03-28 | 浙江工业大学 | 一种基于pso-ga混合优化算法的单目标多生产线调度控制方法 |
| CN103473597A (zh) * | 2013-09-05 | 2013-12-25 | 中国寰球工程公司 | 一种基于遗传算法的天然气液化流程工艺参数优化方法 |
| CN104573820A (zh) * | 2014-12-31 | 2015-04-29 | 中国地质大学（武汉） | 一种用于带约束条件的工程优化问题的遗传算法 |
| CN104770275A (zh) * | 2015-04-14 | 2015-07-15 | 天津大学 | 基于种植结构调整的农业灌溉用水优化配置方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9183742B2 (en) * | 2012-10-26 | 2015-11-10 | Xerox Corporation | Methods, systems and processor-readable media for optimizing intelligent transportation system strategies utilizing systematic genetic algorithms |
-
2016
- 2016-04-06 CN CN201610220442.6A patent/CN105913326B/zh not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102393634A (zh) * | 2011-11-28 | 2012-03-28 | 浙江工业大学 | 一种基于pso-ga混合优化算法的单目标多生产线调度控制方法 |
| CN103473597A (zh) * | 2013-09-05 | 2013-12-25 | 中国寰球工程公司 | 一种基于遗传算法的天然气液化流程工艺参数优化方法 |
| CN104573820A (zh) * | 2014-12-31 | 2015-04-29 | 中国地质大学（武汉） | 一种用于带约束条件的工程优化问题的遗传算法 |
| CN104770275A (zh) * | 2015-04-14 | 2015-07-15 | 天津大学 | 基于种植结构调整的农业灌溉用水优化配置方法 |
Non-Patent Citations (1)
| Title |
|---|
| 基于个体优势遗传算法的水稻生育期模型参数优化;庄嘉祥 等;《中国农业科学》;20131231;第46卷(第11期);第2220-2231页 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105913326A (zh) | 2016-08-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Martre et al. | Model-assisted phenotyping and ideotype design | |
| CN112906298B (zh) | 一种基于机器学习的蓝莓产量预测方法 | |
| US20220075344A1 (en) | A method of finding a target environment suitable for growth of a plant variety | |
| CN105913326B (zh) | 约束性知识与精英个体策略遗传算法融合的作物生育期模型品种参数优化方法 | |
| CN107329511A (zh) | 基于适宜根温区间的水培蔬菜光环境高效调控方法与系统 | |
| CN106650212A (zh) | 基于数据分析的植物智能培育方法及系统 | |
| Lyu et al. | Multi-objective winter wheat irrigation strategies optimization based on coupling AquaCrop-OSPy and NSGA-III: A case study in Yangling, China | |
| CN116595333B (zh) | 土壤-气候智慧型的水稻目标产量与氮肥用量确定方法 | |
| CN117575094B (zh) | 一种基于数字孪生的农作物产量预测与优化方法及设备 | |
| CN105389452A (zh) | 基于神经网络的黄瓜全程光合速率预测模型及建立方法 | |
| Sánchez-Molina et al. | Support system for decision making in the management of the greenhouse environmental based on growth model for sweet pepper | |
| Coussement et al. | Turgor-driven plant growth applied in a soybean functional–structural plant model | |
| CN116258060A (zh) | 基于机器学习的测土配方施肥方法 | |
| CN113553767A (zh) | 一种温室作物光合速率预测模型构建方法及系统 | |
| CN109472320B (zh) | 不确定条件下作物生育期模型品种参数自动校正方法 | |
| CN105654242A (zh) | 一种基于鱼群算法的黄瓜幼苗期二氧化碳优化调控模型及建立与应用 | |
| Zhang et al. | Coupling of environmental factors and growth stages in simulation of maize biomass allocation | |
| KR20240012287A (ko) | Rgb led 광원과 환경 데이터를 이용한 기계 학습 기반 작물 생육 최적화 시스템 | |
| Chen et al. | A water-saving irrigation decision-making model for greenhouse tomatoes based on genetic optimization TS fuzzy neural network | |
| Rani et al. | Analysis of Five Mathematical Models for Crop Yield Prediction. | |
| HOU et al. | Quantitative design of yield components to simulate yield formation for maize in China | |
| Casadebaig et al. | Increased genetic diversity improves crop yield stability under climate variability: a computational study on sunflower | |
| CN105427061A (zh) | 基于改进型鱼群算法的番茄幼苗期光合作用优化调控模型及建立与应用 | |
| Chandrvavanshi et al. | Simulation of crop growth and productivity using simulation model for short duration rice | |
| AYDINHAN et al. | Use of some multicriteria decision-making methods such as grey relational analysis (GRA), the complex proportional assessment (COPRAS), and weighted aggregated sum product assessment (WASPAS) in selection of some Anatolian pine (Pinus nigra Arnold.) origi |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| DD01 | Delivery of document by public notice |
Addressee: Nanjing Agricultural University Document name: Notification of Passing Examination on Formalities |
|
| DD01 | Delivery of document by public notice | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190806 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |